青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

济阳坳陷深层地质构造与演化

Based on geophysical and geological data in Jiyang depression, the paper has identified main unconformity surfaces (main movement surfaces) and tectonic sequences and established tectonic and strata framework for correlation between different sags. Based on different sorts of structural styles and characteristics of typical structures, the paper summarized characteristics and distribution of deep structures, discussed evolution sequence of structure, analyzed the relation between tectonic evolution and generation of petroleum. The major developments are as following: Six tectonic sequences could be divided from bottom to top in the deep zone of Jiyang depression. These tectonic sequences are Cambrian to Ordovician, Carboniferous to Permian, lower to middle Jurassic, upper Jurassic to lower Cretaceous, upper Cretaceous and Kongdian formation to the fourth member of Shahejie formation. The center of sedimentation and subsidence of tectonic sequences distinguished from each other in seismic profiles is controlled by tectonic movements. Six tectonic evolution stages could be summarized in the deep zone in Jiyang depression. Among these stages, Paleozoic stage is croton sedimentation basin; Indosinian stage, open folds of EW direction are controlled by compression of nearly SN direction in early Indosinian (early to middle Triassic) while fold thrust fault of EW – NWW direction and arch protruding to NNE direction are controlled by strong compression in late Indosinian (latter Triassic); early Yanshanian stage (early to middle Jurassic), in relatively weak movement after Indosinian compressional orogeny, fluviolacustrine is deposited in intermontane basins in the beginning of early Yanshanian and then extensively denudated in the main orogenic phase; middle Yanshanian (late Jurassic to early Cretaceous), strike-slipping basins are wide distribution with extension (negative reversion) of NW – SE direction; latter Yanshanian (late Cretaceous), fold and thrust of NE – NNE direction and positive reversion structure of late Jurassic to early Cretaceous strike-slipping basin are formed by strong compression of NW–SE direction; sedimentation stage of Kongdian formation to the fourth member of Shahejie formation of Cenozoic, half graben basins are formed by extension of SN direction early while uplift is resulted from compression of nearly EW direction latterly. Compression system, extension system and strike-slip system are formed in deep zone of Jiyang depression. According to identifying flower structure of seismic profiles and analysis of leveling layer slice of 3D seismic data and tectonic map of deep tectonic interface, strike-slip structures of deep zone in Jiyang depression are distinguished. In the middle of the Jiyang depression, strike-slip structures extend as SN direction, NNW direction in Huimin sag, but NNE in Zhandong area. Based on map of relict strata thickness, main faults activity and regional tectonic setting, dynamic mechanisms of deep structure are preliminary determination. The main reason is the difference of direction and character of the plate’s movement. Development and rework of multi-stage tectonic effects are benefit for favorable reservoir and structural trap. Based on tectonic development, accumulation conditions of deep sub-sags and exploration achievements in recent years, potential zones of oil-gas reservoir are put forward, such as Dongying sag and Bonan sag.

残留盆地分布综合地球物理研究——以大港地区为例

Macro-distribution of residual basins is a basic question in residual basin research，the main object of macro-distribution study is to build strata framework, compute thickness of residual strata and analyze characteristics of residual basins. With the guidance of the theory of integrated geology and geophysical research, the paper assembled series of methods and established the technical chart based on gravity and magnetic data, with restriction of geology, seismic and drilling data. Based on potential field data processing and analysis, forward and inverse computation, region potential field analysis and potential field separation, etc. it computed depth of gravity/magnetic basement and got strata framework. It had got effective results in the research of macro-distribution of residual basin research in the Dagang area. It did the wavelet transform of gravity/magnetic data with multi-kind of wavelet basis using a trou algorithm. From comparison of processing result and their spectral of wavelet analysis, up continuation and filter method, the wavelet approximation is better to fit the regional potential field, and it is an effective method to separate gravity/magnetic effect caused by deep geology bodies. The experiment of matching pursuit shows that te transform domain methods have great advantage in potential data analysis. From the integrated geophysical study of rock property study, gravity/magnetic basement inversion and fault system analysis of the Dagang area, it gets the strata framework and the thickness of pre-Cenozoic residual strata. Comprehensive study with gravity and magnetotelluric profile inversion and interpretation, three prospect plays of macro-distribution of residual basins are fingered out. It has great residual strata thickness in the northern part of Chengning Uplift and there is thrust fault in the deep zone and good up-Paleozoic hydrocarbon source rocks in this area. With integrated analysis, this area will be the most prospective hydrocarbon location of pre-Cenozoic residual basins.

