CN112346129A - 测井曲线提取及合成地震记录制作方法及装置 - Google Patents
测井曲线提取及合成地震记录制作方法及装置 Download PDFInfo
- Publication number
- CN112346129A CN112346129A CN202011181324.1A CN202011181324A CN112346129A CN 112346129 A CN112346129 A CN 112346129A CN 202011181324 A CN202011181324 A CN 202011181324A CN 112346129 A CN112346129 A CN 112346129A
- Authority
- CN
- China
- Prior art keywords
- depth
- amplitude
- curve
- time
- seismic data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 28
- 238000000605 extraction Methods 0.000 title claims description 36
- 238000000034 method Methods 0.000 claims abstract description 46
- 238000009792 diffusion process Methods 0.000 claims abstract description 45
- 230000005284 excitation Effects 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000012937 correction Methods 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 16
- 238000007781 pre-processing Methods 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 11
- 238000009795 derivation Methods 0.000 claims description 8
- 239000013307 optical fiber Substances 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 claims description 5
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 238000012886 linear function Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 230000001629 suppression Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 21
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000005755 formation reaction Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010892 electric spark Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Abstract
本发明公开一种测井曲线提取及合成地震记录制作方法及装置，包括：利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据；对井中地震数据进行初至拾取，得到深度‑时间关系曲线；所述初至拾取的初至位置为初至波峰位置；在所述初至位置拾取振幅值，得到深度‑振幅关系曲线和时间‑振幅关系曲线；根据时间‑振幅关系曲线计算球面扩散因子；根据球面扩散因子、深度‑振幅关系曲线和深度‑时间关系曲线，得到深度‑补偿振幅关系曲线；根据深度‑补偿振幅关系曲线，得到深度‑求导后补偿振幅曲线，作为提取到的测井曲线。利用提取得到的新型测井曲线，结合深度域空变雷克子波，制作出高分辨率的合成地震记录，可解决薄地层识别的难题。
Description
技术领域
本发明涉及地球物理勘探技术领域，尤其涉及测井曲线提取及合成地震记录制作方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
薄储层勘探一直以来就是地震领域的一个难点课题，如何有效提高地震数据的分辨率，是国内外地球物理学家一个重要的追求目标。目前行业内比较认可的极限垂向分辨率的为1/4波长理论，以50Hz地震主频，4000m/s速度计算，可识别的薄层为20m，其精度远远低于测井方法的识别能力，尽管围绕着提高地震分辨率的研究在持续发展，有学者已提出1/8至1/16波长等间接提升极限分辨率的方法；但均在一定的限制条件或借助其它已知地球物理成果才能突破1/4波长。单纯依靠地震数据，想要有效提升纵向分辨能力，需要将两个相互干涉的地震子波的包络清晰分开，即两个子波的频率、振幅、相位等参数完全分离。因此要想分辨小于1/4波长的薄储层仍是十分困难。
发明内容
首先，本发明实施例提供一种测井曲线提取方法，用以准确有效地提取测井曲线，有利于依据该测井曲线实现薄储层的分辨，该方法包括：
利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据；
对井中地震数据进行初至拾取，得到深度-时间关系曲线；所述初至拾取的初至位置为初至波峰位置；
在所述初至位置拾取振幅值，得到深度-振幅关系曲线和时间-振幅关系曲线；
根据时间-振幅关系曲线计算球面扩散因子；
根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线；
根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，作为提取到的测井曲线。
其次，本发明实施例还提供一种合成地震记录制作方法，该合成地震记录制作方法利用上述测井曲线提取方法提取到的测井曲线制作合成地震记录，用以解决薄地层识别的难题，该合成地震记录制作方法包括：
生成时间域高频雷克子波；
利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波；
利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录。
再次，本发明实施例还提供一种测井曲线提取装置，用以准确有效地提取测井曲线，有利于依据该测井曲线实现薄储层的分辨，该测井曲线提取装置包括：
井中地震数据获得模块，用于利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据；
初至拾取模块，用于对井中地震数据进行初至拾取，得到深度-时间关系曲线；所述初至拾取的初至位置为初至波峰位置；
振幅拾取模块，用于在所述初至位置拾取振幅值，得到深度-振幅关系曲线和时间-振幅关系曲线；
球面扩散因子计算模块，用于根据时间-振幅关系曲线计算球面扩散因子；
深度-补偿振幅关系曲线获取模块，用于根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线；
测井曲线获取模块，用于根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，作为提取到的测井曲线。
再次，本发明实施例还提供一种合成地震记录制作装置，该合成地震记录制作装置利用上述测井曲线提取装置提取到的测井曲线制作合成地震记录，用以解决薄地层识别的难题，该合成地震记录制作装置包括：
时间域高频雷克子波生成模块，用于生成时间域高频雷克子波；
深度域空变雷克子波转换模块，用于利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波；
合成地震记录制作模块，用于利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录。
再次，本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述测井曲线提取或合成地震记录制作方法。
最后，本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述测井曲线提取或合成地震记录制作方法的计算机程序。
本发明实施例中，利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据；对井中地震数据进行初至拾取，得到深度-时间关系曲线；所述初至拾取的初至位置为初至波峰位置；在所述初至位置拾取振幅值，得到深度-振幅关系曲线和时间-振幅关系曲线；根据时间-振幅关系曲线计算球面扩散因子；根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线；根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，作为提取到的测井曲线；从而基于井中地震数据，准确有效地提取得到一种新型的测井曲线，有利于依据该测井曲线实现薄储层的分辨。
本发明实施例中，利用测井曲线提取方法提取到的测井曲线制作合成地震记录，其中，生成时间域高频雷克子波；利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波；利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录；从而利用提取得到的新型测井曲线，结合深度域空变雷克子波，制作出高分辨率的合成地震记录，为解决薄地层识别的难题提供了一个有效的手段，对于薄地层识别、油气开发等方面非常具有实用价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中测井曲线提取方法的示意图；
图2为本发明实施例中所采用的速度模型示意图；
图3为本发明实施例中一个深度-补偿振幅关系曲线的示例图；
图4为本发明实施例中一个深度-求导后补偿振幅曲线的示例图；
图5为本发明实施例中合成地震记录制作方法的示意图；
图6为本发明实施例中一个合成高分辨率地震记录的示例图；
图7为本发明实施例中一个多道合成高分辨率地震记录的示例图；
图8为本发明实施例中测井曲线提取装置的示意图；
图9为本发明实施例中测井曲线提取装置的一具体示例图；
图10为本发明实施例中测井曲线提取装置的另一具体示例图；
图11为本发明实施例中测井曲线提取装置的又一具体示例图；
图12为本发明实施例中合成地震记录制作装置的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
发明人考虑到，随着井中地震仪器设备的不断进步，以分布式光纤为代表的许多高密度井中观测系统被应用，使得井中地震方法在薄地层的识别方面作用逐渐突显。基于此，考虑在本发明实施例中基于井中地震数据，提取一种新型的测井曲线，并依据该测井曲线制作高分辨率合成地震记录，为解决薄地层识别的难题提供一个有效的手段，一定程度上拓展井中地震应用领域。
图1为本发明实施例中测井曲线提取方法的示意图，如图1所示，该方法包括：
步骤101、利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据；
步骤102、对井中地震数据进行初至拾取，得到深度-时间关系曲线；所述初至拾取的初至位置为初至波峰位置；
步骤103、在所述初至位置拾取振幅值，得到深度-振幅关系曲线和时间-振幅关系曲线；
步骤104、根据时间-振幅关系曲线计算球面扩散因子；
步骤105、根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线；
步骤106、根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，作为提取到的测井曲线。
由图1所示流程可以得知，本发明实施例基于井中地震数据，可以准确有效地提取得到一种新型的测井曲线，有利于依据该测井曲线实现薄储层的分辨。
具体实施时，先利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据。为了实现更好的激发效果，实施例中，井中地震激发仪器可以包括可控震源仪器。例如可以采用激发一致性较好的可控震源仪器，回避使用炸药震源、气枪震源、电火花震源或重锤震源等。
实施例中，井中地震接收仪器可以包括井中大阵列检波器、井中光栅节点仪器或井中分布式光纤。例如可以采用20级以上井中大阵列检波器、井中光栅节点仪器或井中分布式光纤(DAS)，由于本发明实施例对于记录一致性要求较高，因此不推荐采用数量较少的井中数字检波器或模拟检波器。
实施例中，可以获得高密度井中地震数据，例如可以获得深度采样间距不大于1m，记录时间间隔不大于0.1ms的高密度井中地震数据。
实施例中，获得井中地震数据后，还可以包括：对井中地震数据进行一致性校正。例如，可以包括：对井中地震数据进行激发点间的一致性校正，以及检波器间的一致性校正；其中，检波器布置于在激发点附近地面上，用于监控井中地震激发仪器的输出状态。激发点间的一致性校正，例如可以利用地面监控子波的记录数据，统计激发点间的能量变化，进而完成校正。检波器间的一致性校正，例如可以在检波器下井采集前，进行地面一致性校正；采集完成后，统计多检波点间的能量变化规律，进行一致性校正。
实施例中，获得井中地震数据后，还可以包括：对井中地震数据进行预处理，所述预处理保护数据的原始振幅信息。例如，可以对井中地震数据进行坏道去除处理，和/或，异常噪声压制处理等步骤。不推荐使用同深度叠加、自动增益、道均衡处理等改变振幅相对关系的处理步骤，以保护数据的原始振幅信息。实施例中，可以对一致性校正后的井中地震数据进行预处理。
实施例中，获得井中地震数据后，还可以包括：将井中地震数据处理为相等采集间距的井中地震数据。例如，可以将获得的预处理后数据，插值或抽取成为相等采集间距的井中地震数据。
对井中地震数据进行初至拾取，得到深度-时间关系曲线；所述初至拾取的初至位置为初至波峰位置。实施例中，可以对经一致性校正、预处理并处理为相等采集间距的井中地震数据，进行初至拾取，得到深度-时间关系曲线。例如，拾取井中地震数据的首波时间，得到深度-时间关系曲线t＝f1(h)。其中，t为时间；h为深度。
在所述初至位置拾取振幅值，得到深度-振幅关系曲线a＝f2(h)和时间-振幅关系曲线a＝f3(t)。其中，a为振幅。
实施例中，根据时间-振幅关系曲线计算球面扩散因子(TAR)，可以包括：
对时间-振幅关系曲线a＝f3(t)，对时间t和振幅a均取自然对数，得到曲线ln(a)＝f3(ln(t))，利用一次函数拟合该曲线，读取拟合曲线斜率k，作为球面扩散因子。
实施例中，根据时间-振幅关系曲线计算球面扩散因子，可以包括：在原始频率上，根据时间-振幅关系曲线计算球面扩散因子。需要说明的是，通常井中地震处理过程中，为了避免频率衰减影响，球面扩散因子要在降频后的数据中计算，而在本实施例中，为了使最后提取得到的测井曲线能够更好的实现薄储层分辨的效果，推荐在原始频率上计算。
实施例中，根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线，可以包括：
利用球面扩散因子k和深度-振幅关系曲线a＝f2(h)，得到补偿后的时间-补偿振幅关系曲线a′＝tk·f2(t)；在补偿后的时间-补偿振幅关系曲线a′＝tk·f2(t)中，代入深度-时间关系曲线t＝f1(h)，得到深度-补偿振幅关系曲线a′＝(f1(h))k·f2(f1(h))；其中，a′为补偿后振幅。图2为本发明实施例中所采用的一个速度模型，基于该模型数据采用成熟的数值模拟方法得到深度-补偿振幅关系曲线，同时该模型数据也用于实施成果可靠性的对比检验，图3为本发明实施例中一个深度-补偿振幅关系曲线的示例图。
实施例中，根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，可以包括：
利用深度-补偿振幅关系曲线a′＝(f1(h))k·f2(f1(h))，对不同深度的振幅求取导数计算，得到深度-求导后补偿振幅曲线
其中，a″为求导后补偿振幅。图4为本发明实施例中一个深度-求导后补偿振幅曲线的示例图。
本发明实施例还提供一种合成地震记录制作方法，该合成地震记录制作方法利用上述测井曲线提取方法提取到的测井曲线制作合成地震记录，可以解决薄地层识别的难题。图5为本发明实施例中合成地震记录制作方法的示意图，如图5所示，该合成地震记录制作方法可以包括：
步骤501、生成时间域高频雷克子波；
步骤502、利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波；
步骤503、利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录。
由图5所示流程可以得知，本发明实施例利用提取得到的新型测井曲线，结合深度域空变雷克子波，制作出高分辨率的合成地震记录，为解决薄地层识别的难题提供了一个有效的手段，对于薄地层识别、油气开发等方面非常具有实用价值。
具体实施时，生成时间域高频雷克子波。实施例中，可以生成主频大于100Hz、子波长度为20ms-50ms、采样间隔为与所述井中地震数据相同的、零相位的、正极性的标准雷克子波。
实施例中，利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波，可以包括：
利用所述深度-时间关系曲线t＝f1(h)，将时间域高频雷克子波at＝w(t′)，转换为深度域空变雷克子波
所述深度域空变雷克子波的深度采样间隔与所述测井曲线相同；
其中，t为时间；h为深度；at为时间域高频雷克子波；ah为深度域空变雷克子波；t′为时间域子波采样；h′为深度域子波采样。
实施例中，利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录，可以包括：
利用深度域空变雷克子波ah与所述测井曲线a″进行褶积计算，得到深度域的合成地震记录s；
其中，
a″为求导后补偿振幅；k为球面扩散因子；a＝f2(h)为深度-振幅关系曲线，
图6为本发明实施例中一个合成高分辨率地震记录的示例图。图7为本发明实施例中一个多道合成高分辨率地震记录的示例图。
本发明实施例中还提供了一种测井曲线提取装置和合成地震记录制作装置，如下面的实施例所述。由于这些装置解决问题的原理分别与测井曲线提取方法和合成地震记录制作装置方法相似，因此该装置的实施可以参见相应方法的实施，重复之处不再赘述。
图8为本发明实施例中测井曲线提取装置的示意图，如图8所示，该装置可以包括：
井中地震数据获得模块801，用于利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据；
初至拾取模块802，用于对井中地震数据进行初至拾取，得到深度-时间关系曲线；所述初至拾取的初至位置为初至波峰位置；
振幅拾取模块803，用于在所述初至位置拾取振幅值，得到深度-振幅关系曲线和时间-振幅关系曲线；
球面扩散因子计算模块804，用于根据时间-振幅关系曲线计算球面扩散因子；
深度-补偿振幅关系曲线获取模块805，用于根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线；
测井曲线获取模块806，用于根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，作为提取到的测井曲线。
一个实施例中，井中地震激发仪器可以包括可控震源仪器；井中地震接收仪器可以包括井中大阵列检波器、井中光栅节点仪器或井中分布式光纤。
一个实施例中，井中地震数据获得模块，具体可以用于：
获得深度采样间距不大于1m，记录时间间隔不大于0.1ms的井中地震数据。
如图9所示，一个实施例中，测井曲线提取装置还可以包括：
一致性校正模块901，用于在井中地震数据获得模块获得井中地震数据后，对井中地震数据进行一致性校正。
一个实施例中，一致性校正模块，具体可以用于：
对井中地震数据进行激发点间的一致性校正，以及检波器间的一致性校正；
其中，检波器布置于在激发点附近地面上，用于监控井中地震激发仪器的输出状态。
如图10所示，一个实施例中，测井曲线提取装置还可以包括：
预处理模块1001，用于在井中地震数据获得模块获得井中地震数据后，对井中地震数据进行预处理，所述预处理保护数据的原始振幅信息。预处理模块1001也可以设于图9所示的测井曲线提取装置中。
一个实施例中，预处理模块，具体可以用于：
对井中地震数据进行坏道去除处理，和/或，异常噪声压制处理。
如图11所示，一个实施例中，测井曲线提取装置还可以包括：
井中地震数据处理模块1101，用于在井中地震数据获得模块获得井中地震数据后，将井中地震数据处理为相等采集间距的井中地震数据。井中地震数据处理模块1001还可以设于图9或图10所示的测井曲线提取装置中。
一个实施例中，球面扩散因子计算模块，具体可以用于：
对时间-振幅关系曲线a＝f3(t)，对时间t和振幅a均取自然对数，得到曲线ln(a)＝f3(ln(t))，利用一次函数拟合该曲线，读取拟合曲线斜率k，作为球面扩散因子。
一个实施例中，球面扩散因子计算模块，具体可以用于：
在原始频率上，根据时间-振幅关系曲线计算球面扩散因子。
一个实施例中，深度-补偿振幅关系曲线获取模块，具体可以用于：
利用球面扩散因子k和深度-振幅关系曲线a＝f2(h)，得到补偿后的时间-补偿振幅关系曲线a′＝tk·f2(t)；
在补偿后的时间-补偿振幅关系曲线a′＝tk·f2(t)中，代入深度-时间关系曲线t＝f1(h)，得到深度-补偿振幅关系曲线a′＝(f1(h))k·f2(f1(h))；
其中，a为振幅；t为时间；h为深度；a′为补偿后振幅。
一个实施例中，测井曲线获取模块，具体可以用于：
利用深度-补偿振幅关系曲线a′＝(f1(h))k·f2(f1(h))，对不同深度的振幅求取导数计算，得到深度-求导后补偿振幅曲线
其中，a″为求导后补偿振幅。
图12为本发明实施例中合成地震记录制作装置的示意图，如图12所示，合成地震记录制作装置利用前述测井曲线提取装置提取到的测井曲线制作合成地震记录，合成地震记录制作装置可以包括：
时间域高频雷克子波生成模块1201，用于生成时间域高频雷克子波；
深度域空变雷克子波转换模块1202，用于利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波；
合成地震记录制作模块1203，用于利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录。
一个实施例中，时间域高频雷克子波生成模块，具体可以用于：
生成主频大于100Hz、子波长度为20ms-50ms、采样间隔为与所述井中地震数据相同的、零相位的、正极性的标准雷克子波。
一个实施例中，深度域空变雷克子波转换模块，具体可以用于：
利用所述深度-时间关系曲线t＝f1(h)，将时间域高频雷克子波at＝w(t′)，转换为深度域空变雷克子波
所述深度域空变雷克子波的深度采样间隔与所述测井曲线相同；
其中，t为时间；h为深度；at为时间域高频雷克子波；ah为深度域空变雷克子波；t′为时间域子波采样；h′为深度域子波采样。
一个实施例中，合成地震记录制作模块，具体可以用于：
利用深度域空变雷克子波ah与所述测井曲线a″进行褶积计算，得到深度域的合成地震记录s；
其中，
a″为求导后补偿振幅；k为球面扩散因子；a＝f2(h)为深度-振幅关系曲线，
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述测井曲线提取或合成地震记录制作方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述测井曲线提取或合成地震记录制作方法的计算机程序。
综上所述，本发明实施例中，利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据；对井中地震数据进行初至拾取，得到深度-时间关系曲线；所述初至拾取的初至位置为初至波峰位置；在所述初至位置拾取振幅值，得到深度-振幅关系曲线和时间-振幅关系曲线；根据时间-振幅关系曲线计算球面扩散因子；根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线；根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，作为提取到的测井曲线；从而基于井中地震数据，准确有效地提取得到一种新型的测井曲线，有利于依据该测井曲线实现薄储层的分辨。
本发明实施例中，利用测井曲线提取方法提取到的测井曲线制作合成地震记录，其中，生成时间域高频雷克子波；利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波；利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录；从而利用提取得到的新型测井曲线，结合深度域空变雷克子波，制作出高分辨率的合成地震记录，为解决薄地层识别的难题提供了一个有效的手段，对于薄地层识别、油气开发等方面非常具有实用价值。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (34)
1.一种测井曲线提取方法，其特征在于，包括：
利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据；
对井中地震数据进行初至拾取，得到深度-时间关系曲线；所述初至拾取的初至位置为初至波峰位置；
在所述初至位置拾取振幅值，得到深度-振幅关系曲线和时间-振幅关系曲线；
根据时间-振幅关系曲线计算球面扩散因子；
根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线；
根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，作为提取到的测井曲线。
2.如权利要求1所述的方法，其特征在于，所述井中地震激发仪器包括可控震源仪器；所述井中地震接收仪器包括井中大阵列检波器、井中光栅节点仪器或井中分布式光纤。
3.如权利要求1所述的方法，其特征在于，获得井中地震数据，包括：
获得深度采样间距不大于1m，记录时间间隔不大于0.1ms的井中地震数据。
4.如权利要求1所述的方法，其特征在于，获得井中地震数据后，还包括：
对井中地震数据进行一致性校正。
5.如权利要求4所述的方法，其特征在于，对井中地震数据进行一致性校正，包括：
对井中地震数据进行激发点间的一致性校正，以及检波器间的一致性校正；
其中，检波器布置于在激发点附近地面上，用于监控井中地震激发仪器的输出状态。
6.如权利要求1所述的方法，其特征在于，获得井中地震数据后，还包括：
对井中地震数据进行预处理，所述预处理保护数据的原始振幅信息。
7.如权利要求6所述的方法，其特征在于，对井中地震数据进行预处理，包括：
对井中地震数据进行坏道去除处理，和/或，异常噪声压制处理。
8.如权利要求1所述的方法，其特征在于，获得井中地震数据后，还包括：
将井中地震数据处理为相等采集间距的井中地震数据。
9.如权利要求1所述的方法，其特征在于，根据时间-振幅关系曲线计算球面扩散因子，包括：
对时间-振幅关系曲线a＝f3(t)，对时间t和振幅a均取自然对数，得到曲线ln(a)＝f3(ln(t))，利用一次函数拟合该曲线，读取拟合曲线斜率k，作为球面扩散因子。
10.如权利要求1所述的方法，其特征在于，根据时间-振幅关系曲线计算球面扩散因子，包括：
在原始频率上，根据时间-振幅关系曲线计算球面扩散因子。
11.如权利要求1所述的方法，其特征在于，根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线，包括：
利用球面扩散因子k和深度-振幅关系曲线a＝f2(h)，得到补偿后的时间-补偿振幅关系曲线a′＝tk·f2(t)；
在补偿后的时间-补偿振幅关系曲线a′＝tk·f2(t)中，代入深度-时间关系曲线t＝f1(h)，得到深度-补偿振幅关系曲线a′＝(f1(h))k·f2(f1(h))；
其中，a为振幅；t为时间；h为深度；a′为补偿后振幅。
12.如权利要求11所述的方法，其特征在于，根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，包括：
利用深度-补偿振幅关系曲线a′＝(f1(h))k·f2(f1(h))，对不同深度的振幅求取导数计算，得到深度-求导后补偿振幅曲线
其中，a″为求导后补偿振幅。
13.一种合成地震记录制作方法，其特征在于，所述方法利用权利要求1至12任一所述测井曲线提取方法提取到的测井曲线制作合成地震记录，包括：
生成时间域高频雷克子波；
利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波；
利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录。
14.如权利要求13所述的方法，其特征在于，生成时间域高频雷克子波，包括：
生成主频大于100Hz、子波长度为20ms-50ms、采样间隔为与所述井中地震数据相同的、零相位的、正极性的标准雷克子波。
15.如权利要求13所述的方法，其特征在于，利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波，包括：
利用所述深度-时间关系曲线t＝f1(h)，将时间域高频雷克子波at＝w(t′)，转换为深度域空变雷克子波
所述深度域空变雷克子波的深度采样间隔与所述测井曲线相同；
其中，t为时间；h为深度；at为时间域高频雷克子波；ah为深度域空变雷克子波；t′为时间域子波采样；h′为深度域子波采样。
16.如权利要求15所述的方法，其特征在于，利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录，包括：
利用深度域空变雷克子波ah与所述测井曲线a″进行褶积计算，得到深度域的合成地震记录s；
其中，
a″为求导后补偿振幅；k为球面扩散因子；a＝f2(h)为深度-振幅关系曲线，
17.一种测井曲线提取装置，其特征在于，包括：
井中地震数据获得模块，用于利用井中地震激发仪器和井中地震接收仪器，获得井中地震数据；
初至拾取模块，用于对井中地震数据进行初至拾取，得到深度-时间关系曲线；所述初至拾取的初至位置为初至波峰位置；
振幅拾取模块，用于在所述初至位置拾取振幅值，得到深度-振幅关系曲线和时间-振幅关系曲线；
球面扩散因子计算模块，用于根据时间-振幅关系曲线计算球面扩散因子；
深度-补偿振幅关系曲线获取模块，用于根据球面扩散因子、深度-振幅关系曲线和深度-时间关系曲线，得到深度-补偿振幅关系曲线；
测井曲线获取模块，用于根据深度-补偿振幅关系曲线，得到深度-求导后补偿振幅曲线，作为提取到的测井曲线。
18.如权利要求17所述的测井曲线提取装置，其特征在于，所述井中地震激发仪器包括可控震源仪器；所述井中地震接收仪器包括井中大阵列检波器、井中光栅节点仪器或井中分布式光纤。
19.如权利要求17所述的测井曲线提取装置，其特征在于，井中地震数据获得模块，具体用于：
获得深度采样间距不大于1m，记录时间间隔不大于0.1ms的井中地震数据。
20.如权利要求17所述的测井曲线提取装置，其特征在于，还包括：
一致性校正模块，用于在井中地震数据获得模块获得井中地震数据后，对井中地震数据进行一致性校正。
21.如权利要求20所述的测井曲线提取装置，其特征在于，一致性校正模块，具体用于：
对井中地震数据进行激发点间的一致性校正，以及检波器间的一致性校正；
其中，检波器布置于在激发点附近地面上，用于监控井中地震激发仪器的输出状态。
22.如权利要求17所述的测井曲线提取装置，其特征在于，还包括：
预处理模块，用于在井中地震数据获得模块获得井中地震数据后，对井中地震数据进行预处理，所述预处理保护数据的原始振幅信息。
23.如权利要求22所述的测井曲线提取装置，其特征在于，预处理模块，具体用于：
对井中地震数据进行坏道去除处理，和/或，异常噪声压制处理。
24.如权利要求17所述的测井曲线提取装置，其特征在于，还包括：
井中地震数据处理模块，用于在井中地震数据获得模块获得井中地震数据后，将井中地震数据处理为相等采集间距的井中地震数据。
25.如权利要求17所述的测井曲线提取装置，其特征在于，球面扩散因子计算模块，具体用于：
对时间-振幅关系曲线a＝f3(t)，对时间t和振幅a均取自然对数，得到曲线ln(a)＝f3(ln(t))，利用一次函数拟合该曲线，读取拟合曲线斜率k，作为球面扩散因子。
26.如权利要求17所述的测井曲线提取装置，其特征在于，球面扩散因子计算模块，具体用于：
在原始频率上，根据时间-振幅关系曲线计算球面扩散因子。
27.如权利要求17所述的测井曲线提取装置，其特征在于，深度-补偿振幅关系曲线获取模块，具体用于：
利用球面扩散因子k和深度-振幅关系曲线a＝f2(h)，得到补偿后的时间-补偿振幅关系曲线a′＝tk·f2(t)；
在补偿后的时间-补偿振幅关系曲线a′＝tk·f2(t)中，代入深度-时间关系曲线t＝f1(h)，得到深度-补偿振幅关系曲线a′＝(f1(h))k·f2(f1(h))；
其中，a为振幅；t为时间；h为深度；a′为补偿后振幅。
28.如权利要求27所述的测井曲线提取装置，其特征在于，测井曲线获取模块，具体用于：
利用深度-补偿振幅关系曲线a′＝(f1(h))k·f2(f1(h))，对不同深度的振幅求取导数计算，得到深度-求导后补偿振幅曲线
其中，a″为求导后补偿振幅。
29.一种合成地震记录制作装置，其特征在于，所述合成地震记录制作装置利用权利要求17至28任一所述测井曲线提取装置提取到的测井曲线制作合成地震记录，所述合成地震记录制作装置包括：
时间域高频雷克子波生成模块，用于生成时间域高频雷克子波；
深度域空变雷克子波转换模块，用于利用所述深度-时间关系曲线，将时间域高频雷克子波转换为深度域空变雷克子波；
合成地震记录制作模块，用于利用深度域空变雷克子波与所述测井曲线进行褶积计算，得到深度域的合成地震记录。
30.如权利要求29所述的合成地震记录制作装置，其特征在于，时间域高频雷克子波生成模块，具体用于：
生成主频大于100Hz、子波长度为20ms-50ms、采样间隔为与所述井中地震数据相同的、零相位的、正极性的标准雷克子波。
31.如权利要求29所述的合成地震记录制作装置，其特征在于，深度域空变雷克子波转换模块，具体用于：
利用所述深度-时间关系曲线t＝f1(h)，将时间域高频雷克子波at＝w(t′)，转换为深度域空变雷克子波
所述深度域空变雷克子波的深度采样间隔与所述测井曲线相同；
其中，t为时间；h为深度；at为时间域高频雷克子波；ah为深度域空变雷克子波；t′为时间域子波采样；h′为深度域子波采样。
32.如权利要求31所述的合成地震记录制作装置，其特征在于，合成地震记录制作模块，具体用于：
利用深度域空变雷克子波ah与所述测井曲线a″进行褶积计算，得到深度域的合成地震记录s；
其中，
a″为求导后补偿振幅；k为球面扩散因子；a＝f2(h)为深度-振幅关系曲线，
33.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至16任一所述方法。
34.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至16任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011181324.1A CN112346129A (zh) | 2020-10-29 | 2020-10-29 | 测井曲线提取及合成地震记录制作方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011181324.1A CN112346129A (zh) | 2020-10-29 | 2020-10-29 | 测井曲线提取及合成地震记录制作方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112346129A true CN112346129A (zh) | 2021-02-09 |
Family
ID=74355717
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011181324.1A Pending CN112346129A (zh) | 2020-10-29 | 2020-10-29 | 测井曲线提取及合成地震记录制作方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112346129A (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5642327A (en) * | 1994-06-08 | 1997-06-24 | Exxon Production Research Company | Method for creating a gain function for seismic data and method for processing seismic data |
| US20030030439A1 (en) * | 2001-08-09 | 2003-02-13 | Li Gao | Virtual steering of induction tool for determination of formation DIP angle |
| CN1797037A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种地震波波阻抗反演的方法 |
| US20080027650A1 (en) * | 2005-01-14 | 2008-01-31 | Xiaojun Huang | Method and Apparatus for Estimating Formation Slowness |
| CN104656140A (zh) * | 2014-12-01 | 2015-05-27 | 中国石油天然气集团公司 | 一种基于垂直地震衰减规律的中值滤波方法以及系统 |
| CN104656142A (zh) * | 2013-11-19 | 2015-05-27 | 中国石油天然气集团公司 | 一种利用垂直地震剖面与测井联合的地震层位标定方法 |
| CN105182408A (zh) * | 2015-08-28 | 2015-12-23 | 中国石油天然气集团公司 | 一种合成地震记录的制作方法和装置 |
| CN107831541A (zh) * | 2017-11-17 | 2018-03-23 | 中国石油天然气集团公司 | 基于高密度vsp数据的薄地层识别方法及装置 |
| CN108961352A (zh) * | 2017-05-19 | 2018-12-07 | 中国石油化工股份有限公司 | 一种随钻测井曲线的绘制方法 |
| CN111337977A (zh) * | 2018-12-19 | 2020-06-26 | 中国石油天然气股份有限公司 | 品质因子的确定方法及装置 |
-
2020
- 2020-10-29 CN CN202011181324.1A patent/CN112346129A/zh active Pending
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5642327A (en) * | 1994-06-08 | 1997-06-24 | Exxon Production Research Company | Method for creating a gain function for seismic data and method for processing seismic data |
| US20030030439A1 (en) * | 2001-08-09 | 2003-02-13 | Li Gao | Virtual steering of induction tool for determination of formation DIP angle |
| CN1797037A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种地震波波阻抗反演的方法 |
| US20080027650A1 (en) * | 2005-01-14 | 2008-01-31 | Xiaojun Huang | Method and Apparatus for Estimating Formation Slowness |
| CN104656142A (zh) * | 2013-11-19 | 2015-05-27 | 中国石油天然气集团公司 | 一种利用垂直地震剖面与测井联合的地震层位标定方法 |
| CN104656140A (zh) * | 2014-12-01 | 2015-05-27 | 中国石油天然气集团公司 | 一种基于垂直地震衰减规律的中值滤波方法以及系统 |
| CN105182408A (zh) * | 2015-08-28 | 2015-12-23 | 中国石油天然气集团公司 | 一种合成地震记录的制作方法和装置 |
| CN108961352A (zh) * | 2017-05-19 | 2018-12-07 | 中国石油化工股份有限公司 | 一种随钻测井曲线的绘制方法 |
| CN107831541A (zh) * | 2017-11-17 | 2018-03-23 | 中国石油天然气集团公司 | 基于高密度vsp数据的薄地层识别方法及装置 |
| CN111337977A (zh) * | 2018-12-19 | 2020-06-26 | 中国石油天然气股份有限公司 | 品质因子的确定方法及装置 |
Non-Patent Citations (6)
| Title |
|---|
| SEBASTIAN ROST ET AL: "《Improving Seismic Resolution Through Array Processing Techniques》", SURVEYS IN GEOPHYSICS, vol. 30, 12 May 2009 (2009-05-12), pages 271 - 299 * |
| 张雪建等: "《深度域合成地震记录的制作方法研究》", 石油地球物理勘探, vol. 35, no. 3, 30 June 2000 (2000-06-30), pages 377 - 380 * |
| 李桂花等: "《射线法模拟分析井间地震观测的波场特征》", 地球物理学进展, vol. 22, no. 6, 31 December 2007 (2007-12-31), pages 1847 - 1855 * |
| 江成会: "《测井曲线自动分层》", 勘察科学技术, no. 5, 31 December 1993 (1993-12-31), pages 59 - 62 * |
| 蔡志东等: "《基于VSP的地震地质导向方法与应用试验》", 石油物探, vol. 59, no. 3, 31 May 2020 (2020-05-31), pages 430 - 440 * |
| 陆敬安等: "《综合利用测井和地震数据刻划储层参数横向变化》", 长春科技大学学报, vol. 29, no. 2, 30 April 1999 (1999-04-30), pages 171 - 174 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mousavi et al. | Automatic microseismic denoising and onset detection using the synchrosqueezed continuous wavelet transform | |
| CN102176054B (zh) | 近地表综合信息处理解释方法 | |
| US8737165B2 (en) | Interferometric seismic data processing for a towed marine survey | |
| CN110058305A (zh) | 一种基于卷积神经网络的das地震数据降噪方法 | |
| CN102099705A (zh) | 处理地震数据以供解读的系统和方法 | |
| Zhong et al. | RMCHN: A residual modular cascaded heterogeneous network for noise suppression in DAS-VSP records | |
| CN1837859A (zh) | 三维地震资料处理质量监控技术 | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN110389377B (zh) | 基于波形互相关系数相乘的微震偏移成像定位方法 | |
| CN107678064A (zh) | 一种声波时差实时提取方法 | |
| CN103076626B (zh) | 一种波场净化处理方法 | |
| CN101131435A (zh) | 一种避免动校拉伸的动校正叠加方法 | |
| CN104570087A (zh) | 利用瞬时属性进行地震数据初至及同相轴拾取的方法 | |
| CN113640881B (zh) | 多偏移距二维横向高分辨率瞬态面波探测方法 | |
| US4543632A (en) | Robust estimation method for determining when subsequent data processing can include sign-bit representations of full-waveform seismic traces | |
| CN113219539B (zh) | 光纤声波传感地震数据中井筒波干扰的去除方法及装置 | |
| CN112346129A (zh) | 测井曲线提取及合成地震记录制作方法及装置 | |
| Li et al. | Removing abnormal environmental noise in nodal land seismic data using deep learning | |
| CN110780346A (zh) | 一种隧道超前探测复杂地震波场的分离方法 | |
| CN112099086B (zh) | 一种高分辨率光纤井中地震数据深频分析方法 | |
| CN112764109B (zh) | 偶极横波反射波的分离提取方法及装置 | |
| CA2610997C (en) | Method for coherence-filtering of acoustic array signal | |
| CN112526611A (zh) | 表层地震波品质因子的提取方法及装置 | |
| CN104991273A (zh) | 提取叠前校正道集地震子波的方法 | |
| Ma et al. | A Global and Multi-Scale Denoising Method Based on Generative Adversarial Network for DAS VSP Data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |