南海大陆边缘地壳结构与盆地形成机制--以珠江口－琼东南盆地为例

South China Sea is located in the convergence of Eurasian plate, the Pacific Ocean plate and Indian Ocean-Australia plate. The total area is about 3,500,000 km2, the geologic structure is complicated, and the structure line cut off reciprocal is the marginal sea taking form by that the seafloor spreads during the middle Oligocene. South China Sea continental margin have developed more than 10 large oil-gas bearing basins and a number of medium-small sized basins. These basins contain abundant mineral resources such as oil & gas. The marginal deepwater area in the north part of South China Sea has become our country’s strategic energy prospecting frontier. The deepwater area of Zhujiangkou and Qiongdongnan basins is the research target in this thesis. The thesis studied deep structure and the earth dynamics of the north part of South China Sea margin, and these researches provide scientific basis for oil-gas resources strategic investigation and valuation in deepwater sea area of north part slope of South China Sea. In order to develop the research of rebuilding velocities and density architecture of earth shell in region of interest, in marginal deepwater area in the north part of South China, we adopted 14 long-cable seismic reflection profile data of 3556.41 kilometers in total, the gravity measurement data along profiles (3851.44 kilometers in total), the magnetic observation along profiles (3838.4 kilometers in total) and depth measurement along profile, the logging data of 11 wells in project, the interpreted fault parameter and preexisting geologic and geophysical research achievement. This thesis has carried out concretely studying research as follows: 1. Overlay-velocity data sampling and analysis, interval velocity calculation, time-depth conversion, model building of earth shell velocity and layering character of earth shell are studied on 14 deep sections. Velocity structure in region of interest has revealed: Changchang is the sag with thinnest crust in Qiongdongnan basin; the sedimentary thickness lowers gradually from north to south, and the thickness change from west to east is milder. The sags’ sedimentary velocities in Qiongdongnan basin have obvious demarcation. The velocity of the 8000 meters sedimentary rocks is 4700 m/s in Shunde sag and Baiyun sag, and is the lowest; at that depth, the velocity very different in Liwan sag and Baiyun sag, which is about 800m/s. 2. Extracting gravity data and building of initial crust density model along the section; With Bouguer gravity anomaly data as constraint, revising density distributes of initial model, and building the crust density model. 3. With crust velocity and density as constraint, correcting the effect of thermobaric field and constructing constitution structure of rock in region of interest. By this research, we known that rocks in Zhujiangkou upper crustal layer are chiefly granite-gneiss, quartzite, granodiorite and basalt, however, rocks in Qiongdongnan basin upper earth shell are chiefly composed of granite-gneiss, quartzite, granodiorite, diorite and basalt. 4. Synthetically crust velocity and density structure, gaining expanding factor on crust and entire crust along section. The result is indicated: the expanding factor in every sag rises from northwest to southeast, which have reflected thinning characteristic of crust from continent to ocean. Intra-crustal deformation degree in Changchang and Ledong-Lingshui sag is bigger than that in Songnan-Baodao sag. Entire crust extension factor in Changchang and Songnan-Baodao sag is greater than that in Ledong-Lingshui sag, which can make an explanation of frequently event and longer heating process in middle-east of Qiongdongnan basin. 5. Synthesize multidisciplinary information to discuss the earth dynamics significance of discordogenic seismic profile in deepwater area of Zhujiangkou and Qiongdongnan basins.

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

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.

Pg震相走时的各向异性影响及各向异性结构成像

As the first arrival of seismic phase in deep seismic sounding, Pg is the important data for studying the attributes of the sedimentary layers and the shape of crystalline basement because of its high intensity and reliable detection. Conventionally, the sedimentary cover is expressed as isotropic, linear increasing model in the interpretation of Pg event. Actually, the sedimentary medium should be anisotropic as preferred cracks or fractures and thin layers are common features in the upper crust, so the interpretation of Pg event needs to be taken account of seismic velocity anisotropy. Traveltime calculation is the base of data processing and interpretation. Here, we only study the type of elliptical anisotropy for the poor quality and insufficiency of DSS data. In this thesis, we first investigate the meaning of elliptical anisotropy in the study of crustal structure and attribute, then derive Pg event’s traveltime-offset relationship by assuming a linear increasing velocity model with elliptical anisotropy and present the invert scheme from Pg traveltime-offset dataset to seismic velocity and its anisotropy of shallow crustal structure. We compare the Pg traveltime calculated by our analytic formula with numerical calculating method to test the accuracy. To get the lateral variation of elliptical anisotropy along the profiling, a tomography inversion method with the derived formula is presented, where the profile is divided into rectangles. Anisotropic imaging of crustal structure and attribute is efficient method for crust study. The imaging result can help us interprete the seismic data and discover the attribute of the rock to analyze the interaction between layers. Traveltime calculation is the base of image. Base on the ray tracing equations, the paper present a realization of three dimension of layer model with arbitrary anisotropic type and an example of Pg traveltime calculation in arbitrary anisotropic type is presented. The traveltime calculation method is complex and it only adapts to nonlinear inversion. Perturbation method of travel-time calculation in anisotropy is the linearization approach. It establishes the direct relation between seismic parameters and travetime and it is fit for inversion in anisotropic structural imaging. The thesis presents a P-wave imaging method of layer media for TTI. Southeastern China is an important part of the tectonic framework concerning the continental margin of eastern China and is commonly assumed to comprise the Yangtze block and the Cathaysia block, the two major tectonic units in the region. It’s a typical geological and geophysical zone. In this part, we fit the traveltime of Pg phase by the raytracing numerical method. But the method is not suitable here because the inefficiency of numerical method and the method itself. By the analytic method, we fit the Pg and Sg and get the lateral variation of elliptical anisotropy and then discuss its implication. The northeastern margin of Qinghai-Tibetan plateau is typical because it is the joint area of Eurasian plate and Indian plate and many strong earthquakes have occurred there in recent years．We use the Pg data to get elliptical anisotropic variation and discuss the possible meaning.

拉萨地块中新生代火山作用-从大洋到大陆俯冲转换的构造演化纪录

Late Mesozoic－Cenozoic volcanic rocks are well exposed in Lhasa Terrane, southern Tibet. This research attempts to apply 40Ar/39Ar geochronology, major, trace element and Sr-Nd-O isotopic geochemistry data to constrain the spatio-temporal variations, the composition of source, geodynamic setting. The results indicate that Lhasa Terrane mainly went through three tectonic-magmatic cycle: (1) Phase of Oceanic subduction (140-80Ma). Along with the subducting beneath the Eurasian Plate of Neo-Tethys slab, the oceanic sediment and/or the subducting slab released fluids/melts to metasomatize the subcontinental lithospheric mantle, and induced the mantle wedge partially melt and produced the calc-alkaline continental arc volcanic rocks; (2) Phase of continental-continental collision. Following the subducting of the Neo-Tethys slab, the Indian Plate collided with the Eurasian Plate dragged by the dense Neo-Tethys oceanic lithosphere. The oceanic lithosphere detached from continental lithosphere during roll-back and break-off and the asthenosphere upwelled. The resulting conducted thermal perturbation leads to the melting of the overriding mantle lithosphere and produced the syn-collisional magmatism: the Linzizong Formation and dykes; (3) Following by the detachment of the Tethys oceanic lithosphere, the Indian Lithosphere subducted northward by the drive from the expanding of Indian Ocean. The dense Indian continental lithospheric mantle (±the thickened lower crust) break off, disturb the asthenosphere, and lead to the melting of the overriding mantle lithosphere, which has been metasomatized by the melts/fluids from the subducting oceanic/continental lithosphere and the asthenosphere, and produced the rift-related ultrapotassic rocks.

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

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.

南海张裂大陆边缘数值模拟研究

The South China Sea (SCS) is one of the largest marginal seas in the western Pacific, which is located at the junction of Eurasian plate, Pacific plate and Indian-Australian plate. It was formed by continent breakup and sea-floor spreading in Cenozoic. The complicated interaction among the three major plates made tectonic movement complex and geological phenomena very rich in this area. The SCS is an ideal place to study the formation and evolution of rifted continental margin and sea-floor spreading since it is old enough to have experienced the major stages of the basin evolution but still young enough to have preserved its original nature. As the demand for energy grows day by day in our country, the deep water region of the northern continental margin in the SCS has become a focus of oil and gas exploration because of its huge hydrocarbon potential. Therefore, to study the rifted continental margin of the SCS not only can improve our understanding of the formation and evolution processes of rifted continental margin, but also can provide theoretical support for hydrocarbon exploration in rifted continental margin. This dissertation mainly includes five topics as follows: (1) Various classic lithosphere stretching models are reviewed, and the continuous non-uniform stretching model is modified to make it suitable for the case where the extension of lithopheric mantle exceeds that of the crust. Then simple/pure shear flexural cantilever model is applied to model the basement geometries of SO49-18 profile in the northern continental margin of the SCS. By fitting the basements obtained by using 2DMove software with modeling results, it is found that the reasonable effective elastic thickness is less than 5km in this region. According to this result, it is assumed that there is weak lower crust in the northern continental margin in the SCS. (2) We research on the methods for stretching factor estimation based on various lithosphere stretching models, and apply the method based on multiple finite rifting model to estimate the stretching factors of several wells and profiles in the northern continental margin of the SCS. (3) We improve one-dimension strain rate inversion method with conjugate gradient method, and apply it to invert the strain rate of several wells in the northern continental margin of the SCS. Two-dimension strain rate forward modeling is carried out, and the modeling results show that effective elastic thickness is a key parameter to control basin’s geometry. (4) We simulate divergent upwelling mantle flow model using finite difference method, and apply this newly developed model to examine the formation mechanism of the northwest and central sub-basin in the SCS. (5) We inverse plate thickness and basal temperature of oceanic lithosphere using sea-floor ages and bathymetries of the North Pacific and the North Atlantic based on varied-parameters plate model, in which the heat conductivity, heat capacity and coefficient of thermal expansion depend on temperature or depth. A new empirical formula is put forward based the inversed parameters, which depicts the relation among sea-floor age, bathymetry and heat flow. Then various similar empirical formulae, including the newly developed one, are applied to examine the sea-floor spread issue in the SCS based on the heat flow and bathymetry data of the abyssal sub-basin.