黄海海域前新生代残留盆地综合地质地球物理研究

China’s annual oil import volume has been increasing in recent years, but the oil price in the international market fluctuates and poses a severe threat to China’s economic development and national security. Therefore, it is of great importance to study the gas and oil exploration of Pre-Cenozoic Residual Basins in Yellow Sea. Yellow Sea has widespread and thick Mesozoic and Paleozoic strata that contain multilayer source rock. Hence, Yellow Sea Mesozoic and Paleozoic strata have good conditions of forming Pre-Cenozoic hydrocarbon reservoirs. Pre-Cenozoic Residual Basins are usually buried deep and then transformed many times in its long evolutional history. These characteristics make it difficult to apply a single method in exploring Pre-Cenozoic Residual Basins. On the other hand, it is highly effective to solve key problems of gas and oil exploration of Pre-Cenozoic Residual Basins in Yellow Sea by using integrated geological and geophysical methods which make full use of the advantages of various exploring techniques. Based on the principle of “the region controls the local; the deep restricts the shallow,” this study focuses on Pre-Cenozoic Residual Basins in Yellow Sea to describe the structure frame of its distribution, with gravity, magnetic, seismic, drill-hole and geological data and previous research findings. In addition, the distribution characteristics of Pre-Cenozoic Residual Basins in Yellow Sea are also analyzed. This paper explores the characteristics of error between gravity forward with constant density and gravity forward with variable density through the study on 2-D and 3-D gravity forward in frequency domain. The result shows that there is a linear relationship between error and depth of 2-D geological model but there is a nonlinear relationship between error and depth of 3-D geological model. The error can be removed according to its linear characteristics or statistical nature of nonlinear characteristics. There is also error between gravity inversion with constant density and gravity inversion with variable density due to variable density and edge-effect. Since there are not noticeable rules between the error and the two causes as variable density and edge-effect, this study adopts gravity inversion with variable density and methods to eliminate the edge-effect in basement inversion to improve inversion accuracy. Based on the study on the rock physical properties and strata distribution of Yellow Sea and adjacent regions, this study finds that there is a big density contrast between Cretaceous-Jurassic strata and their substratum. The magnetic basement of south Yellow Sea is regarded as top of Archeozoic-Proterozoic early strata, and there are double magnetic basements in north Yellow Sea. Gravity and magnetic data are used to inverse the gravity basement and magnetic basement of Yellow Sea, with seismic and drill-hole data as constrains. According to data of gravity and magnetic basement distribution, the depth of Cenozoic strata and previous research findings, this paper calculates the thickness of the Mesozoic and Pre-Mesozoic Residual Basins, draws the distribution outline of Pre-Cenozoic Residual Basins in Yellow Sea, and analyzes its macro-distribution characteristics. Gravity inversion is applied on a typical geological profile in Yellow Sea to analyze the characteristics of its fractures and magnetic basements. The characteristics of Pre-Cenozoic Residual Basins distribution outline in Yellow Sea and the fractures and magnetic basements of its typical profile shown by profile inversion provides new geophysical evidence for these structure views such as “the South Yellow Sea and the North Yellow Sea belong to different structural units” and “Sino-Korea and Yangtze blocks combine along Yellow Sea East Fractured Zone in Yellow Sea”.

北天山新生代构造演化与古环境演变

