简介:本文提出了一种矿山微地震震源函数反演的新方法。以球坐标系波动方程为基础,对时空域中的微地震观测信号进行层析成像投影到慢度时间域,我们可以获得微地震波成像能量最大值及其对应的震源位置、发震时间和传播速度以及慢度时间信号。通过分析能量最大值处的慢度时间信号与震源函数之间的关系,推导出了计算矿山微地震震源函数的层析成像公式和利用震源函数计算微地震有效辐射能量的公式。进一步利用最小二乘法,将震源函数的振幅谱拟合成模型震源谱,确定了微地震波的零频极限值和拐角频率,最终计算出该微震事件的震源破裂半径。利用这一方法对理论模型数据和实际资料处理,结果表明:本方法具有既能求出任意一个地震事件的震源位置、发震时间和传播速度,又能同时求出震源函数及其频率特征值和震源破裂半径等震源参数的优点,因此,本方法是一种实时快速计算方法。
简介:阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期(超过100秒),逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。
简介:本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。
简介:基于双相介质理论的AVO正演技术是储层性质描述和流体预测的有效技术手段之一,但是输入参数中基质矿物模量的准确性和双相介质模型的的合理性极大地影响双相介质AVO正演效果。因此,本文采用基于流体因子的基质矿物模量反演方法,自适应反演基质矿物体积模量。引入具有岩石物理意义的多约束条件,使得流体替换技术制作的双相介质模型具有岩石物理意义。保证获得的双相介质AVO特征反映实际地层响应,真实可靠。通过不同岩性岩样的对比分析,说明反演方法的优越性和准确性。同时LH地区实际资料应用,获得孔隙度和流体饱和度等重要岩性参数变化时双相介质AVO特征,特别是不同储层孔隙度在同一入射角对应快纵波和横波反射系数幅值的大小差异和突变角差异是分辨储层孔隙度大小的依据。
简介:中国多年来在岩石圈物理学研究中已取得了一系列的重要成果,并促进了中国地球科学的发展。然而在这21世纪之际,基于国家战略需求和自主创新的方针对岩石圈物理学提出了新的挑战。本文在这样的思维前提下对岩石圈物理学的科学内涵和发展导向进行了较全面的分析、研究和思考,并明确指出当今在这一学科领域应该做些什么,核心问题是什么,又存在哪些关键性的科学问题。研究结果提出,高精度的地球物理场观测与岩石圈内壳、幔精细结构(2维和3维)的刻划;在地球内部力系作用下,深部物质和能量的交换;深部物质运移的物理-力学-化学作用过程及深层动力学响应乃是深化对地球本体的认识和揭示成山、成盆、成岩、成矿和成灾的根本机理所在。为此,在本世纪的上、中叶,在地球科学领域中地球物理学研究必为先导。本文最后对岩石圈物理研究中尚存在的一些问题和困惑进行了分析和讨论。
简介:初至拾取是起伏地表地震资料处理的关键环节,直接影响近地表建模的精度及静校正效果,但起伏地表资料极低的信噪比使传统的自动初至拾取算法几乎失效,而手动拾取要耗费巨大的人工和时间成本。本文研究了改进超虚干涉法,结合多道多域初至质量监控技术,实现起伏地表资料初至波自动拾取。改进超虚干涉法首次将近地表散射波纳入干涉法提高信噪比的范围,并通过折射波和散射波的线性组合使干涉法适用于起伏地表条件下任意复杂射线路径的初至波类型;提出反向干涉和多域干涉的概念,显著增强了所估计的虚源信号;使用波形反褶积滤波器较好地抑制了干涉导致的"假事件"的形成;采用多道多域初至质量监控技术,实现错误"假事件"初至的自动归位,提高了拾取初至的稳定性。本文研究的初至波自动拾取理论与技术具有突出的鲁棒性和稳定批量处理大量实际三维地震数据的能力,在中国西部山区三维地震勘探数据处理的应用效果显著,质量优于某常用商业软件。
简介:在水库现场试验了RISK2型探地雷达探测水库冰厚度的能力,试验时所用天线频率为600MHz;同步钻孔测量雷达探测处的冰厚度;以及在一个点上取样测试分析冰晶体、冰内气泡和冰密度。试验时冰面积雪厚度0.03-0.05m,冰层上部有0.24m粒状冰,其下均为柱状冰;冰内气泡含量呈表层高底层低分布;冰密度随气泡含量变化;冰厚度在平面内不均一。通过探测厚度和实测厚度的对比分析以及气泡含量对介电系数影响的理论分析,建立了积雪、粒状冰和柱状冰三层介质模型,获取雷达波在冰内的理论传递时间。结果发现:能够利用等效介电常数或等效传播速度评价雷达波传递时间,结冰期冰层1/3深度处的对应介电常数或传递速度可以作为等效值;另外因冰内大气泡造成的理论传递时间大于雷达探测时间,其差值随理论传递时间或冰厚的增加呈非线性增加。
简介:ThevariationsintheEarth’srotationareimportanttospacedynamictheoryandnaturaldisastersbecausetheyaffectthelength-of-day(LOD)andconsequentlyhumanlives.Weusemaximumentropymethod(MEM),Lombmethod(LOMB),andphasedispersionminimization(PDM)todeterminethenaturalperiodsofequallyspacedLODtimeseries.WetransformtheobservationalmonthlyLODtimeseries(LODM)tounequallysampledseries(LODMD)byremovingeveryfourth,third,andhalfofthetotalsamples.WealsoapplysplineinterpolationtoLODMDtoyieldequallyspacedtimeseries(LODMDN).Theresultssuggestthatregardlessofthetimeseries,theMEMfrequencyis0.1660month^-1and0.0840month^-1,whereasLOMBandPDMyield0.166month^-1and0.083month^-1,respectively.Furthermore,missingdatathatarelessthanhalfofthetotaldataorsplineinterpolationdonotaffecttheanalysis.Fortheamplitude,neithermissingdatanorsplineinterpolationaffecttheanalysis.
简介:探索地震前兆是一件具有挑战性的研究。本文采用高通滤波的方法,处理分析了来自世界不同位置4个超导重力台站5次大地震前的秒级固体潮重力数据,发现在这些地震发生之前数天内,能够接收到稳定的重力高频扰动信号,其中2008年5月12日的汶川地震前出现一个主频为0.1-0.3Hz的重力扰动,而其他4次地震前出现两个主频为0.12-0.17Hz和0.06-0.085Hz的重力扰动。另一个发现是当地震发生在大陆与海洋板块断裂带时这种重力扰动信号具有远场效应,即重力扰动异常出现在距离震中较远的超导重力仪上,而较近的台站异常较小或没有异常。以上研究结果表明,这种重力扰动信号与地震发生具有较密切的相关性,可作为一种有潜力的临震预报信号;重力扰动信号的远场效应也许能够揭示某些地震发生前地球内部板块之间的相互作用机制。但由于地球上重力固体潮台站分布的非常不均匀,上述结果还需要今后研究更多的样本数据加以证实。
简介:由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。
简介:多分量地震数据中低频缺失是弹性波全波形反演中的一大难题,低频的缺失导致全波形反演无法有效恢复介质的长波长成分进而使反演陷入局部极值。为此,本文提出了反演介质纵横波速度长波长分量的弹性波包络反演方法。该方法利用包络算子具有的解调多分量数据中隐含的低频信息的能力,构造多分量地震数据的包络目标函数进行反演,用以恢复地下介质纵横波速度的长波长成分。一系列数值试验表明,即使在多分量地震数据中缺失低频信息、并且初始模型缺少先验信息的情况下,这种弹性波包络反演方法能够有效降低波形反演的非线性,可以为后续的常规弹性波全波形反演或者深度偏移提供足够精确的初始模型,且该方法对横波速度长波长分量的重建尤为有效。Marmousi-2模型的高精度纵横波速度的反演结果表明,利用该方法反演的纵横波速度作为常规弹性波全波形反演的初始模型,可以显著提高反演结果的精度。此外,本文对弹性波包络反演方法的适用性也进行了初步的研究与讨论。
简介:渗透率评价是储层评价中的一个复杂问题,传统测井方法难以对储层渗透率参数进行直接和有效的评估。储层渗透率对钻井过程中的泥浆滤液侵入有着较大的影响,因此本文提出一种利用泥浆侵入效应来评价储层渗透率的数学模型和方法。首先构造含泥饼的泥浆侵入数值模型,然后以达西渗流理论为基础导出泥浆侵人深度和储层渗透率的近似数学关系,以此可以利用侵入深度测量值来估算储层渗透率。对孔隙度、渗透率和含水饱和度三个主要储层参数的敏感性分析,发现该方法适用于受到泥浆侵入影响的低孔隙度、低渗透率的油层或油水同层。采用现场测井和取芯数据进行数值模拟计算,结果表明估算出来的渗透率曲线和预设的渗透率曲线吻合较好,证明该方法具有一定的可行性。
简介:由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。
简介:地震纹理属性与地震相和储层特征密切相关,广泛地应用在地震资料的解释中。传统的地震纹理属性基于叠后数据提取,受叠加作用的影响,易造成地层特征信息的损失,不利于复杂储层的描述。本文提出叠前纹理属性,其不仅可以精细地刻画不同反射点波形的横向连续性,也能体现AVO、各向异性和介质的均质性。基于叠前纹理属性丰富的地层特征表达能力,结合SOM聚类算法,形成了利用叠前数据进行地震相分析的方法。该方法应用于中国某工区宽方位地震资料,通过对比证实了叠前纹理属性描述地层横向变化的优越性,并能揭示各向异性特征及非均质性特征,基于叠前纹理的分类结果能有效区分不同的地震反射模式,为地震相分析提供了可靠的依据。
简介:为了有效处理人工源的影响,本文开发了带源的CSAMT二维正反演算法,可用于全区(近区、过渡区和远区)资料的反演。引入正则化因子完成磁法二维反演,并且将模型参数调整为磁化率的对数,保证反演过程中磁化率始终为正值。本文基于交叉梯度原理,将CSAMT和磁法进行联合反演,通过搜索交叉梯度项权重的方法,避免了不同异常源引起的两种异常相互干扰的问题。理论模型算例表明基于交叉梯度的联合反演方法优于单独反演。本文开发的带源CSAMT二维正反演算法,有效处理了人工源的影响,保证了最终联合反演算法的可靠性。