共查询到20条相似文献，搜索用时 71 毫秒
1.
流体因子是储层流体识别的重要方法,而叠前地震反演是获得流体因子的有效途径之一.本文从流体因子的构建出发,基于多孔弹性介质岩石物理模型,建立了流体因子与纵横波模量之间的直接关系,避免了流体因子计算所需的密度参数无法准确求取的问题,通过推导基于纵横波模量的Zeoppritz近似方程及弹性阻抗方程,探讨了基于弹性阻抗的纵横波模量直接反演方法,模型与实际应用表明,基于弹性阻抗的纵横波模量直接反演方法合理、可靠,减少了常规方法间接计算纵横波模量带来的累积误差,基于纵横波模量的流体因子计算方法有较好的实际应用效果. 相似文献
2.
本文采用Novak薄层法推导粘弹性地基的扭转动力阻抗,并将其用于考虑桩土相互作用的单桩扭转动力阻抗,又通过传递矩阵法将此公式推广到求取层状地基单桩扭转动力阻抗。而且本文以工程中常用的端承桩为例,推导了层状地基中端承桩扭转动力阻抗的简化计算公式。根据此公式,分析了频率、上覆软土层厚度、上覆软土层刚度等因素对单桩扭转动力阻抗的影响,分析结果表明,随着振动频率提高,层状地基中端承单桩的扭转动力阻抗的实部有缓慢下降的趋势,而虚部则增大;随着上覆土层厚度的增加,层状地基中端承单桩的扭转动力阻抗的实部和虚部均减小;随着上覆土层刚度的减小,层状地基中端承单桩的扭转动力阻抗的实部减小,虚部在低频段减小,而在较高频率段则增大。 相似文献
3.
A method for obtaining the galvanic distortion matrix is presented so that the regional impedance tensor (free of distortion) is recovered. The method is a constrained stochastic heuristic method, which consists in randomly exploring the space of the distortion parameters. Constraints are imposed on the shortest periods of the regional impedance tensor that, at these short periods, tends to be two dimensional (or one dimensional). Depending on the constraints used, two different methods to recover the regional impedance tensor in this 2D/3D case are presented. Method 1 needs to find the strike of the short periods and Method 2 applies to the measurement directions. Twist, shear and anisotropy parameters are obtained. Thus, the regional impedance tensor is recovered with the only exception being the vertical shift due to the gain, which is equal for all the components of the tensor. Examples with synthetic impedance tensors from 2D/3D models perturbed with galvanic distortion are presented to illustrate how the algorithm works. The presence of noise in data is considered and rules for proceeding are provided. The same examples perturbed by Gaussian noise together with experimental data illustrate the capabilities of the algorithm. 相似文献
4.
基于累积法理论的弹性阻抗参数提取 总被引：1，自引：0，他引：1
弹性阻抗中蕴含着丰富的地下介质信息。从弹性阻抗中准确地提取地下介质信息用于储层预测是十分有意义的,传统的提取方法主要基于最小二乘理论的线性方法。本文引入累积法理论,将其用于弹性阻抗参数的提取,用层状模型和部分Marmousi2模型对该方法进行测试。从背景噪声分别为正态分布和均匀分布的弹性阻抗体中提取纵、横波速度和密度,并将其与最小二乘法提取的纵、横波速度和密度做定量对比分析。结果表明两种方法所得结果精度相当,但累积法具有计算量小、简洁直观的优势,因此,将其用于弹性参数的提取是有效可行的。 相似文献
5.
Stability of the time‐domain analysis method including a frequency‐dependent soil–foundation system 
下载免费PDF全文
下载免费PDF全文 A number of methods have been proposed that utilize the time‐domain transformations of frequency‐dependent dynamic impedance functions to perform a time‐history analysis. Though these methods have been available in literature for a number of years, the methods exhibit stability issues depending on how the model parameters are calibrated. In this study, a novel method is proposed with which the stability of a numerical integration scheme combined with time‐domain representation of a frequency‐dependent dynamic impedance function can be evaluated. The method is verified with three independent recursive parameter models. The proposed method is expected to be a useful tool in evaluating the potential stability issue of a time‐domain analysis before running a full‐fledged nonlinear time‐domain analysis of a soil–structure system in which the dynamic impedance of a soil–foundation system is represented with a recursive parameter model. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
6.
V. M. Pavlov 《Izvestiya Physics of the Solid Earth》2009,45(10):850-860
A new method is proposed for calculating synthetic seismograms caused by a force in a plane-parallel medium consisting of
homogeneous elastic isotropic layers. The matrix impedance, i.e., the matrix function of depth, by which motion vector must
be multiplied in order to obtain the stress vector, is introduced for solving a system of ordinary differential equations
with respect to the motion-stress vector, which appears during the separation of variables. An independent nonlinear equation
is obtained for the impedance. The propagator for the motion vector is constructed with the aid of the impedance. The closed
analytical formulas, which do not contain any exponents with positive indices, are obtained both for the impedance and for
the motionvector propagator. The algorithm for the calculation of seismograms, free of limitations on the number and thickness
of layers, as well as on the frequency range of interest, is constructed on the basis of these formulas. The algorithm is
tested with the aid of an analytical solution. 相似文献
7.
针对利用地震道进行相对波阻抗反演中遇到的横向连续性难以保持、初始子波容错度差以及随机噪声干扰影响反演结果等问题,提出了一种基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法.该方法将地震数据剖面的Toeplitz稀疏分解问题分解为两个子反演问题,其一以Toeplitz子波矩阵元素为待反演的参数,用Fused Lasso方法求解,可保证子波具有紧支集且是光滑的;其二以稀疏反射系数矩阵元素为待反演参数,用基于回溯的快速萎缩阈值迭代算法求解,大大降低了目标函数中参数选择的难度.通过交替迭代求解上述两个子反演问题可将地震数据剖面因式分解为一个Toeplitz子波矩阵和一个稀疏反射系数矩阵;然后由反射系数矩阵递推反演可以得到高分辨率的相对波阻抗剖面;利用测井资料加入低频分量后,也可得到高分辨率的绝对波阻抗剖面.Marmousi2模型生成的合成记录算例和实际地震资料算例均表明:本文方法可以从带限地震数据中有效地反演相对波阻抗,反演结果分辨率高并且能够很好地保持地震数据的横向连续性;即使在初始估计子波存在误差和地震数据被随机噪声污染的情况下也能取得较好的效果. 相似文献
8.
Erdal afak 《Soil Dynamics and Earthquake Engineering》2006,26(1):65-70
Foundation impedance functions provide a simple means to account for soil–structure interaction (SSI) when studying seismic response of structures. Impedance functions represent the dynamic stiffness of the soil media surrounding the foundation. The fact that impedance functions are frequency dependent makes it difficult to incorporate SSI in standard time-history analysis software. This paper introduces a simple method to convert frequency-dependent impedance functions into time-domain filters. The method is based on the least-squares approximation of impedance functions by ratios of two complex polynomials. Such ratios are equivalent, in the time-domain, to discrete-time recursive filters, which are simple finite-difference equations giving the relationship between foundation forces and displacements. These filters can easily be incorporated into standard time-history analysis programs. Three examples are presented to show the applications of the method. 相似文献
9.
弹性阻抗是波阻抗的扩展,可以更好的反映储层和流体性质.常规弹性阻抗方程忽略了各向异性的影响,不能正确地描述各向异性介质中弹性阻抗的性质.本文基于Ruger给出的VTI介质反射系数近似公式,结合Connolly推导弹性阻抗的思想,推导出了一种VTI介质中的弹性阻抗方程.通过引入归一化常数,对方程进行了标准化,消除了方程量纲随角度变化而变化的不足.结合常规弹性阻抗反演的流程,得出了VTI介质中弹性阻抗反演的流程.最后,开展了从VTI介质弹性阻抗数据体中提取岩石物性参数的方法研究.数值测试表明,所讨论的参数提取方法具有一定的理论和应用价值. 相似文献
10.
网式大地电磁（Network-MT,N-MT）法采用长数公里至数十公里的电话线为电极线测量电场,很难形成两条笔直且相互垂直的电极线,因此阻抗张量的计算不如大地电磁法中直接.本文依据阻抗张量的旋转规则提出了一种计算N-MT阻抗张量的简便算法.依据该算法计算了中国东北地区5个N-MT测站的阻抗张量,获得了基于阻抗张量的视电阻率、相位曲线和最佳主轴方位角分布图像,为我国东北N-MT资料的进一步处理和解释提供了基础性数据.此外,本文对比分析了朝阳测站中6条N-MT观测电极线上阻抗张量旋转值与观测值之间的差异,重点讨论了产生这种系统性偏差的各种因素,提出其主要因素可能来自“电场等效各向异性”效应,即测站附近的地壳内部存在与观测电极线尺度相比拟的横向非均匀构造,而测站各电极线沿不同方向跨越不均匀构造,此时各电极线上的电场分量不遵循同一电场矢量的分解准则,导致地表观测三角形内阻抗分量不满足统一的阻抗张量旋转规则. 相似文献
11.
用Connolly的弹性阻抗(EI,elastic impedance)公式进行反演只能直接得到纵、横波速度和密度的信息,然后才可间接计算得到纵横波阻抗等其它的参数数据体,这样便增加了一步误差使数据的准确性降低.本文首先针对该方法的这些缺点和不足,提出了以Zoeppritz方程的Fatti近似为基础一种新的弹性阻抗公式,由该公式可得到比用Fatti近似更准确的反射系数,然后对这个公式进行了标准化以实现不同角度的弹性阻抗间量纲的统一,最后用标准化后的公式进行了反演,从反演得到的不同角度的弹性阻抗数据体中可直接提取得到纵横波阻抗数据体.应用实例表明用这种方法提取得到的纵横波阻抗更加稳定、准确,而且能很好地反映储层信息.这种新的方法是对以Connolly公式为基础的传统方法的改进. 相似文献
12.
本文对大地电磁观测阻抗实施一种数学变换--"共轭阻抗变换",发现转换后的观测阻抗和区域阻抗之间存在特定的关系,这种特定关系不受电场局部畸变的影响,而且同样不需要关于地下区域结构维性的假设,在区域结构是三维的情况下也是成立的.对转换后的观测阻抗采用已有的Swift旋转方法即可求得区域主轴方位角,然后采用最优化方法求取区域阻抗相位、振幅以及畸变因子.同时根据转换后的观测阻抗重新定义了不受电场局部畸变影响的构造维性参数.采用合成理论数据验证了新算法的正确性,和Swift、Bahr、GB、相位张量、WAL方法进行了对比分析,并将新方法应用于实测资料的解释,发展了一套消除局部畸变和进行构造维性分析的MT精细资料处理技术.进一步的工作是将新算法推广到多点多频,并采用最优化技术分解,发展一种稳定性好、多测点多频点的阻抗张量分解技术. 相似文献
13.
A direct boundary element procedure is presented to determine the impedance matrix for a three-dimensional foundation supported on an infinitely-long canyon of uniform cross-section cut in a homogeneous half-space. The uniform cross-section of the canyon permits analytical integration along the canyon axis leading to a series of two-dimensional boundary problems involving Fourier transforms of the full-space Green's functions. Solution of these two-dimensional boundary problems leads to a dynamic flexibility influence matrix which is inverted to determine the impedance matrix. The accuracy of the procedure is demonstrated by comparison with previous solutions for a surface-supported, square foundation and results obtained by a three-dimensional boundary element method (BEM) for a foundation of finite-width supported on an infinitely-long canyon. Compared with the three-dimensional BEM, the present method requires less computer storage and is more accurate and efficient. The foundation impedance matrix determined by this procedure can be incorporated in a substructure method for earthquake analysis of arch dams. 相似文献
14.
15.
The Connolly (1999) elastic impedance (EI) equation is a function of P-wave velocity, S-wave velocity, density, and incidence
angle. Conventional inversion methods based on this equation can only extract P-velocity, S-velocity, and density data directly
and the elastic impedance at different incidence angles are not at the same scale, which makes comparison difficult. We propose
a new elastic impedance equation based on the Gray et al. (1999) Zoeppritz approximation using Lamé parameters to address
the conventional inversion method’s deficiencies. This equation has been normalized to unify the elastic impedance dimensions
at different angles and used for inversion. Lamé parameters can be extracted directly from the elastic impedance data obtained
from inversion using the linear relation between Lamé parameters and elastic impedance. The application example shows that
the elastic parameters extracted using this new method are more stable and correct and can recover the reservoir information
very well. The new method is an improvement on the conventional method based on Connolly’s equation.
Wang Baoli graduated with a Bachelor’s degree in Prospecting Information and Engineering from the China University of Petroleum (East
China) in 2004 and earned her Master’s degree from the department of Geophysical Prospecting and Information Technology in
the China University of Petroleum ((East China) in 2006. She now studies for her PhD at the China University of Petroleum
(East China). Her research interest is elastic impedance inversion. 相似文献
16.
密度参数是页岩储层评价的重要敏感参数,但通常认为准确获得密度信息难度较大.本文探索提出了一种稳定的基于弹性阻抗的叠前密度反演方法,以提高优质页岩的识别精度.首先,基于不同的参数化弹性阻抗方程开展密度敏感性分析,优选了敏感性最高的Aki-Richards近似弹性阻抗方程作为反演方程.然后,在弹性阻抗方程线性化的基础上,假设反演方程系数矩阵的逆存在,将密度(取对数)表示为多个角度弹性阻抗(取对数)的加权和,通过井旁道弹性阻抗反演结果与测井数据的回归求取加权系数,避免了常规反演方法对大型系数矩阵的求逆计算,从而提高了密度反演的稳定性,同时由于利用测井资料作为匹配目标,密度反演的精度也得到了提高.模型试算和涪陵页岩气田的实际资料应用实践均表明了该方法能够有效地提高密度反演的精度与稳定性. 相似文献
17.
David L.B. Jupp 《Physics of the Earth and Planetary Interiors》1978,17(2):75-82
The successful interpretation of magnetotelluric data demands well presented, and accurate, estimates of the impedance and tipper functions.In this paper, a method based on the singular value decomposition of the local data matrix, using well known results from multivariate statistical analysis, is proposed as an alternative to current least-squares methods. The advantage obtained is that error on all data channels is recognised, so that impedance estimates are not biased by the choice of “independent noise-free” fields.An important result of the method is an extensive and structured capability for the analysis of the frequency-dependent noise which could be used to improve the resulting impedance estimates. 相似文献
18.
本文在孔隙弹性介质理论的指导下,基于入射角AVO近似方程推导了包含Russell流体项的两项AVO近似方程和相应的弹性阻抗方程,通过分析可知其精度符合反演要求.在贝叶斯理论框架下,建立了包含正则化约束的弹性阻抗反演方法,在此基础上直接提取Russell流体项.该方法可在缺少大角度叠前地震资料的情况下进行叠前直接反演得到流体因子,减少传统方法带来的累积误差.模型试算表明,该方法具有较好的准确度和稳定性.实际工区应用取得了良好的效果,表明该方法有实用性. 相似文献
19.
裂缝储层岩石物理参数的准确获得对地下裂缝预测具有重要意义,而叠前地震反演是获得裂缝岩石物理参数的有效手段.本文从裂缝岩石物理等效模型的构建出发,从测井数据上估测了裂缝岩石物理参数,通过推导含裂缝岩石物理参数的方位各向异性弹性阻抗公式,探讨了基于方位各向异性弹性阻抗的裂缝岩石物理参数地震反演方法.实际工区地震数据应用表明,基于方位各向异性弹性阻抗的裂缝岩石物理参数反演方法合理、可靠,可以降低裂缝岩石物理参数估测的不确定性,为地下裂缝预测提供有力的依据. 相似文献
20.
采用电阻率层析成像二维测量的方法,对河南省三门峡市沿河开发的建筑场地进行了勘探,获取了地下隐伏沙洞的准确位置、埋藏深度、含水情况、尺寸大小等。说明电阻率层析成像二维测量的方法在地下隐伏沙洞的探测中是行之有效的;地下隐伏沙洞在没有水体存在时,电阻率层析成像二维测量成果图中反映的一般都是高阻异常封闭圈。如果有水体存在,其异常则表现为低阻异常封闭圈。 相似文献