井间地震技术研究和应用

Cross well seismic technique is a new type of geophysical method, which observes the seismic wave of the geologic body by placing both the source and receiver in the wells. By applying this method, it averted the absorption to high-frequency component of seismic signal caused by low weathering layers, thus, an extremely high-resolution seismic signal can be acquired. And extremely fine image of cross well formations, structure, and reservoir can be achieved as well. An integrated research is conducted to the high-frequency S-wave and P-wave data and some other data to determine the small faults, small structure and resolving the issues concerning the thin bed and reservoir's connectivity, fluid distribution, steam injection and fracture. This method connects the high-resolution surface seismic, logging and reservoir engineering. In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a research is conducted on cross well seismic technology in general and its important issues in cross well seismic technology in particular. A technological series of integrated field acquisition, data processing and interpretation and its integrated application research were developed and this new method can be applied to oilfield development and optimizing oilfield development scheme. The contents and results in this paper are as listed follows: An overview was given on the status quo and development of the cross well seismic method and problems concerning the cross well seismic technology and the difference in cross well seismic technology between China and international levels; And an analysis and comparison are given on foreign-made field data acquisition systems for cross-well seismic and pointed out the pros and cons of the field systems manufactured by these two foreign companies and this is highly valuable to import foreign-made cross well seismic field acquisition system for China. After analyses were conducted to the geometry design and field data for the cross well seismic method, a common wave field time-depth curve equation was derived and three types of pipe waves were discovered for the first time. Then, a research was conducted on the mechanism for its generation. Based on the wave field separation theory for cross well seismic method, we believe that different type of wave fields in different gather domain has different attributes characteristics, multiple methods (for instance, F-K filtering and median filtering) were applied in eliminating and suppressing the cross well disturbances and successfully separated the upgoing and downgoing waves and a satisfactory result has been achieved. In the area of wave field numerical simulation for cross well seismic method, a analysis was conducted on conventional ray tracing method and its shortcomings and proposed a minimum travel time ray tracing method based on Feraiat theory in this paper. This method is not only has high-speed calculation, but also with no rays enter into "dead end" or "blinded spot" after numerous iterations and it is become more adequate for complex velocity model. This is first time that the travel time interpolation has been brought into consideration, a dynamic ray tracing method with shortest possible path has been developed for the first arrivals of any complex mediums, such as transmission, diffraction and refraction, etc and eliminated the limitation for only traveling from one node to another node and increases the calculation accuracy for minimum travel time and ray tracing path and derives solution and corresponding edge conditions to the fourth-order differential sonic wave equation. The final step is to calculate cross well seismic synthetics for given source and receivers from multiple geological bodies. Thus, real cross-well seismic wave field can be recognized through scientific means and provides important foundation to guide the cross well seismic field geometry designing. A velocity tomographic inversion of the least square conjugated gradient method was developed for cross well seismic velocity tomopgraphic inversion and a modification has been made to object function of the old high frequency ray tracing method and put forward a thin bed oriented model for finite frequency velocity tomographic inversion method. As the theory model and results demonstrates that the method is simple and effective and is very important in seismic ray tomographic imaging for the complex geological body. Based on the characteristics of the cross well seismic algorithm, a processing flow for cross well seismic data processing has been built and optimized and applied to the production, a good section of velocity tomopgrphic inversion and cross well reflection imaging has been acquired. The cross well seismic data is acquired from the depth domain and how to interprets the depth domain data and retrieve the attributes is a brand new subject. After research was conducted on synthetics and trace integration from depth domain for the cross well seismic data interpretation, first of all, a research was conducted on logging constraint wave impedance of cross well seismic data and initially set up cross well seismic data interpretation flows. After it applied and interpreted to the cross well seismic data and a good geological results has been achieved in velocity tomographic inversion and reflection depth imaging and a lot of difficult problems for oilfield development has been resolved. This powerful, new method is good for oilfield development scheme optimization and increasing EOR. Based on conventional reservoir geological model building from logging data, a new method is also discussed on constraining the accuracy of reservoir geological model by applying the high resolution cross well seismic data and it has applied to Fan 124 project and a good results has been achieved which it presents a bight future for the cross well seismic technology.

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

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 seismic survey is the most effective prospecting geophysical method during exploration and development of oil/gas. The structure and the lithology of the geological body become increasingly complex now. So it must assure that the seismic section own upper resolution if we need accurately describe the targets. High signal/noise ratio is the precondition of high-resolution. As one important seismic data processing method, Stacking is an effective means to suppress the records noise. Broadening area of surface stacked is more important to enhance genuine reflection signals and suppressing unwanted energy in the form of coherent and random ambient noise. Common reflection surface stack is a macro-model independent seismic imaging method. Based on the similarity of CRP trace gathers in one coherent zone, CRS stack effectively improves S/N ratio by using more CMP trace gathers to stack. It is regarded as one important method of seismic data processing. Performing CRS stack depends on three attributes. However, the equation of CRS is invalid under condition of great offset. In this thesis, one method based on velocity model in depth domain is put forward. Ray tracing is used to determine the traveltime of CRP in one common reflection surface by the least squares method to regress the equation of CRS. Then we stack in the coherent seismic data set according to the traveltime, and get the zero offset section. In the end of flowchart of implementing CRS stack, one method using the dip angle to enhance the ratio of S/N is used. Application of the method on synthetic examples and field seismic records, the results of this method show an excellent performance of the algorithm both in accuracy and efficiency.

