基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

有源电磁扩散场克希霍夫积分偏移

As active electromagnetic method, field data of CSAMT method follow the equation of diffusion. Propagting in solid earth media, diffusion EM signal has strong attenuation and dispersion, otherwise seismic wave shows weak attenuation and dispersion, therefore the resolution power of CSAMT method is not better than seismic reflection method. However, there is consistence and similarity between EM signal and seismic wave in wave equation, we can apply Kirchhoff integral migration technique, a proven one in seismic method in time domain, to carry out seduo-seismic processing for CSAMT signal in frequency domain so that the attenuation and dispersion could be made compensated in some extent, and the resolution power and interpretation precision of active EM wave could be improved. Satisfying passive homogeneous Helmholtz quation, we proceed with Green theorem and combine the active inhomogenous Helmholtz quation, the Kirchhoff integral formula could be derived. Given practical problems, if we only consider the surface integral value, and assume that the intergral value in other interface is zero, combined with Green theorem in uniform half space, the expression could be simplified, and we can obtain frequency-domain Kirchhoff integral formula in surface, which is also called downward continuation of EM field in frequency domain. With image conditions and energy compensation considered, in order to get image conditions in time domain Fourier inverse transformation in frequency domain can be performed, so we can formulate the active Kirchhoff integral migration expression. At first, we construct relative stratified model, with different frequency series taken into account, then we change the distances between transmitter and reciever, the EM response can be obtained. Analyzing the EM properties, we can clarify near and far zone that can instruct us to carry out transmitter layout in practical application. Combined with field data surveyed in far zone, We perform Kirchhoff integral migration and compare the results with model to interpret. Secondly, with far field EM data, we apply TM mode to get EM response of given 2D model, then apply Kirchhoff integral migration on modelling data and interpret the results.

非自喷层测试及评价解释技术研究与应用

Now low porosity and low permeability reservoir is one of the main targets of exploration for the onshore oilfields of China. Most of the reservoirs are none flowing because of bad formation percolation condition, poor gas oil ratio , low formation pressure coefficient and other factors. In the recent years, a number of domestic oilfields have carried out some research work and achieved some success on oil testing and production technology in such formation. But by now, there is still no systematic and mature technology, particularly testing technology in none flowing formation is still needed further study. Based on study the key problem of well testing and interpretation technology in none flowing formation, solve the important problems in well testing technology, continuously improve and innovate geological information acquisition technology for none flowing reservoir, accurately acquire boundary information and evaluate reservoir flow characteristics. Its wide application remarkable result has shown. The main results and cognitions obtained from research are as follows: 1. This new technology research results help solve the occurrent problems in well testing process for none flowing formations, such as small investigation radius, poor representative of interpretation results from the poor data, low level application of interpretation results. This new technology helps create favorable conditions for early precise reservoir evaluation and reduction of the risk of exploration. 2. The technological difficulties for none flowing well testing are successfully solved by using none flowing formation combined mechanical tool string .This method has been proved by its applications to be able to improve the efficiency of the testing and the quantity of the acquired test data ,and so as to enhance the application of the interpretation results of the test in development of oil fields. 3. The application of the rotary formation tester, selective test valve, well testing string and their allier tools help to resolve problems such as the operation of opening and shutting-in the well under different well conditions, to broaden the scope of well test technology for none flowing formations. 4. Refined Testing Technique for production Wells has greatly shortened the testing dwration and improved the efficiency and accuracy of operation, enriched test results, and at the same time created conditions for conducting multi-well interference well testing.

甘肃临夏盆地龙担剖面晚新生代磁性地层与古环境

