地球物理复杂信号的多尺度熵分析方法

为了研究地球物理观测信号在多尺度空间中的复杂性,本文提出了一种新的广义熵谱（General Entropy Spectral）概念,其在零频率下的值正好给出了传统意义的信息熵。同时,还提出了一种新的局部多尺度熵（Local Multiscale Entropy）分析方法,该方法结合了信息熵和多尺度方法,可用于对数据特征、信息复杂度进行多尺度表达。上述这些新方法适用于分析各种复杂的地球物理信号,以便深入挖掘更多的信号特征和信息。     

A Robust Two-Step Inversion of Complex Magnetotelluric Apparent Resistivity Data

Though two-dimensional inversion is now a standard procedure, interpretation of magnetotelluric (MT) data under the assumption of isotropic and one-dimensional structures is a valuable procedure for a first step interpretation of exploratory and solid earth geophysics investigation data. Because its interpretation requires an efficient inverse modelling, we propose and evaluate an inversion procedure, which consists of two steps. Both steps employ jointly the modulus and the phase of the apparent resistivity function. The first one consists of the use of the asymptotic Bostick-Niblett approach. The second employs the result of the inversion obtained in the first step as a starting model to initialize the linearized inversion performed using a multiple re-weighted least-squares approach. We applied the analysis both to synthetic data and to field data from the Parana Basin in Brazil. The results show that the inversion procedure presents a faster convergence without loss of accuracy, increases the resolving power of the MT technique, and may improve its capability to delineate conductors up to a depth of one hundred kilometers. Therefore a reasonable interpretation of the data employing one-dimensional model can be achieved even in the presence of relatively noisy data, and under conditions that slightly violate the premise of lateral homogeneity.     

Parameter Estimation for Groundwater Models under Uncertain Irrigation Data

The success of modeling groundwater is strongly influenced by the accuracy of the model parameters that are used to characterize the subsurface system. However, the presence of uncertainty and possibly bias in groundwater model source/sink terms may lead to biased estimates of model parameters and model predictions when the standard regression‐based inverse modeling techniques are used. This study first quantifies the levels of bias in groundwater model parameters and predictions due to the presence of errors in irrigation data. Then, a new inverse modeling technique called input uncertainty weighted least‐squares (IUWLS) is presented for unbiased estimation of the parameters when pumping and other source/sink data are uncertain. The approach uses the concept of generalized least‐squares method with the weight of the objective function depending on the level of pumping uncertainty and iteratively adjusted during the parameter optimization process. We have conducted both analytical and numerical experiments, using irrigation pumping data from the Republican River Basin in Nebraska, to evaluate the performance of ordinary least‐squares (OLS) and IUWLS calibration methods under different levels of uncertainty of irrigation data and calibration conditions. The result from the OLS method shows the presence of statistically significant (p < 0.05) bias in estimated parameters and model predictions that persist despite calibrating the models to different calibration data and sample sizes. However, by directly accounting for the irrigation pumping uncertainties during the calibration procedures, the proposed IUWLS is able to minimize the bias effectively without adding significant computational burden to the calibration processes.     

Groundwater Monitoring at the Facilities of Oil and Gas Complex

The environmental aspects of the development of oil and gas deposits are considered. The methodological principles of the formation of a hydrogeological monitoring system in the general structure of comprehensive geoenvironmental monitoring of oil-and-gas facilities are discussed. The structure of comprehensive geoenvironmental monitoring and subsurface-hydrosphere monitoring is presented. It consists of observation, estimate, and forecast subsystems. Requirements are formulated to the source data that forms the basis, on which the observational network is designed; the state of the subsurface hydrosphere is assessed both at the stages of designing the facilities being examined and their monitoring system and during the process of operation of the facilities and systems of monitoring. Changes in the environment and subsurface hydrosphere as its component are predicted.     

Complex Isotopic and Hydrogeochemical Studies of Groundwater in the Kirov Region

Uranium-isotopic and multi-element hydrogeochemical methods are used to assess the natural pollution of fresh groundwater in horizons under development due to a modern intrusion of hypogene water with increased boron and fluorine contents in zones of old tectonic dislocations. The total index of the abovestandard (>MAC) pollution of groundwater in the region of the city of Kirov (the Kirov area) reaches 15. To reduce the effect of natural groundwater pollutants on the health of population, it is recommended that the groundwater intake should be regulated, purifiers should be used, and prospects revealed by the isotopic-and-hydrogeochemical data should be explored. The feasibility of modeling the groundwater formation and tracing the neotectonic dislocations in geological platforms is shown.     

