基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.

新疆萨吾尔金矿带构造-流体地球化学及成矿作用

The Sawuer gold belt is located in the transition belt between Siberian plate and Kazakhstan-Junggar plate. Based on the geological and geochemical studies on the Kuoerzhenkuola and Buerkesidai gold deposits, in Sawuer gold mineralization belt, the time-space structure of mineralization and mineralizing factor are studied, the metallogenic regularity is concluded in thistheses. The ore bodies have the regularity that orebody are of the extensive and compressive in the sallow and depth of volcanic apparatue, respectively, and the vertical extension of orebody is more intensive than the horizontal extension. The gold deposits were controlled by the fractures of volcanic apparatus and regional faults, and featured by the hydrothermal alteration and metasomatism type disseminated mineralization and filling type vein mineralization. By virtue of the geological and geochemical studies on the two deposits that the formation of the two deposits are significantly related to the volcanic activity, we propose new ideas about their origin: (1) the two deposits are located in the same strata, and share the same genesis. (2) both of two deposits are volcanogenic late-stage hydrothermal gold deposits. Based on mapping of volcanic lithofacies and structure for the first time, it is discovered that a volcanic apparatus existed in the study area. Volcanic-intrusive activity can be divided into three cycles and nine lithofacies. where the two deposits are hosted in the same volcanic cycle, in this case, the wall-rock should belong to the same strata. The 40Ar-39Ar age method is employed in this work to analyze the fluid inclusions of quartz in the ore bodies from Kuoerzhenkuola and Buerkesidai gold deposits. The results show that the main mineralization occurred in 332.05 + 2.02-332.59 + 0.5IMa and 335.53 + 0.32Ma~336.78 + 0.50Ma for Kuoerzhenkuola and Buerkesidai gold deposits respectively, indicating that the two deposits are formed almost at the same time, and the metallogenic epoch of the tow deposits are close to those of the hosting rocks formed by volcanic activity of Sawuer gold belt. This geochronological study supplies new evidence for determining the timing of gold mineralization, the geneses of gold deposits? and identifies that in Hercynian period, the Altai developed tectonic-magmatic-hydrothermal mineralization of Early Carboniferous period, except known two metallogenic mineralization periods including tectonic-magmatic-hydrothermal mineralization of Devonian period and Late Carboniferous-Permian period respectively. The study of fluid inclusions indicates that the ore-forming fluid is a type of NaCl-HbO fluid with medium-low temperature and low salinity, Au is transported by the type of auric-sulfur complex (Au (HS)2-), the ore is formed in reduction condition. Hydrogen and oxygen isotopes of fluid inclusions in the major mineralizating stage show that the solutions mainly originated from magmatic water and meteoric water. The fluid mixing and water-rock reaction cause the deposition of Au. The helium and argon isotope compositions of fluid inclusions hosted in pyrite have been measured from Kuoerzhenkuola and Buerkesidai gold deposits in Sawuer gold belt. The results show that the ore-forming fluids of two deposits possessed the same source and is a mixture of mantle- and partial meteoric water-derived fluid, and the reliability of He and Ar isotopic compositions in Hercynian period is discussed. Isotopic studies including H, O, He, C, S, Pb and Sr reveal the same result that the ore-forming fluids of two deposits possessed the same source: the water derived mainly from magmatic water, partially from meteoric water; the mineralizers and ore materials derived mainly from mantle beneath the island arc, and partially from crust. The ore-forming fluids of two deposits are a mixture of mantle-derived fluids being incorporated by crust-derived fluid, and shallow partial meteoric water. Based on these results, it is proposed that the geneses of the two gold deposits are the same, being volcanogenic late-stage hydrothermal gold deposits that the ore-forming fluids filled in fractures of volcanic apparatus and metasomatized the host rocks in the volcanic apparatus. It is the first time we carried out the geophysical exploration, that is, the EH-4 continuous electrical conductivity image system measurement, the results show that relative large-size mineralizing anomalies in underground have been discovered.lt can confirm the law and genesis of the deposits mentioned above, and change the two abandoned mines to current large-size potenial exploration target.

西藏高原草原化嵩草草甸生态系统呼吸及碳平衡

The soil respiration and net ecosystem productivity of Kobresia littledalei meadow ecosystem was investigated at Dangxiong grassland station, one grassland field station of Lhasa Plateau Ecosystem Research Station. Soil respiration and soil heterotrophic respiration were measured at the same time by using Li6400-09 chamber in growing season of year 2004. The response of soil respiration and its components, i.e. microbial heterotrophic respiration and root respiration to biotic and abiotic factors were addressed. We studied the daily and seasonal variation on Net Ecosystem carbon Exchange (NEE) measured by eddy covariance equipments and then the regression models between the NEE and the soil temperature. Based on the researches, we analyzed the seasonal variation in grass biomass and estimated NEE combined the Net Ecosystem Productivity with heterogeneous respiration and then assessed the whether the area is carbon source or carbon sink. 1.Above-ground biomass was accumulated since the grass growth started from May; On early September the biomass reached maximum and then decreased. The aboveground net primary production (ANPP) was 150.88 g m~" in 2004. The under-ground biomass reached maximum when the aboveground start to die back. Over 80% of the grass root distributed at the soil depth from 0 to 20cm. The underground NPP was 1235.04 g m"2.. Therefore annual NPP wasl.385X103kg ha"1, i.e.6236.6 kg C ha"1. 2. The daily variation of soil respiration showed single peak curve with maximum mostly at noon and minimum 4:00-6:00 am. Daily variations were greater in June, July and August than those in September and October. Soil respiration had strong correlation with soil temperature at 5cm depth while had weaker correlation with soil moisture, air temperature, surface soil temperature, and so on. But since early September the soil respiration had a obviously correlation with soil moisture at 5cm depth. Biomass had a obviously linearity correlation with soil respiration at 30th June, 20th August, and the daytime of 27th September except at 23lh October and at nighttime of 27th September. We established the soil respiration responding to the soil temperature and to estimate the respiration variation during monsoon season (from June through August) and dry season (May, September and October). The regression between soil respiration and 5cm soil temperature were: monsoon season (June through August), Y=0.592expfl()932\ By estimating , the soil daily respiration in monsoon season is 7.798gCO2m"2 and total soil respiration is 717.44 gCC^m" , and the value of Cho is 2.54; dry season (May, September and October), Y=0.34exp°'085\ the soil daily respiration is 3.355gCO2m~2 and total soil respiration is 308.61 gCC^m", and the value of Cho is 2.34. So the total soil respiration in the grown season (From May to October) is 1026.1 g CO2IT1"2. 3. Soil heterogeneous respiration had a strong correlation with soil temperature especially with soil temperature at 5cm depth. The variation range in soil heterogeneous respiration was widely. The regression between soil heterogeneous respiration and 5cm soil temperature is: monsoon season, Y=0.106exp ' 3x; dry season, Y=0.18exp°"0833x.By estimating total soil heterotrophic respiration in monsoon season is 219.6 gCC^m"2, and the value of Cho is 3.78; While total soil heterogeneous respiration in dry season is 286.2 gCCbm"2, and the value of Cho is 2.3. The total soil heterotrophic respiration of the year is 1379.4kg C ha"1. 4. We estimated the root respiration through the balance between soil respiration and the soil heterotrophic respiration. The contribution of root respiration to total respiration was different during different period: re-greening period 48%; growing period 69%; die-back period 48%. 5. The Ecosystem respiration was relatively strong from May to October, and of which the proportion in total was 97.4%.The total respiration of Ecosystem was 369.6 g CO2 m" .we got the model of grass respiration respond to the soil temperature at 5cm depth and then estimated the daytime grass respiration, plus the nighttime NEE and daytime soil respiration. But when we estimated the grass respiration, we found the result was negative, so the estimating value in this way was not close. 6. The estimating of carbon pool or carbon sink. The NPP minus the soil heterogeneous respiration was the NEE, and it was 4857.3kg C o ha"1, which indicated that the area was the carbon sink.

陆相砂泥岩层序的高精度检测方法研究及应用