Linxia Basin, situated in the northeast belt of the Tibetan Plateau, is a late Cenozoic depression basin bounded by the Tibetan Plateau and the Chinese Loess Plateau. The Cenozoic deposition, spanning over 30Ma, in which very abundant mammal fossils were discovered, is very suitable for study of uplift processes and geo-morphological evolution of the Tibetan Plateau. The Longdan section (35°31′31.6″N，103°29′0.6″E) is famous for the middle Miocene Platybelodon fauna and the late Miocene Hipparion fauna for a long time and is also one of the earliest known places for wooly rhino, which lies on the east slope of Longdan, a small village of township Nalesi in the south of the Dongxiang Autonomous County, Linxia Hui Nationallity Autonomous Prefecture. The Longdan mammal fauna was discovered at the base of the Early Pleistocene loess deposits at Dongxiang, where the lithology is different from the typical Wucheng Loess on the Chinese Loess Plateau. The rich fossils contain many new species and the major two layers of fossils are in the loess beds. Geologically the fossiliferous area is located in the central part of the Linxia Cenozoic sedimentary basin. Tectonically the Linxia Basin is an intermountain fault basin, bordered by the Leijishan major fault in the south and the north Qinling and Qilianshan major faults in the north. The section is 51.6m thick above the gravel layer, including the 1.6m Late Pleistocene Malan Loess on the top and the other loess-paleosol sequences in the middle of the section. The base of the section is the Jishi Formation, consisting of gravel layer of 13 ~ 17m thick. In this study, 972 bulk samples were collected with an interval of 5cm and other 401 orientied samples were taken with a magnetic compass. In the laboratory, the paleomagnetism, medium grain size, susceptibility, color, micromorphology, anisotropy of magnetic susceptibility were analyzed. From the stratigraphic analysis, the Longdan section from the top 0.3m to the bottom 51.6m, containing 5 normal polarities (N1-N5) and 5 reversal polarities (R1-R5). The paleomagnetic results show N3 is the Olduvai subchron in the middle of the Matuyama chron, and then the chronology of the Longdan mammal fauna is constructed along the section. The Matuyama-Gauss boundary is 45m and N5 enters Gauss chron. The Olduvai subchron with the age of 1.77 ~ 1.95Ma is found just in the upper fossiliferous level of Longdan mammal fauna. Taking the deposit rate of the section into account, the geological age of the upper fossiliferous level of Longdan mammal fauna is estimated to be about 1.9Ma. The lower fossiliferous level is just below the Reunion subchron and its age is estimated to be 2.25Ma. In addition, anisotropy of magnetic susceptibility of the loess-paleosol and other climatic indexes were used for discussing the late Cenozoic paleoenvironmental changes at Longdan, from which the Longdan area should have been an area of predominantly steppe the same as the Longdan mammal fauna.

青藏高原东部电性结构特征及其构造意义

Eastern Himalayan Syntaxis (EHS) and its surroundings (eastern margin of Tibet) is one of the most complicated tectonic areas in the world. As the exhaust opening of the balanced materials of the Tibetan Plateau during the collision of Indan and Eurasian plates, the deep structure beneath EHS surrounding region is referred to as the key to the study of the dynamics of the plateau. EHS3D project, sponsored by NSFC, has been proposed to explore the deep electric features of the area. During the first stage of EHS3D(2006-2008), MT+LMT measurements have been conducted along two lines from Chayu to Qingshuihe (EHS3D-3) and Chayu to Ruoergai (EHS3D-2). This paper will discuss the MT models of EHS3D-3 line. By the data procrssing, including distortion analysis, Robust estimation and strike decomposition, rotated apparent resitivities and phases have been obtained for each station. Then conventional 2-D inversion algorithms (NLCG and RRI) were employed to produce 2-D models. The final preferred 2-D model suggests that the upper crust consists of resistive blocks while in mid-lower crust there are two extensive conductive bodies beneath Lhasa block and Qiangtang terrain respectively. Jinshajiang suture is a gradient belt and Bangong-Nujiang suture appear a conductive belt dipping to the north. . We concluded that the formation of the two conductive bodies attributed to the partial melt and fluids in the lower crust. The regional electric strike derived from decomposition analysis indicates that the crust and upper mantle move in different manners. The upper crust moves like slips of rigid blocks along major slip faults while the lower crust creeps as a flow in the conductive channels.

复杂地质条件下高精度地震勘探技术研究与应用-以泌阳凹陷为例