The Tianshan Mountains is located about 1000-2000 km north of the India-Asia suture and is the most outstanding topography in central Asia, with transmeridional length of nearly 2500 km, north-southern wideness of ~ 300-500 km, peaks exceeding 7000 m above sea level (asl.), and average altitude of over 4000 m asl. Much of the modern relief of the Tianshan Range is a result of contraction driven by the collision of the India subcontinent with the southern margin of Asia, which began in early Tertiary and continues today. Understanding where, when and how the deformation of the Tianshan Mountains occurred is essential to decipher the mechanism of intracontinental tectonics, the process of foreland basin evolution and mountain building, and the history of climate change in central Asia. In order to better constrain the Cenozoic building history of the Tianshan Mountains and the climate change in the southern margin of the Junggar Basin, we carried out multiple studies of magnetostratigraphy, sedimentology, and stable isotopes of paleosol carbonate at the Jingou River section, which is located at the Huoerguosi anticline, the westernest one of the second folds and thrust faults zone in the northern piedmont of the Tianshan Mountains. The Jingou River section with a thickness of about 4160 m is continuous in deposits according to the observed gradual change in sedimentary environments and can be divided into five formations: Anjihaihe, Shawan, Taxihe, Dushanzi and Xiyu in upward sequence. Characteristic remamences were isolated by progressive thermal demagnetization, generally between 300 and 680℃. A total of 1133 out of 1607 samples yielded well-defined ChRMs and were used to establish the magnetostratigraphic column of a 3270-m-thick section from the exposed base of the Anjihaihe Formation to the middle of the Xiyu Formation. Two vertebrate fossil sites and a good correlation with the CK95 geomagnetic polarity time scale suggest that the section was deposited from ~30.5 to ~4.6 Ma and the age of the top of the Xiyu formation is ~2.6 Ma based on an extrapolation of the sedimentation rates. A plot of magnetostratigraphic age vs. height at the Jingou River section shows that significant increases in sedimentation rates as well as notable changes in depositional environments occurred at ~26-22.5 Ma, ~13-11 Ma and ~7 Ma, which represent the initial uplift of the Tianshan Mountains and two subsequent rapid uplift events. In addition, changes in sedimentation rates display characteristic alternations between increases and decreases, which probably indicate that the uplift of the Tianshan Mountains was episodic. We discussed the history of C4 biomass and climatic conditions in the southern margin of the Junggur Basin using the stable carbon and oxygen isotope composition of paleosol carbonates from the Jingou River section during ~17.5-6.5 Ma. The δ13C values indicate that the proportion of C4 biomass was uniform and moderate (15-20 %) during the interval of ~17.5-6.5 Ma. We proposed three hypotheses for this pattern of C4 biomass: (1) counteraction of two opposed factors (global cooling since ~15 Ma and thereafter increased dry and seasonality in central Asia) controlling the growth of C4 grasses, (2) variability in abundance of C3 grasses relative to C3 trees and shrubs if vegetation had ever changed in ecosystems, and (3) the higher latitude of the studied region. The δ18O values show a stepwise negative trend since ~13 Ma which may be attributed to three factors: (1) the temperature decreasing gradually after the middle Miocene (~15 Ma), (2) the increasing contribution of the moistures carried by the polar air masses from the Arctic Ocean to precipitation, and (3) the gradual retreat westward and disappearance of the Paratethys Ocean. Among them, which one played a more important role will need further study of the paleoclimate in central Asia.

中国华北和美国加州地区地震层析成像研究

In recent years seismic tomography has become a powerful tool for studying the three-dimensional crust and mantle structure. In this study, we collected a large number of regional and teleseismic travel-time data and used seismic tomography method to study the relationship between earthquake occurrence and crustal heterogeneity for the 1992 Landers earthquake, heterogeneity and evolution of lithosphere under North China Craton and Southern California, and deep structure and origin of the Changbai intraplate volcano in Northeast China. Our results show a correlation between the seismic rupture zone and crustal heterogeneity. The distribution of the Landers aftershocks is cluster-like and separated or terminated in areas where low-velocity anomalies exist．Most of the large earthquakes with magnitudes >4.0 occurred in or around areas with high P-wave velocity．The possibility is that high-velocity areas are brittle and strong parts which can sustain seismogenic stress，and so can generate earthquakes. Our tomographic images show a very heterogeneous structure in the crust and upper mantle beneath Southern California. Three major anomalies in the upper mantle are revealed clearly beneath the southern Sierra Nevada, Transverse Ranges and Salton Trough. We consider that the high-velocity anomaly beneath the Transverse Ranges was formed through asymmetrical two-side convergence of subcrustal lithosphere and sinking to asthenosphere. Formation of the dense crust root and “drip structure” caused the high-velocity anomaly under the southern Sierra Nevada. The Salton Trough low is the response to the lithospheric extension when the Pacific plate was rifted away from the North American Plate. The tomograpic images beneath the North China Craton show that there exist different lithospheric structures under the different blocks. Complex, prominent low-velocity and high-velocity anomalies are imaged beneath the North China Basin, Trans-North China Orogen (TNCO), and Ordos Block which correspond to rifted, orogenic and cratonic lithospheres, respectively. The thickness of the three-type lithospheres is about 70, 90 and >250 km, respectively. Our results suggest that lithospheric thinning under the eastern part of North China Craton is due to long-term replacement and chemical and thermal erosion of the ancient lithosphere by the hot asthenosphere. The remains of ancient lithosphere exist either in the present upper mantle or have sunk into the mantle transition zone. Our tomographic result of the Changbai volcanic area suggests that the origin of the Changbai volcano is related to the deep dehydration of the subducted Pacific slab and corner flow in the big mantle wedge (BMW) above the stagnant Pacific slab.

西南极南设得兰群岛中新生代火山岩基底沉积学及其大地构造制约

Livingston Island, the second island of South Shetland Island, constains Mesozoic-Cenozoic basement, Mesozoic-Cenozoic volcanic sequences, plutonic intrusions and post-subduction volcanic rocks, which document the history and evolution of an important part of the South Shetland Islands magmatic arc. The sedimentary sequence is named the Miers Bluff Formation (MBF) and is interpreted as turbidite since the first geological study on South Shetland Islands, and is interpreted as turbidite. It base and top are not exposed, but a thickness of more than 3000m has been suggested and seems plausible. The turbidite is overlain by Mid - Cretaceous volcanic rocks and intruded by Eocene tonalites. The age of the Miers Bluff Formation is poorly constrained Late Carboniferous -Early Triassic. Sedimentary Environment, tectonic setting and forming age of sedimentary rocks of the Miers Bluff Formation were discussed by means of the methods of sedimentology, petrography and geochemistry, combinig with the study of trace fossils and microfossil plants. The following conclusions are obstained. A sedimentary geological section of Johnsons Dock is made by outside measuring and watching, and then according the section, the geological map near the Spanish Antarctic station was mapped. Four pebbly mudstone layers are first distinguished, which thickness is about 10m. The pebbly mudstone is the typical rock of debris flow, and the depostional environment of pebbly mudstone may be the channel of mid fan of submarine fan. The sedimentsry structural characteristics and size analysis of sandstones show the typical sedimentary feature of turbidity flow and the Miers Bluff Formation is a deep-water turbidite (include some gravity-flow sediments). The materials of palaeocurrents suggest the continental slope dip to southeast, and indicate the provenance of turbidity sediment in the northwest area. By facies analysis, six main facies which include seven subfacies were recognized, which are formed in mid-fan and lower-fan of submarine, meanwhile, the sedimentary features of each facies and subfacies are summarized. The study of clastic composition, major elements, trace elements and rare earth elements indicates the forming setting of the Miers Bluff Formaton is active continental margin and continental island arc and the provenance is dissected magmatic arc which main composition is felsic gneiss. Many trace fossils of the whole succession were found in the turbidites of the Miers Bluff Formation. All these trace fossils are deep sea ichnofossils. There are fifteen ichnogenus, sixteen ichnospecies. Moreover, a new trace fossil was found and a new ichnogenus and new ichnospecies was proposed - Paleaichnus antarctics ichnogen, et ichnosp, nov.. Except the new ichnogenus and ichnospecies, others had been found in deep-sea flysch turbidites. Some are in mudstone and are preserved in the cast convex of overlying sandstone sole, they formed before turbidity flows occurred and belong to the high-different Graphoglyptida of fiysch mudstone. Others as Fucusopsis and Neonereites are preserved in sandstones and stand for trace assemblages after turbidity sedimentation. These trace fossils are typical members of abyssal "Nereites" ichnofacies, and provide for the depositional environment of the Miers Bluff Formation. Fairly diverse microfossil plants have been recovered from the Miers Bluff Formation, Livingston Island, including spores, pollen, acritarchs, wood fragments and cuticles. Containing a total of about 45 species (forms) of miospores, the palynofiora is quantitatively characterized by the dominance of non-striate bisaccate pollen, but spores of pteridophytes and pollen of gymnosperms are proportionate in diversity. It is somewhat comparable to the subzone C+D of the Alisporites zone of Antarctica, and the upper Craterisporites rotundus zone and the lower Polycingulatisporites crenulatus zone of Australia, suggesting a Late Triassic (possibly Norian-Rhaetian) age, as also evidenced by the sporadic occurrence of Aratrisporites and probable Classopollis as well as the complete absence of bisaccate Striatiti. The parent vegetation and paleoclimate are preliminarily deduced. At last, the paper prooses the provenance of sedimentary rocks of the Miers Bluff Formation locates in the east part to the southern Chile(or Southern South American). In the Triassic period, contrasting with New Zealand, Australia and South American of the Pacific margin of Gondwanaland, the Miers Bluff Formation is deposited in the fore-arc basin or back-arc basin of magmatic arc.

