基于波动方程的层间多次波预测方法研究

In exploration geophysics，velocity analysis and migration methods except reverse time migration are based on ray theory or one-way wave-equation. So multiples are regarded as noise and required to be attenuated. It is very important to attenuate multiples for structure imaging, amplitude preserving migration. So it is an interesting research in theory and application about how to predict and attenuate internal multiples effectively. There are two methods based on wave-equation to predict internal multiples for pre-stack data. One is common focus point method. Another is inverse scattering series method. After comparison of the two methods, we found that there are four problems in common focus point method: 1. dependence of velocity model; 2. only internal multiples related to a layer can be predicted every time; 3. computing procedure is complex; 4. it is difficult to apply it in complex media. In order to overcome these problems, we adopt inverse scattering series method. However, inverse scattering series method also has some problems: 1. computing cost is high; 2. it is difficult to predict internal multiples in the far offset; 3. it is not able to predict internal multiples in complex media. Among those problems, high computing cost is the biggest barrier in field seismic processing. So I present 1D and 1.5D improved algorithms for reducing computing time. In addition, I proposed a new algorithm to solve the problem which exists in subtraction, especially for surface related to multiples. The creative results of my research are following: 1. derived an improved inverse scattering series prediction algorithm for 1D. The algorithm has very high computing efficiency. It is faster than old algorithm about twelve times in theory and faster about eighty times for lower spatial complexity in practice; 2. derived an improved inverse scattering series prediction algorithm for 1.5D. The new algorithm changes the computing domain from pseudo-depth wavenumber domain to TX domain for predicting multiples. The improved algorithm demonstrated that the approach has some merits such as higher computing efficiency, feasibility to many kinds of geometries, lower predictive noise and independence to wavelet； 3. proposed a new subtraction algorithm. The new subtraction algorithm is not used to overcome nonorthogonality, but utilize the nonorthogonality's distribution in TX domain to estimate the true wavelet with filtering method. The method has excellent effectiveness in model testing. Improved 1D and 1.5D inverse scattering series algorithms can predict internal multiples. After filtering and subtracting among seismic traces in a window time, internal multiples can be attenuated in some degree. The proposed 1D and 1.5D algorithms have demonstrated that they are effective to the numerical and field data. In addition, the new subtraction algorithm is effective to the complex theoretic models.

MB倒转在中国黄土中的记录及海陆对比

Chinese loess preserved in northwest and north China are famous for its fine grain size, high accumulate rate and high community and can be good archives for paleoclimate and paleomagnetic variation over the later Cenozoic, with which can be correlated well between marine sediments. Major geomagnetic chrons and long term paleoclimate changes in Quaternary are successfully extracted from Chinese loess-paleosols, as well as short-term geomagnetic excursions and climate instability of high resolution. Magneticstratigraphy based on paleogeomagnetic polarity reversal recorded in Chinese loess is a basic project in loess research since decades ago. True geomagnetic records and exact location of geomagnetic reversal boundary in section is the foundation of magneticstratigraphy. Matuyama-Brunhes (MB) reversal as the youngest one still remains divarication about exact location of its boundary (MBB). L8 and S8 of Luochuan and Xifeng located in the interior of Chinese Loess Plateau(CLP) and Mangshan in southeast part of CLP are chosen to make clear some problems which include magnetic mineral, process of MB reversal, location of MBB, downward displacement scale of magnetic reversal boundary, time lag of paleoclimate record by marine and loess, new correlation between Chinese loess and marine sediments. Rock-magnetic investigations carried on L8 and S8 show that the main mineral are ferrimagnetic assemblage consists of magnetite and maghemite in Luochuan and Xifeng, and magnetite in Mangshan, which all contains little hematite belongs to antiferromagnetic phases. The main carrier of nature remanet magnetism (NRM) is detrital magnetite with pseudo-single domain. Detailed paleomagnetic investigations display that there are several rapid reversals in direction during the process of MB reversal which started at the upper part of S8 and finished at the lower part of L8, and lasted about decades of centimeters to more than 100cm correspond to about 104 years. On the assumption that MBB is located in the middle part of the layer which recorded the very reversal, 11cm is considered as the scale of downward displacement for the MBB in Chinese loess after estimation through correlation between Luochuan and Mangshan records. So this study denies the theory of large scale displacement of MBB and large scale Lock-in depth of NRM acquired from Chinese loess. Time lag of paleoclimate records in terrestrial sediments and marine sediments is considered after reassessment of correlation between low field susceptibility of Chinese loess with marine oxygen isotope in benthic foraminifera. On the basis of traditional correlation between Chinese loess and marine oxygen isotope, this study document a new scheme which correlates L8 and S8 to MIS18 and MIS19, respectively.

滑坡泥石流信息遥感获取与空间预测

On the issue of geological hazard evaluation（GHE）, taking remote sensing and GIS systems as experimental environment, assisting with some programming development, this thesis combines multi-knowledges of geo-hazard mechanism, statistic learning, remote sensing (RS), high-spectral recognition, spatial analysis, digital photogrammetry as well as mineralogy, and selects geo-hazard samples from Hong Kong and Three Parallel River region as experimental data, to study two kinds of core questions of GHE, geo-hazard information acquiring and evaluation model. In the aspect of landslide information acquiring by RS, three detailed topics are presented, image enhance for visual interpretation, automatic recognition of landslide as well as quantitative mineral mapping. As to the evaluation model, the latest and powerful data mining method, support vector machine (SVM), is introduced to GHE field, and a serious of comparing experiments are carried out to verify its feasibility and efficiency. Furthermore, this paper proposes a method to forecast the distribution of landslides if rainfall in future is known baseing on historical rainfall and corresponding landslide susceptibility map. The details are as following: (a) Remote sensing image enhancing methods for geo-hazard visual interpretation. The effect of visual interpretation is determined by RS data and image enhancing method, for which the most effective and regular technique is image merge between high-spatial image and multi-spectral image, but there are few researches concerning the merging methods of geo-hazard recognition. By the comparing experimental of six mainstream merging methods and combination of different remote sensing data source, this thesis presents merits of each method ,and qualitatively analyzes the effect of spatial resolution, spectral resolution and time phase on merging image. (b) Automatic recognition of shallow landslide by RS image. The inventory of landslide is the base of landslide forecast and landslide study. If persistent collecting of landslide events, updating the geo-hazard inventory in time, and promoting prediction model incessantly, the accuracy of forecast would be boosted step by step. RS technique is a feasible method to obtain landslide information, which is determined by the feature of geo-hazard distribution. An automatic hierarchical approach is proposed to identify shallow landslides in vegetable region by the combination of multi-spectral RS imagery and DEM derivatives, and the experiment is also drilled to inspect its efficiency. (c) Hazard-causing factors obtaining. Accurate environmental factors are the key to analyze and predict the risk of regional geological hazard. As to predict huge debris flow, the main challenge is still to determine the startup material and its volume in debris flow source region. Exerting the merits of various RS technique, this thesis presents the methods to obtain two important hazard-causing factors, DEM and alteration mineral, and through spatial analysis, finds the relationship between hydrothermal clay alteration minerals and geo-hazards in the arid-hot valleys of Three Parallel Rivers region. (d) Applying support vector machine (SVM) to landslide susceptibility mapping. Introduce the latest and powerful statistical learning theory, SVM, to RGHE. SVM that proved an efficient statistic learning method can deal with two-class and one-class samples, with feature avoiding produce ‘pseudo’ samples. 55 years historical samples in a natural terrain of Hong Kong are used to assess this method, whose susceptibility maps obtained by one-class SVM and two-class SVM are compared to that obtained by logistic regression method. It can conclude that two-class SVM possesses better prediction efficiency than logistic regression and one-class SVM. However, one-class SVM, only requires failed cases, has an advantage over the other two methods as only "failed" case information is usually available in landslide susceptibility mapping. (e) Predicting the distribution of rainfall-induced landslides by time-series analysis. Rainfall is the most dominating factor to bring in landslides. More than 90% losing and casualty by landslides is introduced by rainfall, so predicting landslide sites under certain rainfall is an important geological evaluating issue. With full considering the contribution of stable factors (landslide susceptibility map) and dynamic factors (rainfall), the time-series linear regression analysis between rainfall and landslide risk mapis presented, and experiments based on true samples prove that this method is perfect in natural region of Hong Kong. The following 4 practicable or original findings are obtained: 1) The RS ways to enhance geo-hazards image, automatic recognize shallow landslides, obtain DEM and mineral are studied, and the detailed operating steps are given through examples. The conclusion is practical strongly. 2) The explorative researching about relationship between geo-hazards and alteration mineral in arid-hot valley of Jinshajiang river is presented. Based on standard USGS mineral spectrum, the distribution of hydrothermal alteration mineral is mapped by SAM method. Through statistic analysis between debris flows and hazard-causing factors, the strong correlation between debris flows and clay minerals is found and validated. 3) Applying SVM theory (especially one-class SVM theory) to the landslide susceptibility mapping and system evaluation for its performance is also carried out, which proves that advantages of SVM in this field. 4) Establishing time-serial prediction method for rainfall induced landslide distribution. In a natural study area, the distribution of landslides induced by a storm is predicted successfully under a real maximum 24h rainfall based on the regression between 4 historical storms and corresponding landslides.

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

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.

电离层胶片频高图数字化转换与分析研究

The ionogram acquired with the ionospheric vertical sounding method is the oldest data in the history of ionospheric research. Using of modern microelectronics and computer technology to digitalize, analyse and preserve the huge amount of historical film ionogram has become more and more important and urgent. This paper introduced the progress of the film ionogram digitalization by using digital image processing technologies to correct and repair film ionogram and convert them in an exchangeable format. An analysis and conversion software, basing on this method, has been developed for the film ionogram analysis, and then it introduces the application of this software by combining the SAO Explorer program for Wuhan film ionogram and pseudo-color ionogram in Yamagawa in Japan. It shows that our method is reliable，and the developed software is used friendly and provides a positive solution in digitalization and analysis of huge amount of historical film ionogram. Firstly, we briefly introduce the film ionogram and the process of its digitalization. By observing a amount of film ionogram, we obtain some common characteristics of the digitalized film ionogram following as: (1) the image rotation are caused by scanning; (2) the vertical axis of a large number of film ionogram exist more or less tilt and bending ; (3) coordinates of the film ionogram appear the non-uniformity phenomena result from the instability of driving motor rotation and the error of altitudinal cursor orientation. Moreover, based on the characteristics of the film ionogram and the SAO Explorer software which is widely used for the digital ionogram analysis in the world, a new method has been developed for film ionogram procession. The method contains the image geometric correction and film ionogram format conversion. The image geometric correction includes such as image rotation correction, vertical correction and coordinates scale correction. After geometric correction, the BMP file format images will be converted to the SBF file format images. Then, we also discuss the data format converting methods, which include two methods of the image data mapping basing on the normalization and logarithm, and the method of the preprocessing of the noise filtering and the threshold setting. Combining with SAO Explorer software, we successfully obtain ionospheric parameters and electron profile from the converted SBF file format digital ionograms. Based on the above method, we developed the software for the film ionogram to realize its correction analysis and conversion of the image format, and then give a introduction for its function and operation. Subsequently, the software are applied into the Wuhan film ionogram which separately observed in the high solar activity year and the low in 1980s last century. The results reveal the converted SBF digital ionogram almost preserve the all echo information of the film ionogram. Furthermore, we expressly discuss the application to the Wuhan film ionogram in 1958 in order to validate the applicability and credibility of the software. And it is showed that the important information of the film ionogram are maintained into the SBF digital ionogram. It is represented that there is credibility for conversion of the software when it applied in the older film ionogram. In sum, this software could apply to the digitalization and analysis of huge amount of historical film ionogram. Last, we extended the function of the software by bring some new conversion method and used it to apply to the pseudo-color ionogram of yamagawa in Japan. The results show that the converted ionogram information basically maintain the importantly ionogram information and the error of scaling of converted SBF file format image is almost acceptable, though there is no preprocessing for the original ionogram. Hence, we could extend the applicable range of the software and apply it to all kinds of simulative ionogram imaging by improving the method and software.

AE指数的可靠性研究

Auroral electrojet index (AE) are usually used to quantitatively describe the activity of the geomagnetic field in the polar region. AE is a means to identify the level of a substorm as well. The auroral electrojet indices (AU, AL, and AE) have served well for more than four decades as measures of magnetospheric substorm activity. However, as substorm studies have progressed considerably during the past several years, the accuracy of the present electrojet indices have become an important issue. Thus it is fortunate for us to reexamine and evaluate the accuracy of the present electrojet indices and improve them if necessary. For a better use of the present indices and for future improvement we examine the limitations of the auroral electrojet indices as an accurate quantitative measure of the auroral electrojets and magnetospheric substorms. Some of the limitations arise from the data availability and also from the present simplified scheme in deriving them, but some of them originate in the definition themselves. In the present paper, we analyze the characteristics of aurora, ionospheric current system and the AE index. It is noted that for the pseudo-breakup events the variations of AE ( > 500 nT) clearly show the sudden increase and slow decrease phases. However, aurora does not expand to poleward or equatorward, and the ionospheric currents presents the features of the magnetic convection. We mainly focus on investigating what time the AE could be used to identify the subtorm.

岩性油藏储层预测的地震理论与方法研究

At present the main object of the exploration and development (E&D) of oil and gas is not the structural oil-gas pools but the subtle lithological oil-gas reservoir. Since the last 90's, the ratio of this kind of pools in newly-added oil reserves is becoming larger and larger, so is the ratio in the eastern oilfields. The third oil-gas resource evaluation indicates the main exploration object of Jiyang depression is the lithological oil-gas pools in future. However, lack of effective methods that are applied to search for this kind of pool makes E&D difficult and the cost high. In view of the urgent demand of E&D, in this paper we deeply study and analyze the theory and application in which the seismic attributes are used to predict and describe lithological oil-gas reservoirs. The great results are obtained by making full use of abundant physics and reservoir information as well as the remarkable lateral continuity involved in seismic data in combination with well logging, drilling-well and geology. ①Based on a great deal of research and different geological features of Shengli oilfield, the great progresses are made some theories and methods of seismic reservoir prediction and description. Three kinds of extrapolation near well seismic wavelet methods-inverse distance interpolation, phase interpolation and pseudo well reflectivity-are improved; particularly, in sparse well area the method of getting pseudo well reflectivity is given by the application of the wavelet theory. The formulae for seismic attributes and coherent volumes are derived theoretically, and the optimal method of seismic attributes and improved algorithms of picking up coherent data volumes are put forward. The method of making sequence analysis on seismic data is put forward and derived in which the wavelet transform is used to analyze not only qualitatively but also quantitatively seismic characteristics of reservoirs.② According to geologic model and seismic forward simulation, from macro to micro, the method of pre- and post-stack data synthetic analysis and application is put forward using seismic in close combination with geology; particularly, based on making full use of post-stack seismic data, "green food"-pre-stack seismic data is as possible as utilized. ③ In this paper, the formative law and distributing characteristic of lithologic oil-gas pools of the Tertiary in Jiyang depression, the knowledge of geological geophysics and the feasibility of all sorts of seismic methods, and the applied knowledge of seismic data and the geophysical mechanism of oil-gas reservoirs are studied. Therefore a series of perfect seismic technique and software are completed that fit to E&D of different categories of lithologic oil-gas reservoirs. ④ This achievement is different from other new seismic methods that are put forward in the recent years, that is multi-wave multi-component seismic, cross hole seismic, vertical seismic, and time-lapse seismic etc. that need the reacquisition of seismic data to predict and describe the oil-gas reservoir. The method in this paper is based on the conventional 2D/3D seismic data, so the cost falls sharply. ⑤ In recent years this technique that predict and describe lithologic oil-gas reservoirs by seismic information has been applied in E&D of lithologic oil-gas reservoirs on glutenite fans in abrupt slop and turbidite fans in front of abrup slop, slump turbidite fans in front of delta, turbidite fans with channel in low slope and channel sanbody, and a encouraging geologic result has been gained. This achievement indicates that the application of seismic information is one of the most effective ways in solving the present problem of E&D. This technique is significant in the application and popularization, and positive on increasing reserves and raising production as well as stable development in Shengli oilfield. And it will be directive to E&D of some similar reservoirs

黄东海及朝鲜海峡冰后期陆架沉积物的环境磁学研究

There are four chapters in this dissertation. The first chapter briefly synthesizes the basic theories, methods and present-day applying situation of environmental magnetism. The second chapter probes into the magnetic mineral diagenesis in the post-glacial muddy sediments from the southeastern South Yellow Sea and its response to marine environmental changes, using the muddy sediment of Core YSDP103 formed in the shelf since about 13 ka BP. The third chapter illustrates the high-resolution early diagenetic processes by investigating the rapidly deposited muddy sediments during the last 6 ka in Cores SSDP-102 and SSDP-103 from the near-shore shelf of Korea Strait. The fourth chapter presents the results of detailed rock magnetic investigation of the surface sediments from the fine-grained depositional area on the outer shelf of the East China Sea in an attempt to provide environmental magnetic evidence for the provenance of the fine-grained deposit. Core YSDP103 was retrieved in the muddy deposit under the cold eddy of the southeastern South Yellow Sea, and the uppermost 29.79 m core represents the muddy sediments formed in the shelf since about 13 ka BP. The lower part from 29.79 to 13.35 m, called Unit A2, was deposited during the period from the post-glacial transgression to the middle Holocene (at about 6 ~(14)C ka BP) when the rising sea level reached its maximum, while the upper part above 13.35 m (called Unit Al) was deposited in a cold eddy associated with the formation of the Yellow Sea Warm Current just after the peak of post-glacial sea level rise. For the the uppermost 29.79 m core, detailed investigation of rock-magnetic properties and analyses of grain sizes and geochemistry were made. The experimental results indicate that the magnetic mineralogy of the core is dominated by magnetite, maghemite and hematite and that, except for the uppermost 2.35 m, the magnetic minerals were subject to reductive diagenesis leading to significant decline of magnetic mineral content and the proportion of low-coercivity component. More importantly, ferrimagnetic iron sulphide (greigite) is found in Unit A2 but absent in Unit Al, suggesting the control of marine environmental conditions on the magnetic mineral diagenesis. Magnetic parameters show abrupt changes across the boundary between the Unit Al and A2, which reflects a co-effect of environmental conditions and primary magnetic components of the sediments on the diagenesis. Alternating zones of high and low magnetic parameters are observed in Unit A2 of Core YSDP103, which is presumably due to periodic changes of the concentration and/or grain size of magnetic minerals carried into the study area. Cores SSDP-102 and SSDP-103, two studied sediment cores from the Korea Strait contain mud sequences (14 m and 32.62 m in thickness) that were deposited during the last 6,000 years. Analyses of grain sizes and geochemistry of the cores have demonstrated that the sediments have uniform lithology and geochemical properties, however, marked down-core changes in magnetic properties suggest that diagenesis has significantly impacted the magnetic properties. An expanded view of early diagenetic reactions that affect magnetic mineral assemblages is evident in these rapidly deposited continental shelf sediments compared to deep-sea sediments. The studied sediments can be divided into four descending intervals, based on magnetic property variations. Interval 1 is least affected by diagenesis and has the highest concentrations of detrital magnetite and hematite, and the lowest solid-phase sulfur contents. Interval 2 is characterized by the presence of paramagnetic pyrite and sharply decreasing magnetite and hematite concentrations, which suggest active reductive dissolution of detrital magnetic minerals, the absolute minimum abundance of magnetite is reached at the end of this interval. Interval 3 is marked by a progressive loss of hematite with depth, and at the base of this interval, 82% to 88% of the hematite component was depleted and the bulk magnetic mineral concentration was reduced to the lowest value in the entire studied mud section. Interval 4 has an increasing down-core enhancement of authigenic greigite, which is interpreted to have formed due to arrested pyritization resulting from consumption of pore water sulfate with depth. This is the first clear demonstration from an active depositional environment for a delay of thousands of years for acquisition of a magnetization carried by greigite. This detailed view of diagenetic processes in continental shelf sediments suggests that studies of geomagnetic field behavior from such sediments should be conducted with care. Paleoceanographic and paleoclimatic studies based on the magnetic properties of shelf sediments with high sedimentation rates like those in the Korea Strait are also unlikely to provide a meaningful signature associated with syn-depositional environmental processes. The rock magnetic properties of the surface sediments from the fine-grained depositional area on the outer shelf of the East China Sea, an area surrounded by sands, are investigated with a view to providing information on the sediment provenance. Multiple magnetic parameters such as magnetic susceptibility (%), anhysteretic remanent magnetization (ARM), saturation rernanent magnetization (SIRM), coercivities of SIRM (Her), and S ratios (relative abundance of low-coercivity magnetic minerals) are measured for all 179 surface samples, and partial representative samples are examined for their magnetic hysteresis parameters, temperature-dependence of magnetic susceptibility and x-ray diffraction spectra. Our research indicates that the magnetic mineralogy is dominated by magnetite with a small amount of hematite and is primarily of pseudo-single domain (PSD) to multidomain (MD) nature with a detrital origin. In the surface sediments, the granulometry of magnetic fractions is basically independent of grain sizes of the sediment containing the magnetic grains, and the composition of magnetic minerals remains almost homogeneous, that is, with a relatively constant ratio of low to high coercivity fraction throughout the area. The magnetic concentration in the study area generally decreases to the east or southeast accompanied by magnetic-particle fining to the east or to the northeast. The geographic pattern of magnetic properties is most reasonably explained by a major source of sediment jointly from the erosion of the old Huanghe River deposit and the discharge of the Changjiang River. The rock magnetic data facilitate understanding of the transport mechanism of fine-grained sediments in the outer shelf of the East China Sea.

矾山杂岩体成因及其与磷（铁）矿的成矿关系

The Fanshan complex consists of layered potassic ultramafic-syenite intrusions. The Fanshan apatite (-magnetite) deposit occurs in the Fanshan complex, and is an important style of phosphorus deposit in China. The Fanshan complex consists of three (First- to Third-) Phases of intrusion, and then the dikes. The First-Phase Intrusive contains ten typical layered rocks: clinopyroxenite, biotite clinopyroxenite, coarse-grained biotite clinopyroxenite, pegmatitic orthoclase-biotite clinopyroxenite, variegated orthoclase clinopyroxenite, interstitial orthoclase clinopyroxenite, biotite rock, biotite-apatite rock, biotite rock and magnetite-apatite rock. This layered intrusive consists of nine rhythmic units. Each rhythmic unit essentially comprises a pair of layers: clinopyroxenite at the bottom and biotite clinopyroxenite at the top. The apatite (-magnetite) deposit is situated near the top of rhythmic Unit no. 6 of the First-Phase Intrusive. The Second-Phase Intrusive contains three typical rocks: coarse-grained orthoclase clinopyroxenite, . coarse-grained salite syenite and schorlomite-salite syenite. The Third-Phase Intrusive includes pseudo-trachytic salite syenite, porphyritic augite syenite, fine-grained orthoclase clinopyroxenite and fine-grained salite syenite. The origin of the Fanshan complex is always paid attention to it in China. Because most layered igneous intrusion in the world not only have important deposit in it, but also carry many useful information for studying the formation of the intrusion and the evolvement of magma. Two sketch maps were drawn through orebodies along no. 25 cross-cut on 425 mL and no. 1 cross-cut on 491 mL in the Fanshan mine. Through this mapping, a small-scaled rhythmic layering (called sub-rhythmic layering in the present study) was newly found at the top of the rhythmic Unit no. 6. The concept of sub-rhythmic layering is defined in this article. The sub-rhythmic layering is recognized throughout this apatite-rich part, except for magnetite-apatite rock. Presence of the layered magnetite-apatite rock is one of the characteristics of the Fanshan apatite (-magnetite) deposit. Thus, from this layer downwards six units of sub-rhythmic layering are recognized in the present study. Each unit consists of biotite clinopyroxenite (or biotite rock and biotite-apatite rock) layer at the bottom and apatite rock layer at the top. To study this feature in detail is an important work for understanding the origin of the Fanshan complex and apatite (-magnetite) deposit. The origin of the Fanshan complex and the relation of the formation of the apatite(-magnetite)deposit will be interpreted by the study of sub-rhythmic layering on the basis of previous research works. The magma formed the Fanshan complex was rich in K2O, early crystallized pyroxene, and after this phase more biotite crystallized, but no amphibole appeared. This indicated that the activity of H2O in the magma was low. Major element compositions of biotite and clinopyroxene (on thin sections) in the sub-rhythmic layering were analyzed using electron microprobe analyzer. The analytical results indicate Mg/(Mg+Fe*+Mn) atomic ratios (Fe*, total iron) of these two minerals rhythmically changed in sub-rhythmic layering. The trends of Mg/(Mg+Fe*+Mn) atomic ratio (Fe*, total iron) of biotite and clinopyroxene indicate that the magma evolved markedly from relatively magnesian bottom layer to less magnesian top layer in each sub-rhythmic unit. A general trend through the sub-rhythmic layering sequence is both minerals becoming relatively magnesian upwards. The formation temperatures for sub-rhythmic layering yield values between 600 and 800 ℃, were calculated using the ratio of Mg/(Mg+Fe+Mn) in the salite and biotite assemblage. The equilibrium pressures in the rhythmic layers calculated using the contents of Al in the salite were plotted in the section map, shown a concave curve. This indicates that the magma formed the First-Phase Intrusive crystallized by two vis-a-vis ways, from its bottom and top to its centre, and the magnetite-apatite rock was crytallized in the latest stage. The values of equilibrium pressures in the sub-rhythmic layering were 3.6-6.8(xlO8) Pa with calculated using the contents of Al in the salite. The characteristics of geochemistry in various intrusive rocks and the rocks or apatite of sub-rhythmic layers indicated that the Fanshan complex formed by the comagmatic crystallization. The contents of immiscible elements and REEs of apatite rock at the top of one sub-rhythmic unit are more than biotite clinopyroxenite at the bottom. The contents of immiscible elements and REEs of apatite of biotite clinopyroxenite at the bottom of one sub-rhythmic unit are higher than apatite rock at the top. The curves of rocks (or apatite) in the upper sub-rhythmic units are between two curves of the below sub-rhythmic unit in the primitive mantle-normalized trace element abundance spider diagram and the primitive mantle-normalized REE pattern. The trend for the contents of immiscible elements and REEs inclines to the same contents from the bottom to the top in sub-rhythmic layering. These characteristics of geochemistry of rocks or apatites from sub-rhythmic layering indicate that the latter sub-rhythmic unit was produced by the residual magma after crystallization of the previous sub-rhythmic unit. The characteristics of petrology, petrochemistry, geochemistry in the Fanshan complex and sub-rhythmic layers and the trends of Mg/(Mg+Fe+Mn) atomic ratio of biotite and clinopyroxene in sub-rhytmic layering rejected the hypotheses, such as magma immiscibility, ravitational settling and multiple and pulse supplement of magma. The hypothesis of differentiation by crystallization lacks of evidences of field and excludes by this study. On the base of the trends of formation temperatures and pressures, the characteristics of petrology, petrochemistry, geochemistry for the Fanshan complex and the characteristics of geochemistry for the rocks (or apatites), the trends of Mg/(Mg+Fe+Mn) atomic ratio of biotite and clinopyroxene in sub-rhytmic layering, and the data of oxygen, hydrogen, strontium and neodymium isotopes, this study suggests that the magma formed the Fanshan complex was formed by low degree partial melting of mantle at a low activity of H2O, and went through the differentiation at the depth of mantle, then multiply intruded and crystallized. The rhythmic layers of the First-Phase Intrusive formed by the magma fractional crystallized in two vis-a-vis ways, from the bottom and top to the centre in-situ fractional crystallization. The apatite (-magnetite) deposit of the Fanshan complex occurs in sub-rhythmic layering sequence. The the origin of the sub-rhythmic layering is substantially the origin of the Fanshan apatite (-magnetite) deposit. The magma formed the rhythmic layers of First-Phase Intrusive was rich in H2O, F and P at the later stage of its in-situ fractional crystallization. The Fanshan apatite (-magnetite) deposit was formed by this residual magma in-situ fractional crystallization. The magnetite-apatite rock was crystallized by two vis-a-vis ways at the latest stage in-situ fractional crystallization in the rhythmic layers. The result was light apatite layer below heavy the magnetite-apatite layer, formed an "inversion" phenomenon.

人工源频率域大地电磁法正反演研究

The CSAMT method is playing an important role in the exploration of geothermal and the pre-exploration in tunnel construction project recently. In order to instruct the interpretation technique for the field data, the forward method from ID to 3D and inversion method in ID and 2D are developed in this paper for the artificial source magnetotelluric in frequency domain. In general, the artificial source data are inverted only after the near field is corrected on the basis of the assumption of half-homogeneous space; however, this method is not suitable for the complex structure because the assumption is not valid any more. Recently the new idea about inversion scheme without near field correction is published in order to avoid the near field correction error. We try to discuss different inversion scheme in ID and 2D using the data without near field correction.The numerical integration method is used to do the forward modeling in ID CSAMT method o The infinite line source is used in the 2D finite-element forward modeling, where the near-field effect is occurred as in the CSAMT method because of using artificial source. The pseudo-delta function is used to modeling the source distribution, which reduces the singularity when solving the finite-element equations. The effect on the exploration area is discussed when anomalous body exists under the source or between the source and exploration area; A series of digital test show the 2D finite element method are correct, the results of modeling has important significant for CSAMT data interpretation. For 3D finite-element forward modeling, the finite-element equation is derived by Galerkin method and the divergence condition is add forcedly to the forward equation, the forward modeling result of the half homogeneous space model is correct.The new inversion idea without near field correction is followed to develop new inversion methods in ID and 2D in the paper. All of the inversion schemes use the data without near field correction, which avoid introducing errors caused by near field correction. The modified grid parameter method and the layer-by-layer inversion method are joined in the ID inversion scheme. The RRI method with artificial source are developed and finite-element inversion method are used in 2D inversion scheme. The inversion results using digital data and the field data are accordant to the model and the known geology data separately, which means the inversion without near field correction is accessible. The feasibility to invert the data only in exploration area is discussed when the anomalous body exists between the source and the exploration area.

Hydroisomerization of long-chain alkane over Pt/SAPO-11 catalysts synthesized from nonaqueous media

SAPO-11 molecular sieves were synthesized from nonaqueous media. The effects of Si and Al sources as well as solvents on the catalytic performance of SAPO-11 were investigated by the hydroisomerization reaction of n-dodecane. The samples were characterized by XRD, XRF, N-2-adsorption, SEM, NH3-TPD, IR-NH3 and Si-29 CP MAS NMR. The SAPO-11 samples synthesized with tetraethoxysilane as the Si source showed higher Si incorporation contents than the SAPO molecular sieves prepared with polymeric Si sources (fumed silica and Si colloidal gel). The reaction results showed that Pt/SAPO-11 catalysts synthesized from ethylene glycol and glycerol media with the monomeric Si and Al sources (tetraethoxysilane, aluminum isopropoxide) exhibited higher catalytic activities than those catalysts with the polymeric Si or Al (pseudo-boehmite) sources, due to the larger external surface area and higher acidity of the former ones. Especially, the catalyst synthesized in an ethylene glycol medium possessed the highest catalytic activity. Over this catalyst, 88% conversion of n-dodecane was achieved at a low temperature of 250 degrees C.