In order to carry out high-precision three-dimensional "integration" for the characteristics of the secondary seismic exploration for Biyang Depression, in the implementation process, through a combination of scientific research and production, summed up high-precision seismic acquisition, processing and interpretation technologies suitable for the eastern part of the old liberated areas, achieved the following results: 1. high-precision complex three-dimensional seismic exploration technology series suitable for shallow depression Biyang block group. To highlight the shallow seismic signal, apply goal-based observing system design, trail from the small panel to receive and protect the shallow treatment of a range of technologies; to explain the use of three-dimensional visualization and coherent combination of full-body three-dimensional fine interpretation identification of the 50-100 m below the unconformity surface and its formation of about 10 meters of the distribution of small faults and improve the small block and stratigraphic unconformity traps recognition. 2. high-precision series of three-dimensional seismic exploration technology suitable for deep depression Biyang low signal to noise ratio of information. Binding model using forward and lighting technology, wide-angle observation system covering the design, multiple suppression and raise the energy of deep seismic reflection processing and interpretation of detailed, comprehensive reservoir description, such as research and technology, identified a number of different types of traps. 3. high-precision seismic exploration technology series for the southern Biyang Depression high steep three-dimensional structure. The use of new technology of seismic wave scattering theory and high-precision velocity model based on pre-stack time migration and depth migration imaging of seismic data and other high-precision processing technology, in order to identify the southern steep slope of the local structure prediction and analysis of sandstone bedrock surface patterns provide a wealth of information.

旋转轴对称介质中的qP波一阶射线追踪

Elastic anisotropy is a very common phenomenon in the Earth’s interior, especial for sedimentary rock as important gas and oil reservoirs. But in the processing and interpretation of seismic data, it is assumption that the media in the Earth’s interior is completely elastic and isotropic, and then the methods based on isotropy are used to deal with anisotropic seismic data, so it makes the seismic resolution lower and the error on images is caused. The research on seismic wave simulation technology can improve our understanding on the rules of seismic wave propagation in anisotropic media, and it can help us to resolve problems caused by anisotropy of media in the processing and interpretation of seismic data. So researching on weakly anisotropic media with rotated axis of symmetry, we study systematically the rules of seismic wave propagation in this kind of media, simulate the process with numerical calculation, and get the better research results. The first-order ray tracing (FORT) formulas of qP wave derived can adapt to every anisotropic media with arbitrary symmetry. The equations are considerably simpler than the exact ray tracing equations. The equations allow qP waves to be treated independently from qS waves, just as in isotropic media. They simplify considerably in media with higher symmetry anisotropy. In isotropic media, they reduce to the exact ray tracing equations. In contrast to other perturbation techniques used to trace rays in weakly anisotropic media, our approach does not require calculation of reference rays in a reference isotropic medium. The FORT-method rays are obtained directly. They are computationally more effective than standard ray tracing equations. Moreover the second-order travel time corrections formula derived can be used to reduce effectively the travel time error, and improve the accuracy of travel time calculation. The tensor transformation equations of weak-anisotropy parameters in media with rotated axis of symmetry derived from the Bond transformation equations resolve effectively the problems of coordinate transformation caused by the difference between global system of coordinate and local system of coordinate. The calculated weak-anisotropy parameters are completely suitable to the first-order ray tracing used in this paper, and their forms are simpler than those from the Bond transformation. In the numerical simulation on ray tracing, we use the travel time table calculation method that the locations of the grids in the ray beam are determined, then the travel times of the grids are obtained by the reversed distance interpolation. We get better calculation efficiency and accuracy by this method. Finally we verify the validity and adaptability of this method used in this paper with numerical simulations for the rotated TI model with anisotropy of about 8% and the rotated ORTHO model with anisotropy of about 20%. The results indicate that this method has better accuracy for both media with different types and different anisotropic strength. Keywords: weak-anisotropy, numerical simulation, ray tracing equation, travel time, inhomogeneity

青藏高原壳幔形变与层圈耦合数值模拟