Ordos Basin is a typical cratonic petroliferous basin with 40 oil-gas bearing bed sets. It is featured as stable multicycle sedimentation, gentle formation, and less structures. The reservoir beds in Upper Paleozoic and Mesozoicare are mainly low density, low permeability, strong lateral change, and strong vertical heterogeneous. The well-known Loess Plateau in the southern area and Maowusu Desert, Kubuqi Desert and Ordos Grasslands in the northern area cover the basin, so seismic data acquisition in this area is very difficult and the data often takes on inadequate precision, strong interference, low signal-noise ratio, and low resolution. Because of the complicated condition of the surface and the underground, it is very difficult to distinguish the thin beds and study the land facies high-resolution lithologic sequence stratigraphy according to routine seismic profile. Therefore, a method, which have clearly physical significance, based on advanced mathematical physics theory and algorithmic and can improve the precision of the detection on the thin sand-peat interbed configurations of land facies, is in demand to put forward.Generalized S Transform (GST) processing method provides a new method of phase space analysis for seismic data. Compared with wavelet transform, both of them have very good localization characteristics; however, directly related to the Fourier spectra, GST has clearer physical significance, moreover, GST adopts a technology to best approach seismic wavelets and transforms the seismic data into time-scale domain, and breaks through the limit of the fixed wavelet in S transform, so GST has extensive adaptability. Based on tracing the development of the ideas and theories from wavelet transform, S transform to GST, we studied how to improve the precision of the detection on the thin stratum by GST.Noise has strong influence on sequence detecting in GST, especially in the low signal-noise ratio data. We studied the distribution rule of colored noise in GST domain, and proposed a technology to distinguish the signal and noise in GST domain. We discussed two types of noises: white noise and red noise, in which noise satisfy statistical autoregression model. For these two model, the noise-signal detection technology based on GST all get good result. It proved that the GST domain noise-signal detection technology could be used to real seismic data, and could effectively avoid noise influence on seismic sequence detecting.On the seismic profile after GST processing, high amplitude energy intensive zone, schollen, strip and lentoid dead zone and disarray zone maybe represent specifically geologic meanings according to given geologic background. Using seismic sequence detection profile and combining other seismic interpretation technologies, we can elaborate depict the shape of palaeo-geomorphology, effectively estimate sand stretch, distinguish sedimentary facies, determine target area, and directly guide oil-gas exploration.In the lateral reservoir prediction in XF oilfield of Ordos Basin, it played very important role in the estimation of sand stretch that the study of palaeo-geomorphology of Triassic System and the partition of inner sequence of the stratum group. According to the high-resolution seismic profile after GST processing, we pointed out that the C8 Member of Yanchang Formation in DZ area and C8 Member in BM area are the same deposit. It provided the foundation for getting 430 million tons predicting reserves and unite building 3 million tons off-take potential.In tackling key problem study for SLG gas-field, according to the high-resolution seismic sequence profile, we determined that the deposit direction of H8 member is approximately N-S or NNE-SS W. Using the seismic sequence profile, combining with layer-level profile, we can interpret the shape of entrenched stream. The sunken lenticle indicates the high-energy stream channel, which has stronger hydropower. By this way we drew out three high-energy stream channels' outline, and determined the target areas for exploitation. Finding high-energy braided river by high-resolution sequence processing is the key technology in SLG area.In ZZ area, we studied the distribution of the main reservoir bed-S23, which is shallow delta thin sand bed, by GST processing. From the seismic sequence profile, we discovered that the schollen thick sand beds are only local distributed, and most of them are distributary channel sand and distributary bar deposit. Then we determined that the S23 sand deposit direction is NW-SE in west, N-S in central and NE-SW in east. The high detecting seismic sequence interpretation profiles have been tested by 14 wells, 2 wells mismatch and the coincidence rate is 85.7%. Based on the profiles we suggested 3 predicted wells, one well (Yu54) completed and the other two is still drilling. The completed on Is coincident with the forecastThe paper testified that GST is a effective technology to get high- resolution seismic sequence profile, compartmentalize deposit microfacies, confirm strike direction of sandstone and make sure of the distribution range of oil-gas bearing sandstone, and is the gordian technique for the exploration of lithologic gas-oil pool in complicated areas.

克拉2异常高压气藏地质特征及开发技术研究

Kela-2 gas field in Tarim Basin is the main supply source for West-to-East Pipeline project, also the largest abnormally-pressured gas field discovered in China currently. The geological characterization, fine geological modeling and field development plan are all the world-class difficult problems. This work includes an integrated geological and gas reservoir engineering study using advanced technology and approaches, the scientific development plan of Kela-2 gas field as well as the optimizations of the drilling, production and surface schemes. Then, it's expected that the Kela-2 gas field can be developed high-efficiently.Kuche depression is one part of the thrust belt of the South Tianshan Mountains, Kela-2 field is located at the Kelasu structural zone in the north of Kuche depression. The field territory is heavily rugged with deeply cut gullies, complex geological underground structure, variable rock types, thrust structure development. Therefore, considerable efforts have been made to develop an integrated technique to acquire, process and interpret the seismic data in complicated mountain region. Consequently a set of seismic-related techniques in the complicated mountain region has been developed and successfully utilized to interpret the structure of Kela-2 gas field.The main reservoir depositional system of Kela 2 gas field is a platform - fan delta - braided river system. The reservoir rocks are medium-fine and extremely fine grained sandstones with high structure maturity and low composition maturity. The pore system structure is featured by medium-small pore, medium-fine throat and medium-low assortment. The reservoir of Kela-2 gas field is characteristic of medium porosity and medium permeability. The pay zone is very thick and its lateral distribution is stable with a good connection of sand body. The overpressure is caused mainly by the strongly tectonic squash activities, and other factors including the later rapid raise and compartment of the high-pressure fluid, the injection of high-pressure fluid into the reservoir.Based on the deliverability tests available, the average binomial deliverability equation is provided applicable for the overall field. The experimental results of rock stress-sensitive tests are employed to analyze the change trend of petrophysical properties against net confining stress, and establish the stress-based average deliverability equation. The results demonstrate the effect of rock deformation on the deliverability is limited to less than 5% in the early period of Kela-2 gas field, indicating the insignificant effect on deliverability of rock deformation.In terms of the well pattern comparisons and development planning optimizations, it is recommended that the producers should be located almost linearly along the structural axis. A total of 9 producers have a stable gas supply volume of 10.76 BCMPY for 17 years. For Kela-2 gas field the total construction investment is estimated at ￥7,697,690,000 RMB with the internal earning rate of 25.02% after taxation, the net present value of ￥7,420,160,000 RMB and the payback period of 5.66 years. The high profits of this field development project are much satisfactory.

龙滩水电站左岸边坡监测信息系统研制及其监测资料分析

Guangxi Longtan Hydropower Station is not only a representative project of West Developing and Power Transmission from West to East in China, but also the second Hydropower Station to Three Gorges Project which is under construction in China. There are 770 X 104m3 creeping rock mass on the left bank slope in upper reaches, in which laid 9 water inlet tunnels and some underground plant buildings. Since the 435m high excavated slope threatens the security of the Dam, its deformation and stability is of great importance to the power station.Based on the Autodesk Map2004, Longtan Hydropower Station Monitoring Information System on Left Bank has been basically finished on the whole. Integrating the hydropower station monitoring information into Geographic Information System(GIS) environment, managers and engineers can dynamically gain the deformation information of the slop by query the symbols. By this means, designers can improve the correctness of analysis, and make a strategic and proper decision. Since the system is beneficial to effectively manage the monitoring-data, equitably save the cost of design and safe construction, and decrease the workload of the engineers, it is a successful application to the combination of hydropower station monitoring information management and computer information system technology.At the same time, on the basis of the geological analysis and rock mass toppling deformation and failure mechanism analysis of Longtan engineering left bank slope, the synthetic space-time analysis and influence factors analysis on the surface monitoring data and deep rock mass monitoring data of A-zone on left bank slope are carried on. It shows that the main intrinsic factor that effects the deformation of Zone A is the argillite limestone interbedding toppling structure, and its main external factors are rain and slope excavation. What's more, Degree of Reinforcement Demand(DRD) has been used to evaluate the slop reinforce effect of Zone A on left bank according to the Engineering Geomechanics-mate-Synthetics(EGMS). The result shows that the slop has been effective reinforced, and it is more stable after reinforce.At last, on the basis of contrasting with several forecast models, a synthetic forecast GRAV model has been presented and used to forecast the deformation of zone A on left bank in generating electricity period. The result indicates that GRAV model has good forecast precision, strong stability, and practical valuable reliability.

