首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 171 毫秒
1.
利用走时反演二维多孔介质渗透率   总被引:2,自引:1,他引:2       下载免费PDF全文
综合利用直接方法和优化方法研究二维多孔介质渗透率反演问题. 问题的提法与前人不同,是由区域内各点流体渗透的走时来反演渗透率. 该反问题的求解可分成两步进行,归结为两个相应的子问题进行研究. 首先由流体走时反演二维多孔介质速度场的数值算法,然后由二维多孔介质速度场反演二维多孔介质渗透率的数值算法,最后给出数值例子. 数值结果表明了所用数值方法的正确性和有效性.  相似文献   

2.
传统的地震逆散射广义Radon变换(GRT)保幅反演方法是建立在散射场一阶Born近似(单散射)的基础上,仅仅适用于弱扰动介质模型.本文从散射场积分方程出发,通过研究二次散射的特征,讨论和验证了基于局部二阶Born近似的GRT非线性保幅反演方法,将传统GRT线性保幅反演算子的适用范围扩展至非均匀强扰动介质.数值测试结果表明:在散射场近似模拟方面,二阶Born近似比一阶Born近似更为准确,二次散射效应主要集中在主散射点周围的局部区域内,超过这一范围,二次散射强度趋于稳定;在保幅反演方面,本文基于局部二阶Born近似的GRT非线性反演算法,明显优于传统的GRT线性反演算法,可以准确重构强扰动介质模型,而计算效率与线性反演方法相当.  相似文献   

3.
逆散射反演是由已知散射波场反演地下介质中散射体各参数了分布的方法,本文从这一思想出发,针对二维大地电磁资料,在BORN近似的前提下,提出了一种新的反地下电导率参数的方法,经模型试算与实际资料的反演验证正确,这一方法特别适于地下无强电性分界面的地区。  相似文献   

4.
用地质雷达数据资料反演二维地下介质的方法   总被引:7,自引:3,他引:7       下载免费PDF全文
从二维麦克斯韦方程组出发推导出反演介电常数和电导率等二维介质物性参数的反演公式.反演的步骤是: 建立初始猜测模型,利用电磁波时间域有限差分法模拟正演数据,用正演数据与观测数据之间的数据残差建立目标函数,通过引入一个由麦克斯韦方程计算的伴随场,将目标函数对介质参数的导数表示成显式形式,应用最优化理论得出对初始猜测模型的修改,用共轭梯度法迭代,最终得到反演结果.用合成数据反演具有粗糙地表的非导电介质的介电常数,用实验数据同时反演介电常数和电导率,并比较了麦克斯韦方程反演结果与声波方程反演结果、波动方程偏移剖面的差异.  相似文献   

5.
变背景速度2.5维反散射技术在地震勘探中的应用问题   总被引:2,自引:0,他引:2  
在三维点源、二维介质构造定义下的Born反散射模型,采用WKBJ近似,建立由地震反射数据求取速度摄动量或反射系数(界面)的反演公式,反演效果与深度偏移类似。本文目的是将所反演的反射系数(界面),用作最优拟合参数反演的反射层位约束信息,也可用作深度偏移处理。  相似文献   

6.
二维介质柱的电磁成像研究   总被引:2,自引:2,他引:0       下载免费PDF全文
从体等效原理出发得到用介质柱介电参数表示的体等效电流,并得到一个积分方程组求解散射问题.利用变分原理,对此方程组求变分,并化为矩阵方程组,消去含介质柱内总场变化的项,得到散射场变化与介质柱介电参数变化之间关系的矩阵方程.即反演方程.用算例检验此法的反演精度、收敛性能、多目标反演能力、抗噪声性能及利用不完备测量散射信息反演目标的能力.结果表明,此法是一种实用性很强的二维介质柱电磁成像方法.  相似文献   

7.
声波方程逆散射反演的近似方法   总被引:7,自引:0,他引:7  
我们在文献[1]里研究了介质参考波速沿某个方向线性变化时的三维声散射理论,导出了声波方程逆散射反演问题解的解析表达式.考虑到应用时的实际条件,本文根据上述反演方法导出2.5维模型的声波方程逆散射反演的波速扰动计算公式,给出该方法在“高频”近似条件下的波速扰动反演计算公式,从而使我们提出的“参考波速线性变化时的声波方程逆散射反演”理论更接近实际应用条件.本文给出的这些反演公式仍然具有原方法的优点,即不但可以使Born近似的假定在大多数情况下能得以满足,而且可以利用快速Fourier变换来快速实现介质波速扰动的反演成象.  相似文献   

