简介：High-qualityseismicgeometryisthekeytoobtainhigh-qualityseismicdata,andcanaffecttheaccuracyofdataprocessingandimaging.Basedontheanalysisoftherelationshipbetweenthequalityofthegeometryandthefouracquisitionparameters(thenumberoftraces,shotlinespacing,andthespaceandnumberofreceiverlines),aqualityevaluationmethodofthegeometrybasedoncomprehensivequalityfactor(CQF)isproposed,andtherelationshipbetweenthegeometryqualityandthefourparametersisgiven.WeusefielddatacollectedinanoilfieldinWesternChinawithcomplexgeology:Firstweuseawideazimuthgeometry.Then,wecalculatetherelationshipcurvebetweengeometryanddataqualitybyvaryingeachparameterwhilekeepingtherestfixed.andtheanalysisresultsaregivenbyusingtheCQFevaluationmethod.Theresultsshowthattheshot-linespacinghasthegreatesteffectonthequalityofthegeometry,andtheincreaseofthereceiverlinespacingcanappropriatelyimprovethequalityofthegeometry,andtheincreaseofthenumberofreceivingtracescanimprovethegeometryquality.Thedifferentacquisitionparametershavedifferenteffectsontheimagingqualityofshallowanddeepevents.Themodelforwardandprestackdepthmigrationareusedtogenerateprestackdepthmigrationprofileswithdifferentacquisitionparameters.Theimagingresultsareconsistentwiththeabovecalculatedresults.Accordingtothedepthofthetargetlayer,thequalityfactorevaluationmethodisappliedtoguidethedesignofthegeometryandoptimizetheacquisitionparameterstoimprovetheimagingaccuracyofseismicdata.

简介：储层物性参数的改变伴随着弹性参数某种程度的改变，而这种改变并非简单的线性关系。加上观测信息的缺失、重叠、噪声破坏及模型理想化等原因，致使这种改变关系的获取具有较大的不确定性。针对传统统计岩石物理物性参数反演方法的不足，以利用弹性参数反演储层物性参数为目的，依据贝叶斯反演框架，我们建立了新的储层物性参数目标反演函数。首先，采用兼具确定性与随机性特点的统计岩石物理模型，考虑到不同弹性参数间的精度存在差异，引入权重系数，建立起储层物性参数与弹性参数间的加权统计关系。其次，基于这种加权统计关系，结合马尔科夫链蒙特卡洛随机模拟技术产生储层物性参数、弹性参数随机联合样本空间作为目标函数求解样本空间。最后，建立解的快速求解准则，求取最大后验概率密度对应的储层物性参数取值作为最终解。实际应用表明，该方法具有较高的反演效率，应用前景较好。

简介：Tikhonov正则化(TR)方法在重磁数据处理中发挥了重要的作用,本文在研究如何利用Tikhonov正则化方法方法解决重力数据3D反演的同时,深入讨论了可进一步提高拟合误差的ExtrapolationTikhonov正则化方法(EXTR)的原理,并就其参数选择方法及各参数对拟合误差、迭代次数及反演结果的影响进行研究。常密度及变密度组合模型试算结果表明,与TR方法相比,EXTR方法不仅可以达到解释人员设定的先验拟合误差水平,在计算时间及迭代次数相应增加的前提下有更高的拟合精度;同时其反演结果也更加紧致,进一步改善了TR反演结果的发散性;并且其反演数据范围更贴近预设模型参数范围,模型特征与预设模型密度分布吻合较好。更多还原

简介：本文提出了一种矿山微地震震源函数反演的新方法。以球坐标系波动方程为基础，对时空域中的微地震观测信号进行层析成像投影到慢度时间域，我们可以获得微地震波成像能量最大值及其对应的震源位置、发震时间和传播速度以及慢度时间信号。通过分析能量最大值处的慢度时间信号与震源函数之间的关系，推导出了计算矿山微地震震源函数的层析成像公式和利用震源函数计算微地震有效辐射能量的公式。进一步利用最小二乘法，将震源函数的振幅谱拟合成模型震源谱，确定了微地震波的零频极限值和拐角频率，最终计算出该微震事件的震源破裂半径。利用这一方法对理论模型数据和实际资料处理，结果表明：本方法具有既能求出任意一个地震事件的震源位置、发震时间和传播速度，又能同时求出震源函数及其频率特征值和震源破裂半径等震源参数的优点，因此，本方法是一种实时快速计算方法。

简介：本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法（FastFourierTransformMovingAverage,FFT-MA）和逐渐变形算法（GradualDeformationMethod,GDM）得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。

