Sequential Gaussian simulation for geosystems modeling: A machine learning approach

Sequential Gaussian Simulation(SGSIM)as a stochastic method has been developed to avoid the smoothing effect produced in deterministic methods by generating various stochastic realizations.One of the main issues of this technique is,however,an intensive computation related to the inverse operation in solving the Kriging system,which significantly limits its application when several realizations need to be produced for uncertainty quantification.In this paper,a physics-informed machine learning(PIML)model is proposed to improve the computational efficiency of the SGSIM.To this end,only a small amount of data produced by SGSIM are used as the training dataset based on which the model can discover the spatial correlations between available data and unsampled points.To achieve this,the governing equations of the SGSIM algorithm are incorporated into our proposed network.The quality of realizations produced by the PIML model is compared for both 2D and 3D cases,visually and quantitatively.Furthermore,computational performance is evaluated on different grid sizes.Our results demonstrate that the proposed PIML model can reduce the computational time of SGSIM by several orders of magnitude while similar results can be produced in a matter of seconds.     

Analysis Method for Predicting Strain in Interior Beds and Sub-Resolution Faults from Area Balance Theory in Extensional Basins

Extensional basins include mainly grabens and half grabens displaced along a lower detachment. Based on area balance theory, there is a linear relationship between a height of regional and the lower detachment h on the outside of the basin and “lost area S“ from the regional in the basin. The pre-growth beds above lower detachment are of the same extensional displacement so that an “S-h diagram“ can be used to determine the depth to lower detachment and to calculate the total extensional displacement of the beds above the lower detachment. The extensional displacement is dominated by the heave of various scale normal faults. The displacement of obvious faults can be immediately figured out from the measured bed-length. The requisite extension calculated by area balance is the layer-parallel strain, which could be accommodated by displacement on sub-resolution faults. Accordingly, the layer-parallel strain can help us predict the magnitude and distribution of sub-resolution faults on the basis of analysis of the structural style and rheological behavior.     

An Effective Method for Free Vibration of Plate on Elastic Half Space

The vibration analysis of a plate on an elastic foundation is an important problem in engineering. It is the interaction of a plate with the three-dimensional half space and the plate is usually loaded from both the upper and lower surfaces. The contact pressure from the soil can not be predefined. According to Lambs solution for a single oscillating force acting on a point on the surface of an elastic half space, and the relevant approximation formulae, a relation between the local pressure and the deflection of the plate has been proposed. Based on this analysis, the reaction of the soil can be represented as the deformation of the plate. Therefore, the plate can be separated from the soil and only needs to be divided by a number of elements in the analysis. The following procedure is the same as the standard finite element method. This is a semi-analytical and semi-numerical method. It has been applied to the dynamic analysis of circular or rectangular plates on the elastic half space, at low or high frequency vibration, and on rigid, soft or flexible foundations. The results show that this method is versatile and highly accurate.     

Bioavailability of crops to environmental hormone Pb in soils and its remediation

Lead （Pb）, as an environmental hormone, acts as an estrogen in animal and human bodies, and its environmental pollution can result in many procreant anomalies, such as maladjustment of incretion, weakening of procreation function and descending of genital immunity. Experiments on pot-growth of paddies and corns in this study indicated that these crops did not show any macro-toxic symptoms in the observed contents of Pb added （0-4000 mg/kg） in soils. However, lead in soils can be easily absorbed by crops and accumulated in the different parts of plants. Statistical analysis showed that Pb contents in the crops are strongly positively correlated to the Pb contents in soils. In the crops, the Pb contents of roots are highest, tens to thousands of times those of straws and seeds. It is suggested that the root is a barrier of Pb absorbed by straws and seeds. Based on this fact, Pb pollution in soils can be gradually reduced by planting crops and subsequently removing the roots from the soil. This biotechnology may be used in environmental remediation, which embraces the quality of sustainable agricultural economy and a healthy environment.     

Thermal Regimes of Lithosphere and Geothermal Resources Potential in Fujian Province, China