8.
声波方程逆散射反演的近似方法   总被引:2,自引:0,他引:2       下载免费PDF全文
我们在文献[1]里研究了介质参考波速沿某个方向线性变化时的三维声散射理论,导出了声波方程逆散射反演问题解的解析表达式.考虑到应用时的实际条件,本文根据上述反演方法导出2.5维模型的声波方程逆散射反演的波速扰动计算公式,给出该方法在“高频”近似条件下的波速扰动反演计算公式,从而使我们提出的“参考波速线性变化时的声波方程逆散射反演”理论更接近实际应用条件.本文给出的这些反演公式仍然具有原方法的优点,即不但可以使Born近似的假定在大多数情况下能得以满足,而且可以利用快速Fourier变换来快速实现介质波速扰动的反演成象.  相似文献   

9.
叠前共炮点道集的奇性反演研究   总被引:1,自引:0,他引:1       下载免费PDF全文
讨论了地震波场的Fourier积分算子奇性反演理论与方法,得出常数背景假设条件下二维共炮点道集奇性反演的解析表达式.构造了两个理论地震模型,分别进行二维有限差分正演与奇性反演的数值计算将二维Fourier积分算子奇性反演方法用于野外的二维地震勘探叠前共炮点记录.在勘探范围内,当介质波速变化相对背景波速符合小扰动假设时,理论模型与实际地震数据的计算结果表明这种奇性反演方法的有效性和实用性  相似文献   

10.
广义散射层析成像反演   总被引:1,自引:1,他引:0       下载免费PDF全文
本文详细地给出了基于非均匀介质的体散射广义散射层析成像反演的基本理论.广义散射层析成像反演可以描述为波场的反向传播和对成像场进行局部波数域滤波的过程.在数值算例中,利用背景速度沿深度方向均匀变化的v(z)介质中的简单的方块作为速度异常体的模型,通过对该模型产生的低频的Born数据和声波的正演数据的测试,在对采集系统进行有限频率带宽和空间孔径的校正来进行局部成像矩阵谱的恢复中,可以看出模型中各点的谱在恢复后的质量无论从覆盖的面积范围还是幅值的均一性上都有着明显的提高;在对速度模型的重建中,广义散射层析成像反演能够很好地恢复速度模型的低频分量,即便是方块速度异常体相对于背景速度的平均速度扰动是23%也能很好地重建模型中的速度,且对于不同的背景速度模型基本上都能很好地恢复Marmousi速度模型的低频分量.所以该方法将基于Born模型的层析成像反演适应范围进行了一定程度的扩展.  相似文献   

11.
计算台站接收函数的最大熵谱反褶积方法   总被引:32,自引:1,他引:32       下载免费PDF全文
提出一种在时间域采用最大熵谱反褶积提取台站接收函数的方法,以最大熵作为自相关函数和互相关函数的递推准则,利用Toeplitz方程及Levinson递推算法,得到预测误差滤波系数的递推公式,从而计算台站接收函数.外推运算过程中,反射系数总是小于1,保证了最大熵谱反褶积的稳定性.时窗外数据熵极大提高了接收函数的分辨率.合成地震图与实测地震图的检验表明, 最大熵谱反褶积是一种在时间域测定台站接收函数的有效方法.   相似文献   

12.
依据瞬变电磁扩散场与虚拟波动场之间的数学表达式,通过波场反变换得到虚拟波场,进而借用地震勘探上成熟的成像方法定位地质界面的位置以及形态.本文首先将时间域波动方程转换为频率域波动方程,将总波动场分为背景场和散射场、总波速分为背景速度和扰动速度,然后运用Green定理得到了散射场关于地质界面速度扰动量的表达式,并使用Bor...  相似文献   

13.
一种新的实时电磁逆散射方法   总被引:2,自引:0,他引:2       下载免费PDF全文
为解决介质圆柱体逆散射问题,提出一种新的在线逆散射方法,通过支持向量机将逆散射问题转化成一个回归估计问题. 该方法可应用于各种逆散射方面, 尤其是目标的几何与电磁参数重构和埋地目标探测. 文中首次将支持向量机方法应用到该领域,设置多个散射场的观测点,通过提取散射场的不同信息作为样本信息训练支持向量机, 建立了介质圆柱体的逆散射模型, 利用该模型重构了介质圆柱体的电磁参数,同时探测了埋地位置. 数值结果显示了该方法的有效性和准确性,为目标的实时逆散射研究提供了一种有效方法.  相似文献   

14.
15.
三维散射算子及其在逆散射中的应用   总被引:1,自引:3,他引:1  
以弹性波为重点,评述了三锥散射算子关系及其在逆散射中的应用。为便于理解,在介绍过程中,还仔细描述一维球对称情况下散射矩阵和逆散射方程推导。本文表明,无论在一维球对称介质情形,还是三维情形(包括声波、电磁波和弹性波),其因果场、反因果场和散射数据之间都存在一个关系,即散射算子关系。如果利用因果场与反因果场之间的时间反转关系,可以得到因果场与散射数据之间的关系或方程。再引进三种场的分解方法,可以得到三种关于因果场的方程,这就是逆散射方程。这些方程分别提出了新的问题。  相似文献   

