建立在小型计算机上的GPSS-F/1000通用仿真软件包

本文阐述了在小型计算机上建立的 GPSS-F/1000通用仿真软件包。该软件包不仅具有原 GPSS-F仿真软件的全部功能,并且使公用数组、系统参数的适配全部活化,通过子程序的形式参量传送赋给,从而,仿真时不必修改子程序的源程序。此外,增加了统计参数估计子程序集和图形绘制子程序集。

基于显式分形插值的高保真地震数据重建方法研究

Based on the fractal theories, contractive mapping principles as well as the fixed point theory, by means of affine transform, this dissertation develops a novel Explicit Fractal Interpolation Function（EFIF）which can be used to reconstruct the seismic data with high fidelity and precision. Spatial trace interpolation is one of the important issues in seismic data processing. Under the ideal circumstances, seismic data should be sampled with a uniform spatial coverage. However, practical constraints such as the complex surface conditions indicate that the sampling density may be sparse or for other reasons some traces may be lost. The wide spacing between receivers can result in sparse sampling along traverse lines, thus result in a spatial aliasing of short-wavelength features. Hence, the method of interpolation is of very importance. It not only needs to make the amplitude information obvious but the phase information, especially that of the point that the phase changes acutely. Many people put forward several interpolation methods, yet this dissertation focuses attention on a special class of fractal interpolation function, referred to as explicit fractal interpolation function to improve the accuracy of the interpolation reconstruction and to make the local information obvious. The traditional fractal interpolation method mainly based on the randomly Fractional Brown Motion (FBM) model, furthermore, the vertical scaling factor which plays a critical role in the implementation of fractal interpolation is assigned the same value during the whole interpolating process, so it can not make the local information obvious. In addition, the maximal defect of the traditional fractal interpolation method is that it cannot obtain the function values on each interpolating nodes, thereby it cannot analyze the node error quantitatively and cannot evaluate the feasibility of this method. Detailed discussions about the applications of fractal interpolation in seismology have not been given by the pioneers, let alone the interpolating processing of the single trace seismogram. On the basis of the previous work and fractal theory this dissertation discusses the fractal interpolation thoroughly and the stability of this special kind of interpolating function is discussed, at the same time the explicit presentation of the vertical scaling factor which controls the precision of the interpolation has been proposed. This novel method develops the traditional fractal interpolation method and converts the fractal interpolation with random algorithms into the interpolation with determined algorithms. The data structure of binary tree method has been applied during the process of interpolation, and it avoids the process of iteration that is inevitable in traditional fractal interpolation and improves the computation efficiency. To illustrate the validity of the novel method, this dissertation develops several theoretical models and synthesizes the common shot gathers and seismograms and reconstructs the traces that were erased from the initial section using the explicit fractal interpolation method. In order to compare the differences between the theoretical traces that were erased in the initial section and the resulting traces after reconstruction on waveform and amplitudes quantitatively, each missing traces are reconstructed and the residuals are analyzed. The numerical experiments demonstrate that the novel fractal interpolation method is not only applicable to reconstruct the seismograms with small offset but to the seismograms with large offset. The seismograms reconstructed by explicit fractal interpolation method resemble the original ones well. The waveform of the missing traces could be estimated very well and also the amplitudes of the interpolated traces are a good approximation of the original ones. The high precision and computational efficiency of the explicit fractal interpolation make it a useful tool to reconstruct the seismic data; it can not only make the local information obvious but preserve the overall characteristics of the object investigated. To illustrate the influence of the explicit fractal interpolation method to the accuracy of the imaging of the structure in the earth’s interior, this dissertation applies the method mentioned above to the reverse-time migration. The imaging sections obtained by using the fractal interpolated reflected data resemble the original ones very well. The numerical experiments demonstrate that even with the sparse sampling we can still obtain the high accurate imaging of the earth’s interior’s structure by means of the explicit fractal interpolation method. So we can obtain the imaging results of the earth’s interior with fine quality by using relatively small number of seismic stations. With the fractal interpolation method we will improve the efficiency and the accuracy of the reverse-time migration under economic conditions. To verify the application effect to real data of the method presented in this paper, we tested the method by using the real data provided by the Broadband Seismic Array Laboratory, IGGCAS. The results demonstrate that the accuracy of explicit fractal interpolation is still very high even with the real data with large epicenter and large offset. The amplitudes and the phase of the reconstructed station data resemble the original ones that were erased in the initial section very well. Altogether, the novel fractal interpolation function provides a new and useful tool to reconstruct the seismic data with high precision and efficiency, and presents an alternative to image the deep structure of the earth accurately.