简介：描述和估计介质含水量、介电常数等属性参数分布是探地雷达探测技术的重要研究内容。雷达波的旅行时间和反射振幅系数与介质含水量、孔隙度与介电常数密切相关。常规通过旅行时间计算波速以估计介质参数的方法,例如透射波法,共中心点速度分析等,对于复杂介质分辨率有限。基于反射振幅的阻抗反演方法可以直接根据反射系数计算雷达波阻抗以估计介质属性参数,从而有效地避开常规方法在计算波速时精度低的问题。本文首先建立了基于高斯型和指数型混合自相关函数的三维多尺度等效随机介质模型刻画地下随机介质参数分布,并在局部加入高斯椭圆方程描述局部随机异常目标体。其次,通过引入锥形函数以降低随机介质模型在离散网格数值计算方法误差。在此基础上,推导了探地雷达递推阻抗反演的基本流程并结合随机介质模型测试了该方法在复杂介质参数估计中的计算精度。最后,对内蒙地区的实测探地雷达数据利用递推阻抗反演方法来估计地下污染物参数,估计介电常数、含水量结果与钻孔实测数据和同期开展的电阻率成像结果有很好的吻合。说明基于递推阻抗反演方法在探地雷达复杂介质属性参数估计中具有很好的应用前景。

简介：针对转换波各向异性叠前时间偏移实践中的四参数速度模型估计中参数间相互影响、难于准确的确定,本文将简化的转换波两参数动校正方程从叠加速度分析拓展到叠前时间偏移速度模型修正,形成了基于两参数动校正方程的新的四参数法转换波叠前时间偏移速度修正新方法和流程。在该方法和流程中,先进行转换波两参数叠加速度分析,获得转换波叠前时间偏移初始的速度和各向异性参数,然后通过分析偏移后的共成像点道集中的剩余时差修正速度和各向异性参数。其中叠前时间偏移速度模型的垂直速度比仍是需要利用纵波和转换波叠加剖面,层位标定后,利用相关的方法获得。有效速度比需结合纵波速度分析,利用Thomsen（1999）公式计算得到初始模型,最终依靠百分比扫描偏移处理获得最终的有效速度比模型。该方法简化了转换波高质量成像的速度估计方法,减小了多个参数估计的不确定性,在实际应用中也取得了较好的成像效果。

简介：河流动力学研究需要河底泥沙物性参数（孔隙率、渗透率和波速等）作为依据,机械取样和浅地层剖面探测是获取物性的重要途径。取样时的机械扰动使测试结果产生偏差,仅利用浅地层剖面数据获取河底特性有限,本文将两者结合起来开展反演研究。将取样测试的级配作为先验信息,进而根据Kozeny-Carman公式确定孔隙率和渗透率的关系。从浅地层剖面数据中提取了水-泥沙界面的声波反射系数。基于等效密度流模型,结合Kozeny-Carman公式和声波反射系数,提出了一种表层淤积泥沙参数反演方法。经过黄河库区试验,得到多个断面的密度和波速等参数,并获取泥沙物性参数的空间变化特征。对比发现,取样点的反演结果与测试结果较为一致,证明了本文提出泥沙参数反演方法的有效性。

简介：基于Cauchy先验分布的贝叶斯AVO反射率反演,可以产生类似于稀疏脉冲反演的结果,增强了对大反射系数的识别能力,反演过程中正则项和参数矩阵的计算都需要用到由上一步迭代所估计出的模型参数信息,因而反问题呈现非线性,此外该过程依赖于模型参数之间的线性统计关系,而且反演结果为反射率而非弹性参数,不利于储层预测和流体检测。贝叶斯AVO波形反演通过将反射率改写为弹性参数差分形式,直接从叠前地震数据中提取弹性参数,无需对弹性参数之间关系做线性假设。本文综合考虑上述两种方法的优点,在研究过程中仍然采用Cauchy先验分布,同时对贝叶斯AVO反射率的反演过程进行了修改,实现了基于Cauchy先验分布的AVO弹性参数弱非线性波形反演。本方法的提出有效避免了模型参数之间的线性假设,同时也能够从地震数据中直接反演得到纵、横波速度和密度。理论合成数据实验证明,此方法可以直接从叠前地震数据中提取弹性参数,而且在含有噪音的情况下也能得到比较准确的反演结果。

简介：弹性参数在甜点区预测和页岩气的开发过程中扮演着重要的角色，因此研究等效弹性参数随页岩气储层属性的变化是一项很有意义的工作。研究中我们用x射线CT扫描技术获得了较为精确的页岩样品微观结构图像。从这些图像中，我们可以获得孔隙度和矿物的详细情况，据此，我们构建了三维数字岩心，并应用有限元法对弹性参数进行了数值模拟，其间深入考察了子样选取、网格划分、求解器类型以及边界条件等，该方法易于区别不同的矿物及其百分含量。本文重点研究孔隙度和干酪根含量对弹性参数的影响，计算结果表明，孔隙度和干酪根含量对弹性性质有较大的影响，当孔隙度和干酪根含量增加时，弹性模量降低，且当孔隙度小于0．75％左右、干酪根含量大于3％左右时弹性参数减小速率较缓。因为孔隙度仅仅为4．5％，孔隙中填充油或气对弹性参数的影响甚微。不同岩心样本具有不同的孔隙度和干酪根含量，传统岩石物理实验不仅昂贵而且费时，而数值模拟是基于数字岩心来计算弹性参数，更加经济、方便。本研究证实了将页岩样品的微观结构图像与弹性模量的计算相结合来预测页岩弹性参数的可行性。

简介：基于曾新吾和MacBeth提出的时域内双源累积旋转方法,建立了频率相关介质中分析多分量VSP数据的横波分裂参数提取算法（DCTF）。该算法可以在频域中针对单个频率提取横波各向异性参数（快横波的极化方向以及快、慢横波间的时间差）,从而避免了目前常用方法中使用窄带通滤波可能带来的误差。通过对地震合成记录的数值分析,确定了该算法的可行性和正确性,并与目前常用方法的应用结果进行了比较。结果显示,频率相关横波分裂参数可以利用DCTF从地震四分量数据中直接提取。在地震频率范围内,含较大尺度裂缝时各向异性参数将表现出频率相关性,这意味着在地震频率范围会出现频散。随着频率的增加,各向异性有降低的趋势。

简介：目前在研究地球物理参数的空间变异特性时常常是基于几何各向异性假设,然而实际中几何各向异性假设常常不成立。本文提出了一种求解地球物理参数各向异性的新方法,首先对各方向上的变异值作归一化处理,然后用拟合的方法得到各方向单位步长上的变异性强度,进而实现对全局各向异性的求解。以上海地区合成孔径雷达干涉测量（InSAR）中差分大气延迟样本数据为例,验证了新方法求解各向异性的优越性：新方法的偏差只有6.4%,而基于几何各向异性假设方法的偏差达到了21.2%,而且新方法的运算效率有非常显著的提高。进一步,通过克里金插值交叉验证实验,证明了基于新算法得到的结构函数的插值效果最好,从而从另外一个角度验证了新算法求解的各向异性更准确,更好地表征了区域化变量的空间结构特征。因此,与几何各向异性假设的方法相比,新方法不仅能更准确的描述地球物理参数的各向异性,而且有更高的运算效率,为更准确地估计出所需的地球物理参数奠定了基础。

简介：在具有垂直对称轴横向各向同性介质中,利用四种参数来确定中间至远偏移距转换波(C-波)动校正.它们是C-波叠加速度VC2,垂直速度比和有效速度比γ0和γeff,以及各向异性参数χeff.我们将这四种参数作为C波叠加速度模型.C-波速度分析的目的就是确定这种叠加速度模型.C-波叠加速度模型VC2,γ0,γeff,和χeff可以由P-波和C-波反射动校正资料获得.然而错误的传播是C-波反射动校正反演中的严重问题.当前短排列叠加速度由于是从双曲线动校正推算而得,因而其精度不足以为各向异性参数提供有意义的反演值.中间偏移非双曲线动校正不再被人们所勿略,而是可以用一个背景γ加以量化.非双曲线分析通过中间偏移距的γ校正量可以产生VC2,若数据不含燥音,其误差小于1%.方法稳健,允许γ启始假定值的误差达20%.该方法也适用垂直非均匀各向异性介质.精度的提高使能够用4分量地震资料计算各向异性参数.为此提出了两种工作流程:双扫描和单扫描流程.理论数据和实际数据的应用表明这两种流程得出的结果其精度相似,但是单扫描流程比双扫描更有效.

简介：致密砂岩具有一定的各向异性。本文通过不同方向样品的声速测量及相关配套实验对致密砂岩进行了孔隙结构及各向异性分析，同时与页岩进行了对比。结果表明，致密砂岩样品所表现出的各向异性主要是由沉积环境差异形成的成分层、薄互层造成，具有典型的TI介质特征。不同方向杨氏模量随围压增加而增大，泊松比变化不明显。各向异性系数均随有效压力增加而减少，同时ε与γ和δ之间存在一定的线性关系。最后总结了不同地区的各向异性规律，为不同地区选择适合的评价模型奠定了基础。本次研究可为致密砂岩的测井评价、地震资料解释及压裂开发提供实验依据。

