CN113552621A - 页岩气地应力确定方法和装置 - Google Patents
页岩气地应力确定方法和装置 Download PDFInfo
- Publication number
- CN113552621A CN113552621A CN202010325413.2A CN202010325413A CN113552621A CN 113552621 A CN113552621 A CN 113552621A CN 202010325413 A CN202010325413 A CN 202010325413A CN 113552621 A CN113552621 A CN 113552621A
- Authority
- CN
- China
- Prior art keywords
- depth
- determining
- stress
- domain
- thickness
- 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
- 238000000034 method Methods 0.000 title claims abstract description 99
- 238000010276 construction Methods 0.000 claims abstract description 35
- 230000008859 change Effects 0.000 claims abstract description 13
- 239000011148 porous material Substances 0.000 claims description 70
- 239000011435 rock Substances 0.000 claims description 56
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 238000013508 migration Methods 0.000 claims description 11
- 230000005012 migration Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 9
- 238000011160 research Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 abstract description 13
- 238000011161 development Methods 0.000 abstract description 9
- 238000005553 drilling Methods 0.000 abstract description 8
- 238000010586 diagram Methods 0.000 description 18
- 238000012986 modification Methods 0.000 description 9
- 230000004048 modification Effects 0.000 description 9
- 230000014509 gene expression Effects 0.000 description 6
- 230000015556 catabolic process Effects 0.000 description 5
- 238000006731 degradation reaction Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000007423 decrease Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000006467 substitution reaction 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/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明供了一种页岩气地应力确定方法和装置，该方法利用了地震数据能够反映地层横向的连续性变化以及地震数据能够获得连续的地应力剖面的特性，提出了一种利用地震数据体进行页岩气地应力的确定方法，能够有效反映全工区内地应力大小的展布情况，特别是横向地应力的展布情况。不仅适用于多井区，也适用于少井或无井区；同时由于地震数据预测具有前瞻性，可进行钻前的地应力预测，对于盆地页岩气勘探开发区构造和成藏都较复杂的情况，建立一套基于地震数据的页岩气地应力预测方法，实现了页岩气地应力的准确预测，为页岩气井位部署和压裂施工进行有效支撑。
Description
技术领域
本发明涉及地球物理信号解释技术领域，尤其涉及一种页岩气地应力确定方法和装置。
背景技术
页岩气储层属于低孔低渗类储层，目前的开发仍依赖于大规模的水力压裂改造才能够获得商业产能。而地应力则是影响压裂改造的关键因素，它将直接影响压裂方案的设计和最终效果。同时地应力不仅是气藏的开发动力之一，且能为井壁稳定性分析和钻井优化设计等提供依据。
目前关于地应力的研究主要集中在地应力测量方法(如水力压裂法和声发射法)和地应力测井计算方法(成像测井、倾角测井等)两个领域。这两种方面更多地关注某一点、或某一井筒处的地应力分布情况。而对于井间的地应力分布，或少井或无井区的地应力分布情况研究关注较少。同时这些方法具有事后性，在钻前应用受到了一定的限制。
盆地页岩气勘探开发区构造和成藏都较复杂，针对页岩气勘探开发，现有技术中没有针对页岩气地应力预测的详细工程解决方案。
发明内容
本发明实施例提供一种页岩气地应力确定方法，实现了页岩气地应力的准确预测，该方法包括：
根据地震数据体，确定深度域解释成果；所述深度域解释成果，包括深度域纵波速度体、深度域横波速度体、深度域密度体、时间厚度数据和深度厚度数据；
根据深度域解释成果，建立地质构造模型；
基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体；
根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体；
根据深度域解释成果和地质构造模型，确定岩石密度体；
根据测井曲线和地质构造模型，确定岩性指示体；
根据岩石密度体，确定垂直主应力体；
根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体；
根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体；
根据垂直主应力体、孔隙压力体、最大水平主应力体和最小水平主应力体，确定页岩气地应力。
本发明实施例还提供一种页岩气地应力确定装置，包括：
深度域解释成果确定模块，用于根据地震数据体，确定深度域解释成果；所述深度域解释成果，包括深度域纵波速度体、深度域横波速度体、深度域密度体、时间厚度数据和深度厚度数据；
地质构造模型建立模块，用于根据深度域解释成果，建立地质构造模型；
时间厚度体和深度厚度体确定模块，用于基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体；
弹性参数体确定模块，用于根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体；
岩石密度体确定模块，用于根据深度域解释成果和地质构造模型，确定岩石密度体；
岩性指示体确定模块，用于根据测井曲线和地质构造模型，确定岩性指示体；
垂直主应力体确定模块，用于根据岩石密度体，确定垂直主应力体；
孔隙压力体确定模块，用于根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体；
最大水平主应力体和最小水平主应力体确定模块，用于根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体；
页岩气地应力确定模块，用于根据垂直主应力体、孔隙压力体、最大水平主应力体和最小水平主应力体，确定页岩气地应力。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述页岩气地应力确定方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现上述页岩气地应力确定方法的计算机程序。
本发明实施例提供的一种岩气地应力确定方法和装置，利用了地震数据能够反映地层横向的连续性变化以及地震数据能够获得连续的地应力剖面的特性，提出了一种利用地震数据体进行页岩气地应力的确定方法，能够有效反映全工区内地应力大小的展布情况，特别是横向地应力的展布情况。不仅适用于多井区，也适用于少井或无井区；同时由于地震数据预测具有前瞻性，可进行钻前的地应力预测，对于盆地页岩气勘探开发区构造和成藏都较复杂的情况，建立一套基于地震数据的页岩气地应力预测方法，实现了页岩气地应力的准确预测，为页岩气井位部署和压裂施工进行有效支撑。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例一种岩气地应力确定方法示意图。
图2为本发明实施例一种岩气地应力确定方法流程图。
图3为本发明实施例一种岩气地应力确定方法的确定深度域解释成果示意图。
图4为运行本发明实施的一种岩气地应力确定方法的计算机装置示意图。
图5为本发明实施例一种岩气地应力确定装置示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
如图1本发明实施例一种岩气地应力确定方法示意图所示，本发明实施例提供一种页岩气地应力确定方法，实现了页岩气地应力的准确预测，该方法包括：
步骤101：根据地震数据体，确定深度域解释成果；所述深度域解释成果，包括深度域纵波速度体、深度域横波速度体、深度域密度体、时间厚度数据和深度厚度数据；
步骤102：根据深度域解释成果，建立地质构造模型；
步骤103：基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体；
步骤104：根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体；
步骤105：根据深度域解释成果和地质构造模型，确定岩石密度体；
步骤106：根据测井曲线和地质构造模型，确定岩性指示体；
步骤107：根据岩石密度体，确定垂直主应力体；
步骤108：根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体；
步骤109：根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体；
步骤110：根据垂直主应力体、孔隙压力体、最大水平主应力体和最小水平主应力体，确定页岩气地应力。
本发明实施例提供的一种岩气地应力确定方法，利用了地震数据能够反映地层横向的连续性变化以及地震数据能够获得连续的地应力剖面的特性，提出了一种利用地震数据体进行页岩气地应力的确定方法，能够有效反映全工区内地应力大小的展布情况，特别是横向地应力的展布情况。不仅适用于多井区，也适用于少井或无井区；同时由于地震数据预测具有前瞻性，可进行钻前的地应力预测，对于盆地页岩气勘探开发区构造和成藏都较复杂的情况，建立一套基于地震数据的页岩气地应力预测方法，实现了页岩气地应力的准确预测，为页岩气井位部署和压裂施工进行有效支撑。
地震数据能够反映地层横向的连续性变化，因此基于地震数据也能获得连续的地应力剖面。利用地震数据，基于岩石力学方法进行地应力的研究和预测。该方法经济易操作，能够有效反映全工区内地应力大小的展布情况，特别是横向地应力的展布情况。不仅适用于多井区，也适用于少井或无井区。同时地震资料预测具有前瞻性，可进行钻前的地应力预测。因此有必要针对威远地区，建立一套基于地震数据的页岩气地应力预测方法。
在进行地应力预测的过程中，需要分别进行垂直主应力的计算、孔隙压力的计算、以及最大最小水平主应力的计算。三个参数的基本计算方式和所需参数如下表1所示：
表1
另外，目前常用的孔隙压力预测方法，多基于页岩气孔隙度随垂直有效应力(即垂直主应力减去孔隙压力)增加呈单调递减的事实。任意深度上的垂直主应力求取较为容易，而孔隙度则可以根据声波速度或电阻率数据估算。但页岩孔隙度变化可能表明有效应力降低，但有效应力的降低也可能是岩性变化导致的。同时，如果发生应力卸载相关的地质过程，如快速抬升和剥蚀，则孔隙度与有效应力之间的关系将不在简单的单调关系。在一个方案中指出，随应力卸载有效应力减小，但沿着应力卸载路径孔隙度几乎不变化。因此，常用的孔隙压力预测方法无法产生理想的预测结果。因此为了更好的获得水平应力场预测结果，需要更高精度的孔隙压力场计算成果。
水平应力场的预测，融合了多种输入参数，包括垂直主应力、孔隙压力、以及力学参数等，因此每一个输入参数都会影响水平应力场的最终预测效果，而其中力学参数的平面构建是一个关键的难点。目前，针对威远地区水平应力场计算中力学参数平面构建的研究较少，本发明实施例就针对该难点进行了研究，并提出了基于地质背景的构造系数构建方法，提高了水平应力场横向预测效果的准确性。
如图1本发明实施例一种岩气地应力确定方法示意图和图2本发明实施例一种岩气地应力确定方法流程图所示，本发明实施例提供的一种页岩气地应力确定方法，可以包括：
根据地震数据体，确定深度域解释成果；所述深度域解释成果，包括深度域纵波速度体、深度域横波速度体、深度域密度体、时间厚度数据和深度厚度数据；根据深度域解释成果，建立地质构造模型；基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体；根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体；根据深度域解释成果和地质构造模型，确定岩石密度体；根据测井曲线和地质构造模型，确定岩性指示体；根据岩石密度体，确定垂直主应力体；根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体；根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体；根据垂直主应力体、孔隙压力体、最大水平主应力体和最小水平主应力体，确定页岩气地应力。
图3为本发明实施例一种岩气地应力确定方法的确定深度域解释成果示意图。如图3所示，本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，根据地震数据体，确定深度域解释成果，包括：
步骤301：根据地震数据体，获得时间域精细层位解释和断层解释，进行断层组合，为各个断层进行分配命名，获得目的层的时间厚度数据；
步骤302：基于时间域精细层位解释、断层解释和道集数据，进行叠前同时反演计算，获得时间域解释成果；所述时间域解释成果，包括：时间域纵波速度体、时间域横波速度体和时间域密度体；
步骤303：构建速度模型，对时间域解释成果进行时深转化，确定深度域的精细层位解释和断层解释；
步骤304：根据深度域的精细层位解释和断层解释，进行断层组合，为各个断层进行分配命名，获得目的层的深度厚度数据；
步骤305：利用深度偏移速度场，将深度域数据转化为时间域数据，进行叠前同时反演，获得时间域反演解释成果；所述深度域数据，包括：深度域的精细层位解释和断层解释，目的层的深度厚度数据；所述时间域反演解释成果，包括时间域反演纵波速度体、时间域反演横波速度体和时间反演域密度体；
步骤306：利用深度偏移速度场，将时间域反演解释成果进行转化，获得深度域纵波速度体、深度域横波速度体和深度域密度体。
在实施例中，综合考虑地区的地质背景认识以及井数据，进行速度场的构建，获得工区的时间域和深度域精细层位解释和断层解释，并获得目的层的时间厚度数据和深度域厚度数据。基于解释成果和道集数据，进行叠前同时反演计算，并利用速度场，获得深度域的纵波速度体、横波速度体以及密度体。
基于时深转化的数据体获得：采用叠前时间偏移+时深转化技术，则首先要获得时间域精细层位解释和断层解释，完成断层组合，并为各个断层进行分配命名，并计算目的层的时间厚度数据。基于时间域精细层位解释、断层解释和道集数据，进行叠前同时反演计算，获得时间域解释成果，前述的时间域解释成果包括：时间域的纵波速度体、横波速度体以及密度体。叠前反演结果与参与反演井点吻合率要达到80％以上。
构建合理的速度模型：在页岩气中，一般采用层速度模型，对时间域解释成果进行时深转化，上述层速度模型在一个实例中可以应用于四川页岩气中。通过时深转化，可以获得深度域的精细层位解释和断层解释、目的层厚度以及深度域的叠前同时反演结果。时深转化后，保证井的深度误差小于10m。在进行建模前进行质控，保证深度层位、反演数据体、井位数据三者的一致性。
基于深度偏移技术的数据体获得：采用深度偏移技术，则通过处理和解释的合作，可获得最终的深度域层位解释和断层解释，完成断层组合，并为各个断层进行分配命名，并计算目的层的深度厚度数据；
利用深度偏移速度场，将深度域数据转化为时间域数据，并进行叠前同时反演，获得时间域反演解释成果；所述深度域数据，包括：深度域的精细层位解释和断层解释，目的层的深度厚度数据；所述时间域反演解释成果，包括时间域反演纵波速度体、时间域反演横波速度体和时间域反演密度体，利用深度偏移速度场，将时间域反演解释成果进行转化，可以获得深度域纵波速度体、深度域横波速度体和深度域密度体；
同样，叠前反演结果与参与反演井点吻合率要达到80％以上，井深度误差要小于10m。同时保证深度层位、反演数据体、井位数据三者的一致性。
在实施例中，根据深度域解释成果，建立地质构造模型，包括：
基于深度域解释成果，建立三维地质网格体。地质构造模型是所有后期属性体的载体，因此在建立时，既要保证合理的简化，也要保证精度需要。
在建地质构造模型的网格时，应综合考虑构造形态、反演成果、工程精度以及模型大小，合理简化模型。在划分地质构造模型的网格时，要注意工程精度和模型计算量之间的平衡。要根据需要调整纵向网格大小，非目的层网格可较粗，目的层段需要较细，同时还要避免相邻网格大小的剧烈突变。
进一步的，基于深度域解释成果，建立三维地质地质构造模型。地质构造模型的建立，分为层位断层预处理、断层模型建立、层序模型建立、骨架模型建立以及网格模型建立五个主要步骤。一般情况下，地震解释数据为散点数据，需要对层位和断层数据进行插值成面，同时，层位一定要有地表层位和目的层下方的层位。在建立地质构造模型前，需要将层位和断层面中的数据奇点进行修正，并进行合理的平滑，保证后期建模的正确性。在断层模型建立时，仅需选择关键断层进行建立即可，同时针对逆断层，需重点处理断层与层位的相切线；层序模型的建立需要根据后期数据的显示和计算进行建立，对于威远地区页岩气，一般将下二叠统底至上奥五峰组底划分为一个层段，便于后期研究。同时，在建立层序模型约束时，一般利用层位和井分层数据进行同时控制，保证井点处层位与分层的完全吻合，确保后期插值的正确性。在骨架模型建立后，要进行模型质控，检查层位和断层的接触关系是否正确，是否出现网格畸点等。在建立网格模型时，应综合考虑构造形态、反演成果、工程精度以及模型大小，合理简化模型。在划分网格时，要注意工程精度和模型计算量之间的平衡。要根据需要调整纵向网格大小，非目的层网格可较粗，目的层段需要较细，同时还要避免相邻网格大小的剧烈突变。针对威远地区，上覆非目的层纵向网格大小从上至下可从20m渐变至5m，而目的层网格纵向尺寸一般为2m左右。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体，包括：
构建覆盖所有目的层段的时间厚度属性segy数据体和深度厚度属性segy数据体；
将时间厚度属性segy数据体和深度厚度属性segy数据体导入地质构造模型中，对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体。
实施例中，上奥五峰组底向下漂移20m生成新层位，利用该层位和地表层位进行时间厚度和深度厚度的segy数据体模型构建。具体构建方法与速度场的构建方法一致，生成一个由线号、道号、X坐标、Y坐标、深度、厚度，组成的文本，并导入landmark速度建模模块中，即可生成segy体。
将所生成的segy体导入建模软件中，并基于前述建立的地质构造模型进行数据采集，即可获得3D网格的时间厚度体和深度厚度体。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体，包括：
根据深度域纵波速度体和深度域横波速度体，基于地质构造模型对地震数据体进行数据采集，确定纵波速度模型和横波速度模型；
根据纵波速度模型和横波速度模型，确定弹性参数体。
在实施例中，如果没有横波速度反演体，可通过工区内井的横波测井数据，根据构造模型进行插值得到。然后再进行弹性参数和力学参数模型的建立。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，所述弹性参数体，包括：泊松比模型和α系数模型；
按如下方式，确定弹性参数体：
其中，Vp为纵波速度；Vs为横波速度；v为泊松比模型；α为α系数模型。
前述提到的设定扫描范围和扫描步长的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，根据深度域解释成果和地质构造模型，确定岩石密度体，包括：
在深度域密度体完整时，将深度域密度体采集到地质构造模型中，获得岩石密度体；
在深度域密度体不完整时，根据密度变化趋势，补齐密度曲线，利用变差函数分析的方式进行插值，建立目的层上方密度体；将目的层上方密度体和原有不完整深度域密度体合并，获得岩石密度体。
实施例中，对于反演所需的密度体，理想的反演密度数据应该是完整的，即包含从地表到目的层的所有反演结果。
如果所得的深度域密度体是完整的，那么可直接将深度域密度体的地震反演数据采集到前述所建立的地质构造模型中，即可获得基于地震的具有横向连续变化的岩石密度体。
如果所得的深度域密度体不是完整的，即仅包含目的层的密度数据，则需要如下三个步骤：(1)根据密度变化趋势，补齐测井密度曲线，使密度曲线包含从地表到目的层上方的所有数据；(2)利用变差函数分析或其他合理的插值方式建立目的层上方的密度体；(3)合并上方密度体和原有不完整深度域密度体，获得完整的岩石密度体。
进一步的，如果没有深度域密度体，则利用密度测井曲线，再补齐曲线后采用合理的插值方式，构建从地表到目的层的岩石密度体。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，根据测井曲线和地质构造模型，确定岩性指示体，包括：
对测井曲线进行归一化处理，利用变差函数的插值方法，基于构造地质模型，建立岩性指示体。
实施例中，在威远地区，自然伽马曲线能够保证页岩气储层段龙马溪的岩性变化，因此利用自然伽马曲线进行归一化处理后，并利用变差函数或合理的插值方法，基于前述建立的地质构造模型，建立目的层的岩性指示体。
由于伽马与弹性参数之间并不存在物理意义上的关系，因此不建议利用伽马反演的方法来获得自然伽马数据体，再采集到网格中的方法来构建。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，根据岩石密度体，确定垂直主应力体，包括：
对岩石密度体进行深度数值积分，确定垂直主应力体。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，按如下方式，确定垂直主应力体：
其中，Sv为垂直主应力体；ρr(h)为随深度变化的上覆岩层密度；h为目的层深度；g为重力加速度。
前述提到的确定垂直主应力体的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体，包括：
根据深度域纵波速度体和岩性指示体，确定有效应力体；
根据时间厚度体，垂直主应力体和有效应力体，确定孔隙压力体。
实施例中，孔隙压力本身是重要的气藏钻井开发参数，同时其准确性也将直接影响水平地应力的计算。孔隙压力的计算方法具有多样化，目前常用的基础方法，除了实验以外，主要有Eaton法、Bower法、phillippone法、地层系数剖面法等。但页岩气地区常存在岩性变化和剥蚀作用。由于岩性、剥蚀作用的影响，常规的方法并不适用于孔隙压力的计算。基于常见纵波速度和有效应力之间的关系形式，利用多因素非线性回归方法，考虑纵波速度、岩性和剥蚀作用的影响，首次提出了上述的页岩气采用的孔隙压力计算方法。
对研究区域内的井进行统计，分别计算每口井目的层段的垂直主应力曲线、孔隙压力曲线、有效应力曲线、纵波速度曲线、归一化自然伽马曲线以及目的层段时间厚度曲线(恒值)。
基于确定有效应力体的公式，利用非线性回归方法，对有效应力拟合参数A，B，C进行求取，可得到回归所得的有效应力曲线，和孔隙压力曲线。
利用时间厚度曲线、垂直主应力曲线和有效应力曲线，以及回归所得的孔隙压力曲线，利用确定孔隙压力体的公式进行线性回归，可确定孔隙压力拟合系数a、b、c，则可以得到最终的孔隙压力回归公式
利用垂直主应力体、纵波速度体、岩性指示体、时间厚度体，根据上述回归关系式，即可计算出孔隙压力体。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，按照如下方式，确定有效应力体：
其中，pe为有效应力体；vp为深度域纵波速度体；Li为岩性指示体；A，B，C为有效应力拟合系数，根据研究区域的井的岩石力学参数进行标定获得。
前述提到的确定有效应力体的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，按照如下方式，确定孔隙压力体：
pp＝a+b(Sv-pe)+cTt
其中，pp为孔隙压力体；Sv为垂直主应力体；pe是有效应力体；Tt是时间厚度体；a，b，c为孔隙压力拟合系数，根据井的岩石力学数据进行标定获得。
前述提到的确定孔隙压力体的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体，包括：
根据时间厚度体和深度厚度体，确定最大水平主应力构造系数和最小水平主应力构造系数；
根据弹性参数体、垂直主应力体、孔隙压力体、最大水平主应力构造系数和最小水平主应力构造系数，确定最大水平主应力体和最小水平主应力体。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，按照如下方式，确定最大水平主应力构造系数和最小水平主应力构造系数：
β1＝X1+X2Tt
β2＝Y1+Y2Dt
其中，β2为最大水平主应力构造系数；β1为最小水平主应力构造系数；Dt为研究段深度厚度体；Tt为研究段时间厚度体，X1、X2、Y1、Y2为区域构造拟合系数，根据井的岩石力学参数数据进行标定获得。
前述提到的确定最大水平主应力构造系数和最小水平主应力构造系数的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
本发明实施例提供的一种页岩气地应力确定方法具体实施时，在一个实施例中，按照如下方式，确定最大水平主应力体和最小水平主应力体：
其中，SHmax为最大水平主应力体；Shmin为最小水平主应力体；v为泊松比模型；α为α系数模型；pp为孔隙压力体；β2为最大水平主应力构造系数；β1为最小水平主应力构造系数；Sv为垂直主应力体。
前述提到的确定最大水平主应力体和最小水平主应力体的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
实施例中，构造系数与区域构造背景紧密相关，一般而言，与埋深存在一定关系。但当地层遭受严重剥蚀，存在应力卸载的情况，构造系数受到剥蚀作用的影响，而储层厚度则能反应剥蚀强度。目前来看，最大最小水平应力构造系数在威远页岩气地区与埋深的相关性一般，但最小水平主应力构造系数β1与时间厚度关系良好，最大水平主应力构造系数β2与深度厚度关系良好。
(1)构造系数体构建
利用测井数据，基于确定最大水平主应力体和最小水平主应力体公式计算各井储层段的构造系数。利用构造系数、目的层时间厚度和深度厚度，基于确定最大水平主应力构造系数和最小水平主应力构造系数公式分别进行拟合，取相关性最强的拟合实验，来确定区域构造拟合系数X1、X2、Y1、Y2。
上奥五峰组底向下漂移20m生成新层位，利用该层位和地表层位进行最大最小主应力构造的segy数据体模型构建。具体构建方法与目的层时间厚度体segy数据构建方法一致，生成一个由线号、道号、X坐标、Y坐标、深度、构造系数，组成的文本，并导入landmark速度建模模块中，即可生成segy体。
将所生成的segy体导入建模软件中，并基于前述构造模型进行数据采集，即可获得3D网格目的层段的最大最下主应力构造数属性体。
(2)水平应力体构建
基于确定最大水平主应力体和最小水平主应力体公式，利用泊松比属性体、α属性体、垂直主应力图、孔隙压力体和最大最小主应力构造系数体，进行计算，即可得目的层段的最大水平主应力体和最小水平主应力体。
图4为运行本发明实施的一种岩气地应力确定方法的计算机装置示意图。如图4所示，本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述页岩气地应力确定方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现上述页岩气地应力确定方法的计算机程序。
本发明实施例中还提供了一种页岩气地应力确定装置，如下面的实施例所述。由于该装置解决问题的原理与一种页岩气地应力确定方法相似，因此该装置的实施可以参见一种页岩气地应力确定方法的实施，重复之处不再赘述。
图5为本发明实施例一种岩气地应力确定装置示意图。如图5所示，本发明实施例还提供一种页岩气地应力确定装置，包括：
深度域解释成果确定模块501，用于根据地震数据体，确定深度域解释成果；所述深度域解释成果，包括深度域纵波速度体、深度域横波速度体、深度域密度体、时间厚度数据和深度厚度数据；
地质构造模型建立模块502，用于根据深度域解释成果，建立地质构造模型；
时间厚度体和深度厚度体确定模块503，用于基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体；
弹性参数体确定模块504，用于根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体；
岩石密度体确定模块505，用于根据深度域解释成果和地质构造模型，确定岩石密度体；
岩性指示体确定模块506，用于根据测井曲线和地质构造模型，确定岩性指示体；
垂直主应力体确定模块507，用于根据岩石密度体，确定垂直主应力体；
孔隙压力体确定模块508，用于根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体；
最大水平主应力体和最小水平主应力体确定模块509，用于根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体；
页岩气地应力确定模块510，用于根据垂直主应力体、孔隙压力体、最大水平主应力体和最小水平主应力体，确定页岩气地应力。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，深度域解释成果确定模块，用于根据地震数据体，确定深度域解释成果，包括：
根据地震数据体，获得时间域精细层位解释和断层解释，进行断层组合，为各个断层进行分配命名，获得目的层的时间厚度数据；
基于时间域精细层位解释、断层解释和道集数据，进行叠前同时反演计算，获得时间域解释成果；所述时间域解释成果，包括：时间域纵波速度体、时间域横波速度体和时间域密度体；
构建速度模型，对时间域解释成果进行时深转化，确定深度域的精细层位解释和断层解释；
根据深度域的精细层位解释和断层解释，进行断层组合，为各个断层进行分配命名，获得目的层的深度厚度数据；
利用深度偏移速度场，将深度域数据转化为时间域数据，进行叠前同时反演，获得时间域反演解释成果；所述深度域数据，包括：深度域的精细层位解释和断层解释，目的层的深度厚度数据；所述时间域反演解释成果，包括时间域反演纵波速度体、时间域反演横波速度体和时间域反演密度体；
利用深度偏移速度场，将时间域反演解释成果进行转化，获得深度域纵波速度体、深度域横波速度体和深度域密度体。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，时间厚度体和深度厚度体确定模块，用于基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体，包括：
构建覆盖所有目的层段的时间厚度属性segy数据体和深度厚度属性segy数据体；
将时间厚度属性segy数据体和深度厚度属性segy数据体导入地质构造模型中，对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，弹性参数体确定模块，用于根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体，包括：
根据深度域纵波速度体和深度域横波速度体，基于地质构造模型对地震数据体进行数据采集，确定纵波速度模型和横波速度模型；
根据纵波速度模型和横波速度模型，确定弹性参数体。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，所述弹性参数体，包括：泊松比模型和α系数模型；
弹性参数体确定模块，用于按如下方式，确定弹性参数体：
其中，Vp为纵波速度；Vs为横波速度；v为泊松比模型；α为α系数模型。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，岩石密度体确定模块，用于根据深度域解释成果和地质构造模型，确定岩石密度体，包括：
在深度域密度体完整时，将深度域密度体采集到地质构造模型中，获得岩石密度体；
在深度域密度体不完整时，根据密度变化趋势，补齐密度曲线，利用变差函数分析的方式进行插值，建立目的层上方密度体；将目的层上方密度体和原有不完整深度域密度体合并，获得岩石密度体。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，岩性指示体确定模块，用于根据测井曲线和地质构造模型，确定岩性指示体，包括：
对测井曲线进行归一化处理，利用变差函数的插值方法，基于构造地质模型，建立岩性指示体。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，垂直主应力体确定模块，用于根据岩石密度体，确定垂直主应力体，包括：
对岩石密度体进行深度数值积分，确定垂直主应力体。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，垂直主应力体确定模块，用于按如下方式，确定垂直主应力体：
其中，Sv为垂直主应力体；ρr(h)为随深度变化的上覆岩层密度；h为目的层深度；g为重力加速度。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，孔隙压力体确定模块，用于根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体，包括：
根据深度域纵波速度体和岩性指示体，确定有效应力体；
根据时间厚度体，垂直主应力体和有效应力体，确定孔隙压力体。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，孔隙压力体确定模块，用于按照如下方式，确定有效应力体：
其中，pe为有效应力体；vp为深度域纵波速度体；Li为岩性指示体；A，B，C为有效应力拟合系数，根据研究区域的井的岩石力学参数进行标定获得。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，孔隙压力体确定模块，用于按照如下方式，确定孔隙压力体：
pp＝a+b(Sv-pe)+cTt
其中，pp为孔隙压力体；Sv为垂直主应力体；pe是有效应力体；Tt是时间厚度体；a，b，c为孔隙压力拟合系数，根据井数据的岩石力学参数进行标定获得。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，最大水平主应力体和最小水平主应力体确定模块，根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体，包括：
根据时间厚度体和深度厚度体，确定最大水平主应力构造系数和最小水平主应力构造系数；
根据弹性参数体、垂直主应力体、孔隙压力体、最大水平主应力构造系数和最小水平主应力构造系数，确定最大水平主应力体和最小水平主应力体。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，最大水平主应力体和最小水平主应力体确定模块，按照如下方式，确定最大水平主应力构造系数和最小水平主应力构造系数：
β1＝X1+X2Tt
β2＝Y1+Y2Dt
其中，β2为最大水平主应力构造系数；β1为最小水平主应力构造系数；Dt为研究段深度厚度体；Tt为研究段时间厚度体，X1、X2、Y1、Y2为区域构造拟合系数，根据井的岩石力学参数数据进行标定获得。
本发明实施例提供的一种页岩气地应力确定装置具体实施时，在一个实施例中，最大水平主应力体和最小水平主应力体确定模块，按照如下方式，确定最大水平主应力体和最小水平主应力体：
其中，SHmax为最大水平主应力体；Shmin为最小水平主应力体；v为泊松比模型；α为α系数模型；pp为孔隙压力体；β2为最大水平主应力构造系数；β1为最小水平主应力构造系数；Sv为垂直主应力体。
综上，本发明实施例提供的一种岩气地应力确定方法和装置，利用了地震数据能够反映地层横向的连续性变化以及地震数据能够获得连续的地应力剖面的特性，提出了一种利用地震数据体进行页岩气地应力的确定方法，能够有效反映全工区内地应力大小的展布情况，特别是横向地应力的展布情况。不仅适用于多井区，也适用于少井或无井区；同时由于地震数据预测具有前瞻性，可进行钻前的地应力预测，对于盆地页岩气勘探开发区构造和成藏都较复杂的情况，建立一套基于地震数据的页岩气地应力预测方法，实现了页岩气地应力的准确预测，为页岩气井位部署和压裂施工进行有效支撑。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (18)
1.一种页岩气地应力确定方法，其特征在于，包括：
根据地震数据体，确定深度域解释成果；所述深度域解释成果，包括：深度域纵波速度体、深度域横波速度体、深度域密度体、时间厚度数据和深度厚度数据；
根据深度域解释成果，建立地质构造模型；
基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体；
根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体；
根据深度域解释成果和地质构造模型，确定岩石密度体；
根据测井曲线和地质构造模型，确定岩性指示体；
根据岩石密度体，确定垂直主应力体；
根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体；
根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体；
根据垂直主应力体、孔隙压力体、最大水平主应力体和最小水平主应力体，确定页岩气地应力。
2.如权利要求1所述的方法，其特征在于，根据地震数据体，确定深度域解释成果，包括：
根据地震数据体，获得时间域精细层位解释和断层解释，进行断层组合，为各个断层进行分配命名，获得目的层的时间厚度数据；
基于时间域精细层位解释、断层解释和道集数据，进行叠前同时反演计算，获得时间域解释成果；所述时间域解释成果，包括：时间域纵波速度体、时间域横波速度体和时间域密度体；
构建速度模型，对时间域解释成果进行时深转化，确定深度域的精细层位解释和断层解释；
根据深度域的精细层位解释和断层解释，进行断层组合，为各个断层进行分配命名，获得目的层的深度厚度数据；
利用深度偏移速度场，将深度域数据转化为时间域数据，进行叠前同时反演，获得时间域反演解释成果；所述深度域数据，包括：深度域的精细层位解释和断层解释，目的层的深度厚度数据；所述时间域反演解释成果，包括时间域反演纵波速度体、时间域反演横波速度体和时间域反演密度体；
利用深度偏移速度场，将时间域反演解释成果进行转化，获得深度域纵波速度体、深度域横波速度体和深度域密度体。
3.如权利要求1所述的方法，其特征在于，基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体，包括：
构建覆盖所有目的层段的时间厚度属性segy数据体和深度厚度属性segy数据体；
将时间厚度属性segy数据体和深度厚度属性segy数据体导入地质构造模型中，对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体。
4.如权利要求1所述的方法，其特征在于，根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体，包括：
根据深度域纵波速度体和深度域横波速度体，基于地质构造模型对地震数据体进行数据采集，确定纵波速度模型和横波速度模型；
根据纵波速度模型和横波速度模型，确定弹性参数体。
5.如权利要求4所述的方法，其特征在于，所述弹性参数体，包括：泊松比模型和α系数模型；
按如下方式，确定弹性参数体：
其中，Vp为纵波速度；Vs为横波速度；v为泊松比模型；α为α系数模型。
6.如权利要求1所述的方法，其特征在于，根据深度域解释成果和地质构造模型，确定岩石密度体，包括：
在深度域密度体完整时，将深度域密度体采集到地质构造模型中，获得岩石密度体；
在深度域密度体不完整时，根据密度变化趋势，补齐密度曲线，利用变差函数分析的方式进行插值，建立目的层上方密度体；将目的层上方密度体和原有不完整深度域密度体合并，获得岩石密度体。
7.如权利要求1所述的方法，其特征在于，根据测井曲线和地质构造模型，确定岩性指示体，包括：
对测井曲线进行归一化处理，利用变差函数的插值方法，基于构造地质模型，建立岩性指示体。
8.如权利要求1所述的方法，其特征在于，根据岩石密度体，确定垂直主应力体，包括：
对岩石密度体进行深度数值积分，确定垂直主应力体。
9.如权利要求8所述的方法，其特征在于，按如下方式，确定垂直主应力体：
其中，Sv为垂直主应力体；ρr(h)为随深度变化的上覆岩层密度；h为目的层深度；g为重力加速度。
10.如权利要求5所述的方法，其特征在于，根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体，包括：
根据深度域纵波速度体和岩性指示体，确定有效应力体；
根据时间厚度体，垂直主应力体和有效应力体，确定孔隙压力体。
11.如权利要求10所述的方法，其特征在于，按照如下方式，确定有效应力体：
其中，pe为有效应力体；vp为深度域纵波速度体；Li为岩性指示体；A，B，C为有效应力拟合系数，根据研究区域的井的岩石力学参数进行标定获得。
12.如权利要求11所述的方法，其特征在于，按照如下方式，确定孔隙压力体：
pp＝a+b(Sv-pe)+cTt
其中，pp为孔隙压力体；Sv为垂直主应力体；pe是有效应力体；Tt是时间厚度体；a，b，c为孔隙压力拟合系数，根据井数据的岩石力学参数进行标定获得。
13.如权利要求1所述的方法，其特征在于，根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体，包括：
根据时间厚度体和深度厚度体，确定最大水平主应力构造系数和最小水平主应力构造系数；
根据弹性参数体、垂直主应力体、孔隙压力体、最大水平主应力构造系数和最小水平主应力构造系数，确定最大水平主应力体和最小水平主应力体。
14.如权利要求13所述的方法，其特征在于，按照如下方式，确定最大水平主应力构造系数和最小水平主应力构造系数：
β1＝X1+X2Tt
β2＝Y1+Y2Dt
其中，β2为最大水平主应力构造系数；β1为最小水平主应力构造系数；Dt为研究段深度厚度体；Tt为研究段时间厚度体，X1、X2、Y1、Y2为区域构造拟合系数，根据井的岩石力学参数数据进行标定获得。
15.如权利要求14所述的方法，其特征在于，按照如下方式，确定最大水平主应力体和最小水平主应力体：
其中，SHmax为最大水平主应力体；Shmin为最小水平主应力体；v为泊松比模型；α为α系数模型；pp为孔隙压力体；β2为最大水平主应力构造系数；β1为最小水平主应力构造系数；Sv为垂直主应力体。
16.一种页岩气地应力确定装置，其特征在于，包括：
深度域解释成果确定模块，用于根据地震数据体，确定深度域解释成果；所述深度域解释成果，包括：深度域纵波速度体、深度域横波速度体、深度域密度体、时间厚度数据和深度厚度数据；
地质构造模型建立模块，用于根据深度域解释成果，建立地质构造模型；
时间厚度体和深度厚度体确定模块，用于基于地质构造模型，通过对深度域解释成果进行数据采集，确定时间厚度体和深度厚度体；
弹性参数体确定模块，用于根据深度域解释成果，基于地质构造模型对地震数据体进行数据采集，确定弹性参数体；
岩石密度体确定模块，用于根据深度域解释成果和地质构造模型，确定岩石密度体；
岩性指示体确定模块，用于根据测井曲线和地质构造模型，确定岩性指示体；
垂直主应力体确定模块，用于根据岩石密度体，确定垂直主应力体；
孔隙压力体确定模块，用于根据深度域解释成果、岩性指示体、时间厚度体和垂直主应力体，确定孔隙压力体；
最大水平主应力体和最小水平主应力体确定模块，用于根据时间厚度体和深度厚度体，结合弹性参数体、垂直主应力体和孔隙压力体，确定最大水平主应力体和最小水平主应力体；
页岩气地应力确定模块，用于根据垂直主应力体、孔隙压力体、最大水平主应力体和最小水平主应力体，确定页岩气地应力。
17.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现权利要求1至15任一项所述页岩气地应力确定方法。
18.一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现权利要求1至15任一项所述页岩气地应力确定方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010325413.2A CN113552621A (zh) | 2020-04-23 | 2020-04-23 | 页岩气地应力确定方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010325413.2A CN113552621A (zh) | 2020-04-23 | 2020-04-23 | 页岩气地应力确定方法和装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113552621A true CN113552621A (zh) | 2021-10-26 |
Family
ID=78100970
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010325413.2A Pending CN113552621A (zh) | 2020-04-23 | 2020-04-23 | 页岩气地应力确定方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113552621A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111948718A (zh) * | 2020-07-21 | 2020-11-17 | 中国石油天然气集团有限公司 | 页岩气储层总有机碳含量预测方法及装置 |
| CN114966855A (zh) * | 2022-07-28 | 2022-08-30 | 煤炭科学研究总院有限公司 | 探测煤岩体高应力区的方法、装置、设备和介质 |
| CN115291279A (zh) * | 2022-07-21 | 2022-11-04 | 中国地质调查局成都地质调查中心（西南地质科技创新中心） | 表征盆外构造变形强度的方法、装置、存储介质和设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030004648A1 (en) * | 2001-03-13 | 2003-01-02 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| CN104977618A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN105467438A (zh) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | 一种基于三模量的泥页岩地应力三维地震表征方法 |
| CN108397130A (zh) * | 2018-02-02 | 2018-08-14 | 中国石油天然气集团有限公司 | 钻井方法和装置 |
-
2020
- 2020-04-23 CN CN202010325413.2A patent/CN113552621A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030004648A1 (en) * | 2001-03-13 | 2003-01-02 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| CN104977618A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN105467438A (zh) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | 一种基于三模量的泥页岩地应力三维地震表征方法 |
| CN108397130A (zh) * | 2018-02-02 | 2018-08-14 | 中国石油天然气集团有限公司 | 钻井方法和装置 |
Non-Patent Citations (1)
| Title |
|---|
| 罗兵;: "涪陵页岩气田一期产建区地应力场预测研究", 江汉石油科技, no. 04, 25 December 2017 (2017-12-25) * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111948718A (zh) * | 2020-07-21 | 2020-11-17 | 中国石油天然气集团有限公司 | 页岩气储层总有机碳含量预测方法及装置 |
| CN111948718B (zh) * | 2020-07-21 | 2023-11-28 | 中国石油天然气集团有限公司 | 页岩气储层总有机碳含量预测方法及装置 |
| CN115291279A (zh) * | 2022-07-21 | 2022-11-04 | 中国地质调查局成都地质调查中心（西南地质科技创新中心） | 表征盆外构造变形强度的方法、装置、存储介质和设备 |
| CN115291279B (zh) * | 2022-07-21 | 2024-05-07 | 中国地质调查局成都地质调查中心（西南地质科技创新中心） | 表征盆外构造变形强度的方法、装置、存储介质和设备 |
| CN114966855A (zh) * | 2022-07-28 | 2022-08-30 | 煤炭科学研究总院有限公司 | 探测煤岩体高应力区的方法、装置、设备和介质 |
| CN114966855B (zh) * | 2022-07-28 | 2022-10-04 | 煤炭科学研究总院有限公司 | 探测煤岩体高应力区的方法、装置、设备和介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2629123B1 (en) | Simulation model optimization | |
| US7373251B2 (en) | Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data | |
| US8359184B2 (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
| AU2013397958B2 (en) | A simulation-to-seismic workflow construed from core based rock typing and enhanced by rock replacement modeling | |
| CN113552621A (zh) | 页岩气地应力确定方法和装置 | |
| US10620339B2 (en) | Static earth model calibration methods and systems using tortuosity evaluations | |
| OA12151A (en) | Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data. | |
| CN110056346B (zh) | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 | |
| CN110244021B (zh) | 一种基于各向异性插值的地层分层方法 | |
| Kamali et al. | 3D geostatistical modeling and uncertainty analysis in a carbonate reservoir, SW Iran | |
| CN113313825A (zh) | 一种油藏储层三维地质建模方法和装置 | |
| CN105242307A (zh) | 复杂碳酸盐储层地震孔隙度获取方法及装置 | |
| Tebo et al. | Use of volume-based 3-D seismic attribute analysis to characterize physical-property distribution: a case study to delineate sedimentologic heterogeneity at the Appleton field, southwestern Alabama, USA | |
| CN104570064A (zh) | 一种砂岩地层横波速度计算方法 | |
| Liu et al. | Geostatistical modeling for fine reservoir description of Wei2 block of Weicheng oilfield, Dongpu depression, China | |
| Lutome et al. | 3D geocellular modeling for reservoir characterization of lacustrine turbidite reservoirs: Submember 3 of the third member of the Eocene Shahejie Formation, Dongying depression, Eastern China | |
| Ali et al. | Applying facies characterization to build 3D rock-type model for Khasib Formation, Amara Oil Field | |
| CN110568501B (zh) | 一种地层中烃源岩纵向非均质特征的定量表征方法 | |
| Masoud et al. | Reservoir Characterization and Geostatistical Model of the Cretaceous and Cambrian-Ordovician Reservoir Intervals, Meghil Field, Sirte Basin, Libya | |
| CN111815769A (zh) | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 | |
| CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
| CN107589457A (zh) | 测井资料尺度粗化的方法和装置 | |
| Papiernik et al. | Hydrocarbon resources assessment for selected area of Lower Paleozoic basin in Poland based on spatial modeling | |
| CN117784245A (zh) | 一种时深转换方法、装置、存储介质及电子设备 | |
| CN114429016A (zh) | 页岩气藏甜点模型构建方法、装置、存储介质以及设备 |
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 |