渤海湾地区末次冰期古环境

The loess-paleosol sequence in China is one of the best archives for studying paleoenvironmental and paleoclimatic processes. The loess deposits in the coast of the Bohai Sea are suitable for the study of aridification in the northern China during glacial periods (fig.2-1). In this paper, stratigraphy was correlated by using magnetic susceptibility, grain size, and thermoluminence (TL) and accelerater mass spectrometer (AMS) ~(14)C ages. Based on the loess records, an interpretation has been made for paleoenvironmental changes on the coast of the Bohai Sea during the last glacials. The results of magnetic susceptibility, grain size, biostratigraphy, TL and ~(14)C dating, suggest that the loess-paleosol sequence in the coast of the Bohai Sea is discontinuous. The loess deposits correlated with the marine δ~(18)O stage 2 are usually absent in some profiles. Also, the thickness of the loess deposits in the same period varied significantly in different sections. In the coast of the Bohai Sea, the higher magnetic susceptibility corresponds to the finer grain size, consistent with the results of the Loess Plateau. It is indicated that the changes in magnetic susceptibility and grain size may record the paleoclimatic fluctuations of the last glacial. Although the loess deposits during the last glacial have been slightly altered by slope runoff, they have still remained main characters as the representative loess deposits of the Loess Plateau. During the last glacial, the less accumulation rate in the coast of the Bohai Sea is similar to that of the desert-loess transition zone in the northwestern Loess Plateau, and the all section contain high concentrations of sand (>60μm), indicate that the aridification in the Bohai Sea occurred during the glacial. But the changes in sand content of loess deposits along a north-south transect of the Bohai Sea and the changes of magnetic susceptibility implicate that desertification might not occur in the shelf of the Bohai Sea during the last glacial. The frequent fluctuations of summer monsoon during the marine δ~(18)O stage 4 are demonstrated by magnetic susceptibility, frequency-dependent susceptibility and the abundance of foraminifera. 46 genera, 71 species of foraminifera were identified from 138 loess samples. Almost all of the foraminifera are present in the last glacial loess, but the distribution patterns of foraminifera show significant temporal changes. The results of magnetic susceptibility, grain size and XRD indicate that not only sea-level changes had influences on foraminifera abundance, but also variations in sediment flux by rivers when sea-level drops might control the abundance of foraminifera. In addition, the diversity and exquitability of the foraminifera suggest that the frequent fluctuations of foraminifera fossil abundance during the marine δ~(18)O stages could be partly attributed to leaching. In summary, the changes in foraminifera of abundance related to the paleoclimatic variations, and the aridification extends to the coast of the eastern China during the last glacial.

大别－苏鲁造山带与华北陆块东部岩石层三维速度结构