利用水位资料反演华北地区构造应力场变化

本文选取华北地区观测条件较好的63口水位观测井资料, 运用小波分析法去除各井水位资料中的短期高频信息, 提取出能反映水位多年动态变化的趋势信息, 并利用各井水位趋势变化数据反演出华北地区多年构造应力场变化图像, 结合其它水文资料, 探讨了华北地区近年来的构造应力场变化特征。     

不同路面条件下地震映像勘探效果的对比研究

讨论了地震映像勘探的主要特点,用实例分析了土质路面、沥青路面等不同介质条件下的勘探效果,明确了这一勘探方法的适用条件.     

A Hierarchical Bayesian Model Averaging Framework for Groundwater Prediction under Uncertainty

Groundwater prediction models are subjected to various sources of uncertainty. This study introduces a hierarchical Bayesian model averaging (HBMA) method to segregate and prioritize sources of uncertainty in a hierarchical structure and conduct BMA for concentration prediction. A BMA tree of models is developed to understand the impact of individual sources of uncertainty and uncertainty propagation to model predictions. HBMA evaluates the relative importance of different modeling propositions at each level in the BMA tree of model weights. The HBMA method is applied to chloride concentration prediction for the “1,500‐foot” sand of the Baton Rouge area, Louisiana from 2005 to 2029. The groundwater head data from 1990 to 2004 is used for model calibration. Four sources of uncertainty are considered and resulted in 180 flow and transport models for concentration prediction. The results show that prediction variances of concentration from uncertain model elements are much higher than the prediction variance from uncertain model parameters. The HBMA method is able to quantify the contributions of individual sources of uncertainty to the total uncertainty.     

Complex Rays and Wave Packets for Decaying Signals in Inhomogeneous,Anisotropic and Anelastic Media

Diffraction and anelasticity problems involving decaying, evanescent or inhomogeneous waves can be studied and modelled using the notion of complex rays. The wavefront or eikonal equation for such waves is in general complex and leads to rays in complex position-slowness space. Initial conditions must be specified in that domain: for example, even for a wave originating in a perfectly elastic region, the ray to a real receiver in a neighbouring anelastic region generally departs from a complex point on the initial-values surface. Complex ray theory is the formal extension of the usual Hamilton equations to complex domains. Liouville's phase-space-incompressibility theorem and Fermat's stationary-time principle are formally unchanged. However, an infinity of paths exists between two fixed points in complex space all of which give the same final slowness, travel time, amplitude, etc. This does not contradict the fact that for a given receiver position there is a unique point on the initial-values surface from which this infinite complex ray family emanates. In perfectly elastic media complex rays are associated with, for example, evanescent waves in the shadow of a caustic. More generally, caustics in anelastic media may lie just outside the real coordinate subspace and one must trace complex rays around the complex caustic in order to obtain accurate waveforms nearby or the turning waves at greater distances into the lit region. The complex extension of the Maslov method for computing such waveforms is described. It uses the complex extension of the Legendre transformation and the extra freedom of complex rays makes pseudocaustics avoidable. There is no need to introduce a Maslov/KMAH index to account for caustics in the geometrical ray approximation, the complex amplitude being generally continuous. Other singular ray problems, such as the strong coupling around acoustic axes in anisotropic media, may also be addressed using complex rays. Complex rays are insightful and practical for simple models (e.g. homogeneous layers). For more complicated numerical work, though, it would be desirable to confine attention to real position coordinates. Furthermore, anelasticity implies dispersion so that complex rays are generally frequency dependent. The concept of group velocity as the velocity of a spatial or temporal maximum of a narrow-band wave packet does lead to real ray/Hamilton equations. However, envelope-maximum tracking does not itself yield enough information to compute synthetic seismograms. For anelasticity which is weak in certain precise senses, one can set up a theory of real, dispersive wave-packet tracking suitable for synthetic seismogram calculations in linearly visco-elastic media. The seismologically-accepiable constant-Q rheology of Liu et al. (1976), for example, satisfies the requirements of this wave-packet theory, which is adapted from electromagnetics and presented as a reasonable physical and mathematical basis for ray modelling in inhomogeneous, anisotropic, anelastic media. Dispersion means that one may need to do more work than for elastic media. However, one can envisage perturbation analyses based on the ray theory presented here, as well as extensions like Maslov's which are based on the Hamiltonian properties.     