From the determinations of surface heat flow based on silica geothermometry, residual heat flow,radiogenic heat production of rocks, thermal conductivities of rocks and so on, the temperatures at differentdepths of the lithosphere in Fujian Province are estimated by using a one-dimensional steady-state model. Ac-cording to these and the research on geological structures. the possibility of existence of various geothermalsystems in Fujian is discussed. The authors consider that in Fujian, geothermal systems related to shallow magmatism do not exist andhot dry rock systems can not be exploited for the time being. However, intermediate-and low-temperaturehydrothermal systems are well developed. They are controlled by active faults and distributed widely insoutheastern Fujian. but do not have an anomalous heat source. The groundwater circulates at great depthsand has a temperature usually less than 150℃ with less discharge but good quality. It can be used directly as avaluable geothermal water resource but not suitable to be developed as an energy source.     

沉积岩有机质的红外光谱及其在石油有机地球化学中的初步应用

Types and degree of muturity of organic matters in sedimentary rocks as index for oil generation and evolution can be determined by infrared spectrometry. In normal cases,oil-generating rocks can be preliminarily identified by infrared spectra of solubable bitumen components. However, the spectra of condensed nuclei of aromatic hydrocarbon should be resorted to whenever the composition and structure of organic matters arc comparative complex and their degree ofmuturity are not all the same. Systematic researches of known samples show that a configuration of condensed aromatic hydrccarbon can be proposed on the basis of the index of aromatic hydrocarbon structure and the ratio of methyl and methylene to aromatic nucleus. It is evidenced by preliminary experience that oil- and coal-generating organisms can be distinguished from each other with the help of this configuration.     

Past and Future of Mathematical Geology

This is a brief review of alternative methods of problem-solving in geoscience with emphasis on the role of mathematical geology. It is desirable to maintain a clear-cut distinction between reliable facts which can be stored in data banks and concepts that can be incorporated in the speciflcations of sta-tistical models designed for specific purposes. If possible, subjective probabilities shclld be incorporated in hypotheses that are to be tested by statistical inference.     

The Econometric Assessment of Losses Caused by Water Pollution in Chongqing, Southwest China

As an important industrial city, Chongqing has experienced rapid economic growth in recent years, but at the same time the growth has been accompanied by serious pollution problems. Water quality monitoring shows that the waters of Chongqing are widely polluted. In this paper, the human capital approach is applied to assessing the economic losses and damages caused by water pollution in Chongqing. Based on previous studies and empirical investigations,water pollution-induced damage can be categorized into two general types: damage to ecosystems, including agriculture; and damage to non-ecosystems, including industry. The economic loss estimation of each type of the damage is conducted by taking into account a group of relevant parameters and methods.     

Tree-Ring Carbon Isotopic Constraints on Carbon-Water Exchanges between Atmosphere and Biosphere in Drought Regions in Northwestern China

The comparison between the carbon isotope and the index of ring width of a pine disc from the Tuomuer Peak region in Xinjiang shows that the effects of climate changes on the tree-ring growth and carbon iso-topic fractionation varies with time. The reason is probably relative to the characters of climate changes and adaptability of the tree-ring growth to climate changes. The relationships between the atmospheric CO2 level and the revised δ13Cair by the tree-ring carbon isotope indicate that the carbon cycle is not in a steady state, but under a stage-change condition in this area. It also can be concluded that the ratio of CO2 from the terrestrial eco-system has increased, and the flux of CO2 exchange between the atmosphere and the biosphere was gradually increasing over the past century. In addition, the results also confirm the validity and superiority of the carbon isotope to the research of the water-use efficiency.     

Capturing and storage of CO<Subscript>2</Subscript> by micron-nano minerals: Evidence from the nature

The increase of CO2 in atmosphere is a main factor leading to "greenhouse effect", which causes more and more serious global environmental problems. The reduction of CO2 is a challenge for the survival of human beings, and it is also a big technical problem. CO2 fluid-rock interaction is a key scientific problem involved in geo-logical storage. The CO2 fluid-rock interaction has a variety of multi-scale changes. Due to great differences in the quantity of surface atoms and surface energy between micron-nano-sized minerals, and ions and crystals, the speed and efficiency of CO2 fluid-rock interaction on a micron-nano scale are much higher than those on other scales. As is known from the natural world, the micron-nano structures of pores and the surface chemical modification of natural porous minerals (zeolite, diatomite, sepiolite, palygorskite, halloysite, etc.) should be further investigated, which can be used as the micron-nano -mineral porous materials with high capacity and high efficiency for capturing CO2. Through simulating the adsorption capacity and process of CO2 by minerals in the natural world, the micron-nano technology is applied to calcium- and magnesium-based minerals (olivine, pyroxene, feldspar, clay, etc.) so as to improve the activity of calcium and magnesium and enlarge the reaction contact area. In this way, the efficiency of capturing and storage of CO2 by calcium- and magnesium-based minerals can be greatly improved. These minerals can also be used as the micron-nano-mineral materials with large capacity and high efficiency for capturing and storing CO2.     

