简介：提高地震资料的分辨率是地震数据处理流程中的重要环节,对后续的精细构造解释起到重要作用。传统的提高分辨率方法大都假设地震资料是稳态的并且噪声水平不随空间发生变化,而实际情况不满足这一假设,导致提高分辨率处理后的效果达不到预期要求。针对这一问题,本文提出了一种基于时频二次谱的提高地震资料分辨率方法。首先,文中提出了基于S变换的时频二次谱,并结合模型论述了时变子波和反射系数在时频二次谱中的特征及其可分离性;其次,依据时变子波和反射系数在时频二次谱中的特征差异,构建了二维滤波器在地震记录的时频二次谱中提取时变子波的振幅谱;再次,文中研究了噪声环境中时变提高分辨率算子设计方法,并提出了依据时频谱能量强弱相对关系自适应确定频带拓宽范围的时变提高分辨率算子设计,进行提高分辨率处理;最后,文中对该方法进行了模型和实际数据的试处理,并与传统谱模拟方法和Q补偿方法的处理结果进行了对比分析,对比结果表明：本方法不需要估计Q值,提高分辨率能力不受震源子波频带的限制,在兼顾信噪比的前提下能够充分提高不同时间局部的地震数据的分辨率。

简介：常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。

简介：地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层（PML）边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件（HABC）具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。

简介：角度域共成像点道集(ADCIGS)是偏移速度分析和振幅随角度变化分析(AVA)的基础数据。传统Kirchhoff叠前深度偏移(KPSDM)按偏移距组织数据,能方便的输出偏移距域共成像点道集(ODCIGS),其高效的角度道集输出是有挑战的。本文提出基于旅行时梯度场的KPSDM角道集输出方法。其核心步骤为:(1)利用任意介质中的动态规划法旅行时计算方法提供炮点和检波点的旅行时场;(2)根据旅行时场的梯度方向计算反射张角;(3)在偏移过程中抽取ADCIGS。由于本文旅行时计算方法没有射线阴影区,也没有对速度光滑性的要求,其角度道集输出在阴影区比传统射线追踪更有优势。基于该角度道集输出方法,本文发展了一种适合大规模三维地震数据的KPSDM及角道集输出的并行实现方案。其基本思想是:(1)按照炮数据来组织输入数据;(2)旅行时场的输入与单炮覆盖范围相联系以节省内存;(3)多炮数据间采用MPI并行处理,单炮深度切片之间采用OpenMp并行处理,可进一步提高内存利用率和并行力度。数值试验结果证明本文角度道集生成方法的优越性和本文实现方案的有效性。

简介：被动源地震数据包含丰富的低频信息,本文有效地提取并利用这些信息对缺失低频的主动源地震数据进行低频重构,提出了基于多正弦窗的被动源多窗谱重构方法,并给出了相应的多震源多道重构公式。与常规互相关法和常规反褶积法重构的被动源记录相比,该方法能重构出更为准确的相对振幅信息。通过分析被动源数据重构前后的频谱特性,发现被动源的低频特性在重构和去噪处理后能更明显的体现出来。并提出了一种用被动源数据重构主动源低频信息的方法,即在功率谱上进行匹配,并在频域进行补偿和平滑。最后进行了数值算例的验证,对低频重构后的数据进行了叠前深度偏移处理。能量匹配方法能够用被动源的低频信息有效地重构主动源缺失的低频信息,低频重构后的记录在偏移成像中能体现更多的细节信息和深部构造。

简介：岩石物理建模及脆性指数构建是影响脆性预测精度的两大重要环节。现有页岩模型对有机质的模拟争议较大，需要寻找合理的岩石物理理论来模拟富有机质页岩。同时，现有脆性公式种类繁多，各公式的适用性值得探究。本文利用Self-ConsistentApproximationandtheDifferentialEffectiveMedium（SCA＋DEM）理论，通过模拟有机质与粘土的耦合性，构建各向异性富有机质岩石物理模型。与前人理论对比，初步验证了本模型的有效性；同时，基于模型构建脆性模板，分析物性参数对各脆性指数公式的影响。结果显示：各脆性公式对不同物性条件下地层的敏感性不同，基于杨氏模量构建的脆性指数对矿物含量的变化较敏感，而基于拉梅系数构建的脆性公式对孔隙度／孔隙流体敏感。应综合各脆性指数公式并结合地层物性信息，以达到最优的预测结果。

简介：河流相储层是一种重要的储层类型，寻找河道砂体是油气勘探的一项重要工作。但是河道变迁比较频繁，叠置比较严重，加之地震资料本身信噪比和分辨率的限制，河流相储层的精细描述至今仍存在一定的困难。本文以胜利油田老河口三维工区为例，首先研究了本区河流相储层的基本特征，发现该类储层剖面上呈强相位、切片上连续性好和频谱上呈低频特征；通过颜色处理和地震精细解释，发现本区河流相储层存在类串珠状特征，为了验证这一观测结果，我们用正演模拟分析了其形成机理，为在类似工区识别此类古河道找到了新的途径；应用多属性融合和RGB显示技术，河道下切特征更加明显，很好地彰显了河道结构特征，提高了河道识别能力；最后，我们研究并开发了多子波检测技术，检测出更多河流相弱反射信息。

简介：标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。

简介：MostcurrentprestackAVAjointinversionmethodsarebasedontheexactZoeppritzequationanditsvariousapproximations.However,theseequationsonlyreflecttherelationbetweenreflectioncoefficients,incidenceangles,andelasticparametersoneithersideoftheinterface,whichmeansthatwave-propagationeffects,suchassphericalspreading,attenuation,transmissionloss,multiples,andeventmismatchingofP-andS-waves,arenotconsideredandcannotaccuratelydescribethetruepropagationcharacteristicsofseismicwaves.ConventionalAVAinversionmethodsrequirethatthesewave-propagationeffectshavebeenfullycorrectedorattenuatedbeforeinversionbuttheserequirementscanhardlybesatisfiedinpractice.Usingaone-dimensional(1D)earthmodel,thereflectivitymethodcansimulatethefullwavefieldresponseofseismicwaves.Therefore,weproposeanonlinearmulticomponentprestackAVAjointinversionmethodbasedonthevectorizedreflectivitymethod,whichusesafastnondominatedsortinggeneticalgorithm(NSGAII)tooptimizethenonlinearmultiobjectivefunctiontoestimatemultipleparameters,suchasP-wavevelocity,S-wavevelocity,anddensity.Thisapproachisrobustbecauseitcansimultaneouslycopewithmorethanoneobjectivefunctionwithoutintroducingweightcoefficients.Modeltestsprovetheeffectivenessoftheproposedinversionmethod.Basedontheinversionresults,wefindthatthenonlinearprestackAVAjointinversionusingthereflectivitymethodyieldsmoreaccurateinversionresultsthantheinversionbyusingtheexactZoeppritzequationwhenthewave-propagationeffectsoftransmissionlossandinternalmultiplesarenotcompletelycorrected.

简介：我们改进了共偏移距-共反射面（COCRS）法，可用以衰减地滚波，即由于低速、低频／高振幅瑞雷波通常产生的相干噪声。COCRS算子是基于双曲线，因此它可以拟合双曲走时的同相轴，如叠前数据中的反射同相轴。相反，地滚波在共中点（CMP）和共炮点道集中是线性的并可以可以利用COCRS算子鉴别与压制。因此，我们在共偏移距剖面之前共炮道集中搜索反射倾斜和曲率。因为这对反射振幅的危害最小化是最理想的，我们只对在地滚波区多次覆盖的数据进行叠加。在CO剖面前搜索CS道聚集是对常规COCRS叠加的另一个改进。我们使用合成和真实数据集测试了所提出的方法，数据采自伊朗西部地区。我们将本方法压制地滚波的结果与f-k滤波和f-k滤波后常规COCRS叠加压制地滚波的结果进行了比较。结果表明，该方法对真伪滚压制效果优于F-K滤波与传统CRS叠加。然而，计算时间高于其他常规的方法，如f滤波。

简介：地震波场正演模拟是地震资料处理、解释中最为重要的技术之一。地震波场正演模拟在大时间步长、长时程的波场延拓中,存在计算不稳定的问题。本文基于声波方程的Hamilton表述,在波动方程求解中用辛差分格式进行时间网格离散,用傅里叶有限差分进行空间网格离散,提出一种新的保结构地震波场正演模拟方法一辛格式傅里叶有限差分法,在保证计算精度的同时提高计算的稳定性。利用声学近似处理空间-波数混合域的积分算子,将该方法推广至各向异性介质。给出各向同性和各向异性条件下的地震正演模拟的计算流程,并将本文方法用于BP盐丘、BPTTI等模型的波场正演模拟。数值算例表明本文开发的方法适用于速度变化剧烈的复杂介质地震波场正演模拟,计算精度高,数值频散小,在各向异性介质正演中能够有效避免qSV波残余,在大时间步长的迭代计算中稳定性好。本文为在辛算法的框架下实现高精度地震正演模拟提供了一种新的选择。