Complex Rays and Wave Packets for Decaying Signals in Inhomogeneous,Anisotropic and Anelastic Media

Diffraction and anelasticity problems involving decaying, “evanescent” or “inhomogeneous” waves can be studied and modelled using the notion of “complex rays”. The wavefront or “eikonal” equation for such waves is in general complex and leads to rays in complex position-slowness space. Initial conditions must be specified in that domain: for example, even for a wave originating in a perfectly elastic region, the ray to a real receiver in a neighbouring anelastic region generally departs from a complex point on the initial-values surface. Complex ray theory is the formal extension of the usual Hamilton equations to complex domains. Liouville's phase-space-incompressibility theorem and Fermat's stationary-time principle are formally unchanged. However, an infinity of paths exists between two fixed points in complex space all of which give the same final slowness, travel time, amplitude, etc. This does not contradict the fact that for a given receiver position there is a unique point on the initial-values surface from which this infinite complex ray family emanates.In perfectly elastic media complex rays are associated with, for example, evanescent waves in the shadow of a caustic. More generally, caustics in anelastic media may lie just outside the real coordinate subspace and one must trace complex rays around the complex caustic in order to obtain accurate waveforms nearby or the turning waves at greater distances into the lit region. The complex extension of the Maslov method for computing such waveforms is described. It uses the complex extension of the Legendre transformation and the extra freedom of complex rays makes pseudocaustics avoidable. There is no need to introduce a Maslov/KMAH index to account for caustics in the geometrical ray approximation, the complex amplitude being generally continuous. Other singular ray problems, such as the strong coupling around acoustic axes in anisotropic media, may also be addressed using complex rays.Complex rays are insightful and practical for simple models (e.g. homogeneous layers). For more complicated numerical work, though, it would be desirable to confine attention to real position coordinates. Furthermore, anelasticity implies dispersion so that complex rays are generally frequency dependent. The concept of group velocity as the velocity of a spatial or temporal maximum of a narrow-band wave packet does lead to real ray/Hamilton equations. However, envelope-maximum tracking does not itself yield enough information to compute synthetic seismogramsFor anelasticity which is weak in certain precise senses, one can set up a theory of real, dispersive wave-packet tracking suitable for synthetic seismogram calculations in linearly visco-elastic media. The seismologically-accepiable constant-Q rheology of Liu et al. (1976), for example, satisfies the requirements of this wave-packet theory, which is adapted from electromagnetics and presented as a reasonable physical and mathematical basis for ray modelling in inhomogeneous, anisotropic, anelastic media. Dispersion means that one may need to do more work than for elastic media. However, one can envisage perturbation analyses based on the ray theory presented here, as well as extensions like Maslov's which are based on the Hamiltonian properties.     

The Geochemical Prerequisites for the Formation of High-Carbonate Inversion Groundwater with Reduced Dissolved Solids Content in Deep Horizons of Oil-and-Gas Bearing Structures

Synthesis of empirical natural materials and thermodynamic computer modeling of geochemical processes in water–rock systems at different boundary conditions (solid-to-liquid ratio, , T) were used to determine the genetic causes of the inverse geochemical zonality that forms in deep horizons of oil-and-gas bearing structures. The geochemical pattern of inversion water was found to form chiefly because of changes in the Eh–pH-conditions of the original groundwater under the effect of organic components of rocks and because of an increase in temperature to 100°C at low values of solid-to-liquid ratios and at no higher than 10^–2 bar.     

不同电极阵列联合反演在古墓探测中的应用

电阻率层析成像是一种广泛应用在水文、考古和地质等浅地表勘探领域的地球物理方法。为了增强电阻率层析成像的分辨率、应对复杂的地质问题,本文提出基于雅可比矩阵的不同电极阵列直流电阻率数据的加权联合反演算法,并以温纳和偶极-偶极电极阵列数据为例,在理论模型和古墓探测的野外实例中测试该算法的有效性。结果表明,加权联合反演结果的横向和纵向分辨率都优于单一电极阵列的反演结果,并在实例中缓解“U形”电极阵列的固有缺陷、减少反演模糊性、更好地约束墓室宽度的反演结果。      