东喜马拉雅构造结新生代碰撞变形与岩石圈演化

The Eastern Himalayan Syntaxis (EHS) is one of the strongest deformation area along the Himalayan belt resulted from the collision between Indian plate and the Eurasian Plate since the 50～60Ma, and has sensitivity tracked and preserved the whole collisional processes. It should depend on the detail geological investigations to establish the deformational accommodate mode, and the uplift history, to elucidate the deep structure and the crust-mantle interaction of the Tibet Plateau of the EHS. The deep-seated (Main Mantle Thrusts) structures were exhumed in the EHS. The MMT juxtapose the Gangdese metamorphic basement and some relic of Gangdese mantle on the high Himalayan crystalline series. The Namjagbawa group which is 1200～1500Ma dated by U/Pb age of zircon and the Namla group which is 550Ma dated by U/Pb age of zircon is belong to High Himalayan crystalline series and Gangdese basement respectively. There is some ophiolitic relic along the MMT, such as metamorphic ocean mantle peridotite and metamorphic tholeiite of the upper part of ocean-crust. The metamorphic ocean mantle peridotites (spinel-orthopyroxene peridotite) show U type REE patterns. The ~(87)Sr/~(86)Sr ratios were, 0.709314～0.720788, and the ~(143)Nd/~(144)Nd ratios were 0.512073～0.512395, plotting in the forth quadrant on the ~(87)Sr/~(86)Sr-~(143)Nd/~(144)Nd isotope diagram. Some metamorphic basalt (garnet amphibolite) enclosures have been found in the HP garnet-kynite granulite. The garnet amphibolites can be divided two groups, the first group is deplete of LREE, and the second group is flat or rich LREE, and their ~(87)Sr/~(86)Sr, ~(143)Nd/~(144)Nd ratios were 0.70563～0.705381 and 0.512468～0.51263 respectively. Trace element and isotopic characteristics of the garnet amphibolites display that they formed in the E-MORB environment. Some phlogolite amphibole harzburgites, which exhibit extensive replacement by Phl, Amp, Tc and Dol etc, were exhumed along the MMT. The Phl-Amp harzburgites are rich in LREE and LILE, such as Rb, K etc, and depletes Eu (Eu~* = 0.36 ～ 0.68) and HFSE, such as Nb, Ta, Zr, Hf, P, Ti etc. The trace element indicate that the Phl-Amp harzburgites have island arc signature. Their ~(87)Sr/~(86)Sr are varied from 0.708912 to 0.879839, ~(143)Nd/~(144)Nd from 0.511993 to 0.512164, ε Nd from- 9.2 to - 12.6. Rb/Sr isochrone age of the phlogolite amphibole harzburgite shows the metasomatism took place at 41Ma, and the Amp ~(40)Ar/~(39)Ar cooling age indcate the Phl-Amp harzburgite raising at 16Ma. There is an intense crust shortening resulted from the thrust faults and folds in the Cayu block which is shortened more 120km than that of the Lasha block in 35～90Ma. With the NE corner of the India plate squash into the Gangdese arc, the sinistral Pai shear fault and the dextral Aniqiao shear fault on the both sides of the Great bent of Yalun Zangbu river come into active in 21～26Ma. On the other hand, the right-lateral Gongrigabu strike-slip faults come into activity at the same period, a lower age bound for the Gongrigabu strike-slip fault is estimated to be 23～24Ma from zircon of ion-probe U/Pb thermochronology. The Gongrigabu strike-slip faults connect with the Lhari strike-slip fault in the northwestern direction and with the Saganing strike-slip at the southeastern direction. Another important structure in the EHS is the Gangdese detachment fault system (GDS) which occurs between the sedimental cover and the metamorphic basement. The lower age of the GDS is to be 16Ma from the preliminary 40Ar/39Ar thermochronology of white mica. The GDS is thought to be related to the reverse of the subducted Indian crust and the fast uplift of the EHS. Structural and thermochronology investigation of the EHS suggest that the eastern Tibet and the western Yunnan rotated clockwise around the EHS in the period of 35～60Ma. Later, the large-scale strike-slip faults (RRD, Gaoligong and Saganing fault) prolongate into the EHS, and connect with the Guyu fault and Gongrigabu fault, which suggest that the Indianchia block escape along these faults. Two kind of magmatic rocks in the EHS have been investigated, one is the mantle-derived amphibole gabbro, dioposide diorite and amphibole diorite, another is crust origin biotit-garnet adamellite, biotit-garnet granodiorite and garnet-amphibole-biotite granite. The amphibole gabbro dioposite diorite and amphibole diorite are rich in LREE, and LILE, such as Ba, Rb, Th, K, Sr etc, depleted in HFSE, such as Nb, Ta, Zr, Hf, Ti etc. The ratio of ~(87)Sr/~(86)Sr are from 0.7044 to 0.7048, ~(143)Nd/~(144)Nd are from 0.5126 to 0.5127. The age of the mantle origin magamatic rocks, which result from the partial melt of the raising and decompression anthenosphere, is 8Ma by ~(40)Ar/~(39)Ar dating of amphibole from the diorite. The later crust origin biotite-garnet adamellite, biotite-garnet granodiorite and garnet-amphibole-biotite granite are characterized by aboudance in LREE, and strong depletion of Eu. The ratios of ~(87)Sr-~(86)Sr are from 0.795035 to 0.812028, ~(143)Nd/~(144)Nd from 0.51187 to 0.511901. The ~(40)Ar/~(39)Ar plateau age of the amphibole from the garnet-amphibole-biotite granite is 17.5±0.3Ma, and the isochrone age is 16.8±0.6Ma. Their geochemical characteristics show that the crust-derived magmatic rocks formed from partial melting of the lower curst in the post-collisional environment. A group of high-pressure kaynite-garnet granulites and enclave of high-pressure garnet-clinopyroxene grnulites and calc-silicate grnulites are outcroped along the MMT. The peak metamorphic condition of the high-pressure granulites yields T=800～960 ℃, P=1.4～1.8Gpa, corresponding the condition of 60km depth. The retrograde assemblages of the high-pressure grnulites occur at the condition of T=772.3～803.3 ℃, P=0.63～0.64Gpa. The age of the peak metamorphic assemblages are 45 ～ 69Ma indicated by the zircon U/Pb ion-plobe thermochronology, and the retrograde assemblage ages are 13～26Ma by U/Pb, ~(40)Ar/~(39)Ar thermochronology. The ITD paths of the high-pressure granulites show that they were generated during the tectonic thickening and more rapid tectonic exhumation caused by the subducting of the Indian plate and subsequent break-off of the subducted slab. A great deal of apatite, zircon and sphene fission-track ages, isotopic thermochronology of the rocks in the EHS show that its rapid raising processes of the EHS can be divided into three main periods. There are 35～60Ma, 13～25Ma, 0～3Ma. 3Ma is a turn in the course of raising in the EHS which is characterized by abruptly acceleration of uplifting. The uplift ratios are lower than 1mm .a~(-1) before 3Ma, and higher than 1mm .a~(-1) with a maximum ratio of 30mm .a~(-1) since 3Ma. The bottom (knick point) of the partial anneal belt is 3.8km above sea level in the EHS, and correspond to age of 3Ma determined by fission-track age of apatite. The average uplift ratio is about 1.4 mm .a~(-1) below the knick point. The EHS has raised 4.3km from the surface of 2.36km above sea level since 3Ma estimated by the fossil partial anneal belt of the EHS. We propose a two-stage subduction model (B+A model) basing on Structural, thermochronological, magmatical, metamorphic and geophysical investigations of the EHS. The first stage is the subduction of the Indian continental margin following after the subduction of the Tethys Ocean crust and subsequent collision with the Gangdese arc, and the second stage is the Indian crust injecting into the lower crust and upper mantle of the Tibet plateau. Slab break-off seems to be occurred between these two stages.