库车坳陷-南天山盆山系统第三纪脆性构造研究

Study on the structural coupling relationship between basin and range is not only helpful to recognize the basin formation and evolution systematically, but also to guide petroleum exploration in the basin. As a late Paleozoic Orogen, the South Tianshan Mountains reactivated and uplifted rapidly during the Cenozoic, and led to the Mesozoic-Cenozoic considerable thick deposits in the Kuqa Depression. The researches of the dissertation were carried out in the Kuqa depression-South Tianshan M ountain s ystem, a nd t he b rittle m icrotectonics w. ere c hosen as t he m ost important object. Based on observations and measurements of the field, we made detailed investigations on the geometry and kinematics of this area, and analyzed the abutting and cutting relationships and relative sequence of many brittle structures, such as joint, shear fractures, faults and some small-scale structures related to them closely. According to those brittle fractures' relationships with stress, the nature and variation of regional palaeostress field during the Cenozoic were studied through inversion of fault slip data and inferring stress state from joint sequences. And the deformation time was estimated primarily via ESR dating of faulting. Results show that the stress field varies as well in times as in space. The maximal principal stress direction shifted from the vertical to the horizontal, and stress regime from weak extension to strong compression from the Paleogene to the Neogene regionally. During the late Neogene, the structural deformation of the South Tianshan and the basin-range boundary was dominated by near N-S extension, while near N-S compressive deformation in the interior of the Kuqa Depression. There exits obvious differential stress state from the north to the south. ESR dating of the faulting during the Cenozoic indicates that, the normal faulting in the north edge of the Kuqa Depression have been active all along from the Miocene to the early Pleistocene, but the thrusting and reverse faulting in the interior only been active from the Pliocene to the early Pleistocene. On the base of those geological data and some geophysical information and theoretical calculation results, we infer that, the different stress regime the basin-range system is ascribed to the vertical uplift of the Tianshan Mountain. It was the vertical uplift that lead to the gravity-driven gliding of thick layers lying on the faulted basement from the South Tianshan Mountain to the Kuqa depression, and to folding and thrusting in the interior and frontal of the Kuqa depression. Combining the structural evolution with petroleum geological conditions of the Kuqa Depression, we think that the strong compressive deformation of the Kuqa Depression during rapid uplifting of the Tianshan Mountains from the Pliocene to the early Pleistocene play crucial role in the structural trap formation and proliferous gas accumulation.

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

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.

电阻率测井反演方法应用研究

