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

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.

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

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

黄骅坳陷高尘头地区储层非均质性及成藏研究

Gaochentou region is located in the southwest direction of Gaochentou village in Huanghua city of Hebei province. In regionally structural position, It lies in Qikou sag In the middle part of Huanghua depression, which belongs to the east part of the south Dagang structure zone in the middle part of Huanghua depression. Its' very beneficial at regional structure in Gaochentou , and It becomes the advantage area for oil and gas gathered and preserved, Sandstone reservoir of Dongying Formation is main bearing bed .Dongying Formation in Gaochentou region of Huanghua depression is consisted of set of mudstone and sandstone interbeds by deposited delta fades . Dongying Formation can be divided into 3 members from above to below: the first member of Dongying Formation (FMDF), the second member of Dongying Formation (SMDF), and third member of Dongying Formation (TMDF). The lithology of the upper part of FMDF was consisted of mostly middle-grained and fine-grained sandstone, and it is small for the oil-bearing area of the sand bodies .The lithology of the lower part is coarse-grained sandstone bodies which are well connected between sandstone bodies of wells, and the lower part was main bed of oil production in Dongying Formation; SMDF and TMDF are consisted of larger scale set of mudstone, in which the sandbodies are lenticular and pinch out quickly, and the lithology was mostly fine sandstone and silt stone, in which there are little oil and gas .Because the reservoirs in this area are largely influenced by the factors such as lithology, fault and others, and the reservoirs have the strong,heterogeneity , there exists the problem of oil-down and water-up for vertical distribution of oil and gas bearing. It is not very clearly for the three dimension distribution of sandstone , and the geology researchs is not enough. So, it can't satisfy the need of further development and production for Gaochentou oilfield.Having the key problem of oil-down and water-up and the mechanism of the reservoir for Gaochentou area, There are as follow study works, the first, is study of the high-resolution correlation of sequence stratigraphy and sedimentary microfacies. Dongying Formation was divided into three parasequence sets and each parasequence set was divided into different amount of parasequences. FMDF, as the main oil and gas producing bed, can be divided into seven parasequences. Oil and gas are discovered in six parasequences except the seventh. On the basis of study of sedimentary microfacies, the sediments of Dongying Formation are considered deposited mainly in delta front subfacies. The microfacies types of Dongying Formation are sub-water distirbutary channel, sub-water natural bank, inter distributary channel bay, distributary channel mouth dam, and delta front mat sand.Seismic facies analysis and logging-constrained inversion technique were applied by Author for transverse prediction of sandstone reservoir. Having 4 modes of interwell single sandbodies correlation technique, Author have described distribution characteristics of sandbodies, and established geological reservoir model of Gaochentou reservoir.Author presented that the reservoirs characteristic have very strong heterogeneity ,and In the section of sandstone interlayed with mudstone,the folium sandstone interlayed with each other, and the wedge shaped sandbodies pinched out in the mudstone. So the pinch-out up sandstone trap and lenticular sandstone trap are easily formed. They are most small scale overlying pinches out in the place of slope. This article applies the concept of deep basin oil to resolve reasonably the problem of which the oil is below the water in Gaochentou area. Combined with the study of sedimentary facies, reservoir and other aspects, the mechanism and patterns of deep basin oil are studied on the basis of characteristics in Gaochentou area.On the basis of the above study, the mechanism of the oil and gas' migration and accumulation in isotropic sandstone and heterogeneous sandstone are thoroughly analyzed through experiments on physical modeling. Experiments on physical modeling show that the discrepancy between sand layers with different permeability and thickness has important influence on the direction, path, and injection layer of oil's migration. At the beginning of the injection of oil and gas in high permeability sand layer, the pressure is low, the migration resistance is small, and the oil and gas are more easily displacing the water in sand. So it can act as good transformation layer or reservoir. But at the beginning of the injection of oil and gas in sand layer with low permeability, the pressure is high, the migration resistance is big, and the oil and gas are more difficultly displacing the water in sand. So it can only act as bad or worse transformation layer or reservoir. Even if it cannot act as transformation layer or reservoir, it can act as water layer or dry layer. The discrepancy between sand layers on permeability and thickness can make discrepancy in injection of oil and gas between different layers. Consequently it leads to small amount of oil and gas injection in sand layers with low permeability. Ultimately it affects the oil's accumulation and distribution in different sand layers.At Last, combining analysis of the structure and pool forming condition, The thesis has established models of reservoir formation to predict the advantage distribution of oil and gas bearing , and put forward the prospective target It is not only of theoretical signification for explosion and importance, but also has realistic value in guiding the progressive petroleum exploration and exploitation.