简介：用基于张量格林函数的体积分方程法对三维异常体进行瞬变电磁响应的正演模拟,首先在频率域内计算电磁场分量的频率域响应,然后利用快速数字滤波技术将计算结果转换到时间域。设计和计算了水平电偶极子源激发下层状水平地层模型背景下的常见地形如山谷、山峰地形的模型,并考察分别把源和接收器放于这些地形中的瞬变电磁场响应,详细分析了这些地形对长偏移距瞬变电磁测深（LOTEM）的影响。结果表明,山谷和山峰地形对LOTEM的结果均有不同程度的影响。当电偶极子源放在山谷谷底时,地形对观测异常场的畸变非常严重;当接收器放在山谷中时,接收器处地形的影响强烈但该影响在空间和时间上只是局部的。总体来讲,不论山峰地形位于何处,其对LOTEM的影响相对较小。当地形处于发射源与接收器之间时,地形对LOTEM的影响非常小,表明在进行LOTEM勘探时,选择发射源的放置比接收器的位置更加重要,野外勘探是尽量把发射源选择在开阔的平坦位置。

简介：为了解决海底起伏变化对地震波场的影响问题,本文提出将（x-z）域中的曲网格映射成（ξ,η）域中的矩形网格,推导出（ξ,η）域中的二维标量声波方程,根据推导出来的波动方程采用逆时有限差分法将海面上采集到的地震波场在（ξ,η）域中向下延拓至海底面,延拓时采用海水的速度,然后采用顺时有限差分法将延拓后的地震波场再反延拓到海面上,延拓时采用海底面以下地层的速度,从而消除了海底起伏带来的负面影响。模型及实际地震资料的计算分析表明该方法不但能够校正由于海底起伏所引起的海底面下地层反射波场的不连续性还能够校正由于海底起伏所引起的地震波的动力学特征的变化。对延拓前后的地震波场进行速度反演,延拓后反演的地层速度比延拓前反演的地层速度的精度提高很多,延拓前后地震波场的叠加剖面对比表明该延拓方法能够明显提高地震波场的成像质量。

简介：在多个区块的测井评价工作发现,孔隙结构直接影响储集层的品质和油气层的电阻率,是测井准确评价流体性质的关键。岩石物理资料表明不同区块内影响储层孔隙结构的因素不同,但效果是一致的,即孔隙结构的复杂程度控制着储层的储集能力和渗透能力。孔隙结构的复杂程度影响了储层中导电流体的分布和含量,从而控制了储层的电阻率。储层出现低阻油气层的内因均为复杂的孔隙结构（骨架导电及工程原因除外）。测井储层评价在分析控制储层孔隙结构复杂程度的地质因素及储层分类的基础上,针对不同类型储层采用不同的模型、参数和标准,可以有效的认识储层品质和识别不同类型储层的流体性质。

简介：在水库现场试验了RISK2型探地雷达探测水库冰厚度的能力,试验时所用天线频率为600MHz;同步钻孔测量雷达探测处的冰厚度;以及在一个点上取样测试分析冰晶体、冰内气泡和冰密度。试验时冰面积雪厚度0.03-0.05m,冰层上部有0.24m粒状冰,其下均为柱状冰;冰内气泡含量呈表层高底层低分布;冰密度随气泡含量变化;冰厚度在平面内不均一。通过探测厚度和实测厚度的对比分析以及气泡含量对介电系数影响的理论分析,建立了积雪、粒状冰和柱状冰三层介质模型,获取雷达波在冰内的理论传递时间。结果发现：能够利用等效介电常数或等效传播速度评价雷达波传递时间,结冰期冰层1/3深度处的对应介电常数或传递速度可以作为等效值;另外因冰内大气泡造成的理论传递时间大于雷达探测时间,其差值随理论传递时间或冰厚的增加呈非线性增加。

简介：在分析Ziolkowski气泡振动模型局限性的基础上,建立了多种实际因素影响下的海上地震勘探单枪子波模型。此模型考虑了气泡壁的热传导作用、枪口节流作用、气泡上浮、液体粘度和枪体本身等对气枪子波的影响。相对于Ziolkowski模型,改进模型所模拟的气枪子波主峰值减小,气泡振动衰减加快,与实测子波吻合性较好。实验分析表明：（1）枪口节流作用控制着气枪脉冲峰值振幅的大小,（2）上浮过程中气泡周围静水压力值减小,气泡振动的周期随之改变,（3）热传导作用和流体粘度是引起气泡振动衰减的主要因素。

简介：在CSEM野外工作中,通常把接收机置于离源3～5倍趋肤深度以外的远区。本文从正演和反演两个角度研究了当源和接收机之间存在导电体时,该导电体对远区中目标体响应的影响。2D有限元正演结果表明导电体主要影响中低频的观测结果,导电体的电导率越低,尺寸越大,影响也越大,并且影响程度随着观测频率的降低而逐渐加大。反演结果表明若不考虑源和勘探区间导电体的存在,反演得到的目标体横向位置会向源的方向移动,带来解释上的误差。在CSEM实际工作中,应该是采用多方位的场源作三维采集和三维反演,以减少源和勘探区间低阻体的影响。更多还原