Formation resistivity is one of the most important parameters to be evaluated in the evaluation of reservoir. In order to acquire the true value of virginal formation, various types of resistivity logging tools have been developed. However, with the increment of the proved reserves, the thickness of interest pay zone is becoming thinner and thinner, especially in the terrestrial deposit oilfield, so that electrical logging tools, limited by the contradictory requirements of resolution and investigation depth of this kinds of tools, can not provide the true value of the formation resistivity. Therefore, resitivity inversion techniques have been popular in the determination of true formation resistivity based on the improving logging data from new tools. In geophysical inverse problems, non-unique solution is inevitable due to the noisy data and deficient measurement information. I address this problem in my dissertation from three aspects, data acquisition, data processing/inversion and applications of the results/ uncertainty evaluation of the non-unique solution. Some other problems in the traditional inversion methods such as slowness speed of the convergence and the initial-correlation results. Firstly, I deal with the uncertainties in the data to be processed. The combination of micro-spherically focused log (MSFL) and dual laterolog(DLL) is the standard program to determine formation resistivity. During the inversion, the readings of MSFL are regarded as the resistivity of invasion zone of the formation after being corrected. However, the errors can be as large as 30 percent due to mud cake influence even if the rugose borehole effects on the readings of MSFL can be ignored. Furthermore, there still are argues about whether the two logs can be quantitatively used to determine formation resisitivities due to the different measurement principles. Thus, anew type of laterolog tool is designed theoretically. The new tool can provide three curves with different investigation depths and the nearly same resolution. The resolution is about 0.4meter. Secondly, because the popular iterative inversion method based on the least-square estimation can not solve problems more than two parameters simultaneously and the new laterolog logging tool is not applied to practice, my work is focused on two parameters inversion (radius of the invasion and the resistivty of virgin information ) of traditional dual laterolog logging data. An unequal weighted damp factors- revised method is developed to instead of the parameter-revised techniques used in the traditional inversion method. In this new method, the parameter is revised not only dependency on the damp its self but also dependency on the difference between the measurement data and the fitting data in different layers. At least 2 iterative numbers are reduced than the older method, the computation cost of inversion is reduced. The damp least-squares inversion method is the realization of Tikhonov's tradeoff theory on the smooth solution and stability of inversion process. This method is realized through linearity of non-linear inversion problem which must lead to the dependency of solution on the initial value of parameters. Thus, severe debates on efficiency of this kinds of methods are getting popular with the developments of non-linear processing methods. The artificial neural net method is proposed in this dissertation. The database of tool's response to formation parameters is built through the modeling of the laterolog tool and then is used to training the neural nets. A unit model is put forward to simplify the dada space and an additional physical limitation is applied to optimize the net after the cross-validation method is done. Results show that the neural net inversion method could replace the traditional inversion method in a single formation and can be used a method to determine the initial value of the traditional method. No matter what method is developed, the non-uniqueness and uncertainties of the solution could be inevitable. Thus, it is wise to evaluate the non-uniqueness and uncertainties of the solution in the application of inversion results. Bayes theorem provides a way to solve such problems. This method is illustrately discussed in a single formation and achieve plausible results. In the end, the traditional least squares inversion method is used to process raw logging data, the calculated oil saturation increased 20 percent than that not be proceed compared to core analysis.

基于MRF模型的油气储层属性随机建模方法研究