中国地下工程氡污染健康影响与防护研究-以内蒙古自治区为例

由于对人体健康的危害，氡越来越受到世界各国的广泛关注。经济的发展，使我国地下公共场所大量涌现，但其空气中氡浓度的大规模调查及系统研究还属空白。本论文采取瞬时法（美国SUN公司的1027型连续测氰仪）与累积法（CR一39SSNTD）相结合的测量方法，对我国23个城市的地下工程进行了为期一年的监测，获得了实测数据700多个。在监测的同时，收集了每个测点的各种相关情况资料。论文以内蒙古自治区的呼和浩特、包头、集宁、鄂尔多斯四个城市为例，从环境地球化学角度出发，在研究了各个测点的各种情况之后，探讨了氧在不同类型工程和不同城市的浓度规律。对长期有人在内工作的地下工程内氛所致工作人员的有效剂量当量进行了估算。最后对氨异常的测点进行了原因分析，并提出了地下工程内氰污染的治理和防范措施。本次研究主要得出以下结论：1、地下工程中氰浓度的季节变化规律为：从总平均值看，夏高冬低。与前人对坑道的调查结果一致，与地上氨浓度的变化规律相反。从两季浓度差来看，两季浓度差从大到小依次为：隧道类工程）单建式铜室类工程＞结建式嗣室类工程。2、总体来看内蒙三类工程的氰浓度值：隧道类工程）结建式铜室类工程）单建式铜室类工程。不管是单季还是年平均值，各城市三类工程的氧浓度都遵循隧道类．工程远大于结建式铜室工程和单建式铜室工程的规律。通过对结建式桐室工程和单建式铜室工程的对比分析，发现在地下工程内工程上部有无建筑物阻隔，并不是影响氨浓度高低的重要因素，即氰的纵向扩散作用小于横向扩散作用。3、呼和浩特、包头、集宁、鄂尔多斯四个城市的氛浓度年平均值分别为：190Bq／m3、700Bq／m3、190Bq／m3、270Bq／m3。去掉仅有一季结果的侧点值，包头、鄂尔多斯的氨浓度年平均值为：62oBq／m3、250Bq／m3。呼和浩特、包头、集宁、鄂尔多斯四个城市的氛浓度单季值的范围分别为：31-830Bq／m3、57-3500Bq／m3、26一65oBq／m，、64一69OBq／m'，；年平均值的范围分别为：39一580Bq／m：'、94一210oBq／m：，、35-460Bq／m3、110-390Bq／m3。呼和浩特、包头、集宁、鄂尔多斯四个城市分别有30.0%、44.4%、33.3%、50.0%的工程年平均氧浓度在200Bq／时以上。四个城市氦浓度的普遍水平：包头）鄂尔多斯）集宁）呼和浩特。4、对所有单建式铜室类地下工程和结建式铜室类工程内工作人员所受剂量当量进行了估算。单建式地下工程内工作人员所受剂量当量的范围是1.247mSv-6．SOgmsv，平均值是2.610Sv。结建式地下工程内工作人员所受剂量当量的范围是1．195Sv-8.531mSv，平均值是2.572mSv。两类地下工程内工作人员所受氨照射剂量当量都很高，有一半的工程剂量当量是地面上公众人员剂量当量0.924mSv的两倍以上。5、不管是从浓度角度还是从剂量估算角度分析，地下工程内的氛水平都比较偏高，需要采取一定防护措施，对现有状况进行改善。6、通风是地下工程防治氡污染最简便易行、见效最快的措施。使用防氡涂料是最经济的措施。