利用宽频带台阵资料研究华北地区地壳及上地幔间断面结构

Intense tectonic renovation has occurred in the eastern continent of china since Mesozoic, as evidenced by the high heat flow, widespread magma extrusion and volcanic activities, and development of large sedimentary basins. To explain the cause and mechanism for the tectonic process in this period, some researchers have put forward various models, such as mantle plume, subduction of the Pacific slab, Yangtze Block-North China Block collision, etc. Their seismological evidence, however, is still scarce..During the period from 2000 to 2003, large temporary seismic arrays were established in North China by the Institute of the Geology and Geophysics, Chinese Academy of Sciences. Total 129 portable seismic stations were linearly emplaced across the western and eastern boundaries of the Bohai Bay Basin, and accumulated a large amount of high-quality data. Moreover, abundant data were also collected at the capital digital seismic network established in the ninth five-year period of national economic and social development. These provide an unprecedented opportunity for us to study the deep structure and associated geodynamic mechanism of lithospheric processes in North China using seismological techniques.Seismology is a kind of observation-based science. The development of seismic observations greatly promotes the improvement of seismologic theory and methodology. At the beginning of this thesis, I review the history of seismic observation progress, and present some routine processing techniques used in the array seismology. I also introduce two popular seismic imaging methods (receiver function method and seismic tomography).Receiver function method has been widely used to study the crustal and upper mantle structures, and many relevant research results have been published. In this thesis I elaborate the theory of this method, including the basic concept of receiver functions and the methodology for data pre-processing, stacking and migration. I also address some problems often encountered in practical applications of receiver function imaging.By using the teleseismic data collected at the temporary seismic arrays in North China, in particular, the traveltime information of P-to-S conversion and multiple reverberations of the Moho discontinuity, I obtain the distributions of the crustal thickness and the poisson ratio at the northwest boundary area of the Bohai Bay Basin and discuss the geological implications of the results.Through detailed intestigations on the crustal structural feature around the middle part of the Tanlu fault, considerable disparity in poisson ratios is found in the western and eastern sides of the Tanlu fault. Moreover, an obvious Moho offset is coincidently observed at the same surface location. A reasonable density model for the Tanlu fault area is also derived by simulating the observed gravity variations. Both receiver function study and gravity anomaly modeling suggest that the crustal difference between the western and eastern sides of the Tanlu fault is mainly resulted from their different compositions.With common conversion point imaging of receiver functions, I estimate the depths of the upper and lower boundaries of the mantle transition zone, i.e., the 410 and 660 km discontinuities, beneath most part of the North China continent The thickness of the transition zone (TTZ) in the study area is calculated by subtracting the depth of .410 km discontinuity from that of the 660km discontinuity. The resultant TTZ is 10-15 km larger in the east than in the west of the study area. Phase transitions at the 410 km and the 660 km discontinuities are known to have different Clapeyron slopes. Therefore, the TTZ is sensitive to the temperature changes in the transition zone. Previous studies have shown that the TTZ would be smaller in the mantle plume areas and become larger when the remnants of subducted slabs are present The hypothesis of mantle plume cannot give a reasonable interpretation to the observed TTZ beneath North China, Instead, the receiver function imaging results favor a dynamic model that correlates the thermal structure of the mantle transition zone and associated upper mantle dynamics of North China to the Pacific plate subduction process.