16.
逆散射级数法预测层间多次波的算法改进及其策略   总被引:12,自引:5,他引:7       下载免费PDF全文
波动方程预测层间多次波有两类算法.一是由Berkhout和Verschuur(1998)提出的基于CFP(共聚焦点)延拓的算法,另一是由Weglein(1997)等人提出的基于逆散射级数法的算法.ISS(逆散射级数)算法具有不依赖速度模型的优点,但是其计算成本很高.本文改进了Weglein提出的1-D预测公式,提高了计算效率.从理论方面,改进算法的计算速度提高了大约12倍.实际结果表明其计算速度更快,提高近80倍左右,主要原因在于后者的空间复杂度小.此外,由于时间-空间域比伪深度-波数域的层间多次波预测有更多的优点,本文推导了1.5-D时间-空间域ISS层间多次波预测算法,该算法有计算速度快、预测噪音小、适应观测系统能力强和不依赖地震子波等优点.在实际应用中,计算成本高是ISS预测层间多次波算法的主要障碍.因此,本文还讨论了预测层间多次波的策略.结合本文提出的改进算法和应用策略,能够在实际地震资料处理中较好地应用ISS预测层间多次波.  相似文献   

17.
Evolutionary computation based on the idea of biologic evolution is one type of global optimization algorithm that uses self-adaptation, self-organization and random searching to solve optimization problems. The evolutionary-simplex algorithm is introduced in this paper. It contains floating encoding which combines the evolutionary computation and the simplex algorithm to ovcrcomc the problems encountered in the genetic algorithm and evolutionary strategy methods.Numerical cxpcrimcnts arc performed using seven typical functions to verify the algorithm. An inverse analysis method to identify structural physical parameters based on incomplete dynamic responses obtained from the analysis in the time domain is prcscntcd by using the evolutionary-simplex algorithm. The modal evolutionary-simplex algorithm converted from the time domain to the modal domain is proposed to improve the inverse efficiency. Numerical calculations for a 50-DOF system show that whcn compared with other methods, the evolutionary-simplex algorithm offers advantages of high precision,cfficient searching ability, strong ability to resist noise, independence of initial value, and good adaptation to incomplete information conditions.  相似文献   

18.
A new method for the 2D inversion of induced polarization (IP) data in the time domain has been developed. The entire IP transients were observed and inverted into 2D Cole-Cole earth models, including resistivity, chargeability, relaxation time and the frequency constant. Firstly, a modified 1D time-domain electromagnetic algorithm was used to calculate the response of a layered polarizable ground. The transient signals were then inverted using the Marquardt method to derive the Cole-Cole parameters of each layer. However, model calculations showed that the EM effects could be neglected for the time range (>1 ms) and for the transmitter–receiver distances (<50 m) used in this study. Therefore, the induction effects were not considered for the solution of the 2D inverse problem and a DC solution was applied. An approximative forward algorithm was introduced in order to calculate the IP transients directly in the time domain and in order to speed up the inverse procedure. The approximation is highly accurate, and this is demonstrated by comparing the approximations with their exact solutions up to 3D. The inverse algorithm presented consists of two steps. The transient voltages of an array data set were inverted separately into a two-dimensional resistivity model for each time channel. The time-dependent resistivity of each cell was then interpreted as the response of a homogeneous half-space. In the 2D inversion algorithm, a 3D DC algorithm was used as a forward operator. The method only requires a standard 2D DC inversion and a homogenous half-space Cole-Cole inversion. The developed algorithm has been successfully applied to synthetic data sets and to a field data set obtained from a waste site situated close to Düren in Germany.  相似文献   

19.
用稳定高效的反Q滤波技术提高地震资料分辨率   总被引:5,自引:2,他引:3       下载免费PDF全文
地震波在地下传播时受到衰减影响,衰减会导致地震波场高频能量的损失和相位畸变.反Q滤波可补偿大地衰减效应.已有的反Q滤波方法存在下列不足:频率域的算法由于算子长度较长,所以计算效率较低;时间域的算法,或者对地震记录上到时较晚的同相轴进行了过度的补偿,或者为防止过度补偿后来的振幅而在最大增益处进行限制,导致振幅的多解性,而且还会影响滤波器的相位效应.本文给出一种通过直接求解时间域的Q模型方程来进行反Q滤波的算法.由于采用带状矩阵解算器,所以具有较高的计算效率,理论数据和实际地震资料的试算结果证明,本方法对地震波的吸收衰减进行了出色的补偿,提高了地震资料的分辨率.  相似文献   

20.
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号