湖南、广西泥盆纪硅岩地球化学与沉积构造背景

During the Devonian, a complicated carbonate platform-basin configuration was created through transtensional rifting in the context of opening of Devonian South China Sea; extensive bedded chert, commonly interbedded with tuffaceous beds, occurred in the narrow, elongate interplatform basins (or troughs) in South China, where they occurred earlier (Early Devonian) in southern Guangxi and later (early Late Devonian) in northern Guangxi-south central Hunan. In order to unravel the origin and distribution of the bedded chert successions, and their relationships to basement faulting activities during the opening of the Devonian South China Sea, studies of element (major, minor and REE) geochemistry and Rb-Sr, Sm-Nd isotopic systematics are carried out upon the chert deposits. These chert deposits commonly have high SiO2 contents and (average 94.01%) and low TFe2O3 (average 0.55%), together with other geochemical parameters, suggestive of both biogenic and hydrothermal origins. However, Fe/Ti ratio are high along the elongate interplatform basins(troughs) to the northwest along Wuxiangling-Zhaisha-Chengbu, and to the southeast along Xiaodong-Mugui-Xinpu, suggesting relatively intense hydrothermal activities there. They generally contain very low total REE contents (∑REE average 31.21ug/g) with mediate negative Ce anomalies (mean Ce/Ce*=0.83) and low Lan/Cen values (average 1.64), indicating an overall continental margin basin where they precipitated. The northward increases in Ce/Ce* values, particularly along the elongate troughs bounded both to the east and west of the Guangxi-Huanan rift basin, suggest a northward enhancement of terrigenous influences, thereby reflecting a gradual northward propagation of open marine setting. Generally low positive Eu anomalies in the chert, except for the apparently high Eu anomalies in the chert from Chengbu (Eu/Eu* up to 4.6), suggest mild hydrothermal venting activities in general, except for those at Chengbu. The initial 87Sr/86Sr (0) ratios of chert generally vary from 0.712000 to 0.73000 , suggesting influences both from terrigenous influx and seawater. The Nd isotopic model ages (tDM or t2DM) and initial εNd (0) values of chert vary mostly from 1.5 to 2.1 Ga, and from –16 to –21, respectively, implying that the silica sources were derived from the provenances of the Palaeoproterozoic crust relics at depth. The high εNd (0) values of chert (-0.22 to 14.7) in some localities, mostly along the elongate troughs, suggest that silica sources may have been derived from deeper-seated mantle, being channeled through the interplate boundary fault zones extending downwards to the mantle. At Wuxiangling, Nanning, chert occurs extensively from the Emsian through the Frasnian strata, both U/Th ratios and tDM ages of chert reached up to a maximum in the early Frasnian corresponding to the extensive development of chert in South China, pointing to a maximum extensional stage of Devonian South China basin, which is supported by the Ce/Ce* values as is opposed to the previous datasets as the coeval minimum values.

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

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 Tarim Basin, extensive carbonates of Lower Paleozoic occur, in which thick Cambrian and Lower Ordovician dolostones are widespread and show a potential perspective in hydrocarbon exploration. So they are viewed as an important target for exploration. Tarim Basin is a poly phase composite basin, which underwent multiphase tectonic modification and volcanic activities; these exerted significant influences on the basin-fills and basin fluid evolution, thereby the diagenetic history, particularly on the deep-buried Lower Paleozoic dolostones. Referring to the classification of dolomite texture proposed by Gregg & Sibley (1984) and Sibley & Gregg (1987). In view of crystal size, crystal shape, crystal surface and contact relation, eight genetic textures of dolomite crystals are identified, based on careful petrographic examinatoins. These textures include: 1) micritic dolomite; 2) relict mimetic dolomite; 3)finely crystalline, planar-e(s), floating dolomite; 4)finely crystalline, planar-e(s) dolomite; 5) finely-coarse crystalline, nonplanar-a dolomite; 6)coarse crystalline, nonplanar saddle dolomite; 7) finely-medium crystalline, planar-e(s) dolomite cement; 8) coarse crystalline, nonplanar saddle dolomite cement, in which the former six textures occurs as in matrix, the latter two in the cements. Detailed geochemistry analysis is carried out on the basis of genetic textures of dolomite and related minerals such as quartz and calcite. The result showed that the calcite has the highest average content in Sr, which can be sorted into two groups; micritic dolomite has the highest average content in Sr among all kinds of dolomites; the REE patterns of all kinds of dolomites is similar to those of marine limestone samples. Saddle dolomite cement has δ13C values from -2.44‰ to 1.27‰ PDB, and δ18O values from -13.01‰ to -5.12‰ PDB, which partially overlap with those of matrix dolomite (δ13C values from -2.83‰ to 2.01‰ PDB, δ18O values from -10.63‰ to -0.85‰ PDB). Saddle dolomite cement has 87Sr/86Sr ratios from 0.7086 to 0.7104, which totally overlap with those of matrix dolomite (0.7084 ~ 0.7116). Compared with saddle dolomite derived from other basins all over the world, the saddle dolomites of Tarim Basin have similar δ13C, δ18O and 87Sr/86Sr ratios values with those of matrix dolomite. This scenario reflects the unusual geological setting and special dolomitizing liquid of Tarim Basin. The values of δ18O, δ13C and 87Sr/86Sr ratios of calcite also can be sorted out two groups, which may been resulted from the one stage of extensive uplift of Tarim Basin from Mesozoic to Cenozoic. Fluid inclusion microthermometry data of the diagenetic mineral indicates that matrix dolomite has relatively low homogenization temperatures (Th) of 80～105oC and salinities of 12.3％ (wt% NaCl equivalent); saddle dolomite has highest Th values, which concentrate in 120～160oC and salinities of 13.5～23.7％ (wt% NaCl equivalent); quartz has relatively low Th of 135～155oC and salinities of 17.8～22.5％ (wt% NaCl equivalent); calcite has relatively low Th of 121～159.5oC and salinities of 1.4～17.5％ (wt% NaCl equivalent). These data suggest that the saddle dolomites could have formed in thermal brine fluids. Based on comprehensive petrographical study, detailed geochemistry and fluid inclusion microthermometry analysis on Lower Paleozoic dolomite of Tarim Basin, three types of dolomitisation mechanism are proposed: Penecontemporaneous dolomitisation (Sabkha dolomitisation & Reflux dolomitisation); Burial dolomitisation (shallow-intermediate burial dolomitisation & Deep burial dolomitisation ); Hydrothermal cannibalized dolomitisation. In view of host-specified occurrences of hydrothermal dolomite, the low abundance of saddle dolomite and high geochemical similarities between saddle dolomite and host dolomite, as well as highest Th and high salinities , the hydrothermal dolomite in Tarim Basin is thus unique, which could have been precipitated in modified fluid in the host dolomite through intraformational thermal fluid cannibalization of Mg ions from the host. This scenario is different from the cases that large scale dolomitizing fluid migration took place along the fluid pathways where abundant saddle dolomite precipitated. Detailed observations on 180 petrographic and 60 casting thin sections show original pores in Lower Paleozoic dolomite were almost died out by complicated diagenetic process after a long time geologic evolution. On the other hand, deep-buried dolomite reservoirs is formed by tectonic and hydrothermal reforming on initial dolomites. Therefore, the distribution of structure-controlled hydrothermal dolomite reservoirs is predicted in Tabei and Tazhong Area of Tarim Basin based on the geophysical data.

东准噶尔和尔赛斑岩铜矿特征及其形成的构造背景

Hersai porphyry copper deposit(PCD) of eastern junggar, newly discovered copper deposit, is located at the eastern segment of the Xiemisitai-Kulankazigan-Zhifang-Qiongheba Paleozoic island arc, Eastern Junggar. The Hersai PCD is developed in a intrusive complex, characterized by intensive and multiform hydrothermal alteration, including potassic alteration, silification, chloritization，sericitization，kaolinitization and carbonatization. Granodiorite, grandiorite porphyry, granite and concealed explosion breccia are hosts of the ore bodies containing veinlet and disseminated ore. Ore-bearing granite (ZK107-1-9), granodiorite (ZK107-1-9) and Ore-barren granodiorite (HES2-1) are selected to date zircon U-Pb age by SHRIMP method, and have an age of 429.4±6.4Ma ,413.0±3.4Ma and 411.1±4.8Ma, respectively, showing that they were emplaced from Late Silurian to Early Devonian. In addition, sample ZK107-1-9 has some hydrothermal zircons with a weighted mean 206Pb/238U age of 404.9±3.7Ma which is interpreted to be related to the granodiorite porphyry. Re-Os dating of five molybdenite samples yielded a weighted average model age of 408.0±2.9Ma, indicating the metallogenic epoch of the Hersai PCD. The ore-forming age is close to the petrogenic time of garnodiorite (411-413Ma), this suggests the ore-forming porphyry is most possiblely granodiorite porphyry. Systematic major - trace elements and Rb-Sr-Sm-Nd-Pb-Hf isotopic characteristics were studied. Analysis results show that these intrusives have some interesting and special characteristics, as following:1) containing both calc-alkaline rocks and high potassium calc-alkaline rocks ; 2) have some characteristics of adakite, but not totally, such as much lower La/Yb ratios and no Eu anomaly or just faint Eu anomaly; 3) have an initial 87Sr/86Sr ratios(0.703852-0.704565) similar to that of BSE, positive εNd(t) values between 6.1 and 7.4, the initial 206Pb/204Pb values (17.576-17.912), 207Pb/204Pb values (15.400-15.453) , 208Pb/204Pb values (37.252-37.466) , and high εHf(t) values (10.2-15.4) close to the value of depleted mantle. These geochemical features suggest that these igneous rocks in the Hersai area not only have some characteristics of island arc, but also some characteristics that only appear in the continental margin arc. It is suggested that Hersai PCD is formed in the subduction setting by the partial melting of young crust. These works and advancements mentioned in the paper are helpful to understand the deposit geology, geochemistry and metallogenesis of Hersai PCD. It is also significant to understand mineralization and tectonic setting in the Qiongheba area.

