共查询到20条相似文献，搜索用时 187 毫秒
1.
2.
饱和储层砂岩中流体粘度引起的超声波速度频散实验研究(英文) 总被引：8，自引：3，他引：5
在实验室对5块储层砂岩进行了模拟地层压力条件下的超声波速度测试。砂岩样品采自WXS凹陷的W地层,覆盖了从低到高的孔隙度和渗透率范围。实验选用了卤水和4种不同密度油作为孔隙流体,结合温度变化,实现了对流体粘度引致的速度频散研究。对实验结果的分析表明:(1)对于高孔隙度和渗透率的样品,无论是哪种流体饱和,观察到的超声波速度测试值和零频率Gassmann预测值的差异较小(约2-3%),基本上可以用Biot模型解释;对于中等孔隙度和渗透率的样品,低粘度流体(<约3mP?S)的频散效应也可以用Biot模型得到合理解释;(2)对于低、中孔隙度和渗透率样品,当流体粘度增加时,喷射流机制起主导作用,导致严重的速度频散(可达8%)。对储层砂岩的微裂隙纵横比进行了估计并用于喷射流特征频率的计算,当高于该特征频率时,Gassmann理论的假设条件受到破坏,实验室测得的高频速度不能直接用于地震低频条件下的W地层砂岩的Gassmann流体替换研究。 相似文献
3.
阵列声波测井波形中包含丰富的地层纵、横波速度信息.如何从这些波形中提取出这些速度信息,得到比较准确的地层纵、横波时差是资料应用人员非常关心的问题.目前,最常用的纵、横时差处理方法是STC———慢度—时间相关法.该方法只有在波形的速度频散不严重时比较有效.而井内声波传播的基本特征是具有频散特征.利用Ma-trix pencil方法可以求出阵列声波测井波形的二维谱.借助于波形的二维谱,本文给出一种求地层的纵、横波时差的方法.本方法充分利用井内声波传播的规律,有效地减小了频散对纵、横波时差计算的影响,得到的纵、横波时差值精度比较高. 相似文献
4.
实验表明,在应变小于10-6范围内,砂岩对地震波的吸收主要由孔隙流体的局部运动引起,而且Q值随频率的变化出现谐振现象。据此,作者认为,地壳上部地震波能量的耗损可以表示为以下两种成分的线性组合:1.由滑动摩擦、热弹性驰豫、位错运动等引起的能量消散,它们主要与岩石的固相成分有关,可近似地用恒Q模型描述;2.由孔隙流体运动引起的能量消散,可以用谐振公式描述。根据这一认识,利用描述因果关系的Kramers-Krnig关系式可推导出表示地壳上部岩石粘弹性（复弹性模量,相速度频散和衰减函数）的公式,它们综合地描述了由各种机制引起的波的频散和吸收,并在谐振Q值等于参考常数Q值时退化为目前常用的Futterman模型。作为这种谐振Q模型的应用,介绍了它用于Q值测量结果外推和频散一吸收研究以及粘弹性介质中反射地震道合成的结果。 相似文献
5.
有限差分是最常用的地震波方程数值模拟方法,但时间和空间离散会产生数值频散.正演是逆时偏移和全波形反演的基本单元,成像和反演的精度很大程度依赖于所采用的数值模拟算法.本文研究了有限差分的时间和空间频散特性及其对逆时偏移和全波形反演的影响.通过时间有限差分+伪谱法、时间频散校正+空间有限差分、时间频散校正+伪谱法获取时间频散、空间频散和无频散数据;发展了抗时间频散、抗空间频散、抗时间+空间频散的逆时偏移和全波形反演方法;采用理论模型和实际资料对提出的方法进行了测试.数值结果表明:逆时偏移同时受时间和空间频散影响,时间频散导致同相轴不聚焦、成像位置偏离,空间频散会产生高频噪声和虚假反射界面;全波形反演在低频大尺度反演中几乎不受时间和空间频散影响,高频精细反演中时间频散引起波形相移、降低反演精度,空间频散增加多解性、导致反演不收敛;抗频散方法可以有效缓解时间和空间频散影响,获得高质量的偏移剖面和反演结果. 相似文献
6.
非弹性的表面波大尺度波形反演的理论方法以如下三方面提出的。（１）计算本征值和本征函数可获得第一阶的频散参数，并把体波的复数形式和波能的吸收引入非弹性中，（２）根据弗拉尔和斯奈德的理论方法，计算了非弹性表面波的激发和合成地震图，球阶非弹性用频散的计算按Ｂｉｓｗａｓ－Ｋｎｏｐｐｏｆ公式被变换成平滑的弹性形式，（３）使用自伴微工子导出关于模型参数的波形Ｆｒｅｃｈｅｔ导数，然后应用非线性反演拟合波形。 相似文献
7.
《应用地球物理》2016,(2)
二氧化碳注入地层驱油过程中会引起地震波速度的频散现象,这是地震波传播过程中相速度随频率发生变化所造成的。所以综合考虑二氧化碳气驱情况下的速度频散效应,提取频散因子,进而动态监测气驱前后的油藏变化,具有方法的可行性。本文从Robinson频散褶积模型出发,首先推导了速度频散因子的定量表征式,构建了含有速度频散因子的反演方程,给出了四个处理步骤的实现过程,其中包括广义S变换谱分解、井旁道子波提取、谱均衡处理和阻尼最小二乘频散因子提取。根据粘弹性介质理论,建立了基于Q值变化的衰减地层模型,对止演后的偏移剖面进行了时频谱分析,在此基础上进行了频散因子提取与对比,验证了方法的有效性。最后,将本文方法应用到胜利油田高89井区注气前后两期资料的目标处理与驱油监测中。从提取的两期地震频散因子剖面和切片来看,方法能较好地指示二氧化碳驱前后的差异性,研究成果与井场实际考察结果非常吻合。 相似文献
8.
9.
由于介观尺度的孔隙流体流动,弹性波传播过孔隙岩层时在地震频段表现出较强的频散和衰减。Johnson理论给出了在任意孔隙形状的条件下,部分气水饱和孔隙介质的理论相速度和品质因子的解析解。本文在Johnson模型的基础上,通过对Q值曲线的低频和高频近似,推导了Q值曲线的近似公式,以及基于孔隙介质基本地球物理参数和孔隙斑块几何形态参数T和比表面积S/V的最大衰减Qmin近似公式。通过与理论值的对比,对Qmin近似公式存在的线性误差进行改正,进一步提高了精度。复杂的斑块形态对最大衰减Qmin和过渡频率ftr的都产生一定影响,且对ftr影响更大。因为数值模拟直接求解介观尺度的Biot孔隙介质方程需要极大的计算量,我们使用Zener模型建立了等效粘弹模型,有效地模拟了地震频带内的衰减和频散现象。 相似文献
10.
本文以Futterman提出的地震波振幅衰减和相速度频散表达式为基础,从井震匹配的角度出发,推导了不进行反Q滤波、仅反Q滤波振幅补偿、仅反Q滤波相位补偿,以及反Q滤波相位与振幅同时补偿四种情况下,地震波速度频散与相速度及地震记录振幅谱间的表达式,并从理论上说明了反Q滤波相位补偿的必要性.通过相同观测系统,炸药震源和可控震源分别激发采集的零偏移距VSP资料实例,验证了本文所推导的速度频散表达式的合理性;通过井震标定实例,进一步说明了反Q滤波相位补偿,可有效消除地震子波速度频散,提升地面地震资料与零偏移距VSP走廊叠加剖面的匹配度,最终提高地震资料成果的可信度. 相似文献
11.
用能量法测定柱状样品的衰减Q值 总被引：1，自引：0，他引：1
对常用的测试衰减Ｑ值的频谱振幅比法和波形反常进行了分析,认为它们对样品几何开头和大小的严格要求是大多数实验室待测样品所不能满足的。因而,提出用波形的能量来计算衰减Ｑ值的方法,阐述了其原理、计算公式和注意事项,讨论了ＱＰ和ＱＥ之间的关系,并对多种样品进行了实验室测量,将本文方法测得的Ｑ值与其它方法测试值进行了比较。认为本方法的Ｑ值介于ＱＰ和ＱＥ之间。 相似文献
12.
Two methods for measuring the attenuation Q value,frequency-amplitude ratio method andwaveform inversion method,were analyzed.The size and shape of most rock samples in thelaboratory can hardly fulfill the rigorous requirements of the above methods.A new testmethod of using the wave energy to calculate the Q value was proposed.The principle,calculation formulae and some points for attention were expounded;the relationship betweenQ_P and Q_E was discussed;and many samples were measured by this method.The result wascompared with those by other methods,and it was considered that the Q value by this methodlies between Q_P and Q_E· 相似文献
13.
频率域正演是频率域波形反演的基础,有效快速的正演差分格式可以保证反演结果的精度和效率.本文以用较小的系数矩阵带宽来高效地压制频域正演频散为目标,综合利用加权平均算子、平均加速度项和优化系数三种方法,提出了优化15点差分格式;并且采用压缩存储方式来存放大型系数矩阵,极大地缩小了内存使用量;进而结合最佳匹配层(PML)边界条件,明显地压制了边界反射;最后,通过与前人方法的对比验证,证实了本方法可以在不明显增加计算量的情况下,较好地压制频散. 相似文献
14.
多震源编码技术可以提高全波形反演的计算效率,但同时会引入串扰噪声使反演结果质量降低. 全变分约束可以有效地压制层内噪声并突出模型界面,其与多震源技术的结合,能在大大提高弹性波全波形反演效率的同时提高反演质量. 本文提出了一种高效的动态多震源全波形反演策略,可以在离散串扰噪声的同时保证照明的均匀性. 根据残留串扰噪声的分布特征,构建了与之匹配的基于各向异性全变分约束的弹性波全波形反演方法. 为了减少周期跳跃效应,将基于稀疏约束的低频重构算法应用于弹性波地震记录,给出了利用快速梯度投影算法求解各向异性全变分约束的全波形反演流程. 模型数据测试结果表明本文的方法不仅能有效地抑制多震源方法导致的串扰噪声,同时也能降低观测数据中的噪声对反演结果的影响. 相似文献
15.
为了提高表层速度反演精度,本文提出了一种新的波形反演方法.该方法只利用初至波波形信息以减少波形反演对初始模型的依赖性,降低反演多解性与稳定性.由于只利用初至波波形信息,所以该方法利用高斯束计算格林函数和正演波场,以减少正演计算量.为了避免庞大核函数的存储,该方法基于Born波路径,利用矩阵分解算法实现方向与步长的累加计算.将此基于Born波路径的初至波波形反演方法应用于理论模型实验,并与声波方程全波形反演和初至波射线走时层析方法相对比,发现该方法的反演效果略低于全波形反演方法,但明显优于传统初至波射线走时层析方法,而计算效率却与射线走时层析相当.同时,相对于全波形反演,本文方法对初始模型的依赖性也有所降低. 相似文献
16.
WTW(Wave equation traveltime+Waveform inversion)反演是基于波动方程的走时反演（WT反演）和波形反演的联合反演方法.WT反演利用波动方程计算走时和走时关于速度的导数，和传统以射线为基础的走时反演相比，具有不必射线追踪、不必拾取初至、不必高频假设以及初始模型和实际模型差别较大时也能较好收敛等优点，但WT反演与波形反演相比其结果分辨率低.与之互补的是，波形反演的反演结果分辨率高，但是当所给初始模型和实际模型相差太大时，波形反演迭代算法容易陷入局部极小点.可见结合两种方法的WTW反演是一种比较好的联合反演方法.常规WTW迭代算法是首先以WT反演为主反演得到地质模型的整体特征，然后再以波形反演为主反演模型细节，该算法耗时和占用计算机存储空间接近WT反演或波形反演的两倍.为了节省运算耗时和计算机存储空间，往往采取首先单独利用WT反演然后再单独利用波形反演的算法.这样做的缺点是不能紧密结合两种反演方法，使得它们的优缺点在每一次迭代中无法得到互补，从而影响了最终的反演结果.针对以上事实，本文提出一种新的方法实现WTW，使得WTW运算速度和存储空间在任何情况下等同于WT反演或波形反演.模型计算表明新的算法具有更好的收敛性. 相似文献
17.
本文针对层状介质并结合梯度法波形反演，提出逐层波形反演的方法. 首先给出介质扰动响应的概念，并在此基础上分析了梯度法波形反演方法. 波形反演实质上是将实测地震记录和预测地震记录的波形残差信息转化为实际地质模型和预测地质模型的模型残差信息. 波形反演的优点是利用大量振幅相位信息得到高分辨率的反演结果， 其缺点是运行耗时大；当初始模型和实际模型相差较大时，迭代算法容易陷入局部极小点，这是因为目标函数和初始模型同实际模型间的差异是非线性的关系. 逐层波形反演方法是使自上而下每一层的目标函数最小，这样总的目标函数也是最小的. 利用二分法速度扫描确定每一层速度不仅提高了运算速度也避免了迭代算法陷入局部极小点的问题. 结合介质扰动响应和目标函数值变化可以更为准确迅速地确定每一层速度和该层界面位置. 相似文献
18.
19.
Feasibility of waveform inversion of Rayleigh waves for shallow shear-wave velocity using a genetic algorithm 总被引：1，自引：0，他引：1
Chong Zeng Jianghai Xia Richard D. MillerGeorgios P. Tsoflias 《Journal of Applied Geophysics》2011,75(4):648-655
Conventional surface wave inversion for shallow shear (S)-wave velocity relies on the generation of dispersion curves of Rayleigh waves. This constrains the method to only laterally homogeneous (or very smooth laterally heterogeneous) earth models. Waveform inversion directly fits waveforms on seismograms, hence, does not have such a limitation. Waveforms of Rayleigh waves are highly related to S-wave velocities. By inverting the waveforms of Rayleigh waves on a near-surface seismogram, shallow S-wave velocities can be estimated for earth models with strong lateral heterogeneity. We employ genetic algorithm (GA) to perform waveform inversion of Rayleigh waves for S-wave velocities. The forward problem is solved by finite-difference modeling in the time domain. The model space is updated by generating offspring models using GA. Final solutions can be found through an iterative waveform-fitting scheme. Inversions based on synthetic records show that the S-wave velocities can be recovered successfully with errors no more than 10% for several typical near-surface earth models. For layered earth models, the proposed method can generate one-dimensional S-wave velocity profiles without the knowledge of initial models. For earth models containing lateral heterogeneity in which case conventional dispersion-curve-based inversion methods are challenging, it is feasible to produce high-resolution S-wave velocity sections by GA waveform inversion with appropriate priori information. The synthetic tests indicate that the GA waveform inversion of Rayleigh waves has the great potential for shallow S-wave velocity imaging with the existence of strong lateral heterogeneity. 相似文献
20.
波动方程深度偏移是解决复杂地质体成像的关键技术,基于波动方程的速度建模为其提供更为精确的速度模型.频率域波形反演是目前研究最为广泛的波动方程速度建模方法之一,它推动了波形反演在勘探尺度下的应用.本文通过对频率域波形反演的实现,分析对比了其有效执行过程中与频率相关的影响因素.介绍了时间域的多尺度反演方法在频率域的一种实现方式,对比分析了输入数据的频点带宽和应用的子波频带范围不同时对反演结果的影响.本文通过设计的山地地质模型对频率域波形反演进行了测试和对比,得到的结论为频率域波形反演的有效计算提供了依据和参考. 相似文献