The 3-D velocity images of the crest and upper mantle beneath the region of 112° -124°E, 28°-39°N including the Dabie-Sulu orogenic belt are reconstructed by using 36405 P-wave arrivals of 3437 regional and 670 distant earthquakes during the period from 1981 to 1996, and gridding the area of 0.5° * 0.5°. The results of tomography demonstrate that: 1. The results of tomographic imaging show a broad heterogeneity in P wave velocity structure for the lithosphere beneath the Dabie-Sulu orogenic belt. 2. In the Dabie orogenic belt, the velocity patterns in the crust are different among various tectonic units. The Dabie and Qinling orogenic belts are remarkable in the tomographic images, and in mm the Hongan and Dabie blocks in the Dabie orogenic belt are also imaged very distinguishably. 3. A velocity (about 5.9～6.0 km/s) layer exists in the Dabie block at depth between 15～25 km, which is coincident with the low-resistance layer at the depth of 12-23 km, being inferred to be the tectonic detachment zone and suggesting that the extension detachment structure was formed in the middle crust. Beneath the southern and northerm Dabie tectonic units, the north-dipping high-velocity (at level of 6.5 ～ 6.6 km/s) block was developed in the crust, which might be correlated with the UHP rockswith low content of the meta-ultramafic rocks. This result is in agreement with the geological observation on the surface. 4. The velocity image at 40 km depth reveals the features at the top of mantle and the configuration of the Moho discontinuity. The depth of the Moho changes slightly along the trend of the orogenic belt. It in Hongan block is less than 40 km, but it is different in the western and eastern parts of the Dabie block, the former is more than 40 km, and the latter less than or equal to 40 km. The remnant of the mountain root exists between the Shangcheng-Macheng fault and the line of Huoshan-Yuexi-Yingshan in the Dabie orogenic belt, and beneath the southern and northern Dabie tectonic units. However, the thickness of the Moho is about 40 km and there is no obvious changes, which suggest that the Dabie orogenic belt has been experienced quite in the gravity equilibration. The Moho's depth in the Sulu is less than 40 km. 5. There is a dipping slab-like high-velocity body in the uppermost mantle. It is sandwiched by slow velocities and exists beneath the Dabie-Sulu orogenic belt in the range of depths between the Moho discontinuity and 110 km at least. This high-velocity body outlines a picture of the slab interpreted as the remnant of the Triassic subducted YZ. 6. The Sulu orogenic belt displays "crocodilian" velocity structure, the upper crust of the Yangtze thrusted over the Huabei crest, and the Huabei crust indented into the Yangtze crust, where the ancient subduction zone of the Yangtze lithosphere located. Based on the previous geological data, this structure is not related with the collision between the Yangtze and Sino-Korean Blocks, but caused by the sinistral offset of the Tan-Lu Fault. Studied on the velocity structure of the eastern Huabei lithosphere indicates: 1. The 'present-day' lithosphere of the eastern Huabei is between 40-100 km thick with greatly thinned lithosphere around the Bohai Sea. Generally, thickness of the lithosphere in this region decreased eastwards. 2. The attenuation of the lithosphere is attributed to the strongly uplift of the asthenosphere. In the area between the Taihang Mountains and the Tan-Lu Fault, there is a 'lever' with red low velocity belt, it is clearly defined, transverse continuity, depth between 100-150 km, local variations visible, and an upwards trend towards the Bohai Sea. Generally, the velocity structure in the mantle beneath the lithosphere displays irregular column-shape consisting of alternating high and low velocities, and when cold high velocity ancient lithosphere connects with the hot low velocity mantle materials forming precipitous compact structure. More heat pathways from the mantle occur towards the Tan-Lu Fault. 3. The strongly irregular characteristics of the contact between the asthenosphere and the lithosphere is induced by the long-term hot, chemical erosion and alteration on the contact. 4. There are still preserved high velocity lithosphedc root beneath Huabei with 'block-shape' distribution and surrounded by hot materials. Results of our studies indicate that the evolution models of the eastern China mantle are characterized by the direct contact between the uplifted lithosphere and the Huabei Craton accompanying the upwelling of the deep mantle materials. At the contact betwen the lithosphere and the asthenosphere, the upwelled mantle materials replaced and altered the lower lithosphere forming the metasome through the hot and chemical modifications impacted on the Craton lithosphere, and changed it into the lithosphere gradually, resulting in the lithospheric thinning. Thus, the lithospheric thinning is the result of the upwelling of the asthenosphere.

沾化凹陷泥质岩油气藏成藏特征研究