Stochastic reservoir modeling is a technique used in reservoir describing. Through this technique, multiple data sources with different scales can be integrated into the reservoir model and its uncertainty can be conveyed to researchers and supervisors. Stochastic reservoir modeling, for its digital models, its changeable scales, its honoring known information and data and its conveying uncertainty in models, provides a mathematical framework or platform for researchers to integrate multiple data sources and information with different scales into their prediction models. As a fresher method, stochastic reservoir modeling is on the upswing. Based on related works, this paper, starting with Markov property in reservoir, illustrates how to constitute spatial models for catalogued variables and continuum variables by use of Markov random fields. In order to explore reservoir properties, researchers should study the properties of rocks embedded in reservoirs. Apart from methods used in laboratories, geophysical means and subsequent interpretations may be the main sources for information and data used in petroleum exploration and exploitation. How to build a model for flow simulations based on incomplete information is to predict the spatial distributions of different reservoir variables. Considering data source, digital extent and methods, reservoir modeling can be catalogued into four sorts: reservoir sedimentology based method, reservoir seismic prediction, kriging and stochastic reservoir modeling. The application of Markov chain models in the analogue of sedimentary strata is introduced in the third of the paper. The concept of Markov chain model, N-step transition probability matrix, stationary distribution, the estimation of transition probability matrix, the testing of Markov property, 2 means for organizing sections-method based on equal intervals and based on rock facies, embedded Markov matrix, semi-Markov chain model, hidden Markov chain model, etc, are presented in this part. Based on 1-D Markov chain model, conditional 1-D Markov chain model is discussed in the fourth part. By extending 1-D Markov chain model to 2-D, 3-D situations, conditional 2-D, 3-D Markov chain models are presented. This part also discusses the estimation of vertical transition probability, lateral transition probability and the initialization of the top boundary. Corresponding digital models are used to specify, or testify related discussions. The fifth part, based on the fourth part and the application of MRF in image analysis, discusses MRF based method to simulate the spatial distribution of catalogued reservoir variables. In the part, the probability of a special catalogued variable mass, the definition of energy function for catalogued variable mass as a Markov random field, Strauss model, estimation of components in energy function are presented. Corresponding digital models are used to specify, or testify, related discussions. As for the simulation of the spatial distribution of continuum reservoir variables, the sixth part mainly explores 2 methods. The first is pure GMRF based method. Related contents include GMRF model and its neighborhood, parameters estimation, and MCMC iteration method. A digital example illustrates the corresponding method. The second is two-stage models method. Based on the results of catalogued variables distribution simulation, this method, taking GMRF as the prior distribution for continuum variables, taking the relationship between catalogued variables such as rock facies, continuum variables such as porosity, permeability, fluid saturation, can bring a series of stochastic images for the spatial distribution of continuum variables. Integrating multiple data sources into the reservoir model is one of the merits of stochastic reservoir modeling. After discussing how to model spatial distributions of catalogued reservoir variables, continuum reservoir variables, the paper explores how to combine conceptual depositional models, well logs, cores, seismic attributes production history.

地球化学新方法在隐状矿床找矿预测中的应用研究