Groundwater contaminant source identification based on QS-ILUES

When groundwater pollution occurs,to come up with an efficient remediation plan,it is particularly important to collect information of contaminant source(location and source strength)and hydraulic conductivity field of the site accurately and quickly.However,the information can not be obtained by direct observation,and can only be derived from limited measurement data.Data assimilation of observations such as head and concentration is often used to estimate parameters of contaminant source.As for hydraulic conductivity field,especially for complex non-Gaussian field,it can be directly estimated by geostatistics method based on limited hard data,while the accuracy is often not high.Better estimation of hydraulic conductivity can be achieved by solving inverse groundwater problem.Therefore,in this study,the multi-point geostatistics method Quick Sampling(QS)is proposed and introduced for the first time and combined with the iterative local updating ensemble smoother(ILUES)to develop a new data assimilation framework QS-ILUES.It helps to solve the contaminant source parameters and non-Gaussian hydraulic conductivity field simultaneously by assimilating hydraulic head and pollutant concentration data.While the pilot points are utilized to reduce the dimension of hydraulic conductivity field,the influence of pilot points’layout and the ensemble size of ILUES algorithm on the inverse simulation results are further explored.     

Zonal Metamorphic Complex in the Tongulack Mountain Ridge, Altai, Russia, and Explanation of Its Origin with the Help of Thermal Modelling

A moderate pressure / high temperature zonal metamorphic complex in the Tongulack Mountain Ridge, Altai, Russia, is described, and the applicability of the models of magmatic intrusion and fluid flow to explanation of its origin discussed. The Precambrian complex was formed at 500-700℃ and 3.0-5.5 kbars; it is a linear, 25-30 km wide, thermal anticline with a curved axis showing symmetric metamorphic zoning. The metamorphism was isochemical by its nature, as is corroborated by the chemical compositions of the rocks. Four zones can be recognized within the metamorphic complex: chloritic (on the peripheries), cordieritic, sillimanitic and staurolite-out (in the centre). The zones are separated by successive isograds: cordierite, staurolite-in or sillimanite and staurolite-out. It is argued that the origin of the metamorphic zoning can be explained best by a combined fluid-magmatic model; conductive heat flow from the intrusion predominated considerably over the fluid flux in heat transfer: the fluid flow     

断裂带里的粘土岩化问题

Fault clay is the combined results of 10 ig-term dynamic-, thermal-metamorphism and tectonic geochemical process. Approximately similar fault clay can be produced from different rocks as long as the abovemeutioned actions develope to a certain extent.Fault clay is characterized by a great variety of clay minerals as a result of the depletion in silica and the relative enrichment in Ah Different types of clayization can be recognized depending on the environmen and the extent to which alkalis are leached. The process that results in the formation Of fault clay, i.e., the disintegration of rocks, dissolution, transportation, deposition mstasormatism, re-crystallization and the alteration of clay minerals already formed goes Under the heading of tectonic geochemistry. It is induced by the pressure and temperature that prevail under dynamic-, the rmal-metamorphism and is controlled entirely by tectonic movement. Therefore, it differs from weathering and leaching in terms of either the mechanism or the characteristies of material exchange or the texture and structure of the resulted clays. Clayization is an indieator that the fault activity has reached to a stage of matoration. On the other hand, clayization is also a process of mineralization which gives rise under favorite conditions elay deposit of economic slginficance.     

Advanced prediction of tunnel boring machine performance based on big data

Predicting the performance of a tunneling boring machine is vitally important to avoid any possible accidents during tunneling boring.The prediction is not straightforward due to the uncertain geological conditions and the complex rock-machine interactions.Based on the big data obtained from the 72.1 km long tunnel in the Yin-Song Diversion Project in China,this study developed a machine learning model to predict the TBM performance in a real-time manner.The total thrust and the cutterhead torque during a stable period in a boring cycle was predicted in advance by using the machine-returned parameters in the rising period.A long short-term memory model was developed and its accuracy was evaluated.The results show that the variation in the total thrust and cutterhead torque with various geological conditions can be well reflected by the proposed model.This real-time predication shows superior performance than the classical theoretical model in which only a single value can be obtained based on the single measurement of the rock properties.To improve the accuracy of the model a filtering process was proposed.Results indicate that filtering the unnecessary parameters can enhance both the accuracy and the computational efficiency.Finally,the data deficiency was discussed by assuming a parameter was missing.It is found that the missing of a key parameter can significantly reduce the accuracy of the model,while the supplement of a parameter that highly-correlated with the missing one can improve the prediction.     

碳酸盐岩生油岩的有机地球化学、岩石学特征

Features of oil-generating carbonate rocks are studied in terms of petrology and organic geochemistry of fluorescent materials in rocks. It is postulated by the authors that the organic carbon content in oil-generatlng carbonate rocks is necessarily lower than that in shales, and the possible cause of this observation is discussed. The quality of oil-ganerating carbonate rocks may be judged by aromatic structure index,amino acid, pigment index, and DTA data on organic matter. However, the application of presently available organic geochemical criteria to non-oil-generating reservoir rocks may always give misleading results. To circumvent this difficulty, it is necessary first to make distinction between primary and secondary organic matter by cxamining the fluorescent materials under a fluorescence microscope. Measurementof radical content of the organic substances and DTA are useful techniques to determine the degree of maturation. The maturation of organic substances proceeds at a lower rate in carbonate rocks than in shales. It is likely, therefore, that some older carbonate strata can also be promising areas for petroleum exploration. Oil-generating carbonate rocks consist mainly of mierite limestones rich in organic matter and are thought to be formed under low-encrgy environment. Consequently, petrological(especially mierolithofacies) studies are very important in evaluating carbonate source rocks.     

Analysis of the Achievments in Tectonic Dynamopetrologenic and Dynamo—metallogenic Simulating Experiments and Their Geological Significance

Under the action of tectonic stress ore fluids carrying ore-forming materials can migrate from a higher stress district to a lower one.This is a fact which has been widely accepted by geological cir-cles.However,can the components migrate in solid rocks under the action of stress?This problem has long attracted the attention of many geologists.The author has made a series of simulating experi-tal system.The results show that mylonite was produced as a consequence of shear flowage deformation in dolomite .The contents of the elements in porphyryclast and flowage deformation matrix were analyzed by means of electron microprobe,and the results show that Pb and Zn are obvi-ously concentrated in the shear-flowage deformation zone.Furthermore,Pb is concentrated in its cen-tre while Zn is concentrated on its margins .Newly-formed micro-grained galena can be detected locally.     

Effects of temperature and pressure differences on water seepage in breccia

Effects of temperature and pressure differences on water seepage in breccia were investigated by using the physicochemical seepage instrument.The results show that the relationship of flow and pressure differences can be expressed by a linear equation,and the seepage coefficient is linearly correlated with temperature.The relation-ship between seepage flow and temperature could be described with the linear equation.The constant and tempera-ture seepage coefficient showed a linear relation with pressure.Binary quantitative equation for the seepage flow,temperature and pressure was obtained,and explained with experimental data and theoretical analysis.     

http://www.sciencedirect.com/science/article/pii/S1674987110000058

<正>The thermo-electric coefficients of twenty-six magnetite samples,formed either by magmatism or metamorphism,were tested by the thermo-electric instrument BHET—06.Results showed that the coefficient is of a constant value of about -0.05 mV/℃.It is emphasized that because every magnetite grain was tested randomly,the coefficient is independent of the crystallographic direction.This fact means the thermal voltage generated from a single magnetite crystal can be accumulated,and as a result a new thermo-electric field can arise when a gradient thermal field exists and is active within the earth's crust.Because magnetite is widespread in the earth's crust(generally appearing more in the middle-lower crust),there is more-than-random probability that the additional thermo-electric field can be generated when certain thermal conditions are fulfilled.We,therefore,used the thermo-electric effect of magnetite to study the mechanism responsible for the presence of abnormal geo-electric fields during earthquake formation and occurrence, because gradient thermal fields always exist before earthquakes.The possible presence of additional thermo-electric fields was calculated under theoretical seismological conditions,using the following calcu-lation formula:E= - 0.159(σ×△T×φ×ρ_2×[(h~2-2x~2)cosα+ 3hxsinα]/ρ_1(h~2 +x~2)~(5/2)).In the above formula,σis thermo-electric coefficient of magnetite,△T is the temperature difference acting on it,φis a sectional area on a block of magnetite vertically perpendicular to the direction of the thermal current.ρ_1 andρ_2 are the respective resistivities of magnetite and the crust,and h,α,and x,respectively,h is the depth of embedded magnetite block,αmeans the angle created by the horizontal line and ligature of the two poles of magnetite block,and x is the distance from observation point to projective center point of the magnetite block on earth surface.According to simulations calculated with this formula,additional thermo-electric field intensity may reach as high as n to n×10~2 mV/km.This field is strong enough to cause obvious anomalies in the background geo-electric field,and can be easy probed by earthquake monitoring equipment. Therefore,we hypothesize that geo-electric abnormalities which occur during earthquakes may be caused by the thermo-electric effect of magnetite.     

Importance of measurement of rare-earth elements from the Soil Information and Monitoring System： A normalization method for detection of chromium contamination in soils

The Soil Information and Monitoring System （TIM） is an independent subsystem of the integrated Environmental Information and Monitoring System in Hungary. Based on physiographical-soil-ecological units 1202 representative observation points were selected and exactly defined by geographical coordinates using GPS. The first sampling was done in 1992. Some soil parameters are measured every year, some others every 3 years or every 6 years depending on their changeability （stability）. Toxic element, heavy metal and rare earth element contents in soil samples were determined after hot and strong acid digestion by ICP-OES equipment. The rare earth elements were measured in TIM soil samples and some of them can be found in a relatively high amount, e.g. cerium and neodymium concentration are 3-4 times higher than the very often-measured copper concentration. Close correlation was found among some of the rare earth element concentrations. The Y-Yb correlation is the strongest, its r-value is higher than 0.95. The reason of difference among the yttrium concentrations of samples is not the deviation of ICP-OES measurement that is much smaller than the difference among the individual soil samples. Similar, close correlation was obtained between the soil yttrium and chromium concentration. The amount of chromium contamination can be calculated as the distance from the Y-Cr trend line. We call this method yttrium normalization. The principle of this method is described as follows： such elements than chromium and yttrium insist in the soil in a strongly bonded form. Therefore, in a non-contaminated area there is a close correlation between the yttrium and chromium concentrations. If the measured chromium concentration in a sample is not on the Cr-Y trend line than the distance it indicates the amount of antropogenic chromium. Contamination can be detected in that case of course if the contaminant does not contain yttrium.     

Modeling the effect of wetting on the mechanical behavior of crushable granular materials

It is well known that the compressibility of crushable granular materials increases with the moisture content,due to the decrease of particle strength in a humid environment.An existing approach to take into account the effect of grain breakage in constitutive modeling consists in linking the evolution of the grain size distribution to the plastic work.But how the material humidity can affect this relationship is not clear,and experimental evidence is quite scarce.Based on compression tests on dry and saturated crushable sand recently reported by the present authors,a new non-linear relationship is proposed between the amount of particle breakage and the plastic work.The expression contains two parameters:(1)a material constant dependent on the grain characteristics and(2)a constant depending on the wetting condition(in this study,dry or saturated).A key finding is that the relationship does not depend on the stress path and,for a given wetting condition,only one set of parameters is necessary to reproduce the results of isotropic,oedometric,and triaxial compression tests.The relationship has been introduced into an elastoplastic constitutive model based on the critical state concept with a double yield surface for plastic sliding and compression.The breakage ratio is introduced into the expression of the elastic stiffness,the critical state line and the hardening compression pressure.Incremental stress-strain computations with the model allow the plastic work to be calculated and,therefore,the evolution of particle crushing can be predicted through the proposed non-linear relationship and reintroduced into the constitutive equations.Accurate predictions of the experimental results in terms of both stress-strain relationships and breakage ratio were obtained.