Application of altimetry data assimilation on mesoscale eddies simulation

Mesoscale eddy plays an important role in the ocean circulation. In order to improve the simulation accuracy of the mesoscale eddies, a three-dimensional variation (3DVAR) data assimilation system called Ocean Variational Analysis System (OVALS) is coupled with a POM model to simulate the mesoscale eddies in the Northwest Pacific Ocean. In this system, the sea surface height anomaly (SSHA) data by satellite altimeters are assimilated and translated into pseudo temperature and salinity (T-S) profile data. Then, these profile data are taken as observation data to be assimilated again and produce the three-dimensional analysis T-S field. According to the characteristics of mesoscale eddy, the most appropriate assimilation parameters are set up and testified in this system. A ten years mesoscale eddies simulation and comparison experiment is made, which includes two schemes: assimilation and non-assimilation. The results of comparison between two schemes and the observation show that the simulation accuracy of the assimilation scheme is much better than that of non-assimilation, which verified that the altimetry data assimilation method can improve the simulation accuracy of the mesoscale dramatically and indicates that it is possible to use this system on the forecast of mesoscale eddies in the future.

大气季节内振荡对东亚夏季风活动的影响研究

利用ERA40逐日再分析资料、NCEP/NCAR2逐日再分析资料、中国740个测站日降水资料、上海台风研究所提供的西太平洋热带气旋资料、Kaplan等重建的月平均SSTA资料、NOAA逐日长波辐射（OLR）等资料，应用离散功率谱分析、带通滤波、EOF分析等统计方法，研究了东亚夏季风(EASM）的移动特征、东亚地区季节内振荡(ISO)的基本特征、季节内振荡对东亚夏季风活动的影响、季节内振荡对东亚夏季风异常活动的影响机理。主要结论如下： （1）综合动力和热力因素定义了可动态描述东亚夏季风移动和强度的指数，并利用该指数研究了东亚夏季风的爆发和移动的季节内变化及其年际和年代际变化特征。研究发现，气候平均东亚夏季风前沿分别在28候、33候、36候、38候、40候、44候出现了明显的跳跃。东亚夏季风活动具有显著的年际变率，主要由于季风前沿在某些区域异常停滞和突然跨越北跳或南撤引起，造成中国东部旱涝灾害频繁发生。东亚夏季风的活动具有明显的年代际变化，在1965年、1980年、1994年发生了突变，造成中国东部降水由“南旱北涝”向“南涝北旱”的转变。 （2）东亚季风区季节内变化具有10~25d和30~60d两个波段的季节内振荡周期，以30-60d为主。存在三个主要低频模态，第一模态主要表征了EASM在长江中下游和华北地区活动期间的低频形势；第二模态印度洋-菲律宾由低频气旋式环流控制，主要表现了ISO在EASM爆发期间的低频形势；第三模态主要出现在EASM在华南和淮河活动期间的低频形势。第一模态和第三模态是代表东亚夏季风活动异常的主要低频形势。 （3）热带和副热带地区ISO总是沿垂直切变风的垂直方向传播。因此，在南海-菲律宾东北风垂直切变和副热带西太平洋北风垂直切变下，大气热源激发菲律宾附近交替出现的低频气旋和低频反气旋不断向西北传播，副热带西太平洋ISO以向西传播为主。中高纬度地区，乌拉尔山附近ISO以向东、向南移动或局地振荡为主；北太平洋中部ISO在某些情况下向南、向西传播。 （4）季风爆发期，伴随着热带东印度洋到菲律宾一系列低频气旋和低频反气旋， 冷空气向南输送，10~25天和30~60天季节内振荡低频气旋同时传入南海加快了南海夏季风的爆发。在气候态下，ISO活动表现的欧亚- 太平洋（EAP）以及太平洋-北美（PNA）低频波列分布特征(本文提出的EAP和PNA低频波列与传统意义上的二维定点相关得到的波列不同)。这种低频分布形式使得欧亚和太平洋中高纬度的槽、脊及太平洋副热带高压稳定、加强，东亚地区的低频波列则成为热带和中高纬度ISO相互作用影响东亚夏季风活动的纽带。不同的阶段表现不同的低频模态，30~60d低频模态的转变加快了EASM推进过程中跳跃性；30-60d低频模态的维持使得EASM前沿相对停滞。 （5）30-60d滤波场，菲律宾海域交替出现的低频气旋和低频反气旋不断向西北传播到南海-西太平洋一带。当南海-西太平洋地区低频气旋活跃时，季风槽加强、东伸，季风槽内热带气旋（TC）频数增加；当南海-西太平洋低频反气旋活跃时，季风槽减弱、西退，TC处于间歇期，生成位置不集中。 （6）在El Nino态下，大气季节内振荡偏弱，北传特征不明显，但ISO由中高纬度北太平洋中部向南和副热带西太平洋向西的传播特征显著，东亚地区ISO活动以第三模态为主，EASM集中停滞在华南和淮河流域，常伴随着持续性区域暴雨的出现，易造成华南和江淮流域洪涝灾害，长江和华北持续干旱。在La Nina态下，大气季节内振荡活跃，且具有明显的向北传播特征，PNA低频波列显著，东亚地区ISO活动以第一模态单峰为主；EASM主要停滞在长江中下游和华北地区，这些地区出现异常持续强降水，华南和淮河流域多干旱；在El Nino态向La Nina态转换期，ISO活动以第一模态双峰为主，长江中下游常常出现二度梅。

热带印度洋－太平洋热力异常联合模及其对我国旱涝的影响研究

热带太平洋、印度洋的热力异常对我国夏季降水都有重要影响。然而，由于两大洋的海温变化并非完全孤立，而是相互联系的，因此将两大洋作为一个整体来考虑，研究其热力异常的时空变化，建立一个能够较好表征两大洋热状况变化的联合模指数，进而探讨其对我国旱涝的影响既有重要的理论意义又有明显的实用价值。 本文根据1955－2003年的全球海洋上层热含量资料、SODA数据、我国降水资料、850hPa月平均纬向风资料及副热带高压指数等资料，利用经验正交函数（EOF）分解法和典型相关分析等方法，分析了热带印度洋－太平洋上层热含量距平场的时空变化，提出了一个热带印度洋－太平洋热力异常联合模指数，并研究了该联合模的变化原因及其对我国夏季降水的影响。得到的研究结果主要如下。 1．热带印度洋-太平洋热力异常联合模具有显著的年际和年代际变化，并在1976年前后发生了一次由冷到暖的气候跃变。该联合模的变化特征能够较好地反映出太平洋ENSO和印度洋偶极子的基本信息。 2. 热带印度洋－太平洋热力异常联合模的年代际变化主要是由西太平洋暖池的纬向变异和赤道中东太平洋纬向风异常引起的，而该联合模的年际变化则是西太平洋暖池和东印度洋暖池的纬向变异以及赤道中太平洋和赤道东印度洋纬向风异常共同作用的结果。 3．热带印度洋－太平洋热力异常联合模对我国夏季降水有一年到半年的超前影响，其中以冬季联合模对次年夏季降水的影响最显著。冬季联合模对我国夏季降水的影响可能是通过夏季西太平洋副热带高压的变动来实现的。

油气勘探中的弱磁异常提取与区域场分析

Petroleum and natural gas is an important strategic resources. The short of the reserves will block the development of economy and threaten the safety of nation, along with the main oil fields of our country coming to the height of power and splendor of the exploitation and exploration. Therefore, it makes a great sense to inaugurate new explorative field and increase the reserves of petroleum and natural gas. Magnetic exploration is a main method of geophysics exploration. the developing observation apparatus and the perfect processing method provide wide space for magnetic exploration in these years. The method of magnetic bright spot is an application of magnetic exploration. The vertical migration of the hydrocarbon changes physical and chemical environment above the hydrocarbon reservoir, the new environment make tervalent iron translate into bivalent iron, that produce small scale magnetic anomaly, that is magnetic bright spot. The method of magnetic bright spot explores oil and gas field by the relation between the hydrocarbon and magnetic anomaly. This paper systemically research to pick-up and identify magnetic bright spot combining an oil field item, then point out advantaged area. In order to test the result, the author use the seismic information to superpose the magnetic bright spot, that prove the magnetic bright spot is reliable. then, the author complete a software to pick and identify the magnetic bright spot. The magnetic basement is very important to research forming and evolvement of the basin, especially, it is a crucial parameter of exploring residual basin in the research on pre-Cenozoic residual. This paper put forward a new method to inverse the interface of the magnetic layer on the basis of previous work, that is the method of separation of magnetic field step by step. The theory of this method is to translate the result of magnetic layer fluctuation to the result of magnetization density change, and the magnetic layer is flat, the paper choose thickness of magnetic layer as unit thickness, and define magnetic layer as a unit-thickness layer in order to convenient calculation, at the same time, define the variational magnetization density as equivalent magnetic density. Then we translate the relation between magnetic field and layer fluctuation to the relation between magnetic field and equivalent magnetic density, then, we can obtain the layer fluctuation through calculating equivalent magnetic density. Contrast to conventional parker method, model experimentation and example checkout prove this method is effective. The merit of this method is to avoid flat result in a strongly fluctuant area because of using a uniform average depth, the result of this method is closer to the fact, and this method is to inverse equivalent magnetic density, then translate equivalent magnetic density to layer fluctuation, this lays a foundation to inverse variational magnetic density in the landscape orientation and portrait.

基于随机介质理论的多尺度储层非均质性定量描述

Application of long-term exploration for oil and gas shows that the reservoir technology of prediction is one of the most valuable methods. Quantitative analysis of reservoir complexity is also a key technology of reservoir prediction. The current reservoir technologies of prediction are based on the linear assumption of various physical relationships. Therefore, these technologies cannot handle complex reservoirs with thin sands, high heterogeneities in lithological composition and strong varieties in petrophysical properties. Based on the above-mentioned complex reservoir, this paper conducts a series of researches. Both the comprehending and the quantitative analysis of reservoir heterogeneities have been implemented using statistical and non-linear theories of geophysics. At the beginning, the research of random media theories about reservoir heterogeneities was researched in this thesis. One-dimensional (1-D) and two-dimensional (2-D) random medium models were constructed. The autocorrelation lengths of random medium described the mean scale of heterogeneous anomaly in horizontal and deep directions, respectively. The characteristic of random medium models were analyzed. We also studied the corresponding relationship between the reservoir heterogeneities and autocorrelation lengths. Because heterogeneity of reservoir has fractal nature, we described heterogeneity of reservoir by fractal theory based on analyzing of the one-dimensional (1-D) and two-dimensional (2-D) random medium models. We simulated two-dimensional (2-D) random fluctuation medium in different parameters. From the simulated results, we can know that the main features of the two-dimensional (2-D) random medium mode. With autocorrelation lengths becoming larger, scales of heterogeneous geologic bodies in models became bigger. In addition, with the autocorrelation lengths becoming very larger, the layer characteristic of the models is very obvious. It would be difficult to identify sandstone such as gritstone, clay, dense sandstone and gas sandstone and so on in the reservoir with traditional impedance inversion. According to the obvious difference between different lithologic and petrophysical impedance, we studied multi-scale reservoir heterogeneities and developed new technologies. The distribution features of reservoir lithological and petrophysical heterogeneities along vertical and transverse directions were described quantitatively using multi-scale power spectrum and heterogeneity spectrum methods in this paper. Power spectrum (P spectrum) describes the manner of the vertical distribution of reservoir lithologic and petrophysical parameters and the large-scale and small-scale heterogeneities along vertical direction. Heterogeneity spectrum (H spectrum) describes the structure of the reservoir lithologic and petrophysical parameters mainly, that is to say, proportional composition of each lithological and petrophysical heterogeneities are calculated in this formation. The method is more reasonable to describe the degree of transverse multi-scale heterogeneities in reservoir lithological and petrophysical parameters. Using information of sonic logs in Sulige oil field, two spectral methods have been applied to the oil field, and good analytic results have been obtained. In order to contrast the former researches, the last part is the multi-scale character analysis of reservoir based on the transmission character of wave using the wavelet transform. We discussed the method applied to demarcate sequence stratigraphy and also analyzed the reservoir interlayer heterogeneity.

东天山红山高硫型浅成低温铜-金矿床氧化带的形成机理和硫酸盐矿物的系列发现及意义

Located in the Paleozoic uplift along the southern margin of Tu-Ha basin in eastern Xinjiang, the newly discovered Hongshan Cu-Au deposit occurs in the superimposed Mesozoic volcanic basin upon the north section of later Paleozoic Dananhu-Tousuquan accretionary arc. Kalatage Cu-Au orebelt is controlled by NWW-trend faults, and includes Hongshan and Meiling Cu-Au deposits. The host rocks of Hongshan ore district are mainly rhyolitic-dacitic ignimbrites, whereas Cu-Au mineralization is closely related to quartz porphyry, rhyolitic porphyry and granitic porphyry. Mineralization styles are dominantly veinlet-disseminated and veinlet, occasionally stockwork. The mineral association is chalcopyrite, pyrite, bornite, chalcocite and sphalerite. The hydrothermal alteration consists of silicfication, sericitization, alunitization, pyrophylitization, illitization, hydromuscovitization, and chloritization. Hongshan Cu-Au deposit, on the edge of the desert, is one of the driest areas in eastrn Tianshan. Moreover， the highest temperature has been up to 60℃, and the average rainfall receives only 34.1mm/y. The light rainfall and rapid evaporation in the vicinity of this deposit have allowed the formation of a great variety of water-soluble sulfates. Oxidization zone of this deposit lies on the upper part of primary sulfide orebodies appearing with a depth of 50-60m, which is dominant in sulfate minerals. 1. Based on the field observation, the volcanic and sub-volcanic rock composition, hydrothermal alteration, ore structure and mineralization characteristics, this paper proposed that the Hongshan Cu-Au deposit belongs to a transitional type from high-sulfide epithermal to porphyry Cu-Au deposit, which corresponds with the typical HS-epithermal deposit such as Zijinshan Au-Cu deposit in Fujian Province, SE-China. 2. The Hongshan copper-gold deposit was controlled by the tectonic, stratum, magma activity and volcanic apparatus, whereas Au mineralization is closely related to quartz porphyry, rhyolitic porphyry and fine grained pyritization in hydrothermal activity, and Cu mineralization is closely related to quartz porphyry and hydrothermal explosive breccia. 3. Oxidation zone of Hongshan Cu-Au deposit lies on the upper part of primary sulfide orebodies deposit. 23 sulfate minerals were identified in this work. The results of samples XRD and chemical analysis were furthermore confirmed through thermal, infrared spectrum and mössbauer spectrum analysis. Among those, nine minerals as Ferricopiapite, Cuprocopiapite, Rhomboclase, Parabutlerite, Krausite, Yavapaiite, Metasideronatrite Kroehnkite and Paracoquimbite were founded in China for the first time. And Paracoquimbite was secondly reported in the world (first case reported at 1938 in Chile). 4. EPMA analysis shows that Al impurity in crystal lattice is important to polytype formation of paracoquimbite and coquimbite besides stack fault. 5. Compared with Meiling Cu-Au deposit in the same Kalatage ore belt from the characteristics of δ34S of barite, lithofacies, hydrothermal alteration and homogeneous temperature, Hongshan Cu-Au deposit belongs to the same metallogenic system of HS-epithermal type as Meiling Cu-Au deposit. But Hongshan Cu-Au deposit has less extensive alteration and shallower denudation. 6. Sulfur isotope analyses show that δ34S values of pyrites vary in the range of +1.86‰~+5.69‰, with an average of 3.70‰, mostly in the range of +1.86‰~+3.20‰, and δ34Scp<δ34Spy. Therefore ore-forming fluid of porphyry comes from mantle and was contaminated by the earth’s crust. Sulfur isotope has reached balance in ore-forming process. 7. Sulfur isotope analyses show that δ34S values of sulfates vary in the range of +2.15‰~+6.73‰, with an average of +3.74‰, mostly equals as δ34S values of primary sulfides in Hongshan Cu-Au deposit. So supergene sulfates inherit sulfur of primary sulfide. δ34S values are mostly same in different sulfates. As well as pyrite and chalcopyrite, volcanic hot spring and associated native sulfur underground also provide water medium and sulfur during the formation process of sulfate. 8. According to the EPMA of sample chalcopyrite and pyrite in Hongshan Cu-Au, the value of Cu/Ni is 0.98-34.72, mostly close to the value of 5, which shows that Hongshan deposit is a typical volcanogenic magmaic hypothermal deposit. Au and Ag, Zn, Te and Bi are positive correlation, Cu and Hg, Se, Sb are positive correlation, indicates Au and Cu don’t locate in the factor of mineralization of same mineralization groups. The reasons of gold concentration in the oxidation zone are: 1). Change of redox potential (Eh) makes gold to deposit from the liquid of mineralization zone; 2). PH is one of the most factors of gold’s deposition; 3). Soluble complex and colloid of gold can be adsorbed easily. 9. The biotite and hornblende K-Ar isotopic ages from the wall rock－quartz diorite, biotite granite and monzonite granite are 231.99±3.45Ma, 237.97±2.36Ma and 296.53±6.69Ma respectively. The ore-bearing rhyolitic breccia lava contains breccia of the biotite granite which indicates the volcanism and related Cu-Au mineralization occurred later than the granite, possibly in Mesozoic. K-Ar ages of granitoids in Sanya, Baishiquan and Hongliugou area and Molybdenite Re-Os age of Baishan Mo deposit all are in Triassic. Besides late Paleozoic magmatism, igneous magmatic event of Mesozoic was widespread in eastern Tianshan. 10. The K-Ar age dating indicates that the K-Ar age of Voltaite occurred below surface 1m is 56.02±3.98Ma, K-Ar age of Ferricopiapite occurred below surface 1.5m is 8.62±1.12Ma, K-Ar age of Yavapaiite occurred below surface 14 m is 4.07±0.39Ma, and K-Ar age of Voltaite occurred below surface 10 m is 14.73±1.73Ma. So the age interval of oxidation zone of Hongshan copper-golden bed is between 60 -3.38Ma. Oxidization occurred at Caenozoic era (from 65Ma), which can be identified through comparing with different deposits oxidation zone in other countries. The coupling between global tectonic event and climatic change event which occur from Caenozoic era has some effect on epigeosphere system, which can act on the surface of bed oxidation zone similarly. It induces that the age mentioned above coincide with collision of India-Asia and multistage uplifting of Qinhai-Tibet Plateau happened subsequently. Bed oxidation zone is the effect and record of collision and uplifting of Tibet Plateau. The strong chemical weathering of surface accumulation to which was leaded by PETM event occurred Paleocene and Eocene is the reason of Voltaite sharply rises. On the contrary, Ferricopiapite formed due to the global cold weather. The predecessor did much research through biota, isotopes, susceptibility, but this paper try to use different sulfate mineral instead of climatic change. So the research of sulfate minerals not only indicates a great deal of oxidized zone feature, but also the intergrowth of sulfate minerals may be used to trace paleoenviroment and paleoclimate of oxidation zone. 11. Analysis of the information of alteration and mineralization features of four bore cores, induced activity polarization well logging and Eh-4 geophysical section, deep mineralization anomaly objects of Hongshan ore districts shows low resistance, middle and high polarization, measurements of Eh-4 consecutive conductance section show the existing of concealed porphyry ore body deeper than 450m, on the top of and around rock body there are low resistance body ranged from 100-300Ω•m, this area may be the ore-bearing part. In a word, Hongshan Cu-Au deposit deposit is a combine of upper HS-style epithermal Au deposit and deeper porphyry mineralization system. It has great potential to find large HS-style epithermal-porphyry Au-Cu deposits. This paper consists of seven chapters and twenty seven sections. The geological character of deposit is basic condition in this work. Constitute of oxidation zone, research of sulfate mineral, relation between oxidation and primary zone, K-Ar ages of potassic sulfate are key parts of thesis. Genesis of ore deposit is the further expansion of this research. Analysis of ore-controlling factors is the penetration above basic. Analysis of potential is application of exploration.