Based on multi-principle (such as structures, tectonics and kinematics) exploratory data and related results of continental dynamics in the Tibetan plateau, the author reconstructed the geological-geophysical model of lithospherical structure and tectonic deformation, and the kinetics boundary conditions for the model. Then, the author used the numerical scheme of Fast Lagrangian Analysis of Continua (FLAC), to stimulate the possible process of the stress field and deformational field in the Tibetan plateau and its adjacent area, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. With the above-mentioned results, the author discussed the relationship between crustal movement in shallow layer and the deformational process in interior layers, and its possible dynamic constraints in deep. At the end of the paper, an integrative model has been put forward to explain the outline images of crust-mantle deformation and coupling in the Tibetan Plateau. (1) The characteristics of crust-mantle structure of the Tibetan plateau have been shown to be very complex, and vertical and horizontal difference is significant. The general characteristics of crust-mantle of the Tibetan plateau may be that it's layering in depth direction, and shows blocking from south to north and belting from east to west, mainly according to the results of about 20 seismic sections, such as wide-angle seismic profiles, CMP, seismic tomography and so on. (2) The crust had shortened about 2200km, while the shortening is different for different block from south to north in the Tibetan plateau. It is about 11.5mm/a in Himalayan block, about 9.0mm/a in Lhas-Gangdese block, about 7.0mm/a in Qiangtang block and Songpan-Ganzi-Kekexili block, about 8.0mm/a in Kunlun-Qaidam, and about ll.Omm/a in Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. Which - in demonstrates the shortening rate decreases from south to north, but this rate increases near the north edge of the Tibetan plateau. The crust thickening rate is about 0.4mm/a in the whole Tibetan plateau; and this rate is about 0.5mm/a in Himalayan block, about 0.4mm/a in Lhas-Gangdese block, about 0.3mm/a in Qiangtang block, about 0.2mm/a in Songpan-Ganzi-Kekexili block and about O.lmm/a in Kunlun-Qaidam-Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. This implies that the thickening rate decreases in the blocks of the Tibetan plateau. From south to north, the displacement of eastern boundary in the Tibetan plateau is about 37mm/a in Himalayan block, about 45mm/a in Lhas-Gangdese block, about 47mm/a in Qiangtang block, about 43mm/a in Songpan-Ganzi-Kekexili block, and about 35mm/a in Kunlun-Qaidam-Qilian block, since the collision-matching between the Indian continent and Eurasia continent had happened about 50Ma ago. This implies that the rate of eastward displacement is biggest in the middle of plateau, and decreases to both sides. The transition of S-N compression stress field in Tibetan Plateau, since about 28Ma+ ago, may be caused by two reasons: On one hand, the movement direction of Eurasia continent changed from northward to southward about 28Ma± ago in the northern plateau. On the other hand, the front belt that is located between India continent's and Eurasia continent's convergence-collision, had moved southward to high Himalayan from Indus-Brahmaputra suture almost at the same time in southern plateau. Affected by the stress field, the earlier tectonics rotated clockwise, NE and NW conjugate strike-slip faults developed, and the SN rift formed. This indicated that the EW movement started. The ratio between upper crust and lower crust of different blocks from south to north in the Tibetan plateau during the process of deformation are as following: about 3.5~5:1 in Himalayan block, about 1~5: 3-4 (which is about 1:3o--4 in south and about 4~5:3 in north) in Lhas-Gangdese block, about 1:3~447mm/a in these blocks: Which is located to the north of Banggong-nujiang suture.

黄土高原粉尘堆积源区环境及风化历史初步研究

In this study, we examined the surface features of quartz grains, the quartz oxygen isotopic ratios and the mineralogical compositions of the loess - paleosol - red clay sediments systematically. The surface features of quartz grains do not show significant changes of the dust deposits through the past seven million years. The particles were mainly created in the process of glacial and frost weathering of high mountains. Then the surfaces were altered in some degree by the flood and wind abrasion. The surface features registered all these processes. The assemblages of surface features changed for four times in the past seven million years, the occurrence ages are: 5.0~4.2MaBP, about 3.6MaBP, about 2.6MaBP and about 0.9MaBP, respectively. This may indicate that there were uplift events of the Tibetan Plateau during those times. The oxygen isotopic compositions of quartz in the sediments represent the oxygen isotopic compositions of the initial dusts because of the stable properties of quartz both physically and chemically. The oxygen isotopic compositions of 4~16um quartz changed significantly at about 2.6MaBP, decreasing from about 19.5%o to about 18.5%o. This decrease of quartz oxygen isotopic ratio suggests that the environments of the dust source areas changed at that time, or the range of dust source area changed at that time. The environmental change may result from the structural evolution of the Tibetan Plateau and global cooling at that time. The coarse fractions (>30μm) of the dust deposits were examined using the EDXA device for mineral identification. The quartz content has a decrease trend during 7~2MaBP, then increase rapidly at about 2MaBP. After 2MaBP, quartz content continues to decrease. The Ca-plagioclase content / quartz content ratio increase at about 3.6MaBP. The ratio shows a peak of 3-6 fold values at about 2.5~1.8MaBP, the cause of this is still unknown. The Ca-plagioclase content / quartz content ratio continues to increase after 1 MaBP. The flowing can be regarded as the conclusion remarks of this study: Some of the red clay sediment of the Chinese Loess Plateau (at least Lingtai and Jingchuan red clays) is eolian in origin. The quartz grains from dust deposits throughout the past seven million yeas showed the clues of glacial and frost processes. This indicates that the high mountains of western China reached a certain altitude to favor the glacial and/or frost processes at least seven millions years before. The weathering intensities of the past seven nnillion yeas have a decreasing trend. In about 5~4.5MaBP, the weathering is relatively weak, and the dust supply is relatively low. At about 3.6MaBP and 2.6MaBP, the dust supply increased significantly. The mineralogical composition, the quartz surface feature and the quartz oxygen isotope composition were influenced by the uplift of the Tibetan Plateau. The Plateau may have reached a certain altitude to generate the arid regions of inland China and favor the glacial and frost weathering. And it underwent a phased uplift, which have uplift events at about 3.6MaBP and 2.6MaBP.

2200万年来风尘沉积的粒度变化及古气候意义

Three eolian deposit formations, including Quaternary loess (QL, Liu et al.3 1985), Hipparion red earth (HRE, also called red clay, Liu et al., 1985) and Miocene loess (ML, Guo et al., 2002) constitute a set of unique paleoclimatic archives in northern China dated back to at least 22Ma ago. The Miocene loess is a recently discovered loess-soil sequence. Detailed investigation has been made on its origin, chronology and paleoclimatic significance (Guo et al., 2002), but further work is still needed to obtain detailed paleoclimate information, and mechanical links behind paleoclimatic changes. In this study, grain size analysis of QL, HRE and ML has been conducted on two sections: Qinan and Xifeng. The objective is focused on comparison of the grain size distribution characteristics (GSDC) among different eolian deposit formations, and reconstruction of the Asian monsoon circulation in the past 22 Ma. Results show that GDSC of ML resembles that of QL and HRE, and GDSC of ML is especially similar to HRE. Both ML and HRE contain a significant proportion of fine fraction, however, QL has a large amount of coarse sediments. This is mainly due to that the wind system transported aeolian dust was weaker in the late Tertiary than that in the Quaternary. Grain size difference between loess and paleosol in ML is much smaller than that in QL, indicating that the climatic fluctuations during the late Tertiary were much smaller than that happened in the Quaternary The grain size records of the past 22 Ma reveal several evolution phases of the Asia winter monsoon. -2.7 Ma BP is the most important boundary in the process of the winter monsoon evolution: the wind strength have significantly enhanced since 2.7 Ma ago. During a period between -22.0 and -3.6 Ma, three periods with relatively stronger winter monsoon are recorded in the QA-I section, between 21.2 and 19.9, and 16.0 and 13.3, and 8.7 to 6.9 Ma, respectively. From 3.6 to 2.7, the winter monsoon was enhanced gradually. In the Miocene time, the intensified winter monsoon phases (between 21.2 and 19.9, and 16.0 and 13.3, and 8.7 and 6.9 Ma) seemed to have a close relationship with the uplift of the Tibetan Plateau and/or the ongoing global cooling, but the forcing mechanism behind the Asia winter monsoon evolution need to be further investigated. During the Pliocene-Pleistocene time, the Asia winter monsoon strengthened at 3.6 and 2.7Ma ago are in good agreement with the ongoing global cooling and the Arctic ice sheet development. In the mean time, much evidence suggests that an intense uplift of the Tibetan Plateau occurred at ~3.6 Ma, which is synchronous with a major increase in Asia winter monsoon. Therefore, two major factors may be invoked to explain the winter monsoon enhancement: Arctic ice sheet development and Tibetan uplift. We propose that changes in location and intensity of the Siberian-Mongolian high that were caused by the Tibetan uplift and Arctic ice sheet development might be an important factor for Asian winter monsoon evolution in the Pliocene-Pleistocene.

青藏高原东部电性结构特征及地球动力学意义

As the largest and highest plateau on the Earth, the Tibetan Plateau has been a key location for understanding the processes of mountain building and plateau formation during India-Asia continent-continent collision. As the front-end of the collision, the geological structure of eastern Tibetan Plateau is very complex. It is ideal as a natural laboratory for investigating the formation and evolution of the Tibetan Plateau. Institute of Geophysics, Chinese Academy of Sciences (CAS) carried out MT survey from XiaZayii to Qingshuihe in the east part of the plateau in 1998. After error analysis and distortion analysis, the Non-linear Conjugate Gradient inversion(NLCG), Rapid Relaxation Inversin (RRI) and 2D OCCAM Inversion algorithms were used to invert the data. The three models obtained from 3 algorithms provided similar electrical structure and the NLCG model fit the observed data better than the other two models. According to the analysis of skin depth, the exploration depth of MT in Tibet is much more shallow than in stable continent. For example, the Schmucker depth at period 100s is less than 50km in Tibet, but more than 100km in Canadian Shield. There is a high conductivity layer at the depth of several kilometers beneath middle Qiangtang terrane, and almost 30 kilometers beneath northern Qiangtang terrane. The sensitivity analysis of the data predicates that the depth and resistivity of the crustal high conductivity layer are reliable. The MT results provide a high conductivity layer at 20~40km depth, where the seismic data show a low velocity zone. The experiments show that the rock will dehydrate and partially melt in the relative temperature and pressure. Fluids originated from dehydration and partial melting will seriously change rheological characteristics of rock. Therefore, This layer with low velocity and high conductivity layer in the crust is a weak layer. There is a low velocity path at the depth of 90-110 km beneath southeastern Tibetan Plateau and adjacent areas from seismology results. The analysis on the temperature and rheological property of the lithosphere show that the low velocity path is also weak. GPS measurements and the numerical simulation of the crust-mantle deformation show that the movement rate is different for different terranes. The regional strike derived from decomposition analysis for different frequency band and seismic anisotropy indicate that the crust and upper mantle move separately instead of as a whole. There are material flow in the eastern and southeastern Tibetan Plateau. Therefore, the faults, the crustal and upper mantle weak layers are three different boundaries for relatively movement. Those results support the "two layer wedge plates" geodynamic model on Tibetan formation and evolution.

井间地震技术研究和应用

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.

西南天山托云盆地白垩-新生代古地磁结果及其构造意义

Western China is regarded as an assemblage of blocks or microplates. The India/Asia postcollisional kinematics of these blocks has attracted many geologists to pay attentions, especially on the geodynamics and intracontinental deformation of Tibetan and adjoining parts of central Asia. So far there are still many debates on the amount of continental shortening and extrusion within Western China blocks. Paleomagnetism plays a very important role in the paleogeographic reconstruction and depiction of kinematics of the blocks, however the unequilibrium of paleomagentic data obtained from Western China prevents paleomagnetists from studying the kinematics and intracontinental deformation on the Tibetan plateau and the central Asia. Moreover, shallower inclinations observed in the Cretaceous and Cenozoic terrestrial red sediments in central Asia makes it difficult to precisely estimate the northward convergence of Tibetan plateau and its adjacent areas since the onset of the Indian/Asian collision. In this thesis, detailed rock magnetic, chronological and paleomagnetic studies have been carried out on the Tuoyun Basin in the southwestern Tianshan to discuss the possible continental shortening and tectonic movements since the Cretaceous-Tertiary. Ar-Ar geochronological study has been conducted on the upper and lower basalt series from the Tuoyun Basin, yielding that the lower and upper basalt series were extruded during 115-113 Ma and 61.8-56.9 Ma, respectively. Both the age spectrum and inverse isochron show that the samples from the upper and lower basalt series have experienced no significant thermal events since extrusion of the baslts. Rock magnetic studies including temperature dependence of magnetization and susceptibility during a heating-cooling cycle from temperature up to 600 ℃ suggest that the baslt samples from the lower and upper basalt series are ferromagnetically predominant of magnetite and a subordinate hematite with a few sites of titanomagnetite. The predominant magnetic mineral of the intercalated red beds is magnetite and hematite. Anisotropy of magnetic susceptibility shows that both the baslts and the intercalated red beds are unlikely to have undergone significant strain due to compaction or tectonic stress since formation of the rocks. The stable characteristic remanent magnetization (ChRM) isolated from the most samples of the upper and lower basalt series and intercalated red beds, passes fold test at the 99% confidence level. Together with the geochronological results, we interpret the characteristic component as a primary magnetization acquired in the formation of rocks. Some sites from both the upper and lower basalts yielded shallower inclinations than the reference field computed from the Eurasia APW, we prefer to argue that these shallow inclinations might be related to geomagnetic secular variation, whereas the shallow inclination in the intercalated red beds is likely to be related to detrital remanent magnetization. Paleomagnetic results from the early Cretaceous-Paleogene basalts indicate that no significant N-S convergence has taken place between the Tuoyun Basin and the south margin of Siberia. Furthermore, the Cretaceous and Tertiary paleomagnetic results suggest that the Tuoyun Basin was subjected to a local clockwise rotation of 20°-30° with respect to Eurasia since the Paleocene time, which is probably subsequent to the Cenozoic northward compression of the Pamir arc.