利用华北地区承压井水位动态反演含水层体应变

The change of the confined aquifer level reflects the pore pressure change, and the pore pressure change of the aquifer is closely related to the aquifer pressure. This paper uses the tidal response of the well water level data in the North China region to calculate the tidal factor of each well and extract the effective water trend information. Then, the volumetric strain of an existing confined aquifer well in the North China region is inverted, and the contour maps are plotted on a half-year scale from 2009 to 2012. Results show that it can reflect the state of stress and strain in deep crust to a certain extent in the North China region.     

Effect of Different Sampling Methodologies on Measured Methane Concentrations in Groundwater Samples

Analysis of dissolved light hydrocarbon gas concentrations (primarily methane and ethane) in water supply wells is commonly used to establish conditions before and after drilling in areas of shale gas and oil extraction. Several methods are currently used to collect samples for dissolved gas analysis from water supply wells; however, the reliability of results obtained from these methods has not been quantified. This study compares dissolved methane and ethane concentrations measured in groundwater samples collected using three sampling methods employed in pre‐ and post‐drill sampling programs in the Appalachian Basin. These include an open‐system collection method where 40 mL volatile organic analysis (VOA) vials are filled directly while in contact with the atmosphere (Direct‐Fill VOA) and two alternative methods: (1) a semi‐closed system method whereby 40 mL VOA vials are filled while inverted under a head of water (Inverted VOA) and (2) a relatively new (2013) closed system method in which the sample is collected without direct contact with purge water or the atmosphere (IsoFlask^®). This study reveals that, in the absence of effervescence, the difference in methane concentrations between the three sampling methods was relatively small. However, when methane concentrations equaled or exceeded 20 mg/L (the approximate concentration at which effervescence occurs in the study area), IsoFlask^® (closed system) samples yielded significantly higher methane concentrations than Direct‐Fill VOA (open system) samples, and Inverted VOA (semi‐closed system) samples yielded lower concentrations. These results suggest that open and semi‐closed system sample collection methods are adequate for non‐effervescing samples. However, the use of a closed system collection method provides the most accurate means for the measurement of dissolved hydrocarbon gases under all conditions.     

复杂地基条件下桩-土-核岛结构相互作用模型研究

合理有效地模拟桩-土-结构动力相互作用是软土地基条件下核岛厂房结构抗震适应性分析及地基处理的关键环节。以某拟建核岛厂房实际工程为研究背景,结合SuperFLUSH软件平台,以Goodman单元模拟桩与桩周土间的接触效应,采用等价线性法描述近场软土地基非线性特性,并在模型底部和侧面引入黏性边界模拟半无限地基辐射阻尼效应,从而建立土质地基条件下桩-土-核岛结构相互作用分析模型。进而,通过对原状地基和嵌岩桩处理地基条件下核岛厂房的楼层反应谱、结构节点相对位移（绝对值）的对比分析,探讨考虑桩-土间接触效应的嵌岩桩基对核岛厂房结构的影响规律。研究成果可为实际工程中类似土质地基条件下核岛厂房结构的地基处理提供参考。     

Advantages of the Optimum Pulse Moment in Surface NMR and Application in Groundwater Exploration

The surface nuclear magnetic resonance (SNMR) method is widely used in groundwater detection because of its sensitivity to hydrogen in water and direct water detection. However, low signal-to-noise ratios (SNRs) restrict the development of this technique. An optimum pulse sequence is designed according to correspondence between the pulse moment strength and its best detection depth. Because only selection of the pulse intensity distribution according to the target aquifer depth is required and the “on-resonance” pulse pattern is still employed, this pulse sequence emission can be easily achieved using existing SNMR instrumentation. Numerical simulation results and field experiments show that, compared with traditional exponential growth pulses, the optimum pulse sequence effectively improves the SNR of the SNMR method. The aquifer boundary, water content, and pore characteristics of the inversion result are thus more consistent with characteristics of underground structures. Additionally, because the optimum pulse sequence focuses most of the pulse moments in the target depth range, in situations where two aquifers are separated by a relatively narrow aquitard, it is better able to resolve the individual aquifers than the traditional pulses. Optimum pulse moments improve the SNR by enhancing the signal amplitude, compared with various filtering methods, and obtain a better detection effect. This kind of pulse sequence can be used as an alternative pulse sequence form of the SNMR method.