With the progress of prospecting, the need for the discovery of blind ore deposits become more and more urgent. To study and find out the method and technology for the discovery of blind and buried ores is now a priority task. New geochemical methods are key technology to discover blind ores. Information of mobile components related to blind ores were extracted using this new methods. These methods were tested and applied based on element' s mobile components migrating and enriched in geophysical-geochemical process. Several kinds of partial extraction techniques have tested based on element' s occurrence in hypergenic zone. Middle-large scale geochemical methods for exploration in forest and swamp have been tested. A serious of methods were tested and applied effetely about evaluation of regional geochemical anomaly, 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system instead of the normal net. 1. Element related with ores can be mobiled to migrate upwards and be absorpted by surface soil. These abnomal components can be concentrated by natural or artificial methods. These trace metalic ions partially exist in dissovlvable ion forms of active state, and partially have been absorbed by Fe-Mn oxide, soil and organic matter in the soil so that a series of reaction such as complex reaction have take place. Employing various partial extraction techniques, metallic ions related with the phase of the blind ores can be extracted, such as the technique of organic complex extraction, Fe-Mn oxide extraction and the extraction technique of metallic ions of various absorption phases. 2.1:200000 regional geochemical evaluation anomaly methods: Advantageous ore-forming areas were selected firstly. Center, concentration, morphological feature, belt of anomaly were choosed then. Geological and geochemical anomalies were combined. And geological and geochemical background information were restrained. Xilekuduke area in Fuyun sheet , Zhaheba area in Qiakuerte sheet, the west-north part in Ertai sheet and Hongshanzui anomaly in Daqiao sheet were selected as target areas, in Alertai, in the north of Xinjiang. in Xilekuduke area, 1:25000 soil geochemical methods sampling based on the net in dendritic water system was carried out. Cu anomaly and copper mineralization were determined in the center area. Au , Cu anomalies and high polarization anomaly were determined in the south part. Prospecting by primary halo and organic complex extraction were used to prognosis blind ore in widely rang outcrop of bedrock. 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system were used in transported overburden outside of mining area. Shallow seismic method and primary halo found a new blind orebody in mining area. A mineralization site was fou and outside of Puziwan gold mine, in the north of Shanxi province. Developing middle-large scale geochemical exploration method is a key technique based 1:200000 regional geochemical exploration. Some conditions were tested as Sampling density , distribution sites of sample, grain size of sample and occurrence of element for exploration. 1:50000 exploration method was advanced to sample clast sediment supplement clast sediment in valley. 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system was applied to sample residual material in A or C horizon. 1:2000 primary or soil halo methods used to check anomalies and determine mineralization. Daliang gold mineralization in the northern Moerdaoga was found appling these methods. Thermomagnetic method was tested in miniqi copper-polymetallic ore. Process methods such as grain size of sample, heated temperature, magnetic separating technique were tested. A suite of Thermomagnetic geochemical method was formed. This method was applied in Xiangshan Cu~Ni deposit which is cover by clast or Gobi in the eastern Xinjiang. Element's content and contrast of anomaly with Thermomagnetic geochemical method were higher than soil anomaly. Susceptibility after samples were heated could be as a assessment conference for anomaly. In some sectors thermo-magnetic Cu, Ni, Ti anomalious were found outside deposits area. There were strong anomal ies response up ore tested by several kind of partial extraction methods include Thermomagnetic, enzyme leach and other partial extractions in Kalatongke Cu-Ni deposit in hungriness area in the northern of Xinjiang. Element's anomalies of meobile were mainly in Fe-Mn oxide and salt. A Copper mineralization site in Xilekuduke anomaly area had been determined. A blind ore was foung by shallow seismic and geochemical method and a mineralization site was found outside this mining area in Puziwan gold deposit in shanxi province. A Gold mineralization site was found by 1:50000 geochemical exploration in Daliang, Inner Mongolia.

油气藏自然电位异常与油气赋存的相关性研究