阿尔泰南缘侧向增生过程中的VMS矿床：地质、地球化学特征与成矿机制

The Chinese Altai is one of the most important volcanogenic massive sulfide (VMS) deposit districts in China. All orebodies were lenticular or bedded and stratabounded by a suite of early Devonian volcanic-sedimentary rocks. Hydrothermal feeder zones developed under some of the orebodies. All the ores are massive or laminated, and show typical characteristics of VMS deposit. Based on the mineralizing time and the metal assembles, we divide 3 metallogenic stages: 1, Fe orefroming stage associated with basaltic and sedimentary rocks during very early Devonian; 2, Cu-Pb-Zn oreforming stage associated with rhyolitic and sedimentary rocks during early Devonian; 3， Cu-Zn oreforming stage in the dacitic and basaltic rocks during mid. Devonian. The hosting rocks for all orebodies are different, but they show very similar geochemical and isotopic characteristics. All the felsic rocks show enriched lighted rare earth elements (REE) patterns (La/Yb>5), and with an obvious Eu negative anomalies (Eu/Eu*<0.6). In the meanwhile, all the mafic rocks show flat REE pattern and no Eu anomalies. The Ashele basalt show an apparent Ce negative anomalies (Ce/Ce* <0.76), All the volcanic roks in Chinese Altai show the decoupled property between the high field strength elements (HFSE) and large ion lithophile elements (LILE). The negative Nb, Ta characteristics with respect to adjacent elements indicate that subduction-modified source. The Nd(t) of the hosting rocks for all orebodies changed in a small range (-1.5~5), and the (87Sr/86Sr)i change in a big range. The initial Sr value of the hosting rocks in Mengku and Tiemuerte are obviously affected by the seawater (0.705~0.710), and initial Sr values of hosting rocks Ashele change in a small range (0.704~0.706). All Sr-Nd isotopes of ores have the same range with the hosting rocks, indicating that both the ores and volcanic rocks have the same island arc source. The mean sulfur isotopes of sulfides from Ashele and Mengku are 6.2‰ and 3.4‰, respectively, indicating a deep magmatic source. However, the sulfur isotopes of sulfides from Keketale, Tiemuerte and Keyinbulake changed in -15.8‰~9.9‰, -23.5‰~1.87‰, -8.3‰~1.6‰, respectively. And the big sulfur isotope range indicated that the sulfur of the ores was a combination biogenic and magmatic source. All volcanic rocks from the VMS deposits in the southern Chinese Altai show a typical subduction related environments. Based on the regional and locally geological evidence, here we propose that the southern Chinese Altai is an island arc system, and all VMS deposits formed during the lateral accretion process. No VMS deposit formed during the formation of the island arc during Silurian; Fe VMS deposit formed during the beginning of the opening of the backarc basin in very early Devonian; Cu-Pb-Zn VMS deposits formed during the mature stage of the backarc basin in early Devonian; at last the Cu-Zn VMS deposit formed during the rifted stage of the island arc itself.

湘西地区上震旦统上部－下寒武统底部黑色岩系沉积学及地球化学

Widespread black chert-shales occur in the Ediacaran-Cambrian(E-C) boundary successions along the flank of Yangtze Platform, South China, remarkable changes in sedimentology, geochemistry and biology were recorded. Although extensive studies were carried out upon this boundary succession, the origin of black chert-shales still remain controversial. This paper focuses on the E-C black chert-shales in western Hunan, South China, upon which detailed depositional and geochemical changes are documented, accordingly a depositional model for black chert-shales is proposed. Stratigraphic anatomy across the depositional strike demonstrates that the shallow-water Dengying dolostone along the platform margin sharply pass basinward into the Liuchapo chert successions, which indicate syndepositional extensional faulting at depth could have occurred along the platform margin. The deep-water Niutitang phosphorite-rich black shales are either underlain by the Dengying dolostones on the platform margin toward platform interior or directly by the Liuchaopo chert successions farther basinwards. By detailed investigation, silica chimneys are firsly identified approximately in the chert along platform margin; two types of silica chimneys, including mounded and splayed/funnelized chert(generally brecciated) bodies are further sorted out. The mounded chert are exitbited by domed or hummocky surfaces on the top and irregular spongy to digitiform internal fabrics; within the silica mounds, abundant original vesicles/voids and/or channels were mostly plugged by initial chalcedony, quartze crystals with minor dolomite and bladed barite crystals. Splayed/funnelized brecciated chert “intrusion” cross-cut the uppermost dolostones capping to the horizon underneath, and are directly overlain by the Niutitang phosphorite-rich black shales. Their similarities to the silica chimneys reported from the oceanic spreading centres suggest a similar origin responsible for these unique silica bodies which is also supported by the microthermonmetric data and element geochemistry. High P, Ba, Fe contents and positive correlation between Fe and TOC concentrations in the Niutitang black shales indicate a high palaeo-productivity in the Early Cambrian ocean. The low Th/U and the high V/Cr, V/Sc, V/(V+Ni) ratios in the black shales suggest an anoxic water condition during this interval. Furthermore, Positive Eu anomalies and high Ba contents in the sediments also imply a hydrothermal influence on the formation of Niutitang black shales. To better constrain the placement of deep-water successions straddling the E-C boundary and the timing of hydrothermal silica chimneys, sensitive high-resoluton ion microprobe(SHRIMP) U-Pb dating of zircon grains from tuffs within the chert succession of Liuchapo Formation at Ganziping was conducted and yields a weighted-mean 206Pb/238Pb age at 536.6±5.5Ma, younger than E-C boundary age(542.0±0.3Ma). This age combined with carbon isotopic data is then proposed to correspond to the U-Pb age of zircons(538.2±1.5Ma) from the Zhongyicun member of Meishucun Formation at Meishucun in eastern Yunna, thus, the E-C boundary in Gazngziping was placed between the Dengying formations and Liuchapo formatioms. therefore, the silica chimneys took place at the beginning of the Cambrian period. The temporal coincidence of silica chimneys and negative excursions of δ13C and δ34Spy pairs suggest hydrothermal activities were likely responsible for the isotopic changes. Under such a circumstance, vast amounts of greenhouse gases(CO2, CH4, H2S), with highly 13C-depleted carbon and 34S-depleted sulfur would be released into the ocean and atmosphere. A positive shift in δ34Scas and Δ34S values from the late Ediacaran to the Early Cambrian could be a reflection of enhanced bacterial sulfate reduction(BSR), strengthened by the intensified oceanic anoxia stimulated by hydrothermal activities. Based on the analyses of sedimentology and geochemistry, a model- “oceanic anoxia induced by hydrothermal–volcanic activies” was proposed to responsible for the formation of black chert-shales during this E-C transition. Under this case, hydrothermal-volcanic activies could release large large amount of greenhouse into atmosphere and metal micronutrients into the ocean, which may lead to global warming, stratified ocean, thereby a high palaeoproductivity; on the other hand, the massive releasing of reduced hydrothermal fluids with abundant H2S, could have in turn enhanced the ocean anoxia. All of these were favourable the for preservation of organic matter, and subsequent extensive deposition of black silica-shales.

特提斯喜马拉雅带晚白垩世-始新世沉积源区分析及其构造意义

Tethyan Himalayan Sequence (THS) is located at the frontier of the India-Asia collision zone, which can preserve critical information about collision. This paper reports detailed petrology, geochemistry, spinels electron microprobe data, and in situ U-Pb ages and Lu-Hf isotopic data on detrital zircons from the late Cretaceous to early Eocene strata in Gyantze and Gamba area, south Tibet that provide important constraints on the early tectonic evolution of the India-Asia collision. In Gyantze, the lithic arkose in Zongzhuo mélange is characterized by, SiO2 =80.4%, Al2O3=8.6%, Na2O=1.6%, K2O=1.1%, LaN/YbN=8.90, and εNd (0) =-10.27. Spinels compositions are characterized by low TiO2 (generally <0.1%) and a Cr number mainly between 70 and 80. The largest population of detrital zircons is within the 73-169Ma range with high εHf (t) and > 500 Ma with complex εHf (t) values. The lithic arkose in Rilang conglomerate is characterized by, SiO2 =56.5%, Al2O3=15.6%, Na2O=4.7%, K2O=0.6%, LaN/YbN=5.00-5.29, and εNd (0) =1.92. Spinels of 2006T98 display high TiO2 (generally >0.2%) and a Cr number mainly between 70 and 85, other spinels are characterized by low TiO2 (generally <0.2%) and a Cr number mainly between 60 and 90. The largest population of detrital zircons is within 90-146 Ma range with high εHf (t). The lithic arkose in Jiachala formation is characterized by, SiO2 =64.6%, Al2O3=12.1%, Na2O=1.9%, K2O=1.8%, LaN/YbN=7.73-9.13, and εNd (0) =-5.52~-8.43. Spinels in the Jiachala formation have low TiO2 (generally <0.2%) and a Cr number between 39 and 88. Detrital zircons have a wide range of age distribution of 82-3165Ma with complex εHf (t). In Gamba, The quartze sandstone in Jidula formation is characterized by, SiO2=97.4%, Al2O3=0.9%, Na2O=0.03%, K2O=0.18%, LaN/YbN=18.70-21.684, and εNd (0) between -13.1~-7.4. While the lithic arkose in Zhepure formation is characterized by, SiO2=68.4%, Al2O3=7.3%, Na2O=1.15%, K2O=0.52%, LaN/YbN=6.09-8.99, and εNd(0)=-5.8~-6.3. Based on our geochemical analysis, spinles electron microprobe data, U–Pb ages and Hf isotope data for detrital zircons of the late Cretaceous-Eocene strata in Gyantze and Gamba, southern Tibet, the following major conclusions can be drawn: 1. In Gyantze, the Zongzhuo mélange was mainly derived from accretionary prism/THS of continental slop and Gangdese arc. Rilang conglomerate was totally from Gangdese arc. The Jiachala formation was derived from THS, suture zone and Gangdese arc. 2. In Gamba, the Jidula formation was from India craton, while the Zhepure formation was derived from THS, suture zone and Gangdese arc. 3. The deposite of Zongzhuo mélange and Rilang conglomerate (73-55Ma) marks the collision between India and Asia. 4. Late Paleocene-Eocene tectonic evolution is consistent with foreland basin system.

地震波场延拓+菲涅尔体叠加方法及成像效果

In the increasingly enlarged exploration target, deep target layer(especially for the reservoir of lava) is a potential exploration area. As well known, the reflective energy becomes weak because the seismic signals of reflection in deep layer are absorbed and attenuate by upper layer. Caustics and multi-values traveltime in wavefield are aroused by the complexity of stratum. The ratio of signal to noise is not high and the fold numbers are finite(no more than 30). All the factors above affect the validity of conventional processing methods. So the high S/N section of stack can't always be got with the conventional stack methods even if the prestack depth migration is used. So it is inevitable to develop another kind of stack method instead. In the last a few years, the differential solution of wave equation was hold up by the condition of computation. Kirchhoff integral method rose in the initial stages of the ninetieth decade of last century. But there exist severe problems in it, which is are too difficult to resolve, so new method of stack is required for the oil and gas exploration. It is natural to think about upgrading the traditionally physic base of seismic exploration methods and improving those widely used techniques of stack. On the other hand, great progress is depended on the improvement in the wave differential equation prestack depth migration. The algorithm of wavefield continuation in it is utilized. In combination with the wavefield extrapolation and the Fresnel zone stack, new stack method is carried out It is well known that the seismic wavefield observed on surface comes from Fresnel zone physically, and doesn't comes from the same reflection points only. As to the more complex reflection in deep layer, it is difficult to describe the relationship between the reflective interface and the travel time. Extrapolation is used to eliminate caustic and simplify the expression of travel time. So the image quality is enhanced by Fresnel zone stack in target. Based on wave equation, high-frequency ray solution and its character are given to clarify theoretical foundation of the method. The hyperbolic and parabolic travel time of the reflection in layer media are presented in expression of matrix with paraxial ray theory. Because the reflective wave field mainly comes from the Fresnel Zone, thereby the conception of Fresnel Zone is explained. The matrix expression of Fresnel zone and projected Fresnel zone are given in sequence. With geometrical optics, the relationship between object point in model and image point in image space is built for the complex subsurface. The travel time formula of reflective point in the nonuniform media is deduced. Also the formula of reflective segment of zero-offset and nonzero offset section is provided. For convenient application, the interface model of subsurface and curve surface derived from conventional stacks DMO stack and prestack depth migration are analyzed, and the problem of these methods was pointed out in aspects of using data. Arc was put forward to describe the subsurface, thereby the amount of data to stack enlarged in Fresnel Zone. Based on the formula of hyperbolic travel time, the steps of implementation and the flow of Fresnel Zone stack were provided. The computation of three model data shows that the method of Fresnel Zone stack can enhance the signal energy and the ratio of signal to noise effectively. Practical data in Xui Jia Wei Zhi, a area in Daqing oilfield, was processed with this method. The processing results showed that the ability in increasing S/N ratio and enhancing the continuity of weak events as well as confirming the deep configuration of volcanic reservoir is better than others. In deeper target layer, there exists caustic caused by the complex media overburden and the great variation of velocity. Travel time of reflection can't be exactly described by the formula of travel time. Extrapolation is bring forward to resolve the questions above. With the combination of the phase operator and differential operator, extrapolating operator adaptable to the variation of lateral velocity is provided. With this method, seismic records were extrapolated from surface to any different deptlis below. Wave aberration and caustic caused by the inhomogenous layer overburden were eliminated and multi-value curve was transformed into the curve.of single value. The computation of Marmousi shows that it is feasible. Wave field continuation extends the Fresnel Zone stack's application.

辽河盆地东部凹陷热历史及其与油气关系

This dissertation focuses on the basin geothermal history, tectonothermal evolution and the relationship between geothermal field evolution and hydrocarbon generation. Based on the research of present-day geothermal field, geothermal history of Eastern Subdepression of Liaohe Basin was reconstructed with available data from drillings, loggings, seismic cross-sections, BHTs and thermal indicators. 12 heat flow density values were calculated. Ranging from 45.7 mW/m~2 to 70.0 mW/m~2, the mean value of these determinations exhibits 58.0(±5.83mW/m~22). The heat flow density in the uplift and ramp is greater than that in the sag. The main factors affecting the heat flow density are the morphostructure of basement and thickness of sedimentary cover. The Subdepression experienced a much higher heat flow period from 43 Ma to about 25 Ma. The heat flow increased gradually from Sha3 stage to Dongying stage, and reached the peak at the late of Dongying stage, then cooled down. Structural subsidence analysis shows that the subsidence of Eastern Subdepression can be divided into two phases: earlier (25-43Ma) initial (rift) and late (since 25Ma) thermal subsidence. The lower present-day heat flow and the higher palaeo-heat flow corresponding to structural subsidence stage as well as the typical rift subsidence style in Eastern Subdepression provide with some insights to the tectonic-thermal evolution of the basin. The source rocks of Sha3 member began to generate oil in the Shal stage, and entered oil-window at the late of Dongying stage. The source rocks of Shal began to generate oil at the late of Dongying stage, and being at the stage of lower maturation-maturation now. Most of Dongying source rocks are not mature now. The late of Dongying stage is the critical time for the oil system.

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

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

中国东海陆架盆地构造-热演化研究

East China Sea Shelf Basin (ECSSB), as a basin with prospect of oil & gas resource and due to its special geological location on the west margin of the west Pacific, attracts a lot of attention from many geologists in the world.Based on systematic temperature measurements, bottom hole temperature (BHT) and the oil temperature data, the geothermal gradients in the ECS SB are calculated and vary from 25 to 43°C/km, with a mean of 32.7°C/km. The geothermal gradient in Fuzhou Sag has the higher value(40.6°C/km) in Taibei Depression than that in others. The lower value (27.2 °C/km) occurs in in Xihu Depression. The middle values occurs in Jiaojiang and Lishui sags in Taibei Depression with a mean value of 34.8 °C/km. Incorporated with the measured thermal conductivity, heat flow values show that the ECSSB is characterized by present-day heat flow around 70.6mW/m2, varying between 55 and 88 mW/m2. No significant difference in heat flow is observed between the Xihu and the Taibei Depressions. These heat flow data suggest that the ECSSB is geothermally not a modem back-arc basin.Applying the paleogeothermal gradient based method, thermal history is reconstructed using vitrinite reflectance (VR) and apatite fission track (AFT) data. The results suggest that the thermal history was different in the Taibei and the Xihu depressions. Paleo-heat flow values when the pre-Tertiary formations experienced their maximum temperature at the end of the Paleocene reached a mean of 81 mW/m2 in the Taibei Depression, much higher than the present-day value. The lower Tertiary sediments in the Xihu Depression experienced maximum temperatures at the end of Oligocene and reached a mean paleo-heat flow value of 83.4 mW/m2. The time, when the paleo-heat flow reached the maximum value, suggests that the ECSSB rifted eastward.Tectonic subsidence analysis shows that the timing of the major rifting episode was different across the ECSSB. The rifting occurred from the Late Cretaceous to the early Eocene in the Taibei Depression, followed by thermal subsidence from the late Eocene to the end of Miocene. In contrast, in the Xihu Depression the initial subsidence lasted until the early Miocene and thermal subsidence to the end of Miocene. From Pliocene to the present, an accelerated subsidence took place all along the West Pacific margin of the east Asia.The thermal lithosphere thickness is determined by temperature profile in the lithosphere, the mantle adiabat or the dry basalt solidus. It indicates that the thermal lithosphere reached the thinnest thickness at the end of Eocene in the Taibei Depression and the end of Oligocene in the Xihu Depression, respectively, corresponding with a value of 57-66km and 56-64km. In Taibei Depression, the lithosphere thickness decreased 16-22km from the end of Mesozoic to Paleocene. After Paleocene, the thickness increased 13-16km and reached 71-79 km at present-day. In Xihu Depression, From the end of Oligocene to present-day, the thickness increased 10-13km and reached 69-76km at present-day. The evolution of the lithosphere thickness is associated closely with the lithosphere stretching.Combining the reconstructed thermal history and the burial history, the maturation of the Jurassic oil-source rock shows that the main hydrocarbon generation phase was in the mid-Jurassic and a secondary hydrocarbon generation occurred at the end of Paleocene. The secondary generation was controlled mainly by the tectono-thermal background during the Paleocene.