西天山造山带古生代构造演化与地壳增生—来自花岗岩和蛇绿岩的证据

The Central Asian Orogen Belt (CAOB), which is different from the subductional orogen and the collisional orogen, is known as the most important site of crustal growth in the Phanerozoic, and it has been a ‘hot spot’ for studying the orogenic belts. The Chinese West Tianshan Orogen is occupying the west-southern part of the CAOB and is of great importances to understand the orogenic processes and the continental growth in the Central Asia. The West Tianshan Orogen had undergone complex tectonic evolutional processes in Paleozoic times and large volumes granitic rocks have recorded important information about these processes. Litter is known about Phanerozoic continental growth in the Western Tianshan area so far, compared with the other areas of the CAOB, such as eastern Junggar, western Junggar, Altai and Alakol. The aim of this dissertation is to set up the chronology frame of granitoids in western Tianshan, provide new evidence for the tectonic evolution and discuss the Paleozoic continental growth in this area, on the basis of the studies on the isotopic chronology, major element, trace element and Nd-Sr isotopic geochemistry of granitoids and the isotopic chronology and geochemistry of the ophiolites in this area, especially the Kule Lake ophiolites. 25 precise SHRIMP U-Pb zircon and LA-ICPMS U-Pb zircon ages have been obtained in this dissertation. The granitic rocks in western Tianshan had been formed during two periods: the granitic gneiss with an age of 896Ma, possibly representing the forming age of the Precambrian basement; the granitic rocks with ages varying from 479Ma to 247Ma, recording the Paleozoic orogenic process of western Tianshan. The granitoids in western Tianshan are composed of intermediate-basic rocks, intermediate rocks, intermediate-acid rocks and acid rocks, mainly intermediate-acid rocks and acid rocks. They are mostly granite, granodiorite, quartz syenite and monzodiorite. Different types of granitic rocks are exposed in different tectonic units. The granitoids on the northern margin of the Yili Plate mainly formed in late Paleozoic (413Ma ~ 281Ma), those with ages varying from 413Ma to 297Ma show continental arc affinities and the magnesian calc-alkalic metaluminous diorite of 281Ma display the geochemical characteristics similar to those of granites formed during the post-orogenic period. The granitiods on the southern margin of the Yili Plate include the adakite diorite of 470Ma which was formd by partial melting of thickened lower crust, the post-collisional alkali-feldspar granite of 430Ma, the volcanic arc granite of 348Ma and the Triassic post-collisional granite. The granitoids in the Central Tianshan Plate formed in 479Ma ~ 247Ma, mainly in 433Ma ~ 321Ma. The granitic rocks with ages of 479Ma ~ 321Ma are magnesian calc-alkalic to alkalic rocks with continental arc affinities. A few post-collisional granitoids of 276Ma ~ 247Ma may have inherited the geochemical characteristics of pre-existing arc magma. The granitic rocks in Southern Tianshan (northern margin of the Tarim plate) formed two stages, 420Ma ~ 411Ma and ca. 285Ma. The magnesian calcic to alkalic granites of 420Ma ~ 411Ma may formed during the extension process of the continental margin. The granite of 285Ma includes mostly ferroan calc-alkalic to alkali-calcic rocks with high SiO2 and high alkaline contents, and obviously negative anomaly of Eu, Ba, Sr, P, Ti, similar to the geochemical characteristics of the A-type granite which is formed during post-collisional extension. The Kule Lake ophiolite in southern Tianshan shows the affinity of N-MORB. A SHRIMP zircon U-Pb age of 425±8Ma has obtained for gabbros. Some zircons have given another group of 206Pb/238U age 918Ma, which may indicate the information of the pre-exist old basement rock. The small oceanic basin represented by Kule Lake ophiolite probably developed on the split northern margin of Tarim block. A model for Paleozoic tectonic evolution of the West Tianshan Orogen has been proposed here on the basis of the new results obtained in this dissertation and the previous published data. In Early Cambrian, the Terskey Ocean occurred along the North Nalati fault (NNF), and it separated the Yili plate from the Central Tianshan plate which was probably connected with the Tarim plate. The Terskey Ocean probably subducted towards south under the Central Tianshan plate and towards north under the Yili plate simultaneously. In the early stage of Late Ordovician, the Terskey Ocean had been closed, and the Yili and Central Tianshan plates collided. Meanwhile, extension happened within the joint Central Tianshan and Tarim plates gradually and the Paleo-South Tianshan Ocean had been formed. In Early Silurian, the Paleo-South Tianshan Ocean began to subduct beneath the composite Yili-Central Tianshan plate, which was intruded by volcanic arc granitoids. In Middle Silurian, the Paleo-South Tianshan Ocean, which had reached a certain width, was subducting strongly. And this subduction may have produced voluminous granitoids in the Central Tianshan plate. In the latest stage of Carboniferous, the Paleo-South Tianshan ocean closed, and the Yili-Central Tianshan plate and Tarim plate collided. In Late Cambrian, Paleo-Junggar Ocean occurred to north of the Yili plate; and started to subduct towards south under the Yili plate in Ordovician. This subduction may have produced a magma arc on the northern margin of the Yili plate. In Late Carboniferous, the Paleo-Junggar Ocean had been closed. The Yili-Central and Junggar plates amalgamated together. The West Tianhan Orogen may undergo a post-collisional collapse since Permian. And the magmatic activities may continue to early Triassic. The initial 87Sr/86Sr ration of the granitic rocks in the western Tianshan Mountains varies from 0.703226 to 0.716343, and Nd（t）from -6.50 to 2.03. The characteristics of Sr-Nd isotope indicate that the source of granitic material is not a sole source, which may be produced by mantle-crust magma mixing. In Paleozoic time, lateral growth of the continental crust along active continental margins was dominant, whereas the vertical growth of continental crust resulted from post- collisional mantle derived magmas was not obvious.

