基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

班公湖带多不杂超大型富金斑岩铜矿床的成岩成矿年代学、岩石学及高氧化岩浆¬流体-成矿作用

Duobuza copper deposit, newly discovered typical gold-rich porphyry copper deposit with superlarge potential, is located in the Tiegelong Mesozoic tectonic -magmatic arc of the southern edge of Qiangtang block and the northern margin of Bangonghu-Nujiang suture. Quartz diorite porphyrite and grandiorite porphyry, occurred in stock, are the main ore-bearing porphyries. As the emplacement of porphyry stock, a wide range of hydrothermal alteration has developed. Within the framework of the ore district, abundant hydrothermal magnetite developed, and the relationship between precipitation of copper and gold and hydrothermal magnetite seems much close. Correspondingly, a series of veinlets and network veinlets occurred in all alteration zones. Therefore, systematic research on such a superlarge high-grade Duobuza gold-rich porphyry copper deposit can fully revealed the metallogenic characteristics of gold-rich porphyry copper deposits in this region, establish metallogenetic model and prospecting criteria, and has important practical significance on the promotion of regional exploration. In addition, this research on it can enrich metallogenic theory of strong oxidation magma-fluid to gold-rich porphyry copper deposit, and will be helpful to understand the metallogenic characteristics in early of subduction of Gangdese arc stages and its entire evolution history of the Qinghai-Tibet Plateau, the temporal and spatial distribution of ore deposits and their geodynamics settings. Northern ore body of Duobuza copper deposit have been controlled with width (north-south) about 100 ~ 400 m, length (east-west) about 1400 m, dip of 200 °, angle of dip 65 °~ 80 °. And controlled resource amount is of 2.7 million tons Cu with grade 0.94% and 13 tons Au with 0.21g/tAu. Overall features of ore body are large scale, higher grade copper, gold-rich. Ore occurred in the body of granodiotite porphyry and quartz diorite porphyrite and its contact zone with wall rock. Through the detailed mapping and field work studies, some typies of alteration are identificated as follows: albitization, biotititation, sericitization, silication, epidotization, chloritization, carbonatization, illitization, kaolinization and so on. The range of alteration is more than 10km2. Wall alteration zone can be divided into potassic alteration, moderate argillization alteration, argillization, illite-hydromuscovite or propylitization from ore-bearing porphyry center outwards, but phyllic alteration has not well developed and only sericite-quartz veins occurred in local area. Moreover, micro-fracture is development in ore district , and correspondingly a series of veinlets are development as follows: biotite vein (EB type), K-feldspar-biotite-chalcopyrite-quartz vein, magnetite-antinolite-K-feldspar vein, quartz-chalcopyrite-magnetite veins (A-type), quartz-magnetite-biotite-K-feldspar vein, chalcopyrite veinlets in potassic alteration zone; (2) chalcopyrite occurring in the center vein–quartz vein (B type), chalcopyrite veinlets, chalcopyrite-gypsum vein in intermediate argillization alteration; (3) chalcopyrite- pyrite-quartz vein, pyrite-quartz vein, chalcopyrite-gypsum veins, quartz-gypsum- molybdenite-chalcopyrite vein in argillization alteration; (4) gypsum veins, quartz-(molybdenite)-chalcopyrite vein, quartz-pyrite vein, gypsum- chalcopyrite vein, potassium feldspar veinlets, Carbonate veins, quartz-magnetite veins in the wall rock. In short, various veins are very abundant within the framework of the ore district. The results of electronic probe microscopy analysis (EMPA) indicate that Albite (Ab 91.5~99.7%) occurred along the rim of plagioclase phenocryst and fracture, and respresents the earliest stages of alteration. K-feldspar (Or 75.1~96.9%) altered plagioclase phenocryst and matrix or formed secondary potassium feldspar veinlets. Secondary biotite occurred mainly in phenocryst, matrix and veinlets, belong to magnesium-rich biotite formed under the conditions of high-oxidation magma- hydrothermal. Chloritization developed in all alteration zones and alterd iron- magnesium minerals such as biotite and hornblende and then formed chlorite veinlets. As the temperature rises, Si in the tetrahedral site of chlorite decreased, and chlorite component evolved from diabantite to ripiolite. The consistent 280℃~360℃ of formation temperature hinted that chlorite formed on the same temperature range in all alteration zones. However, formation temperature range of chlorite from the gypsum-carbonate-chlorite vein was 190℃~220℃, and it may be the product of the latest stage of hydrothermal activity. The closely relationship between biotite and rutile indicate that most of rutiles are precipitated in the process of biotite alteration and recrystallization. In addition, the V2O3 concentration of rutile from ore body in Duobuza gold-rich porphyry copper deposit is >0.4％, indicate that V concentration in rutile has important significance on marking main ore body of porphyry copper deposit. Apatites from Duobuza deposit all are F-rich. And apatite in the wall rock contained low MnO content and relatively high FeO content, which may due to the basaltic composition of the wall rocks. The MnO in apatite from altered porphyry show a strong positive correlation with FeO. In addition, Cl/F ratio of apatite from wall rock was highest, followed by the potassic alteration zone and potassic alteration zone overprinted by moderate argillization alteration was the lowest. SO2 in Apatite are in the scope of 0 to 0.66%, biotite in the apatite has the highest SO2, followed by the potassic alteration zone, potassic alteration zone overprinted by moderate argillization alteration, and the lowest in the surrounding rocks, which may be caused by the decrease of oxygen fugacity of hydrothermal fluid and S exhaust by sulfide precipitation in potassic alteration. Magnetite in the wall rock have higher Cr2O3 and lower Al2O3 features compared with altered porphyry, this may be due to basalt wall rock generally has high Cr content. And magnetites have higher TiO2 content in potassic alteration than moderate argillization alteration overprinted by potassic alteration, argillization and wall rock, suggested that its formation temperature in potassic alteration was the highest among them. The ore minerals mainly are chalcopyrite and bornite, and Au contents of chalcopyrite, bornite, and pyrite are similar with chalcopyrite slightly higher. The Eu* negative anomaly of disseminated chalcopyrite was relatively lower than chalcopyrite in veinlets. Within a drill hole, the Eu* negative anomaly of disseminated chalcopyrite was gradually larger from bottom to top. Magnetite has the same distribution model, with obvious negative Eu* abnormal, and ΣREE in great changes. The gypsum has the highest ΣREE content and the obvious negative anomaly, and biotite obviously has the Eu* abnormal. Based on the petrographic and geochemical characteristics, five series of magmatic rocks can be broadly classified; they are volcanic rocks of the normal island arc, high-Nb basaltic rocks, adakites, altered porphyry and diorite. The Sr, Nd, Hf isotopes and geochemistry of various series of magmatic rock show that they may be the result of mixing between basic magma and various degrees of acid magma coming from lower crust melted by high temperature basic underplating from partial melting of the subduction sediment melt metasomatic mantle wedge. Furthermore S isotope and Pb isotope of the sulfide, ore-bearing porphyries and volcanic rocks indicated ore-forming source is the mantle wedge metasomatied by subduction sediment melt. Oxygen fugacity of magma estimated by Fe2O3/FeO of whole rock and zircon Ce4+/Ce3+ indicated that the oxidation of basalt-andesitic rocks is higher than ore-forming porphyry, and might imply high-oxidation characteristics of underplated basic magma. Its high oxidative mechanism is likely mantle sources metasomatied by subduction sediment magma, including water and Fe3+. And such high oxidation of basaltic magma is conducive to the mantle of sulfides in the effective access to melt. And the An component of dark part within plagioclase phenocryst zoning belong to bytownite (An 74%), and its may be a result of magma composition changes refreshment by basaltic magma injection. SHRIMP zircon U-Pb and LA-ICP-MS zircon U-Pb geochronology study showed that the intrusions and volcanic rocks from Duobuza porphyry copper deposit belong to early Cretaceous magma series (126~105Ma). The magma evolution series are as follows: the earliest diorite and diorite porphyrite → ore-bearing porphyry and barren grandiorite porphyry →basaltic andesite → diorite porphyrite → andesite → basaltic andesite, and magma component shows a evolution trend from intermediate to intermediate-acid to basic. Based on the field evidences, the formation age of high-Nb basalt may be the latest. The Ar-Ar geochronology of altered secondary biotite, K-feldspar and sericite shows that the main mineralization lasting a interval of about 4 Ma, the duration limit of whole magma-hydrothermal evolution of about 6 Ma, and possibly such a long duration limit may result in the formation of Duobuza super-large copper deposit. Moreover, tectonic diagram and trace element geochemistry of volcanic rocks and diorite from Duobuza porphyry copper deposit confirm that it formed in a continental margin arc environment. Zircon U-Pb age of volcanic rocks and porphyry fall in the range of 105~121Ma, and Duobuza porphyry copper deposit locating in the north of the Bangonghu- Nujiang suture zone, suggested that Neo-Tethys ocean still subducted northward at least early Cretaceous, and its closure time should be later than 105 Ma. Three major inclusion types and ten subtypes are distinguished from quartz phenocrysts and various quartz veins. Vapor generally coexisting with brine inclusions, suggest that fluid boiling may be the main ore-forming mechanism. Raman spectrums of fluid inclusions display that the content of vapor and liquid inclusion mainly contain water, and vapor occasionally contain a little CO2. In addition, the component of liquid inclusions mainly include Cl-, SO42-, Na+, K+, a small amount of Ca2+, F-; and Cl- and Na+ show good correlation. Vapor mainly contains water, a small amount of CO2, CH4 and C2H6 and so on. The daughter minerals identified by Laman spectroscopy and SEM include gypsum, chalcopyrite, halite, sylvite, rutile, potassium feldspar, Fe-Mn-chloride and other minerals, and ore-forming fluid belong to a complex hydrothermal system containing H2O-NaCl-KClFeCl2CaCl2. H and O isotopic analysis of quartz phenocryst, vein quartz, magnetite, chlorite and gypsum from all alteration zones show that the ore-forming fluid of Duobuza gold-rich porphyry copper deposit consisted mainly of magmatic water, without addition of meteric water. Duobuza gold-rich porphyry copper deposit formed by the primary magmatic fluid (600-950C), which has high oxidation, ultra-high salinity and metallogenic element-rich, exsolution direct from the magma, and it is representative of the typical orthomagmatic end member of the porphyry continuum. Moreover, the fluid evolution model of Duobuza gold-rich porphyry copper deposit has been established. Furthermore, two key factors for formation of large Au-rich porphyry copper deposit have been summed up, which are ore-forming fluids earlier separated from magma and high oxidation magma-mineralization fluid system.

