CN112180443B - 页岩气二维地震甜点区优选方法及装置 - Google Patents
页岩气二维地震甜点区优选方法及装置 Download PDFInfo
- Publication number
- CN112180443B CN112180443B CN201910598604.3A CN201910598604A CN112180443B CN 112180443 B CN112180443 B CN 112180443B CN 201910598604 A CN201910598604 A CN 201910598604A CN 112180443 B CN112180443 B CN 112180443B
- Authority
- CN
- China
- Prior art keywords
- data
- faults
- area
- stratum
- seismic
- 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.)
- Active
Links
- 235000021185 dessert Nutrition 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005457 optimization Methods 0.000 title claims abstract description 31
- 238000010586 diagram Methods 0.000 claims abstract description 32
- 238000009933 burial Methods 0.000 claims abstract description 25
- 230000002349 favourable effect Effects 0.000 claims abstract description 6
- 238000006243 chemical reaction Methods 0.000 claims description 27
- 238000004590 computer program Methods 0.000 claims description 13
- 238000002679 ablation Methods 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 9
- 238000011157 data evaluation Methods 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 6
- 235000009508 confectionery Nutrition 0.000 claims description 4
- 238000005553 drilling Methods 0.000 claims description 4
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 11
- 238000009826 distribution Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 1
- 238000004519 manufacturing process Methods 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
- 239000005416 organic matter Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Mining & Mineral Resources (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种页岩气二维地震甜点区优选方法及装置，该方法包括获取页岩气工区的二维地震资料和地质测井资料；二维地震资料至少包括偏移数据体资料和偏移速度体资料；根据二维地震资料和地质测井资料生成优选参数，优选参数包括：一级剖面资料、地震反射构造图、地震反射埋深图、优质页岩厚度、一类断层信息和二类断层信息；根据优选参数确定甜点区优选结果。本发明实施例根据二维地震资料对甜点区进行优选，得到的甜点区优选结果指导三维地震资料采集范围的部署，便于在更有利的页岩气区块进行三维地震资料采集，从而进行批量页岩气开发，合理利用了现有的二维地震资料、降低了成本，也降低页岩气勘探开发选区工作困难的问题。
Description
技术领域
本发明涉及地球物理技术领域，尤其是涉及一种页岩气二维地震甜点区优选方法及装置。
背景技术
甜点区是指页岩气相对富足的、具有经济开采价值的页岩区域。目前，现有页岩气勘探开发都是直接利用三维地震资料进行选区，而三维地震资料的采集成本较高，尚不可能大面积铺开，因此，大部分页岩气可开发地区，初期无充足的三维地震资料，给甜点区优选工作造成困难。
发明内容
本发明提供了一种页岩气二维地震甜点区优选方法及装置，可以利用二维地震资料进行页岩气甜点区的优选，优选结果可用于指导三维地震资料采集范围的部署，降低了成本，从而降低页岩气勘探开发选区工作的困难。
第一方面，本发明实施例提供了一种页岩气二维地震甜点区优选方法，获取页岩气工区的二维地震资料和地质测井资料；所述二维地震资料至少包括偏移数据体资料和偏移速度体资料；根据所述二维地震资料和所述地质测井资料生成优选参数，优选参数包括：一级剖面资料、地震反射构造图、地震反射埋深图、优质页岩厚度、一类断层信息和二类断层信息；根据所述优选参数确定甜点区优选结果。
第二方面，本发明实施例还提供一种页岩气二维地震甜点区优选装置，获取模块，用于获取页岩气工区的二维地震资料和地质测井资料；所述二维地震资料至少包括偏移数据体资料和偏移速度体资料；参数生成模块，用于根据所述二维地震资料和所述地质测井资料生成优选参数，优选参数包括：一级剖面资料；地震反射构造图；地震反射埋深图；优质页岩厚度；一类断层信息；二类断层信息；确定模块，用于根据所述优选参数确定甜点区优选结果。
第三方面，本发明实施例还提供一种电子设备，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述页岩气二维地震甜点区优选方法。
第四方面，本发明实施例还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质，所述程序代码使所述处理器执行上述页岩气二维地震甜点区优选方法。
本发明实施例带来了以下有益效果：本发明实施例提供了一种页岩气二维地震甜点区优选方法及装置，该方法首先获取页岩气工区的二维地震资料和地质测井资料，根据二维地震资料中的偏移数据体资料、偏移速度体资料以及地质测井资料生成优选参数，优选参数包括一级剖面资料、地震反射构造图、地震反射埋深图、优质页岩厚度、一类断层和二类断层，最后，根据优选参数在页岩气工区中确定甜点区优选结果。本发明实施例根据二维地震资料对甜点区进行优选，得到的甜点区优选结果指导三维地震资料采集范围的部署，便于在更有利的页岩气区块进行三维地震资料采集，从而进行批量页岩气开发，合理利用了现有的二维地震资料、降低了成本，也降低页岩气勘探开发选区工作困难的问题。
本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的页岩气二维地震甜点区优选方法流程示意图；
图2为本发明实施例提供的获取工区的一级剖面的资料和范围的步骤的流程图；
图3为本发明实施例提供的获取时间域的地层反射层位和断层的步骤的流程图；
图4为本发明实施例提供的得到深度域的地震反射层位和断层的步骤的流程图；
图5为本发明实施例提供的得到目的层的埋深图的步骤的流程图；
图6为本发明实施例提供的获得优质页岩厚度平面预测分布预测图的步骤的流程图；
图7为本发明实施例提供的获得优质页岩气二维地震甜点区优选成果图步骤的流程图；
图8为本发明实施例提供的页岩气二维地震甜点区优选方法的流程图；
图9为本发明实施例提供的页岩气二维地震甜点区优选装置结构框图；
图10为本发明实施例提供的电子设备结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前，中国页岩气勘探开发处于起步阶段，在页岩气勘探初期，利用钻井、露头地质调查资料及样品分析结果，从富有机质页岩区域分布、岩性、有机质特征、储层特征、含气性、地层压力等初步优选出甜点区，但要优选出可开发的页岩气区块，必须借助于地震解释成果(目的层埋深、地震资料品质、断层分级评价成果、构造成果图件等)进行优选。大部分页岩气可开发地区无三维地震资料，只有一些二维地震资料，必须利用现有的二维地震成果进行初步的页岩气甜点区的选区，指导三维地震的部署及开发井位的部署。目前，国内外尚未有采用二维地震针对页岩气二维地震甜点区优选的相关专利和文献。
基于此，本发明实施例提供的一种页岩气二维地震甜点区优选方法及装置，可以利用二维地震资料对页岩气甜点区进行优选，优选结果可以用于部署三维地震勘探和页岩气的开发。
为便于对本实施例进行理解，首先对本发明实施例所公开的一种页岩气二维地震甜点区优选方法进行详细介绍。需要说明的是，在执行本发明实施例提供的方法之前，需要人为对页岩气工区选定目的层，目的层是要研究的地层，本发明实施例是基于目的层执行，对于不同的目的层可以分别执行本发明实施例的方法，得到不同的甜点区优选结果。
本发明实施例提供了一种页岩气二维地震甜点区优选方法，参见图1所示的一种页岩气二维地震甜点区优选方法流程图，该方法包括以下步骤：
步骤S102，获取页岩气工区的二维地震资料和地质测井资料。
在本发明实施例中，页岩气工区可以是已经得到了初步优选甜点区的，待进一步优选甜点区的地面区域。对二维地震资料进行预处理可以得到多种类别的资料，包括偏移数据体资料和偏移速度体资料等，其中，偏移数据体资料是二维地震资料中用于从时间域描述地震位置的资料，偏移速度体是二维地震资料中用于描述地震速度的资料。地质测井资料可以包括相关钻井资料、测井资料和地质资料等。
步骤S104，根据二维地震资料和地质测井资料生成优选参数。优选参数包括：一级剖面资料、地震反射构造图、地震反射埋深图、优质页岩厚度、一类断层信息和二类断层信息。
在本发明实施例中，优选参数可以根据二维地震资料和地质测井资料生成。一级剖面资料是指二维地震资料中选出的较优质的且可用性较高的数据，一级剖面资料包括一级剖面资料中与本发明实施例页岩气工区对应的资料和范围。一类断层和二类断层是对二维地震资料中的断层进行评价得到的较优质的且可用性较高的数据。地震反射构造图用于描述页岩气工区对应地下区域的结构情况，地震反射埋深图用于描述页岩气工区对应地下区域结构的埋深情况。
步骤S106，根据优选参数确定甜点区优选结果。
在本发明实施例中，根据优选参数可以在页岩气工区内选定或删除指定区域，经过选定及删除处理后的页岩气工区作为甜点区优选结果。甜点区优选结果是从页岩气工区内最终确定的区域，实现了从初步优选的甜点区中进一步选择更优质甜点区。
本发明实施例提供了一种页岩气二维地震甜点区优选方法，该方法首先获取页岩气工区的二维地震资料和地质测井资料，根据二维地震资料中的偏移数据体资料、偏移速度体资料以及地质测井资料生成优选参数，优选参数包括：一级剖面资料、地震反射构造图、地震反射埋深图、优质页岩厚度、一类断层和二类断层，最后，根据优选参数在页岩气工区中确定甜点区优选结果。本发明实施例根据二维地震资料对甜点区进行优选，得到的甜点区优选结果可以用于部署三维地震资料采集，便于有目的地进行三维地震资料采集，降低成本，从而降低页岩气勘探开发选区工作的困难。
可以根据优选参数从多个维度对甜点区进行筛选，以得到更好的甜点区优选结果，因此，根据二维地震资料和地质测井资料生成优选参数，具体可以包括以下步骤：
(1)根据二维地震资料评价标准对二维地震资料进行评价，得到一级剖面资料。
参见图2所示的取工区的一级剖面的资料和范围的步骤的流程图，在本发明实施例中，二维地震资料评价标准可以根据实际需求进行选择和调整，本发明实施例不作具体限定。对二维地震资料进行评价可以得到多级剖面资料，例如，一级剖面资料和二级剖面资料等，获取其中的一级剖面资料的资料和范围。
(2)对页岩气工区内的二维地层剖面层位和断层进行解释，获取时间域内的地层反射层位和断层。
地震资料解释是指依据时间剖面的波形特征和地质规律赋予地震反射层明确的地质意义。在本发明实施例中，对二维地震资料中的二维地层剖面层位和断层进行解释，得到时间域内的地层反射层位和断层。
(3)建立页岩气工区的时深转换速度体，并根据时深转换速度体对时间域内的地层反射层位和断层进行时深转换，得到深度域的地层反射层位和断层。
在本发明实施例中，时深转换速度体是用于将时间域数据转换成深度域数据的参数。深度域的地层反射层位和断层可以用于描述页岩气工区地面以下区域的地层反射层位和断层的位置深度。
(4)对时间域的断层或深度域的断层进行平面组合，得到断层平面成果，并根据断层平面成果和深度域的地层反射层位生成地震反射构造图和地震反射埋深图。
将多个时间域的断层或多个深度域的断层进行组合连接，可以得到断层平面成果，断层平面成果包括组合连接后的数据，获取断层平面成果中的数据，按照断层断距大小，断层倾向，走向和该区地质背景进行组合，可以生成地震反射构造图，基于地震反射构造图和地面海拔，可以得到目的层的地震反射埋深图。
(5)根据断层断开的地层和断距对深度域内的断层进行分类评价，得到一类断层信息和二类断层信息。
在本发明实施例中，对深度域内的断层进行分类评价是指根据断层断开的地层和断距确定断层的种类信息，种类信息包括一类断层信息和二类断层信息等，其中，一类断层比二类断层的规模更大，断距更大，断开地层更多。
(6)根据地质测井资料和时间域内的地层反射层和断层进行地震反演，得到优质页岩厚度。
把常规的界面型反射剖面转换成岩层型的测井剖面，将地震资料变成可与测井资料直接对比的形式，实现这种转换的处理过程叫地震反演。优质页岩厚度是资源评价的重要因素，通过地震反演以及多种地震信息融合可以得到优质页岩厚度的分布信息。在本发明实施例中，根据地质测井资料和时间域内的地层反射层和断层进行地震反演，可以得到优质页岩厚度的分布情况。
需要说明的是，在实际执行中，上述得到一级剖面资料的资料和范围、得到地震反射构造图和地震反射埋深图以及一类断层信息和二类断层信息以及得到优质页岩厚度的顺序可以调整，本发明实施例不作具体限定。
参见图7所示的获得优质页岩气二维地震甜点区优选成果图步骤的流程图，在得到优选参数后，根据各项优选参数确定甜点区优选结果的具体步骤，包括：
根据一级剖面资料的资料和范围，删除甜点区中非一级资料区域；根据优质页岩厚度，删除优质页岩厚度小于目标厚度的区域；根据地震反射构造图，删除距离剥蚀线目标距离范围内的区域；根据一类断层信息和二类断层信息，分别删除距离第一类断层第一距离范围内的区域和距离第二断层第二距离范围内的区域；第一距离和第二距离的值不相同；根据地层压力系数成果，删除压力系数小于目标系数的区域；将删除甜点区中非一级资料区域、删除页岩厚度小于目标厚度的区域、删除距离剥蚀线目标距离范围内的区域、删除距离第一类断层第一距离范围内的区域、距离第二类断层第二距离范围内的区域和删除压力系数小于目标系数的区域后剩余的区域作为目标甜点区；根据地震反射埋深图，对埋深区域进行分类，得到有利区分类结果；将有利区分类结果和目标甜点区共同作为甜点区优选结果。
在本发明实施例中，根据一级剖面资料的资料和范围可以在页岩气工区初步选定的甜点区中确定对应的地面区域，将非一级剖面资料对应的地面区域删除。删除是指确定为非选择区域。目标厚度是根据经验或者测试设置的厚度，例如，可以为30米，将优质页岩厚度小于30米区域对应的地面区域删除。地震反射构造图中包括剥蚀线，目标距离可以根据经验或者测试设置，例如，可以为8公里，将距离剥蚀线8公里范围内的区域对应的地面区域删除。第一距离和第二距离均可以根据经验或者测试设置，且二者的值不相同，例如，可以确定第一距离为1.5公里，第二距离为0.7公里，则可以删除距离第一类断层1.5公里范围内的区域对应的地面区域和删除距离第二类断层0.7公里范围内的区域对应的地面区域。地层压力系数成果可以使用现有成果数据，目标系数可以根据经验或者测试设置，例如，可以为1.2，将压力系数小于1.2区域对应的地面区域删除。
对页岩气工区初步选定的甜点区完成上述各项删除操作后，得到目标甜点区，最后根据地震反射埋深图信息对目标甜点区的埋深区域进行分类，得到有力区分类结果，其中，埋深图信息包括埋深位置。例如，可以将埋深1500～3500米划分为一类有利区，将埋深3500～4000米划分为二类有利区，将埋深4000～4500米划分为三类有利区。不同类别的有利区表明有利区的可开发程度不同。将目标甜点区及其对应的有利区分类结果共同作为甜点区优选结果。
为了得到更准确的时间域内的信息，参见图3所示的获取时间域的地层反射层位和断层的步骤的流程图，对页岩气工区内的二维地震剖面层位和断层进行解释，获取时间域内的地层反射层位和断层，具体包括以下步骤：
(1)获取地质测井资料中的声波测井曲线和密度曲线，并根据声波测井曲线和密度曲线生成声波合成记录。
在本发明实施例中，声波测井曲线是振幅关于时间的曲线，密度曲线是密度关于时间的曲线，对声波测井曲线和密度曲线进行合成。
需要说明的是，可以根据可通过声波测井曲线和密度曲线联合生成合成记录，也可通过声波测井曲线和密度(恒值)生成合成记录。
(2)根据声波合成记录对偏移数据体资料进行标定，得到页岩气工区中钻井的时深曲线和时间地层反射层位。
将合成记录道与相邻井的井旁地震道的波组关系、波形特征和层间时差相似的特征的道进行对齐，以得到时深曲线。时深曲线是深度关于时间的曲线。将测井上的地质分层标定到偏移体数据中的二维地震剖面上，在地震剖面上确定相应的地层反射特征，可以得到页岩气工区中钻井的时间地层反射层位。
(3)根据时间地层反射层位对二维地震剖面层位和断层进行解释，得到时间域内的地层反射层位和断层。
根据时间地层反射层位，并采用地震解释的相关方法结合页岩气工区的地质情况对二维地震剖面层位和断层进行解释，得到时间域的地层反射层位和断层。
为了提高时间域到深度域的转换精度，参见图4所示的得到深度域的地震反射层位和断层的步骤的流程图，建立页岩气工区的时深转换速度体，并根据时深转换速度体对时间域内的地层反射层位和断层进行时深转换，得到深度域的地层反射层位和断层，包括以下步骤：
(1)对时深曲线进行取点平滑处理，得到包括有均匀的纵向上的速度控制点的时深曲线。
对时深曲线进行取点平滑处理是指间隔选取时深曲线的横纵坐标，得到新的时深曲线。本发明实施例中，可以按照相同的间隔选取时深曲线的横纵坐标，得到的时深曲线中包括均匀的在纵向上的速度控制点。
(2)对时间域内的地层反射层位进行取点平滑处理，得到均匀的横向上的地层反射层位。
对时间域内的地层反射层位进行取点平滑处理是指间隔选取地层反射层位，在本发明实施例中，可以按照相同的间隔选取地层反射层位，得到均匀的横向上的地层反射层位。
(3)根据均匀的横向上的地层反射层位对地层格架进行约束，并根据均匀的纵向上的速度控制点的速度，利用偏移速度体对横向的速度的变化趋势进行约束，得到各地层段的速度，以得到页岩气工区的时深转换速度体。
在本发明实施例中，利用均匀的横向上的地层反射层位对地层格架进行约束，根据均匀的纵向上的速度控制点的速度，利用偏移速度体对横向的速度的变化趋势进行约束，纵向上则是通过平滑取点来确定精度，得到各地层段的速度，从而得到页岩气工区的时深转换速度体。
(4)根据时深转换速度体对时间域内的地层反射层位和断层进行时深转换，得到深度域的地层反射层位和断层。
在本发明实施例中，将时间域内的地层反射层位和断层作为时深转换速度体的输入，经转换可以得到深度域的地层反射层位和断层。时深转换速度体用于将时间域的数据转换为深度域的数据，以用于描述地层反射层位和断层的位置深度。
参见图5所示的得到目的层的埋深图的步骤的流程图，根据断层平面成果和深度域的地层反射层位生成地震反射构造图和地震反射埋深图，包括：
根据时间域内的地层反射层位确定剥蚀线，并根据深度域的层位进行网格化加上剥蚀线和断层平面成果生成地震反射构造图；根据地面海拔加上地震反射构造图生成地震反射埋深图。
在本发明实施例中，根据解释的层位，画出相应地层的剥蚀线，加上组合好的断层平面成果，将等值线进行网格化，得到地震反射构造图。利用目的层的构造图等值线值加上地面相应的地面海拔，再将剥蚀线和组合好的断层加上，得到地震反射埋深图。
参见图6所示的获得优质页岩厚度平面预测分布预测图的步骤的流程图，根据地质测井资料和时间域内的地层反射层和断层进行地震反演，得到优质页岩厚度，包括：
根据地质测井资料确定时间域内的地层反射层和断层中优质页岩段的响应模式；根据响应模式对时间域内的地层反射层和断层进行地震反演，得到优质页岩厚度。
在本发明实施例中，根据工区测井解释的优质页岩段，利用测井响应特征和地震响应特征得到优质页岩的响应模式。可以采用如密度曲线、声波曲线等进行拟合和分析，得到相应区块的与地震相关度最好的曲线，得到该曲线对应的各类参数即得到响应模式。根据优质页岩的响应模式，对时间域内的地层反射层和断层利用速度反演和伽马反演等，可以得到优质页岩厚度。
本发明实施例提供了一种页岩气二维地震甜点区优选方法及装置，参见图8所示的页岩气二维地震甜点区优选方法的流程图，本发明实施例可以利用二维地震资料，根据地震资料品质评价结果、利用二维地震资料评价标准对地震资料进行评价，获得一级剖面的资料和范围；然后对工区内的二维地震剖面层位和断层进行解释，获取时间域的地层反射层位和断层；再建立工区的时深转换速度场，对工区内解释的时间地震层位和断层进行时深转换，获得深度域的地震反射层位和断层；同时根据二维地震剖面上解释的断层进行平面组合，获得断层平面成果；在进行构造图的编制，加入地面海拔进行埋深图的编制，获得地震反射构造图和埋深图；并对解释的断层根据断层断开的地层和断距进行断层分类分级评价，获得I类断层和II类断层；同时利用地震反演方法对进行地震地震反演，获得优质页岩厚度；最后根据地震资料评价结果、预测的优质页岩埋深、优质页岩厚度、构造成果中的剥蚀线、断层的分类分级评价标准、以及利用的地层压力系数进行页岩气二维地震甜点区优选，获得页岩气二维地震甜点区优选结果，优选结果可以用于三维地震勘探和页岩气的开发。
本发明实施例还提供一种页岩气二维地震甜点区优选装置，参见图9所示的页岩气二维地震甜点区优选装置结构框图，该装置包括：获取模块81，用于获取页岩气工区的二维地震资料和地质测井资料；二维地震资料至少包括偏移数据体资料和偏移速度体资料；参数生成模块82，用于根据二维地震资料和地质测井资料生成优选参数，优选参数包括：一级剖面资料；地震反射构造图；地震反射埋深图；优质页岩厚度；一类断层信息；二类断层信息；确定模块83，用于根据优选参数确定甜点区优选结果。
本发明实施例还提供一种电子设备，参见图10所示的电子设备结构示意框图，该电子设备包括存储器91、处理器92，存储器中存储有可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述任一种方法的步骤。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的电子设备的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
本发明实施例还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质，程序代码使处理器执行上述任一种方法的步骤。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (9)
1.一种页岩气二维地震甜点区优选方法，其特征在于，包括：
获取页岩气工区的二维地震资料和地质测井资料；所述二维地震资料至少包括偏移数据体资料和偏移速度体资料；
根据所述二维地震资料和所述地质测井资料生成优选参数，所述优选参数包括：一级剖面资料；地震反射构造图；地震反射埋深图；优质页岩厚度；一类断层信息；二类断层信息；一类断层信息和二类断层信息是对二维地震资料中的断层进行评价得到的数据；一类断层的断距大于二类断层的断距；一级剖面资料根据二维地震资料评价标准对二维地震资料进行评价得到；
根据所述优选参数确定甜点区优选结果；
根据所述优选参数确定甜点区优选结果，包括：
根据所述一级剖面资料，删除所述甜点区中非一级剖面资料区域；
根据所述优质页岩厚度，删除所述优质页岩厚度小于目标厚度的区域；
根据所述地震反射构造图，删除距离剥蚀线目标距离范围内的区域；
根据所述一类断层信息和二类断层信息，分别删除距离一类断层第一距离范围内的区域和距离二类断层第二距离范围内的区域；所述第一距离和所述第二距离的值不相同；
根据地层压力系数成果，删除压力系数小于目标系数的区域；
将删除所述甜点区中非一级剖面资料区域、删除所述页岩厚度小于目标厚度的区域、删除距离剥蚀线目标距离范围内的区域、删除距离一类断层第一距离范围内的区域、距离二类断层第二距离范围内的区域和删除压力系数小于目标系数的区域后剩余的区域作为目标甜点区；
根据所述地震反射埋深图，对埋深区域进行分类，得到有利区分类结果；
将所述有利区分类结果和所述目标甜点区共同作为甜点区优选结果。
2.根据权利要求1所述的方法，其特征在于，根据所述二维地震资料和所述地质测井资料生成优选参数，包括：
根据二维地震资料评价标准对所述二维地震资料进行评价，得到一级剖面资料；
对所述页岩气工区内的二维地层剖面层位和断层进行解释，获取时间域内的地层反射层位和断层；
建立所述页岩气工区的时深转换速度体，并根据所述时深转换速度体对所述时间域内的地层反射层位和断层进行时深转换，得到深度域的地层反射层位和断层；
对所述时间域的断层或所述深度域的断层进行平面组合，得到断层平面成果，并根据所述断层平面成果和深度域的地层反射层位生成地震反射构造图和地震反射埋深图；
根据断层断开的地层和断距对深度域内的断层进行分类评价，得到一类断层信息和二类断层信息；
根据所述地质测井资料和所述时间域内的地层反射层和断层进行地震反演，得到优质页岩厚度。
3.根据权利要求2所述的方法，其特征在于，对所述页岩气工区内的二维地震剖面层位和断层进行解释，获取时间域内的地层反射层位和断层，包括：
获取所述地质测井资料中的声波测井曲线和密度曲线，并根据所述声波测井曲线和所述密度曲线生成声波合成记录；
根据所述声波合成记录对所述偏移数据体资料进行标定，得到所述页岩气工区中钻井的时深曲线和时间地层反射层位；
根据所述时间地层反射层位对所述二维地震剖面层位和断层进行解释，得到时间域内的地层反射层位和断层。
4.根据权利要求3所述的方法，其特征在于，建立所述页岩气工区的时深转换速度体，并根据所述时深转换速度体对所述时间域内的地层反射层位和断层进行时深转换，得到深度域的地层反射层位和断层，包括：
对所述时深曲线进行取点平滑处理，得到包括有均匀的纵向上的速度控制点的时深曲线；
对所述时间域内的地层反射层位进行取点平滑处理，得到均匀的横向上的地层反射层位；
根据所述均匀的横向上的地层反射层位对地层格架进行约束，并根据均匀的纵向上的速度控制点的速度，利用偏移速度体对横向的速度的变化趋势进行约束，得到各地层段的速度，以得到所述页岩气工区的时深转换速度体；
根据所述时深转换速度体对所述时间域内的地层反射层位和断层进行时深转换，得到深度域的地层反射层位和断层。
5.根据权利要求2所述的方法，其特征在于，根据所述断层平面成果和深度域的地层反射层位生成地震反射构造图和地震反射埋深图，包括：
根据所述时间域内的地层反射层位和断层确定剥蚀线，并根据所述剥蚀线和所述断层平面成果生成地震反射构造图；
根据所述地震反射构造图和所述深度域的地层反射层位生成地震反射埋深图。
6.根据权利要求2所述的方法，其特征在于，根据所述地质测井资料和所述时间域内的地层反射层和断层进行地震反演，得到优质页岩厚度，包括：
根据所述地质测井资料确定所述时间域内的地层反射层和断层中优质页岩段的响应模式；
根据所述响应模式对所述时间域内的地层反射层和断层进行地震反演，得到优质页岩厚度。
7.一种页岩气二维地震甜点区优选装置，其特征在于，包括：
获取模块，用于获取页岩气工区的二维地震资料和地质测井资料；所述二维地震资料至少包括偏移数据体资料和偏移速度体资料；
参数生成模块，用于根据所述二维地震资料和所述地质测井资料生成优选参数，所述优选参数包括：一级剖面资料；地震反射构造图；地震反射埋深图；优质页岩厚度；一类断层信息；二类断层信息；一类断层信息和二类断层信息是对二维地震资料中的断层进行评价得到的数据；一类断层的断距大于二类断层的断距；一级剖面资料根据二维地震资料评价标准对二维地震资料进行评价得到；
确定模块，用于根据所述优选参数确定甜点区优选结果；
确定模块，具体用于：
根据所述一级剖面资料，删除所述甜点区中非一级剖面资料区域；
根据所述优质页岩厚度，删除所述优质页岩厚度小于目标厚度的区域；
根据所述地震反射构造图，删除距离剥蚀线目标距离范围内的区域；
根据所述一类断层信息和二类断层信息，分别删除距离一类断层第一距离范围内的区域和距离二类断层第二距离范围内的区域；所述第一距离和所述第二距离的值不相同；
根据地层压力系数成果，删除压力系数小于目标系数的区域；
将删除所述甜点区中非一级剖面资料区域、删除所述页岩厚度小于目标厚度的区域、删除距离剥蚀线目标距离范围内的区域、删除距离一类断层第一距离范围内的区域、距离二类断层第二距离范围内的区域和删除压力系数小于目标系数的区域后剩余的区域作为目标甜点区；
根据所述地震反射埋深图，对埋深区域进行分类，得到有利区分类结果；
将所述有利区分类结果和所述目标甜点区共同作为甜点区优选结果。
8.一种电子设备，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至6任一项所述的方法的步骤。
9.一种具有处理器可执行的非易失的程序代码的计算机可读介质，其特征在于，所述程序代码使所述处理器执行上述权利要求1至6任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910598604.3A CN112180443B (zh) | 2019-07-04 | 2019-07-04 | 页岩气二维地震甜点区优选方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910598604.3A CN112180443B (zh) | 2019-07-04 | 2019-07-04 | 页岩气二维地震甜点区优选方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112180443A CN112180443A (zh) | 2021-01-05 |
| CN112180443B true CN112180443B (zh) | 2024-03-01 |
Family
ID=73914535
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910598604.3A Active CN112180443B (zh) | 2019-07-04 | 2019-07-04 | 页岩气二维地震甜点区优选方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112180443B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114060022B (zh) * | 2021-11-25 | 2022-08-05 | 西南石油大学 | 一种基于多尺度断层发育特征的页岩气产能预测方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104853822A (zh) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN104977618A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN105700013A (zh) * | 2016-04-26 | 2016-06-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 构建页岩气甜点产能模型的方法 |
| CN105986815A (zh) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | 一种用于识别页岩地层地质甜点的方法 |
| CN106199727A (zh) * | 2016-08-30 | 2016-12-07 | 中国华能集团清洁能源技术研究院有限公司 | 一种应用地震层速度识别页岩气甜点的方法 |
| CN106443770A (zh) * | 2015-08-07 | 2017-02-22 | 中国石油化工股份有限公司 | 一种页岩气地质甜点的预测方法 |
| CN107102354A (zh) * | 2016-12-21 | 2017-08-29 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | 一种页岩甜点地震综合评价方法 |
| CN107545512A (zh) * | 2016-06-29 | 2018-01-05 | 中国石油化工股份有限公司 | 基于动态富集的页岩油甜点综合评价方法 |
| CN108363100A (zh) * | 2018-01-16 | 2018-08-03 | 太原理工大学 | 基于封存单元与粗糙集的煤层气甜点区地震地质识别方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8626447B2 (en) * | 2010-09-13 | 2014-01-07 | Chevron U.S.A. Inc. | System and method for sweet zone identification in shale gas reservoirs |
| CN111441758B (zh) * | 2018-12-29 | 2021-03-30 | 中国石油天然气股份有限公司 | 页岩油气甜点区的预测方法及装置 |
-
2019
- 2019-07-04 CN CN201910598604.3A patent/CN112180443B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104977618A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN104853822A (zh) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN105986815A (zh) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | 一种用于识别页岩地层地质甜点的方法 |
| CN106443770A (zh) * | 2015-08-07 | 2017-02-22 | 中国石油化工股份有限公司 | 一种页岩气地质甜点的预测方法 |
| CN105700013A (zh) * | 2016-04-26 | 2016-06-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 构建页岩气甜点产能模型的方法 |
| CN107545512A (zh) * | 2016-06-29 | 2018-01-05 | 中国石油化工股份有限公司 | 基于动态富集的页岩油甜点综合评价方法 |
| CN106199727A (zh) * | 2016-08-30 | 2016-12-07 | 中国华能集团清洁能源技术研究院有限公司 | 一种应用地震层速度识别页岩气甜点的方法 |
| CN107102354A (zh) * | 2016-12-21 | 2017-08-29 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | 一种页岩甜点地震综合评价方法 |
| CN108363100A (zh) * | 2018-01-16 | 2018-08-03 | 太原理工大学 | 基于封存单元与粗糙集的煤层气甜点区地震地质识别方法 |
Non-Patent Citations (4)
| Title |
|---|
| Sweet Spots Discrimination in Shale Gas Reservoirs Using Seismic and Well-logs Data.A Case Study from the Worth Basin in the Barnett Shale;Leila Aliouane 等;Energy Procedia;第59卷;第22-27页 * |
| 利用地球物理综合预测方法识别页岩气储层甜点――以四川盆地长宁区块下志留统龙马溪组为例;陈胜 等;天然气工业;第37卷(第05期);第20-30页 * |
| 四川盆地W地区龙马溪组页岩气压裂效果评估和产能预测研究;李大军 等;石油物探;第56卷(第05期);第735-745页 * |
| 四川盆地昭通区块龙马溪组页岩气"甜点"预测方法及应用;刘伟 等;石油地球物理勘探;第53卷(第S2期);第211-217、223页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112180443A (zh) | 2021-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9097821B2 (en) | Integrated workflow or method for petrophysical rock typing in carbonates | |
| JP3390445B2 (ja) | 地質学的特性を推定するために統計学的較正技術を使用する地震トレース解析方法 | |
| AU2010332157B2 (en) | Workflow for petrophysical and geophysical formation evaluation of wireline and LWD log data | |
| CA2931435A1 (en) | Method for developing a geomechanical model based on seismic data, well logs and sem analysis of horizontal and vertical drill cuttings | |
| Sams et al. | Stochastic inversion for high resolution reservoir characterisation in the Central Sumatra Basin | |
| KR102111207B1 (ko) | 셰일가스 스윗 스팟 도출 방법 | |
| Allo et al. | Characterization of a carbonate geothermal reservoir using rock-physics-guided deep neural networks | |
| Faraji et al. | Integration of core data, well logs and seismic attributes for identification of the low reservoir quality units with unswept gas in the carbonate rocks of the world’s largest gas field | |
| Arata et al. | Look ahead geosteering via real time integration of logging while drilling measurements with surface seismic | |
| CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
| CN112180443B (zh) | 页岩气二维地震甜点区优选方法及装置 | |
| CN112528106A (zh) | 火山岩岩性识别方法 | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| Price et al. | Dynamic calibration of the Shaikan Jurassic full-field fractured reservoir model through single-well DST and multi-well interference discrete fracture network simulation | |
| CN106291748A (zh) | 一种获取烃源岩有机碳含量的方法 | |
| Lanzarone et al. | The value of VSP data through early phases of field appraisal and development: A modeling and acquisition case study in the Gulf of Mexico | |
| Lavenu et al. | Modeling Fractures in a Heterogeneous Carbonate Reservoir Onshore UAE: A Case Study | |
| CN104345337B (zh) | 一种用于地震反演的时控储层参数建模方法 | |
| Gauthier et al. | Static and dynamic characterization of fracture pattern in the Upper Jurassic reservoirs of an offshore Abu Dhabi field: from well data to full field modeling | |
| Fretwell et al. | A new approach to 3-D geological modeling of complex sand injectite reservoirs: The Alba field, United Kingdom central North Sea | |
| Webb et al. | A time-lapse case study in West Africa: Integrating disciplines for a complete reservoir study and field management | |
| Dickmann et al. | Correlating rock support and ground treatment means with in-tunnel seismic data | |
| Berkhout | Seismic resolution: A key to detailed geologic information | |
| CN111352154B (zh) | 一种基于宽方位地震的储层预测方法 | |
| Al Badi et al. | Using Zero-Offset VSP to Identify Interbed Multiples Generators, A Case Study from South Oman |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |