改进的约束变密度界面反演策略及其应用

计算密度分界面的起伏变化在区域地质构造研究和石油矿产资源勘探中具有重要意义.已有密度界面反演方法更多侧重约束变密度界面反演算法,而对约束信息的准确性、研究区横向密度变化往往考虑不足,影响了最终反演结果的可信度.本文在变密度界面正反演算法基础之上,结合实际需求,提出已知深度信息约束下的变密度界面反演策略.该策略主要包括变密度约束反演算法、已知深度约束信息校验和分区变密度模型三个方面.其中反演算法提供了带已知深度约束信息的密度界面迭代反演方法;约束信息校验用于评估约束信息精度,通过调差降低约束信息的系统误差;在反演过程中引入水平密度分区以应对不同构造背景密度界面模型,提高反演结果的可信度.最后将本策略应用于南海莫霍面深度反演计算中,结果显示借助已知约束信息,利用分区密度模型能够获得更为可信的深度反演结果,验证了该策略的正确性.     

约束变密度界面反演法在渤海深部结构研究中的应用

在分析研究约束变密度界面反演基础理论及其方法应用的基础上,本文以我国渤海深部结构为研究对象,利用变密度反演法,对渤海海域主要坳陷进行沉积厚度的反演及其与断层等的综合解释工作,而后对地壳与地幔之间存在的重要分界面-莫霍面,基于已有的对莫霍面信息有较为系统、整体性反应的研究区重力数据,对其进行变密度反演计算,获得了研究区莫霍面的深度分布情况.通过与已知的地质地球物理资料,主要地震OBS测深等前人研究成果相互对比验证,取得了较为良好的效果.在此基础上,进一步证明约束变密度界面反演法在沉积界面反演、深部构造研究等方面,可以为其他勘探区域的研究提供较为准确可靠的方法参考,从而丰富发展地球物理反演工作.     

唐山震区结晶基底磁性界面的初步探讨

扼要地叙述了频率域磁性单界面反演方法的基本原理,运用该程序对唐山震区航磁资料进行了数字处理,获得了震区结晶基底磁性界面的定量特征。结果表明唐山处在华北平原北部与燕山隆起毗邻的褶皱拗陷中,震源区开平向斜的底部基底磁质界面为一拗陷,周围受3条深断裂切割围限。文章最后,结合人工地壳测深、地形变、地应力等成果进行了分析研究,探索了唐山地震的发生及孕育过程     

重力密度界面反演方法研究进展

重力位场的界面反演是位场处理解释中的重要问题.本文针对重力密度界面的反演问题,较全面的介绍了一些有代表性的方法,包括空间域反演方法和频率域反演方法.重点阐述了当前占主流的频率域中的Parker-Oldenburg迭代方法,在此基础上讨论了重力密度界面反演目前存在的问题和以后的发展方向.     

利用地震回折波资料反演界面位置与速度分布

论述了利用地震回折波资料反演界面位置与速度分布的方法，推导了地震波走时对于界面位置偏导数的计算公式。数值模拟和实测资料的计算结果表明了该方法的有效性和编制的计算程序的实用性。该方法最突出的特点是充分地利用了透射波资料中所含的界面位置的信息。界面位置的分辨率与界面两边的速度反差有关，速度差别越大，则分辨率越高。     

基于深度学习U-net网络的重力数据界面反演方法

重力数据的密度界面反演是位场数据解释中的一项主要工作,在区域构造演化、深部莫霍面确定等领域的研究中发挥重要作用.近年来,数据驱动的深度学习方法广泛地应用在地球物理数据处理与反演中,本文提出一种基于深度学习U-net网络的重力数据密度界面反演方法.首先,对半椭球体界面模型进行随机抽取和组合进而形成地下起伏界面数据集,并基于Parker正演理论对界面数据集进行重力异常正演计算,为深度学习网络模型的训练提供特征完备的数据源;其次,设计了基于U-net网络模型的深度学习界面反演算法,在传统的损失函数基础上增加光滑损失项和过拟合抑制项,提高重力界面反演结果的光滑性和收敛效率;最后通过测试样本集进行反演预测,验证建立深度学习网络模型的泛化性.本文通过理论模型和实际数据试验分析了本文方法在密度界面反演中的有效性和实用性,基于改进损失函数约束的深度学习界面反演方法有效地提高了密度界面反演的收敛效率和计算稳定性.     

反演磁性地质界面的意义与方法

磁性地质界面的反演是地球物理勘探的重要课题之一,对于提高磁测资料的地质解释效果有重要意义.通过磁测资料反演磁性地质界面主要分两步：第一步,从总异常中分离有磁性界面起伏引起的磁异常;第二步,由分离的磁异常反演磁性界面的深度.本文分别从这两个方面回顾了提取异常和反演界面的较有代表性的技术.     

三维密度界面的正反演研究和应用

重力位场的界面反演是位场处理解释中的重要问题.本文将基于快速傅里叶变换的频率域界面反演方法Parker-Oldenburg公式推广到物性可随深度变化的三维情况,得出了密度可以横向、纵向任意变化的重力界面正反演公式.该方法在计算时可以合理地选取地面下某一深度作为基准面以减小界面起伏,使迭代易于收敛.理论模型试验表明该方法反演精度高,收敛速度快,在密度界面反演中具有广泛的实用价值.最后利用该方法反演华北地区莫霍面的深度,反演结果得到了地震测深数据的验证.     

具有已知深度点的二维单一密度界面的反演

密度界面的反演是地球物理勘探的重要课题之一,对研究盆地、区域构造和深部构造有重要意义.本文采用一种密度界面上已知若干控制点深度的迭代反演法,对二维单一密度界面的起伏模型进行计算及讨论,并对大港地区的实测剖面进行反演,取得了较为理想的效果.     

三维密度界面的正反演研究和应用

重力位场的界面反演是位场处理解释中的重要问题.本文将基于快速傅里叶变换的频率域界面反演方法Parker-Oldenburg公式推广到物性可随深度变化的三维情况,得出了密度可以横向、纵向任意变化的重力界面正反演公式.该方法在计算时可以合理地选取地面下某一深度作为基准面以减小界面起伏,使迭代易于收敛.理论模型试验表明该方法反演精度高,收敛速度快,在密度界面反演中具有广泛的实用价值.最后利用该方法反演华北地区莫霍面的深度,反演结果得到了地震测深数据的验证.     

Joint solution of leveling and gravity data in crustal deformation monitoring

ＪｏｉｎｔｓｏｌｕｔｉｏｎｏｆｌｅｖｅｌｉｎｇａｎｄｇｒａｖｉｔｙｄａｔａｉｎｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｍｏｎｉｔｏｒｉｎｇＺＵＳＨＥＮＧＺＨＡＮＧ１）（张祖胜）ＹＵＡＮＸＩＹＡＮＧ２）（杨元喜）ＨＡＮＲＯＮＧＳＵＮ１）（孙汉荣）Ｚ...     

Asymptotic formulae for correcting finite element predicted natural frequencies of gravity and embankment dams

Extraction of natural frequencies of a gravity dam or an embankment dam plays an important role in the seismic design of the dam because the seismic response of a dam is dependent largely on the dynamic characteristics of the dam. Owing to the lack of exact solutions and the geometry of a dam, numerical methods such as finite element methods have been often used to extract the natural frequencies of the dam. Since the finite element method is an approximate one, the resulting finite element solution to the natural frequency of a dam cannot be safely used unless its accuracy is evaluated within the acceptable range for the seismic design of the dam. To solve this problem, some asymptotic formulae for correcting the finite element predicted natural frequencies of a gravity dam and an embankment dam have been developed in this paper. Since the present asymptotic formulae are derived from the fact that the finite element solution tends to the exact one if the finite element size used approaches zero, they provide a corrected solution of higher accuracy for the natural frequency of a dam so that the accuracy of a finite element solution can be evaluated against this corrected solution. After the correctness and usefulness of the present formulae are assessed, two practical examples have been given to show how the asymptotic formulae can be used to correct and evaluate the discretization error for the finite element predicted natural frequencies of gravity dams and embankment dams.     

Numerical modeling of wave equation by a truncated high-order finite-difference method

Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations.However, the computation cost generally increases linearly with increased order of accuracy.Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note th...     

Analytical and numerical solutions for wave propagation in water-saturated porous layered half-space

A fluid-saturated one-layer continuum underlain by a rigid half-space is considered. An exact solution is developed in frequency domain for analyzing disturbance induced by a strip footing located at the surface with vertical harmonic excitation. Since the analytical solution can be used only for very simple conditions, a finite element model has been developed also and compared with the exact solution. It is shown that finite element results are in close agreement with the results which have been obtained by a transformation technique. The proposed finite element scheme can take into account the complex geometry and inhomogeneity for practical problems. Besides this, the analytical results exhibit the overall characteristic of wave propagation in porous media and will provide a representative test problem which can be used for a quantitative evaluation of the accuracy of various numerical solution methods.     

A Gravimetric Geoid Computation and Comparison with GPS Results in Northern Andalusia (Spain)

Two new GPS surveys have been carried out to check the accuracy of an existing gravimetric geoid in a test area located in northern Andalusia (Spain). The fast collocation method and the remove-restore procedure have been used for the computation of the quasigeoid model. The Spanish height system is based on orthometric heights, so the gravimetrically determined quasigeoid has been transformed to a geoid model and then compared to geoid undulations provided by GPS and levelling at benchmarks belonging to the Spanish first-order levelling network. The discrepancies between the gravimetric solution and GPS/levelling undulations amount to ±2 cm for one survey and ±5 cm for another after fitting a plane to the geoid model.     

Improvements and Corrections to AT123D Code