山东平度大庄子金矿床地质地球化学特征及成因——兼论层间滑动角砾岩型金矿成矿模式

The author selected the geological and geochemical characteristics and the genesis of the Dazhuangzi gold deposit in the Pingdu City as the central content of the study. The author summarized geological features of the other gold deposits formed within the same geological setting along the margin of the Jiaolai pull-apart basin and compared these gold deposits with the Dazhuangzi gold deposit. On the basis of the first-hand data obtained from field investigation and from mining production reports, ore-controlling structures, geological characteristics and mineralization regularities of the Dazhuangzi gold deposit are studied in detail. According to the analyzing results of petrochemistry, trace element, rare earth element and fluid inclusion etc., the geochemical characteristics, the genesis and the ore-forming material source of the Dazhuangzi gold deposit and that of the other similar gold deposits along the margin of the Jiaolai Basin are proposed. The study results suggest that the Dazhuangzi gold deposit belongs to the typical interstratified glide breccia type gold deposit, which is controlled by the interstratified glide fault structure located along the margin of the Mesozoic pull-apart Jiaolai basin. The interstratified glide fault structure is in the outer part of unconformity belt between the overlying strata and the basement of the pull-apart basin, being along the marble strata of the Jingshan group. The formation of the ore-controlling structure is related closely with the evolution of the Jiaolai Basin in the Mesozoic. The ore-controlling structure underwent the structural stress changes from compressive to tensional and then to compressive stress with strike slipping features sequentially, which were coincided with the regional tectonic stress evolution. The interstratified glide breccia type gold mineralization mainly occurs in the siliceous-marble breccias and cataclastic rocks within the interstratified glide fault structure. The gold minerogenetic epoch is later than 120Ma when the ore-controlling structure was tensioning and strike-slipping. The occurrences of the ore controlling structure and the gold ore bodies are the same as that of the unconformity belt. The geological and geochemical studies show that the source of the ore-forming material is alike with that of the volcanic rocks of the Qingshan formation, which is widespread in the Jiaolai Basin. Both of them came from the deep crust or even the upper mantle. Based on the geological characteristics and the minerogenetic regularities of the Dazhuangzi gold deposit, a genetic model of the deposit is constructed. In addition, the author used the remote sensing image and exploration results of geochemical and geophysical methods to point out several prospecting areas for further exploration. Through comprehensive study on the interstratified glide fault structure and on the interstratified glide breccia type gold deposits along the Jiaolai pull-apart basin, three types of interstratified glide structures and related gold mineralization are set up according to evolution and distribution of main fault as well as related secondary faults in time and space. They are named as Penjiakuang type, Dazhuangzi type and Fayunkuang type. The author summarized up the minerogenetic characteristics and regularities controlled by these three different types of interstratified glide structures respectively, and set up a general minerogenetic model of the interstratified glide breccia type gold deposit.

山西堡子湾金矿角砾岩带控矿规律、成矿时代及隐伏矿体定位预测

The Puziwan gold deposit is located at the northern boundary of the North China platform where there are very favorable ore-forming conditions. The deposit is a medium-sized gold deposit associated with silver and copper polymetallic minerallization which were closely related to cryptoexplosive processes and fluidization and underwent multi-epoch superposition of mineralization. The mineralization consists of cryptoexplosive breccia-type, fluidizing-type, porphyry-type, quartz vein-type, etc. The ore-controlling law of ore-bearing brecciated zone is systematically studied in the paper. The shape, scale and attitude of ore bodies are dominat in the upper, the ore bodies of fluidizing breccia-type are dominant in the upper and the middle part, the porphyry-type ore bodies are dominant in the bottom. The quartz vein-type ore bodies are impenetrated in all the brecciated zone. The metallogenic epoch of Puziwan gold deposit is collated and stipulated. The Rb-Sr isochron age of quartz porphyry (wall rock) is 233 Ma, refering to the Ar-Ar age of gold-bearing quartz, excluded the former conclusion that the so-called metallogenic epoch (245.9Ma) is in Indosinian epoch. The nonage metallogenic age of Buziwan gold deposit shoule be in the Yanshanian epoch (142.5Ma). By applying the sub-specimen sampling technique, the Rb-Sr isochron age of gold-bearing pyrites in late mineralization epoch is dated to be 64 Ma. In conclusion, the metallogenetic epoch of Puziwan gold deposit is in late Yanshanian-early Himalayan epoch. On the above basis, the metallogenic model of the "train type" and new texture model of ore deposit are established. by applying the methods of geology, shallow seismic exploration, remote sensing, partial geochemical extraction and the study of inclusions in the late granite porphyry, the author has made the metallogenic progosis in the depth and the prephery of Puziwan gold deposit and eight prospecting targets are proposed.

地球化学新方法在隐状矿床找矿预测中的应用研究

With the progress of prospecting, the need for the discovery of blind ore deposits become more and more urgent. To study and find out the method and technology for the discovery of blind and buried ores is now a priority task. New geochemical methods are key technology to discover blind ores. Information of mobile components related to blind ores were extracted using this new methods. These methods were tested and applied based on element' s mobile components migrating and enriched in geophysical-geochemical process. Several kinds of partial extraction techniques have tested based on element' s occurrence in hypergenic zone. Middle-large scale geochemical methods for exploration in forest and swamp have been tested. A serious of methods were tested and applied effetely about evaluation of regional geochemical anomaly, 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system instead of the normal net. 1. Element related with ores can be mobiled to migrate upwards and be absorpted by surface soil. These abnomal components can be concentrated by natural or artificial methods. These trace metalic ions partially exist in dissovlvable ion forms of active state, and partially have been absorbed by Fe-Mn oxide, soil and organic matter in the soil so that a series of reaction such as complex reaction have take place. Employing various partial extraction techniques, metallic ions related with the phase of the blind ores can be extracted, such as the technique of organic complex extraction, Fe-Mn oxide extraction and the extraction technique of metallic ions of various absorption phases. 2.1:200000 regional geochemical evaluation anomaly methods: Advantageous ore-forming areas were selected firstly. Center, concentration, morphological feature, belt of anomaly were choosed then. Geological and geochemical anomalies were combined. And geological and geochemical background information were restrained. Xilekuduke area in Fuyun sheet , Zhaheba area in Qiakuerte sheet, the west-north part in Ertai sheet and Hongshanzui anomaly in Daqiao sheet were selected as target areas, in Alertai, in the north of Xinjiang. in Xilekuduke area, 1:25000 soil geochemical methods sampling based on the net in dendritic water system was carried out. Cu anomaly and copper mineralization were determined in the center area. Au , Cu anomalies and high polarization anomaly were determined in the south part. Prospecting by primary halo and organic complex extraction were used to prognosis blind ore in widely rang outcrop of bedrock. 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system were used in transported overburden outside of mining area. Shallow seismic method and primary halo found a new blind orebody in mining area. A mineralization site was fou and outside of Puziwan gold mine, in the north of Shanxi province. Developing middle-large scale geochemical exploration method is a key technique based 1:200000 regional geochemical exploration. Some conditions were tested as Sampling density , distribution sites of sample, grain size of sample and occurrence of element for exploration. 1:50000 exploration method was advanced to sample clast sediment supplement clast sediment in valley. 1:25000 bedrock or soil geochemical methods sampling based on the net in dendritic water system was applied to sample residual material in A or C horizon. 1:2000 primary or soil halo methods used to check anomalies and determine mineralization. Daliang gold mineralization in the northern Moerdaoga was found appling these methods. Thermomagnetic method was tested in miniqi copper-polymetallic ore. Process methods such as grain size of sample, heated temperature, magnetic separating technique were tested. A suite of Thermomagnetic geochemical method was formed. This method was applied in Xiangshan Cu~Ni deposit which is cover by clast or Gobi in the eastern Xinjiang. Element's content and contrast of anomaly with Thermomagnetic geochemical method were higher than soil anomaly. Susceptibility after samples were heated could be as a assessment conference for anomaly. In some sectors thermo-magnetic Cu, Ni, Ti anomalious were found outside deposits area. There were strong anomal ies response up ore tested by several kind of partial extraction methods include Thermomagnetic, enzyme leach and other partial extractions in Kalatongke Cu-Ni deposit in hungriness area in the northern of Xinjiang. Element's anomalies of meobile were mainly in Fe-Mn oxide and salt. A Copper mineralization site in Xilekuduke anomaly area had been determined. A blind ore was foung by shallow seismic and geochemical method and a mineralization site was found outside this mining area in Puziwan gold deposit in shanxi province. A Gold mineralization site was found by 1:50000 geochemical exploration in Daliang, Inner Mongolia.