多波分量地震资料处理解释方法

To deal with some key problems in multi-component seismic exploration, some methods are introduced in this thesis based on reading amounts of papers about multi-component seismic theories and methods. First, to find a solution for the detection of the fracture density and orientation in igneous, carbonate and shale reservoirs, a large amount of which exist in domestic oil fields with low exploration and development degree, a new fast and slow shear waves separation method called Ratio Method based on S-wave splitting theory is discussed in this thesis, through which the anisotropy coefficient as well as fracture parameters such as density and azimuthal angle can be acquired. Another main point in this thesis involves the application of seismic velocity ratio (Vp/Vs) to predict the Hthological parameters of subsurface medium. To deal with the unfeasibility of velocity ratio calculation method based on time ratio due to the usually low single-noise ratio of S-wave seismic data acquired on land, a new method based on detailed velocity analysis is introduced. Third, pre-stack Kirchhoff integral migration is a new method developed in recent years, through which both S and P component seismic data as well as amplitude ratio of P/S waves can be acquired. In this thesis, the research on untilizing the P and S wave sections as well as amplitude ratio sections to interpret low-amplitude structures and lithological traps is carried out. The fast and slow shear wave separation method is then be applied respectively to detect the density and azimuthal angle of fractures in an igneous rock gas reservoir and the coal formation in a coal field. Two velocity ratio-calculating methods are applied respectively in the lithological prediction at the gas and coal field after summarizing a large amount of experimental results draw domestically and abroad. P and S wave sections as well as amplitude ratio sections are used to identify low-amplitude structures and lithological traps in the slope area of a oil-bearing sedimentary basin. The calculated data concerning fracture density and azimuthal angle through the introduced method matches well with the regional stress and actual drilling data. The predicted lithological data reflects the actual drilling data. Some of the low-amplitude and lithological traps determined by Kirchhoff migration method are verified by the actual drilling data. These results indicate that these methods are very meaningful when dealing with complex oil and gas reservoir, and can be applied in other areas.

基于可复制性的独立成分分析——应用于功能磁共振成像数据

This dissertation systematically depicted and improved the application of Independent Component Analysis (ICA) to Functional Magnetic Resonance Imaging (fMRI), following the logic of verification, improvement, extension, and application. The concept of “reproducibility” was the philosophy throughout its four concluded studies. In the “verification” study, ICA was applied to the resting-state fMRI data, verified the resultant components with reproducibility, and examined the consistency of the results from ICA and traditional “seed voxel” method. At the meantime, the limitation of ICA application on fMRI data analysis was presented. In the “improvement” study, an improved ICA algorithm based on reproducibility, RAICAR, was developed to aid some of the limitations of ICA application. RAICAR was able to rank ICA components by reproducibility, determine the number of reliable components, and obtain more stable results. RAICAR provided useful tools for validation and interpretation of ICA results. In the “extension” study, RAICAR as well as the concept of “reproducibility” was extended to multi-subject ICA analysis, and gRAICAR algorithm was developed. gRAICAR allows some variation across subjects, examining common components among subjects. gRAICAR is also capable to detect potential subject grouping on some components. It is a new way for exploratory group analysis on fMRI. In the “application” study, two newly developed methods, RAICAR and gRAICAR, were used to investigate the effect of early music training on the brain mechanism of memory and learning. The results showed brain mechanism difference in memory retrieval and learning process between two groups of subjects. This study also verified the usefulness and importance of the new methods.