简介：Tilt—depth法可以用于快速反演磁源的上顶埋深，但其无法反演磁源底部深度，同时该方法的反演解过于单一。针对于此，本文在有限厚度台阶磁异常导数公式基础上，推导出了可以同时反演地质体上、下界面埋深的改进型Tilt—depth法计算公式，并采用在Tilt梯度图上选取多特征点进行联合反演的计算模式来提高反演解的可靠性。二维及三维组合模型试验均证实了改进型Tilt—depth法可以有效地反演出磁性体的上顶与下底深度，且上顶埋深的反演精度明显高于常规方法的。最后将改进型Tilt—depth法用于松辽盆地长岭断陷航磁数据反演，其上顶埋深的反演结果与钻孔钻遇火山岩的深度基本一致，证实改进型Tilt—depth法相对于常规方法的反演结果更加准确。

简介：坐标变换法通过将物理空间的曲网格映射为计算空间的矩形网格,将起伏地表转化为水平地表,同时将物理空间的波动方程转化为计算空间的波动方程,在计算空间完成数值模拟,坐标变换的方法对处理起伏自由边界具有较好的适应性和应用效果。本文在传统坐标变换方法的基础上,根据计算区域速度差异采用不同的网格大小和采样时间步长,提出了一种基于时空双变网格的起伏地表坐标变换正演模拟方法。在编程实现算法的基础上,通过典型模型波场模拟试算结果分析可知：（1）变网格方法与常规方法波场模拟误差在0.5%左右;（2）变网格方法计算效率视不同的变网格区域面积及变网格大小可提高几倍量级,在本文模型和计算参数下提高约5倍。（3）在满足模拟精度及频散条件要求下,变网格方法较全局细网格算法能显著节约计算内存。为此,针对起伏地表数值模拟,本文方法具有较高的模拟计算精度和一定的适应性。

简介：JointPP–PSinversionoffersbetteraccuracyandresolutionthanconventionalP-waveinversion.P-andS-waveelasticmodulideterminedthroughdatainversionsarekeyparametersforreservoirevaluationandfluidcharacterization.Inthispaper,startingwiththeexactZoeppritzequationthatrelatesP-andS-wavemoduli,acoefficientthatdescribesthereflectionsofP-andconvertedwavesisestablished.Thismethodeffectivelyavoidserrorintroducedbyapproximationsorindirectcalculations,thusimprovingtheaccuracyoftheinversionresults.Consideringthattheinversionproblemisill-posedandthattheforwardoperatorisnonlinear,priorconstraintsonthemodelparametersandmodifiedlow-frequencyconstraintsarealsointroducedtotheobjectivefunctiontomaketheproblemmoretractable.Thismodifiedobjectivefunctionissolvedovermanyiterationstocontinuouslyoptimizethebackgroundvaluesofthevelocityratio,whichincreasesthestabilityoftheinversionprocess.TestsofvariousmodelsshowthatthemethodeffectivelyimprovestheaccuracyandstabilityofextractingPandS-wavemodulifromunderdetermineddata.Thismethodcanbeappliedtoprovideinferencesforreservoirexplorationandfluidextraction.

简介：特征值相干属性是一种高效、成熟的裂缝预测技术。但即使在连续、稳定的目的层段,部分裂缝隐没于相近高灰度值区域而难以识别;在破碎带发育区域,由于成像能量难以完美聚焦,中、小尺度的裂缝区域呈现相近低灰度值的云雾状模糊而无法区分。本文提出一种以直方图均衡化处理为核心技术的新型裂缝增强方法。该方法能够强化相干图像中不连续性信息与背景的差异性,突出裂缝的线性结构:采用相干图像阈值逐级调节的方式,实现对不同尺度裂缝的预测。该方法同时提升了显性、隐性裂缝识别能力和精度。更多还原

简介：边界识别是重磁数据解释中的常用方法之一,依据其结果可划分出地质体的水平范围。边界识别结果受地质体埋深及导数计算误差的影响所识别边界与真实边界之间存在一定的差距,且边界识别法无法直观地给出地质体的深度信息。为了获得异常体的水平位置和深度信息,本文提出空间归一化边界识别方法,其对不同深度的边界识别函数进行归一化计算,空间归一化边界识别法的最大值对应于异常体的水平位置和深度。常规边界识别结果的误差随理深的减小而减小,而空间归一化边界识别法是通过最大值来判断地质体的位置,最大值是在地质体处获得,因此归一化边界识别方法所获得的结果是准确的。通过理论模型试验证明归一化边界识别方法能有效地完成异常体的水平位置和深度的计算,所获得的水平位置和深度信息与理论值相一致,为下一步的勘探计划提供了更加可靠的依据。将其应用于实际航磁数据的解释,获得了断裂的具体分布形式。更多还原

简介：超宽带信号的数字化在很大程度上取决于采样时刻的准确性,因此是否具有高精度的可编程定时器是实现超宽带信号数据采集的关键.本文基于斜波发生器原理实现了高精度可编程定时电路,可编程计数范围为16bit,实际最小定时单位达8ps,折合采样频率高达50GHz,迄今为止尚未有关于同样指标产品的报道.该定时器已应用于作者研制的地质雷达信号数据采集系统中,取得了很好的实际应用效果.

简介：在实验室对5块储层砂岩进行了模拟地层压力条件下的超声波速度测试。砂岩样品采自WXS凹陷的W地层,覆盖了从低到高的孔隙度和渗透率范围。实验选用了卤水和4种不同密度油作为孔隙流体,结合温度变化,实现了对流体粘度引致的速度频散研究。对实验结果的分析表明：（1）对于高孔隙度和渗透率的样品,无论是哪种流体饱和,观察到的超声波速度测试值和零频率Gassmann预测值的差异较小（约2-3%）,基本上可以用Biot模型解释;对于中等孔隙度和渗透率的样品,低粘度流体（〈约3mP·S）的频散效应也可以用Biot模型得到合理解释;（2）对于低、中孔隙度和渗透率样品,当流体粘度增加时,喷射流机制起主导作用,导致严重的速度频散（可达8%）。对储层砂岩的微裂隙纵横比进行了估计并用于喷射流特征频率的计算,当高于该特征频率时,Gassmann理论的假设条件受到破坏,实验室测得的高频速度不能直接用于地震低频条件下的W地层砂岩的Gassmann流体替换研究。

简介：天然气水合物沉积层在常规叠后地震剖面上，一般表现为似海底反射（BSR）特征，常作为水合物识别的重要标志，但由于地震数据分辨率的限制与多解性的存在，类似的反射特征不一定是水合物的表现。本文先从南中国海神狐海域地质条件分析人手，综述研究区域天然气水合物发育的构造、沉积与运移环境；其次以加权方程与Zoeppritz公式的正演模拟为基础，讨论天然气水合物与BSR特征的关系；然后选取经过三口已钻获天然气水合物站位的地震测线，提取角度域共成像道集用于BSR的AVA特征研究；最后依据叠前地震弹性参数同步反演思路，对水合物重点赋存区域开展精细预测研究，获取了研究区域含水合物沉积层稳定带的分布特征与饱和度数。该方法可作为采用多道地震数据估算天然气水合物含量的一种有效技术。

简介：由于孔隙型碳酸盐岩储层的复杂孔隙结构和强非均质性,岩石导电效率与含水孔隙度之间理想的线性关系相对复杂或并不存在,故前人基于岩石导电效率理论建立含水饱和度的计算公式并不完全适用于碳酸盐岩储层。基于岩石导电效率理论,推导了岩石导电效率的计算公式,阐明了岩石导电效率与含水孔隙度之间线性关系相对复杂或并不存在的根本原因,发现了岩石导电效率与电阻率之间的幂函数关系,分析了所建立含水饱和度计算公式误差的主控因素,得出了岩石导电效率的计算精度是该方法是否能推广应用的关键。与Archie公式相比,基于岩石导电效率理论建立的含水饱和度计算公式能更准确地计算研究靶区孔隙型碳酸盐岩储层的含水饱和度。在伊拉克某油田和印尼某气田3口井碳酸盐岩储层中的应用,表明当计算岩石导电效率的相对误差不大于0.1时,所计算储层含水饱和度的绝对误差不大于0.1,基本满足孔隙型碳酸盐岩储层精细评价的需求。