The Stack Spontaneous Potential (SSP) is a direct hydrocarbon location technology and a new hydrocarbon detection method with independent intellectual property. A subsurface hydrocarbon accumulation associated with the upward hydrocarbon micro seepage induces a relatively strong negative potential abnormal zone, of which the anomaly can be measured on the surface with specially designed instruments through careful field measuring procedures. With special software programmed according to a unique geochemical and geophysical model, the original data are analyzed, processed and interpreted on the computer, and then on a series of resulting anomaly distribution maps and/or profiles, the favorable surface locations of the hydrocarbon accumulations can be easily identified. The study of the SSP has been started since 1989, and especially from 1996 to 1997, both profile and area tests were conducted in the Daqing Oilfield. On the testing line of 15kms, there are 6 wells in total, among which some are oil-producing wells, and some are water-producing wells. The final matching ratio of the favorable oil well locations and the possible water well locations predicted by the SSP to those of known wells was up to 83 percent. In the area test, of which the acreage is 800 km2, the matching ratio compared with the existing wells was 87 percent; furthermore, regarding to wells subsequently drilled after the test, the matching ratio was 85 percent. The matching ratio in the development area is more about 10 percent than those of in exploration area. The reason is that, comparing the exploration area, the development area acreage is less and the container rocks are more simplex. In development area there is not so much interference of SSP also. Since 1997 the SSP has been tested and applied all over China to a number of hydrocarbon bearing basins and known oil fields, including the Daqing, Jiangsu, Changqing, Shengli, Nanyang, Jianghan and Zhongyuan Oilfields, only to name a few. The SSP surveys in total areas of over 10,000km2 in more than 30 regions in China so far have been completed in various exploration and development stages, the satisfactory outcomes of which have further evidenced that the dependence between the SP anomaly and abundance of hydrocarbon. Up to date, a substantial amount of successful tests and actual surveys finished in exploration and development practices have evidenced that the SSP is significantly more reliable in comparison with any other similar direct hydrocarbon indication technique generally known to the oil industry, such as the Redox. The SSP can be applied to search for almost all kinds of hydrocarbon accumulations, regardless of the type of traps, such as structural, stratigraphic, buried hill traps, and so on; however, it is interesting to be noted that the SSP seems to be particularly effective in detecting the stratigraphic oil traps according to our practices. On the other hand, there is virtually no surface geographical constrains in terms of field data acquisition, except for those water covered areas, because of the inherent characteristics of the technology itself. Furthermore, utilizing the SSP requires no special considerations to subsurface geological conditions in regard to formation resistivity, since the SSP measurements will not be influenced by either overly high or overly low resistivity of formations lying above the hydrocarbon accumulations. There are two kind of theory, of which, as we know one is called hypbyssal theory such as "Redox"[61 the other is call plutonic theory such as cracking of hydrocarbon [8][9] and natural polarization [3], to describe the mechanism of SP anomaly of oil reservoir and to indicate that the dependence between the SP abnormality and abundance of hydrocarbon has be existed theoretically/The quantitative dependence, which has not been founded due to the complicity of container rocks, be discovered during the exploration and development practices is the crux to the quantitative analysis of SP Anomaly processing. Based on the thorough study of the complex of collector rocks, every kind of thickness of collector rock can be conversed to be a standard effective thickness; the thickness is called apparent effective thickness (AET). The conversation coefficient (ai, 1=1,2,3) could be determined by the variety of every collector rock storability (CRS). The discoveration of quantitative: dependence between AET and the amplitude of SSP, in the practices of exploration and development, is a promotion for the SSP supplied in the oil exploration, and make the data analysis forward to the quantitative stage.

环渤海地区地壳和上地幔地震层析成像研究

In Circum-Bohai region (112°~124°E, 34°~42°N ), there exists rich gas-petroleum while inner-plate seismic activity is robust. Although the tectonic structure of this region is very complicated, plenty of geological, geophysical and geochemical researches have been carried out.In this paper, guided by the ideas of "One, Two, Three and Many" and "The depth controls the shallow, the regional constrains the local", I fully take advantage of previous results so as to establish a general image of this region. After collecting the arrival-time of P-wave phases of local events and tele-seismic events recorded by the stations within this region from 1966 to 2004, I process all these data and build an initial model. Then, a tomography image of crust and upper-mantle of this region is obtained. With reference to previous results, we compare the image of various depths and five cross-profiles traverse this region along different direction. And finally, a discussion and conclusion is made.The principle contents is listed as below: 1) in the first chapter, the purpose and meaning of this thesis, the advance in seismic tomography, and the research contents and blue-print is stated; 2) in the second chapter, I introduce the regional geological setting of Circum-Bohai region, describe the tectonic and evolutionary characteristics of principle tectonic units, including Bohai Bay Basin, Yanshan Fold Zone, Taihangshan Uplifted Zone, Jiao-Niao Uplifted Zone and Luxi Uplifted Zone, and primary deep faults; 3) In the third chapter, the previous geophysical researches, i.e., gravity and geomagnetic characters, geothermal flow, seismic activity, physical character of rocks, deep seismic sounding, and previous seismic tomography, are discussed; 4) in the fourth chapter, the fundamental theory and approach of seismic tomography is introduced; 5) in the fifth chapter, the technology and approaches used in this thesis, including collecting and pre-processing of data, the establishment of initial velocity model and relocation of all events; 6) in the sixth chapter, I discuss and analyze the tomography image of various depth and five cross-sections; 7)in the seventh chapter, I make a conclusion of the results, state the existing problems and possible solutions.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

济阳地区深部结构与含油气性讨论

Oil and gas exploration of marine strata in China's Pre-Cenozoic residual basins is regarded as a worldwide puzzle because of existent problems and cruxes. Objectively speaking, the subsurface geologic structure is complicated, and the surface conditions of some areas are tough. On the other hand, there are still many problems to be solved in oil and gas exploration technologies of Pre-Cenozoic marine fades, and theoretic cognition about petroleum geology is not profound yet. Therefore, it is principal to explore integrated geophysical research ways of Pre-Cenozoic residual basins. Seismic prospecting and geophysical integrated interpretation technologies aimed at middle Paleozoic marine facies with deeper burial and complicated geologic conditions have not formed due to bad quality of deep strata data. Pre-Cenozoic strata, and especially extension, thickness and internal structure of Paleozoic strata can not be recognized from seismic profiles, so it is hard to systematically cognize structural features and oil-gas resources prospect of Pre-Cenozoic basins. To further investigation of fabric and structural features, basin prototype, formation and evolution pattern of Pre-Cenozoic basins, and also their control over formation, migration and aggregation of oil and gas, will play a guiding and promotive role in developing new surveying areas, selecting advantageous zones and predicting oil-gas resources.This paper follows the modem macrocontrol theory of "Region constrains local, deep strata controls shallow ones", and uses the integrated geophysical method of "One guide, two hinges, three combinations and multi feedbacks'*. Based on several years of geological and geophysical results of the Shengli Oilfield, and 14 newly-joint regional seismic profiles, deep structure and oil-gas bearing capacity of the Jiyang area are discussed and new cognitions are drawn as below.Seismic identification marks Tr, Tg, Tgl and Tg2 are established for importantPre-Cenozoic geological interfaces, and promoted to the whole Jiyang area.Through area-wide tracking and clogging of important seismic reflection marker,the isochronic framework of pre-Tertiary basin is set up in the Jiyang area for the firsttime, which is vital for basin research.Integrated with geological and geophysical research results, the Jiyang area isdivided into four first-order tectonic sequences- basement, lower tectonic layer,upper tectonic layer, and top tectonic layer. The basement and lower tectonic sequence which are related to Pre-Cenozoic are studied with emphasis.Through the research of regional seismic profiles, the point of view is given thatthe Kongdian Formation of Jiyang is structural transition period. The top-bottomunconformable interface of the Kongdian Formation is found out for the first time,and the basin model is determined primarily, which lay a basis for prototype basinresearch of the Jiyang Kongdian Formation.The distribution status of Middle-Paleozoic is delineated in the Jiyang area.The maximum thickness of Paleozoic lies in the top of the south declivity of half-graben. The thickness gets thinner towards the center of Mesozoic and Cenozoic half-graben basin, and even disappears. Structural action in the west-north affects the distribution of Paleozoic residual strata.6. The features of second-order tectonic sequence of the Jiyang depression isstudied and its evolution history of is rebuilt.Combined with the 5-stage evolution history of the China continent and structure evolution features of the Jiyang area, the structure sedimentary process since Paleozoic is divided into 5 periods - basement forming , Indosinian orogenic, Yanshan negative reversal, Himalayan extension and Neogene subsidence period.Combined with the research results of gravity, magnetic surveying and regionalprofiles, this paper brings forward the idea for the first time that the western boundaryof the Jiyang depression is the Ningjin-Yangpan fracture zone, and forms aside-column assemblage with the Wudi fracture zone.The opinion that under Middle-Cenozoic basins in the middle Jiyang area theremight superimpose an old residual basin is given for the first time. And if it is provedto be true, a new exploration space will be pioneered for Jiyang and even north China.There exists many types of tectonic-stratigraphic traps formed under piezotropy,extension and compound action in Pre-Cenozoic Jiyang. Therein all kinds of burialhills are the most important oil-gas trap type of Pre-Cenozoic, which should besurveyed layeredly according to the layout of oil sources.As such a new challenging project and field, the paper systematically analyses different geophysical responses of the Jiyang area, frames the deep structure of the area, and preliminarily recognizes the Pre-Cenozoic residual basins. It breaks through to a certain extent in both theory and practice, and is expected to provide new geophysical and geotectonic clues for deep exploration in Shengli.