均衡重力异常计算与综合解释——以黄海地区为例

Petroleum and Natural Gas is an important strategic resources. The reserves of Petroleum and Natural Gas can’t meet the need of our country, which also blocks the development of economy and threatens the safety of national. Therefore, it makes a great sense to bring “the second round of oil & gas exploration” into effect and study the exploration of oil and gas of Pre-Cenozoic residual basins in China. The integrated geophysical exploration is the main way to research the Pre-Cenozoic residual basins. Gravity exploration is one of the most important exploration methods, which has played an important role in oil and gas prospecting, such as compartmentalizing geotectonic elements, delineating the distribution range of sedimentary basins, searching oil and gas structure, abstracting oil and gas information, and so on, from its naissance. The isostatic gravity anomalies is significant for exploration, which can help us research deep crustal structure, the equilibrium state of earth, the geologic structure of shallow crust, the basement shape of sedimentary basins and the genetic evolution of sedimentary basins. In the paper, we stress the implication and physical meanings systemically, and discuss the calculation theory. On the basis of previous work, we test different isostatic compensation models and parameters to find out their influences to the result of isostatic gravity anomalies. In addition, we improve the method of isostatic gravity anomalies calculation and give a system of isostatic gravity anomalies calculation which is proved has satisfying effect. From the research above, we find that the results of Platt model and Airy model are consistent, which have similar form and almost the same value. However, by contrast, the Airy model is proved has better adaptability than Platt model. The two main parameters——crust thickness and density difference of crust and mantle, both have influence to the isostatic gravity anomalies, but the latter have more. Finally, we adopt the regional field extending edge method to make the result more of actual geologic condition. On the methods above, we calculate the isostatic gravity anomalies field in Yellow Sea area from the Bouguer gravity anomalies and the water depth and altitude data. And then the isostatic gravity anomalies character is analyzed and the integrated geological-geophysical interpretation is made on the basis of summarizing the previous research result systemically and analyzing other geophysical data and geological information. From the research, we find that the Yellow Sea area belongs to continental type crust equilibrium regions, where the isostatic gravity anomalies field is placid and has less fluctuation values, which implies that the area is in equilibrium state to different extends.