Mudstone reservoir is a subtle reservoir with extremely inhomogeneous, whose formation is greatly related to the existence of fracture. For this kind of reservoir, mudstone is oil source rock, cover rock and reservoir strata, reservoir type is various, attitude of oil layer changes greatly, and the distribution of oil and gas is different from igneous or clastic rock reservoir as well as from carbonate reservoir of self-producing and self-containing of oil and gas. No mature experience has been obtained in the description, exploration and development of the reservoir by far. Taking Zhanhua depression as an example, we studied in this thesis the tectonic evolution, deposit characteristics, diagenesis, hydrocarbon formation, abnormal formation pressure, forming of fissure in mudstone reservoir, etc. on the basis of core analysis, physical simulation, numerical simulation, integrated study of well logging and geophysical data, and systematically analyzed the developing and distributing of mudstone fissure reservoir and set up a geological model for the formation of mudstone fissure reservoir, and predicted possible fractural zone in studied area. Mudstone reservoir mainly distributed on the thrown side of sedimentary fault along the sloping area of the petroleum generatiion depression in Zhanhua depression. Growing fault controlled subsidence and sedimentation. Both the rate of subsidence and thickness of mudstone are great on the thrown side of growing fault, which result in the formation of surpressure in the area. The unlocking of fault which leads to the pressure discharges and the upward conduct of below stratum, also makes for the surpressure in mudstone. In Zhanhua depression, mudstone reservior mainly developed in sub-compacted stratum in the third segment of Shahejie formation, which is the best oil source rock because of its wide spread in distribution, great in thickness, and rich in organic matter, and rock types of which are oil source mudstone and shale of deep water or semi-deep water sediment in lacustrine facies. It revealed from core analysis that the stratum is rich in limestone, and consists of lamina of dark mudstone and that of light grey limestone alternately, such rock assemblage is in favor of high pressure and fracture in the process of hydrocarbon generation. Fracture of mudstone in the third segment of Shahejie formation was divided into structure fracture, hydrocarbon generation fracture and compound fracture and six secondary types of fracture for the fist time according to the cause of their formation in the thesis. Structural fracture is formed by tectonic movement such as fold or fault, which develops mainly near the faults, especially in the protrude area and the edge of faults, such fracture has obvious directivity, and tend to have more width and extension in length and obvious direction, and was developed periodically, discontinuously in time and successively as the result of multi-tectonic movement in studied area. Hydrocarbon generation fracture was formed in the process of hydrocarbon generation, the fracture is numerous in number and extensively in distribution, but the scale of it is always small and belongs to microfracture. The compound fracture is the result of both tectonic movement and hydrocarbon forming process. The combination of above fractures in time and space forms the three dimension reservoir space network of mudstone, which satellites with abnormal pressure zone in plane distribution and relates to sedimentary faces, rock combination, organic content, structural evolution, and high pressure, etc.. In Zhanhua depression, the mudstone of third segment in shahejie formation corresponds with a set of seismic reflection with better continuous. When mudstone containing oil and gas of abnormal high pressure, the seismic waveform would change as a result of absorb of oil and gas to the high-frequency composition of seismic reflection, and decrease of seismic reflection frequency resulted from the breakage of mudstone structure. The author solved the problem of mudstone reservoir predicting to some degree through the use of coherent data analysis in Zhanhua depression. Numerical modeling of basin has been used to simulate the ancient liquid pressure field in Zhanhua depression, to quantitative analysis the main controlling factor (such as uncompaction, tectonic movement, hydrocarbon generation) to surpressure in mudstone. Combined with factual geologic information and references, we analyzed the characteristic of basin evolution and factors influence the pressure field, and employed numerical modeling of liquid pressure evolution in 1-D and 2-D section, modeled and analyzed the forming and evolution of pressure in plane for main position in different periods, and made a conclusion that the main factors for surpressure in studied area are tectonic movement, uncompaction and hydrocarbon generation process. In Zhanhua depression, the valid fracture zone in mudstone was mainly formed in the last stage of Dongying movement, the mudstone in the third segment of Shahejie formation turn into fastigium for oil generation and migration in Guantao stage, and oil and gas were preserved since the end of the stage. Tectonic movement was weak after oil and gas to be preserved, and such made for the preserve of oil and gas. The forming of fractured mudstone reservoir can be divided into four different stages, i.e. deposition of muddy oil source rock, draining off water by compacting to producing hydrocarbon, forming of valid fracture and collecting of oil, forming of fracture reservoir. Combined with other regional geologic information, we predicted four prior mudstone fracture reservoirs, which measured 18km2 in area and 1200 X 104t in geological reserves.

南海天然气水合物地热学研究

Through generalizing the thermal field characteristics in gas hydrates distribution area in the world, the favorable thermal conditions for gas hydrates in the South China Sea are analyzed firstly. On the basis of above analysis, focused on the gas hydrates stability zone (GHSZ), the dissertation initiated the gas hydrates studies with geothermal methods in the South China Sea which will provide useful information for gas hydrates resource exploration and evaluation in the future. On the basis of study on hydrates phase equilibrium and the GHSZ affecting factors, the potential planar distribution of gas hydrates is determined by studying the temperature and pressure conditions in the sea bottom with different water depth, and the thickness of GHSZ is attained by solving the hydrates phase boundary curve equation and geothermal gradient curve equation. The result shows that, if the chemical composition of hydrocarbons contains methane only and the salt content of water is 3.5%, hydrates can form and keep stable at sea bottom at water depth not less than 550m, and the thickness of GHSZ is more than 200m in Xisha Through, Southeastern area of Dongsha Islands, Southwestern basin of Taiwan Island, northern area of Nansha Trough. The GHSZ is thicker with heat flow, geothermal gradient, and thermal conductivity decreasing, and with water depth increasing. Geothermal field simulating also attains the base of GHSZ in Xisha through, which is less than the depth of BSR. Although the present T-P conditions is not the most favorable for gas hydrates through 6Ma history, gas hydrates are still profitable in Xisha Through, Southeastern area of Dongsha Islands, Southwestern basin of Taiwan Island, Luzon Trough and northern area of Nansha Trough by systemic study on the sedimentary and structural characteristics, the conditions of T-P and natural gas source, considering geochemical and geophysical indications found in the South China Sea.

天然地震转换波偏移成像方法研究

Butovskaya, a scholar of Former Soviet Union, first determined the depth of basalt layer in Tashkent Zone by using converted waves on seismogram in 1952. From then on, more and more scholars developed the comprehensive research that imaged the earth interior structures by applying converted waves information. With the digitalization of earthquake observation, The inversion imaging of complete or partial waveform record can efficiently improve inversion quality and widen its usage scope, therefore great progress is made in converted wave imaging. This paper makes a certain study in converted wave imaging on that basis. Transmitted PP waves and converted PS waves are generated when a P-wave propagates through an interface separating two media with large impedance contracts. A PS converted wave is a seismic body wave, which result from the conversion of an incident parent P wave at a boundary within the crust to a refracted S wave. The thickness of a single crustal layer can theoretically be determined by observing, with three-componented seismometer at a single station, the difference in time of the arrival of the parent P wave and the arrival of the PS converted wave. For a multilayered media, PS converted wave arrivals corresponding to each of the layers can theoretically be observed, provided the station is sufficiently from the source of the parent P wave to allow initial penetration of the P wave beneath the deepest layer considered. To avoid the difficulty of picking up transmitted P-wave and converted wave phases, this paper proposed a converted wave migration method by estimating the travel time difference between PS converted wave and PP transmitted wave. To verify its validity, we apply the converted wave PS migration algorithm to synthetic data generated by three forward modeling. The migration results indicate that PS converted wave may be migrated to reconstruct the transmitting interface. This technique is helpful to investigate the deep earth structures by using earthquake data.

濮城油田沙三中6-10高分辨率层序地层学研究

The foundation of reservoir model and residual oil prediction have been the core of reservoir detailed description for improved oil production and enhanced oil recovery. The traditional way of sandstone correlation based on the geometrical similarity of well-logs which emphasizes "based on the cycle and correlating from larger to smaller" has shown its theoretical limits when explaining the correlating and the scale, geometry, continuity, connectivity of sandstones and the law of the reservoir property. It has been an urgent and difficult subject to find new theory and methods to solve the reservoir correlation and property prediction. It's a new way to correlate strata and found framework of reservoir through the process-response analysis in the base-level cycles. And it is also possible to analyze the reservoir property in reservoir framework. Taking the reservoir of zonation 6-10 in S3~2 of Pucheng Oil Field in Henan Province as an example, we founded the detailed reservoir stratigraphic framework through base-level correlation. In the strata frame, sediment distribution and its development are discussed based on sediment volume partitioning and facies differentiation analysis. Reservoir heterogeneities and its relation to base-level are also discussed. The analysis of primary oil distribution shows the base-level controlled oil distribution in reservoir. In this paper, subjects as following are discussed in detail. Based on the analysis of sedimentary structure and sedimentary energy, the facies model was founded. Founding stratigraphy framework through base level analysis In the studying zone, one long term cycle, 6 middle term cycles and 27 short term cycles was identified and correlated. 3 Predicting the property of reservoir for improving oil development The base level controlled the property of sandbody. The short and very short term cycle controlled the pattern of heterogeneities in sandbody, and the middle and long term cycle controlled the area and inter-layer heterogeneities. On the lower location of the middle and long term base level, the sandbody is well developed, with a wide area and large thickness, while on the high location of base level, there is an opposite reservoir character. 4 The studying of reservoir development response and oil distribution making a solid base for development adjustment Primary oil distribution is controlled by base level location. It tells that the sandbody on the high base level location was poor developed for its difficulty to develop. While on the low location of the base level, the sandbody is well developed for its relative easy to develop and dominant role in the development, but high residual oil for its high original oil content.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

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

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.