The utilization of the double focal transformation for sparse data representation and data reconstruction

In many cases, seismic measurements are coarsely sampled in at least one dimension. This leads to aliasing artefacts and therefore to problems in the subsequent processing steps. To avoid this, seismic data reconstruction can be applied in advance. The success and reliability of reconstruction methods are dependent on the assumptions they make on the data. In many cases, wavefields are assumed to (locally) have a linear space–time behaviour. However, field data are usually complex, with strongly curved events. Therefore, in this paper, we propose the double focal transformation as an efficient way for complex data reconstruction. Hereby, wavefield propagation is formulated as a transformation, where one‐way propagation operators are used as its basis functions. These wavefield operators can be based on a macro velocity model, which allows our method to use prior information in order to make the data decomposition more effective. The basic principle of the double focal transformation is to focus seismic energy along source and receiver coordinates simultaneously. The seismic data are represented by a number of localized events in the focal domain, whereas aliasing noise spreads out. By imposing a sparse solution in the focal domain, aliasing noise is suppressed, and data reconstruction beyond aliasing is achieved. To facilitate the process, only a few effective depth levels need to be included, preferably along the major boundaries in the data, from which the propagation operators can be calculated. Results on 2D and 3D synthetic data illustrate the method's virtues. Furthermore, seismic data reconstruction on a 2D field dataset with gaps and aliased source spacing demonstrates the strength of the double focal transformation, particularly for near‐offset reflections with strong curvature and for diffractions.     

基于截断牛顿法的VTI介质声波多参数全波形反演

不同类别参数间的相互耦合使多参数地震全波形反演的非线性程度显著增加,地震波速度与各向异性参数取值数量级的巨大差异也会使反演问题的性态变差.合理使用Hessian逆算子可以减弱这两类问题对反演的影响,提高多参数反演的精度,而截断牛顿法是一种可以比较准确地估计Hessian逆算子的优化方法.本文采用截断牛顿法在时间域进行了VTI介质的声波双参数同时反演的研究.不同模型的反演试验表明,在VTI介质声波双参数同时反演中,截断牛顿法比有限内存BFGS(Limited-memory Broyden-Fletcher-Goldfarb-Shanno,L-BFGS)法能更准确地估计Hessian逆算子,进而较好地平衡两类不同参数的同时更新,得到了比较精确的反演结果.     

顾及多方向观测值权比反演地球重力场的动力积分法

考虑到不同坐标系下各个方向观测值对反演地球重力场的频谱贡献不同,建立了顾及多方向观测值权比的动力积分法,并利用该方法反演了高精度的GOCE HL-SST卫星重力场模型.首先,分析了不同坐标系下各个方向观测值与地球重力场信息的响应关系,其中惯性系(IRF)下X、Z方向的观测值分别对扇谐系数、带谐系数最为敏感,Z方向的解算精度在全频段均略高于X、Y方向;地固系(EFRF)下各个方向的独立解算精度均与能量守恒法的解算精度相当;局部指北坐标系(LNOF)下X、Z和Y三个方向的解算精度依次递减,且Y方向在47阶附近有明显"驼峰"现象.其次,比较了不同坐标系下顾及三个方向观测值权比的加权解算模型,其中加权联合解算模型精度在20至70阶次均明显优于等权解算模型,在带谐项和共振阶次精度提升明显,且LNOF下的加权联合解算精度要优于IRF和EFRF.最后,比较了GOCE和CHAMP卫星的模型解算精度,采用本文计算方法,仅利用2个月GOCE轨道观测值解算的模型精度优于包含更长观测时段信息的AIUB-CHAMP01S和EIGEN-CHAMP03S模型,且略优于ASU-GOCE-2months模型.     

地震波频散效应与反Q滤波相位补偿

本文以Futterman提出的地震波振幅衰减和相速度频散表达式为基础,从井震匹配的角度出发,推导了不进行反Q滤波、仅反Q滤波振幅补偿、仅反Q滤波相位补偿,以及反Q滤波相位与振幅同时补偿四种情况下,地震波速度频散与相速度及地震记录振幅谱间的表达式,并从理论上说明了反Q滤波相位补偿的必要性.通过相同观测系统,炸药震源和可控震源分别激发采集的零偏移距VSP资料实例,验证了本文所推导的速度频散表达式的合理性;通过井震标定实例,进一步说明了反Q滤波相位补偿,可有效消除地震子波速度频散,提升地面地震资料与零偏移距VSP走廊叠加剖面的匹配度,最终提高地震资料成果的可信度.     

Learning of pore pressure response and dynamic soil behavior from downhole array measurements

Downhole arrays are deployed to measure motions at the ground surface and within the soil profile, with some arrays instrumented to also record the pore pressure response within soft soil profiles during excitation. The measurements from these arrays have typically been used in conjunction with parametric and nonparametric inverse analysis approaches to identify soil constitutive model parameters for use in site response analysis or to identify averaged soil behavior between locations of measurement. The self-learning simulations (SelfSim) inverse analysis framework, previously developed and applied under total stress conditions, is extended to effective stress considerations and is employed to reproduce the measured motions and pore pressures from downhole arrays while extracting the underlying soil behavior and pore pressure response of individual soil layers. SelfSim is applied to the 1987 recordings from the Imperial Valley Wildlife Liquefaction Array. The extracted soil behavior suggests a new functional form for modeling the degradation of the shear modulus with respect to excess pore pressures. The extracted pore pressure response is dependent on the number and amplitude of shear strain cycles and has a functional form similar to current strain-based pore pressure generation models.     

Two-dimensional physical-based inversion of confined and unconfined aquifers under unknown boundary conditions

An inverse method is developed to simultaneously estimate multiple hydraulic conductivities, source/sink strengths, and boundary conditions, for two-dimensional confined and unconfined aquifers under non-pumping or pumping conditions. The method incorporates noisy observed data (hydraulic heads, groundwater fluxes, or well rates) at measurement locations. With a set of hybrid formulations, given sufficient measurement data, the method yields well-posed systems of equations that can be solved efficiently via nonlinear optimization. The solution is stable when measurement errors are increased. The method is successfully tested on problems with regular and irregular geometries, different heterogeneity patterns and variances (maximum K_max/K_min tested is 10,000), and error magnitudes. Under non-pumping conditions, when error-free observed data are used, the estimated conductivities and recharge rates are accurate within 8% of the true values. When data contain increasing errors, the estimated parameters become less accurate, as expected. For problems where the underlying parameter variation is unknown, equivalent conductivities and average recharge rates can be estimated. Under pumping (and/or injection) conditions, a hybrid formulation is developed to address these local source/sink effects, while different types of boundary conditions can also exert significant influences on drawdowns. Local grid refinement near wells is not needed to obtain accurate results, thus inversion is successful with coarse inverse grids, leading to high computation efficiency. Furthermore, flux measurements are not needed for the inversion to succeed; data requirement of the method is thus not much different from that of interpreting classic well tests. Finally, inversion accuracy is not sensitive to the degree of nonlinearity of the flow equations. Performance of the inverse method for confined and unconfined aquifer problems is similar in terms of the accuracy of the estimated parameters, the recovered head fields, and the solver speed.     

声波方程频率域有限元参数反演

推导出频率域有限元声波正演方程,为了消除边界反射,将Clayton-Engquist旁轴波动方程吸收边界条件引入频率域,并对有限元刚度矩阵和质量矩阵进行压缩存储,利用广义共轭梯度法求解有限元方程获得正演解.在此基础上,推导出在某一频率下波场数据残差δU与单元物性参数修改量δλ之间关系的Jacobi矩阵,反演方法允许利用地面二维炮集全波场资料与给出初始模型参数的正演值的差值δU,迭代求得δλ.由于计算机内存的限制,方法计算不允许有过多数目的未知数个数,因此还提出了对同一介质物性单元的Jacobi矩阵元素进行压缩组装的措施,从而使反演的未知量个数减少,结合采用共轭梯度迭代法,使得只需利用有效波频段的少数一些频率即可进行迭代反演.正演和反演理论模型的数值模拟结果表明方法是有效的.     

Solution of the inverse geodetic problem for navigation

Abstract

Simplified geodetic formulas for very long and medium geodesies used in the sea and air navigation as well as methodology and formulas for investigating errors involved in these simplifictions for any hyperbolic system are presented.     

CMT data inversion using a Bayesian information criterion to estimate seismogenic stress fields

We developed an inversion method to estimate the stress fields related to earthquake generation (seismogenic stress fields) from the centroid moment tensors (CMT) of seismic events by using Akaike's Bayesian information criterion (ABIC). On the idea that the occurrence of an earthquake releases some part of the seismogenic stress field around its hypocentre, we define the CMT of a seismic event by a weighted volume integral of the true but unknown seismogenic stress field. Representing each component of the seismogenic stress field by the superposition of a finite number of 3-D basis functions (tri-cubic B-splines), we obtain a set of linear observation equations to be solved for the expansion coefficients (model parameters). We introduce prior constraint on the roughness of the seismogenic stress field and combine it with observed data to construct a Bayesian model with hierarchic, highly flexible structure controlled by hyper-parameters. The optimum values of the hyper-parameters are objectively determined form observed data by using ABIC. Given the optimum values of the hyper-parameters, we can obtain the best estimates of model parameters by using a maximum likelihood algorithm. We tested the validity of the inversion method through numerical experiments on two synthetic CMT data sets, assuming the distribution of fault orientations to be aligned with the maximum shear stress plane in one case and to be random in the other case. Then we applied the inversion method to actual CMT data in northeast Japan, and obtained the pattern of the seismogenic stress field consistent with geophysical and geological observations.     

Numerical inversion of deformation caused by pressure sources: application to Mount Etna (Italy)

The interpretation of geodetic data in volcanic areas is usually based on analytical deformation models. Although numerical finite element (FE) modelling allows realistic features such as topography and crustal heterogeneities to be included, the technique is not computationally convenient for solving inverse problems using classical methods. In this paper, we develop a general tool to perform inversions of geodetic data by means of 3-D FE models. The forward model is a library of numerical displacement solutions, where each entry of the library is the surface displacement due to a single stress component applied to an element of the grid. The final solution is a weighted combination of the six stress components applied to a single element-source. The pre-computed forward models are implemented in a global search algorithm, followed by an appraisal of the sampled solutions. After providing extended testing, we apply the method to model the 1993–1997 inflation phase at Mt Etna, documented by GPS and EDM measurements. We consider four different forward libraries, computed in models characterized by homogeneous/heterogeneous medium and flat/topographic free surface. Our results suggest that the elastic heterogeneities of the medium can significantly alter the position of the inferred source, while the topography has minor effect.