Although based on exact analytical solutions, semi‐analytical solute transport models can have significant numerical error in applications with high frequency oscillatory source terms and when parameter value combinations cause series solution approximations to converge slowly. Methods for correcting these numerical errors are presented and implemented in the AT123D code, which employs Green's functions to represent point, linear, and rectangular prismatic source zones. In order to increase its computational accuracy, a Romberg numerical integration scheme was added to AT123D with prespecified error criteria, variable time stepping, and partitioning of the integral to handle rapidly changing source terms. More rapidly converging series solution approximations for the Green's functions were also incorporated to improve both accuracy and computational efficiency for finite‐depth aquifers. AT123D also has been modified to eliminate redundant calculations at points where approximate steady‐state conditions have been reached to improve computational efficiency during numerical integration. These modifications help to decrease computer run times that can be excessive for three‐dimensional problems with large numbers of computational points, small time steps, and/or long simulation time periods. Errors in the original AT123D code also were corrected in this modified version, AT123D‐AT, in order to accurately simulate finite‐duration (pulse) source releases.     

Analytical study on accelerating falling of non-spherical particle in viscous fluid

Unsteady motion of a vertically falling non-spherical particle has attracted considerable attention due to its frequent applications in nature and industry.A series of semi-analytical methods have been used to raise the results＇ accuracy as well as widening the region of convergence.The current study pursued a new analytical solution for the unsteady motion of a rigid non-spherical particle in a quiescent Newtonian fluid,based on the Optimal Homotopy Analysis Method.With a view towards obtaining the highest level of accuracy and ensuring the convergence of the analytical results,the averaged residual errors were obtained and minimized.In addition to flexibility,it was also proven that the proposed method can lead to completely reliable and precisely accurate results.Based on the series solution,the effects of physical parameters on the terminal settling velocity（i.e.the greatest velocity that a falling body may reach）and the acceleration time（i.e.the time that a particle reaches the settling velocity） are investigated.     

瞬变电磁场时域格林函数解

近源时间域电磁场具有信号强、探测深度大、精度高等优点,但传统勘探电磁场理论中偶极子近似在近源会引起较大误差,导致这一优势的发挥受到了制约.开展直接时间域电磁场解析式研究,是解决这一问题的途径之一.本文提出在点电荷微元假设下,引入时域格林函数,求取瞬变电磁场时间域解析解.采用积分运算法,把电磁场阻尼波动方程的求解问题转化为求其格林函数积分形式解的问题;建立辅助路径解决奇点问题,利用复分析中的约当引理、留数定理和广义函数等理论和方法,推导计算出时间域格林函数的时空四重广义积分.得到达朗贝尔方程的直接时域格林函数精确解析式,与传统方法“比拟”出的公式具有相同的形式,验证了本文推导的时域格林函数解析公式的正确性;推导出扩散方程的直接时间域解析解.通过与时变点电荷源时间域的电磁响应近似表达式进行对比,得出本文所推导的公式计算精度较高的结论;建立了全空间回线源瞬变电磁场问题的直接时间域求解公式.为解决全场区瞬变电磁场精细探测直接时域解析问题提供了基础理论.     

Application of a nonlinear slug test model

Knowledge of the hydraulic conductivity distribution is of utmost importance in understanding the dynamics of an aquifer and in planning the consequences of any action taken upon that aquifer. Slug tests have been used extensively to measure hydraulic conductivity in the last 50 years since Hvorslev's (1951) work. A general nonlinear model based on the Navier-Stokes equation, nonlinear frictional loss, non-Darcian flow, acceleration effects, radius changes in the wellbore, and a Hvorslev model for the aquifer has been implemented in this work. The nonlinear model has three parameters: beta, which is related primarily to radius changes in the water column; A, which is related to the nonlinear head losses; and K, the hydraulic conductivity. An additional parameter has been added representing the initial velocity of the water column at slug initiation and is incorporated into an analytical solution to generate the first time step before a sequential numerical solution generates the remainder of the time solution. Corrections are made to the model output for acceleration before it is compared to the experimental data. Sensitivity analysis and least squares fitting are used to estimate the aquifer parameters and produce some diagnostic results, which indicate the accuracy of the fit. Finally, an example of field data has been presented to illustrate the application of the model to data sets that exhibit nonlinear behavior. Multiple slug tests should be taken at a given location to test for nonlinear effects and to determine repeatability.     

Influence of soil to structure stiffness on the accuracy of the pushover method for underground structures

The pushover method for underground structures is a seismic analysis method featured by high calculation accuracy and a simple implementation process. The method has been widely used in seismic design and other related scientific research; however, the influence of different soil-structure flexibility ratios on the accuracy of this method is still not well understood. In this study, we select the cross-section structures beneath the Daikai subway station as the research object and establish 12 finite element analysis models with different soil-structure flexibility ratios using ABAQUS. All models are computed by the dynamic time-history method or the pushover method. Furthermore, the dynamic time-history solution result is taken as the standard solution, and the precision and application of the pushover analysis method are discussed based on the parameters of peak interlayer displacement and peak internal force of the middle column section. The results show that the soil-structure flexibility ratio has a significant influence on the calculation accuracy of the pushover method, and the calculation accuracy of this method is the most ideal when the soil-structure flexibility is equal to 1. The research results can provide significant references for the seismic design of underground structures or the improvement of simplified seismic analysis methods.