地球化学场的多重分形分析与分解

Geochemical and Geophysical anomaly play an important role in mineral exploration，their spatial structure character include singularity and self－similar。The singularity of an anomaly reflects the enrichment characters of the geochemical element ,The anomaly separation by multifractal model is useful in mineral anomaly assessment。In recent years, The continuous multifractal mode of the geochemical fields was studied ,it can be separated into the simple continuous multifractal mode and the high concentration multifractal mode, and the C－A、C－D、 S－A、MSDV、W－A method to decompose the anomaly were presented。Those are succeeded in interpretation of Geochemical and Geophysical anomaly。 This study makes a summarization to these method， we present a multifractal method based on wavelet transform to analyze the multifractal fields 。The singularity and spectrum are calculated through tracing the wavelet maximum modulus in different measure，and then the fields can be decomposed by the characters of the singularity。 It is demonstrated to be useful in interpretation of Au anomaly in Gekou－Shicheng region Rushang Shandong Province 。 Based on the multifractal theory , Using the concentration—area（C—A）method ,We study two geochemical fields in Chifeng area , Inner Mongolia。The results show that the geochemical fields have three different multifractal modes。Based on these ,we discuss the enrichment mode of the geochemical elements and their distributions in space and get the anomaly lower limit ,then the geochemical backgrounds、regional anomalies and local anomalies are distinguished

基于三重震相波形模拟研究亚洲东北部上地幔过渡带速度结构

Both the global and regional P wave tomographic studies have revealed significant deep structural heterogeneities in subduction zone regions. In particular, low-velocity anomalies have been observed beneath the descending high-velocity slabs in a number of subduction zones. The limited resolution at large depths and possible trade-off between the high and low velocities, however, make it difficult to substantiate this feature and evaluate the vertical extent of the low-velocity structure. From broadband waveform modeling of triplicated phases near the 660-km discontinuity for three deep events, we constrained both the P and SH wave velocity structures around the base of the upper mantle in northeast Asia. For the two events beneath the southern Kurile, the rays traveled through the lowermost transition zone and uppermost lower mantle under the descending Pacific slab. Our preferred models consistently suggest normal-to-lower P and significantly low SH velocities above and below the 660-km discontinuity extending to about 760-km depth compared with the global IASP91 model, corroborating previous observations for a slow structure underneath the slab. In contrast, both high P and SH velocity anomalies are shown in our preferred models for the Japan subduction zone region, likely reflecting the structural feature of a slab stagnant above the 660-km discontinuity. The velocity jumps across the 660-km discontinuity were found to be on average 4.5% and 7% for P and S waves under the south Kurile, and 3% and 6% under the Japan subduction zone. The respective velocity contrasts in the two regions are consistent with mineralogical models for colder slab interior and hotter under-slab areas. Based on mineral physics data, the depth-averaged ~1.5% P and ~2.5% SH velocity differences in the depth range of 560-760 km between the two regions could be primarily explained by a 350~450K temperature variation, although the presence of about 0.5wt%~1wt% water might also contribute to the subtle velocity variations near the base of the transition zone in the southern Kurile. From our modeling results, we speculate that the slow structure in the southern Kurile may be correlated to the low velocity zone observed previously around the 410-km discontinuity under Northern Honshu. Both are probably associated with a thermal anomaly rooted in the lower mantle beneath the subduction zone in northeast 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.

白垩纪超静磁带地球磁场古强度测定

The long-term variations in the strength of the geomagnetic field provide important constraints on the chemical-physical processes of the Earth’s interior. Especially, the intensity of the geomagnetic field during the Cretaceous normal superchron (CNS) is crucial to understand the geodynamo. But a paucity of paleointensity further limits to obtain essential knowledge interior process in the deep earth. In order to improve the experimental efficiency, this study tried to apply two new rock magnetic methods (FORC diagram and low-temperature demagnetization technique) to determine the paleointensity. First, some problems in the theory and technique in paleointensity experiments were introduced. A combined palaeomagnetic and geochronologic study was further conducted on a basaltic lava sequence at Jianchang in Liaoning Province, northeastern China. Radiometric 40Ar/39Ar dating indicates that the volcanism occurred at about 119 Ma within the marine anomaly C34n in Cretaceous normal superchron (CNS). Rock magnetic investigations show that pseudo-single domain (PSD) titanium-poor titanomagnetite is dominant in the studied lava flows. Both stepwise thermal and alternating field demagnetizations isolate the well-defined normal characteristic remanent magnetization (ChRM) in three independent lava flows with a mean direction of D/I = 6.0/51.9 degree(a95 = 12.3degree). Palaeointensity was determined using the modified Thellier method with systematic partial thermoremanent magnetization (pTRM) checks on total 72 samples, but only 10 samples exhibit ideal linear behavior on the Arai-plot in the temperature interval of 300-560 C and yield an average paleointensity of (25.8+/-1.4)uT. In addition, slopes of the line defined by the initial and the final points on the Arai-plot for the other 18 samples with characteristic PSD features give an average paleointensity estimation value of (24.8+/-1.9)uT. The consistency of these two approaches confidently demonstrates the fidelity of our results. The overall mean field strength determined using both approaches are thus estimated to be (25.2+/-0.7 )uT. This value corresponds to the virtual dipole magnetic moments (VDM) of (4.5+/-0.1)E22 Am^2, which is about half of the value of present field. This finding suggests that palaeointensity just at the onset of the CNS is characterized by a weak magnetic field.

地震资料叠前储层预测综合研究--以辽中凹陷东二下段浊积扇气藏为例

Reservoir prediction techniques from prestack seismic are among the most important ones for exploration of lithologic hydrocarbon reservoir. In this paper, we set the turbidite fan sandstone reservoir in Liao-Zhong depress as our researching target, and aims to solve the apllication difficulties on pre-stack inversion in the area, where the drilling data is scarce and the reservoir is lateral varied. Meanwhile, AVO analysis and pre-stack inversion for gas-bearing detection is systematically researched. The seismic reflection characters of gas-bearing sandstone in turbidite fan with different fluid content are defined, after analyzing results from AVO seismic simulation and porous fluid replacement of real log data, and under the guides of the seismic characters from classical gas-bearing sandstone reservoir and numerical simulation for complicate gas-bearing sandstone. It is confirmed that detecting gas-bearing sandstone in turbidite fan via AVO technologies is feasible. In terms of AVO analysis, two AVO characters, fluid detection factor and product of intercept and gradient, can effectively identify top and bottom boundaries and lateral range of tuibidite gas sand by comparing real drilling data. Cross-plotting of near and far angle stack data could avoid the correlation existing in P-G analysis. After comparing the acoustic impedance inversions with routine stacked data and AVO intercept, impedance derived from AVO intercept attribute could reduce the acoustic impedance estimating error which is caused by AVO. On the aspect of elastic impedance inversion, the AVO information in the pre-stack gathers is properly reserved by creating partial angle stack data. By the far angle elastic impedance alone, the gas sand, with abnormally low range of values, can be identified from the background rocks. The boundary of gas sand can also be clearly determined by cross-plotting of near and far angle elastic impedances. The accuracy of far angle elastic impedance is very sensitive to the parameter K, and by taking the statistical average of Vp/Vs on the targeted section in key wells, the accuracy of low frequency trends is gurranteed; the intensive absorsion within the area of the gas sand, which tends to push the spectral of seismic data to the lower end, will cause errors on the inversion result of elastic impedance. The solution is to confine the inversion on the interested area by improving the wavelet. On the aspect of prestack AVA simultaneous inversion, the constraint of local rock-physical trends between velocities of P-wave、S-wave and density successfully removes the instability of inversion, thus improves the precision of the resulting elastic parameters. Plenty of data on rock properties are derived via AVO analysis and prestack seismic data inversion. Based on them, the fluid anomaly is analysized and lithological interpretation are conducted. The distribution of gas sand can be consistently determined via various of ways, such as cross-plotting of P and G attributes, near and far partial angle stack data, near and far angle elastic impedances, λρ and Vp/Vs, etc. The shear modulo and density are also reliable enough to be used for lithological interpretation. We successfully applied the AVO analysis and pre-stack inversion techniques to gas detecting for turbidite fan sand reservoir in Liao-Zhong depression.

含激电信息的人工源频率域电磁法正演研究——大功率人工固定源极低频电磁法激电效应预研究

ABSTRACT Base on Dias-model, this paper study the forward modeling to IP response with MT souece and CSAMT source with the author’s program，that is to provide the theoretical basis for extracting the IP information from the far field, near field and transition field in extremely low frequency electromagnetic sounding with artificial strong fixed source signal. The outline of the paper is as follows. Firstly, review the history of the complex-resistivity-model, and analyses the bases for choosing Dias model; meanwhile, effects and responses of each parameter in Dias model are analyses. Afterward, study the forward modeling to 1D layered model with MT source, numerically simulate the IP effects of classic geoelectric objects; Find the clear anomaly and relationship between the peak value of the amplitude anomaly ratio and phase anomaly difference with and without IP parameters to be considered within a frequency range. On the basis of the modeling of MT, Study the CSAMT modeling with a dipole-source, obtain the anomaly responses and the relationship between anomaly and characters of object. Base on the infinite line source, study the 2D IP effect of geoelectric objects with 2D modeling. Calculate the response with different source distance, object depth and different wall rock resistivity systemically. Finally, conclude the viewpoints and give the discussion of the result, and point out the lack of this research also.

太阳活动变化对中低纬电离层的影响

With the variations of solar activity, solar EUV and X-ray radiations change over different timescales (e.g., from solar cycle variation to solar flare burst). Since solar EUV and X-ray radiations are the primary energy sources for the ionosphere, theirs variations undoubtedly produce significant and complicated effects on the ionosphere. So the variations of solar activity significantly affect the ionosphere. It is essential for both ionospheric theory and applications to study solar activity effects on the ionosphere. The study about solar activity variations of the ionosphere is an important part of the ionospheric climatology. It can enhance the understanding for the basic processes in the ionosphere, ionospheric structure and its change, ionosphere/thermosphere coupling, and so on. As for applications, people need sufficient knowledges about solar activity variations of the ionosphere in order to improve ionospheric models so that more accurate forecast for the ionospheric environments can be made. Presently, the whole image about the modalities of ionospheric solar activity variations is still unknown, and related mechanisms still cannot be well understood. This paper is about the effects of the 11-year change in solar activity to the low- and mid-latitude ionosphere. We use multi-type ionospheric observations and model to investigate solar activity effects on the electron density and ionospheric spatial structure, and we focus on discussing some related mechanisms. The main works are as follows: Firstly, solar activity variations of ionospheric peak electron density (NmF2) around 1400 LT were investigated using ionosonde observations in the 120°E sector. The result shows that the variation trend of NmF2 with F107 depends on latitudes and seasons. There is obvious saturation trend in low latitudes in all seasons; while in middle latitudes, NmF2 increases linearly with F107 in winter but saturates with F107 at higher solar activity levels in the other seasons. We calculated the photochemical equilibrium electron density to discuss the effects induced by the changes of neutral atmosphere and dynamics processes on the solar activity variations of NmF2. We found that: (1) Seasonal variation of neutral atmosphere plays an important role in the seasonal difference of the solar activity variations of NmF2 in middle latitudes. (2) Less [O]/[N2] and higher neutral temperature are important for the saturation effect in summer, and the increase of vibrational excited N2 is also important for the saturation effect. (3) Dynamics processes can significantly weaken the increase of NmF2 when solar activity enhances, which is also a necessary factor for the saturation effect. Secondly, solar activity variations of nighttime NmF2 were investigated using ionosonde observations in the 120°E sector. The result shows that the variation trends of NmF2 with F107 in nighttime are different from that in daytime in some cases, and the nighttime variation trends depend on seasons. There is linear increase trend in equinox nighttime, and saturation trend in summer nighttime, while the increase rate of NmF2 with F107 increases when solar activity enhances in winter nighttime (we term it with “amplification trend”). We discussed the possible mechanisms which affect the solar activity variations of nighttime NmF2. The primary conclusions are as follows: (1) In the equatorial ionization anomaly (EIA) crest region, the plasma influx induced by the pre-reversal enhancement (PRE) results in the change of the variation trend between NmF2 and F107 from “saturation” to “linear” after sunset in equinoxes and winter; while the recombination process at the F2-peak is the primary factor that affects the variation trend of NmF2 with F107 in middle latitudes. (2) The recombination coefficient at the F2-peak height reaches its maximum at moderate solar activity level in winter nighttime, which induces NmF2 attenuates more quickly at moderate solar activity level. This is the main reason for the amplification trend. (3) The change of the recombination process at the F2-peak with solar activity depends on the increases of neutral parameters (temperature, density et al.) and the F2-peak height (hmF2). The seasonal differences in the changes of neutral atmosphere and hmF2 with solar activity are the primary reasons for the seasonal difference in the variation trend of nighttime NmF2 with F107. Finally, we investigated the solar activity dependence of the topside ionosphere in low latitudes using ROCSAT-1 satellite (at 600 km altitude) observations. The primary results and conclusions are as follows: (1) Latitudinal distribution of the plasma density is local time, seasonal, and solar activity dependent. In daytime, there is a plasma density peak at the dip equator. The peak is obviously enhanced at high solar activity level, and the strength of the peak strongly depends on seasons. While at sunset, two profound plasma density peaks (double-peak structure) are found in solar maximum equinox months. (2) Local time dependence of the latitudinal distribution is due to the local time variation of the equatorial dynamics processes. Double-peak structure is attributed to the fountain effect induced by strong PRE. Daytime peak enhances with solar activity since the plasma density increases with solar activity more strongly at the dip equator due to the equatorial vertical drift, and its seasonal dependence is mainly due to the seasonal variations of neutral density and the equatorial vertical drift. In the sunset sector, seasonal and solar activity dependences of the latitudinal distribution are related to the seasonal and solar activity variations of PRE. (3) The variation trend of the plasma density with solar activity shows local time, seasonal, and latitudinal differences. That is different from the changeless amplification trend at the DMSP altitude (840 km). Profound saturation effect is found in the dip equator region at equinox sunset. This saturation effect in the topside ionosphere is realated to the increase of PRE with solar activity. Solar activity variation trend of the topside plasma density was discussed quantitatively by Chapman-α function. The result shows that the effect induced by the change of the scale height is dominant at high altitudes; while the variation trend of ROCSAT-1 plasma density with solar activity is suggested to be related to the changes of the peak height, the scale height, and the peak electron density with solar activity.