CN113376692B - 致密砂岩气水平井压裂改造方案优化方法及装置 - Google Patents
致密砂岩气水平井压裂改造方案优化方法及装置 Download PDFInfo
- Publication number
- CN113376692B CN113376692B CN202010161440.0A CN202010161440A CN113376692B CN 113376692 B CN113376692 B CN 113376692B CN 202010161440 A CN202010161440 A CN 202010161440A CN 113376692 B CN113376692 B CN 113376692B
- Authority
- CN
- China
- Prior art keywords
- curve
- avo
- gas
- transformation
- parameter curve
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000004048 modification Effects 0.000 title claims abstract description 27
- 238000012986 modification Methods 0.000 title claims abstract description 27
- 230000009466 transformation Effects 0.000 claims abstract description 68
- 239000011435 rock Substances 0.000 claims abstract description 40
- 230000007480 spreading Effects 0.000 claims abstract description 15
- 238000003892 spreading Methods 0.000 claims abstract description 15
- 230000005012 migration Effects 0.000 claims abstract description 13
- 238000013508 migration Methods 0.000 claims abstract description 13
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 230000000704 physical effect Effects 0.000 claims abstract description 6
- 230000002159 abnormal effect Effects 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 abstract description 4
- 239000007789 gas Substances 0.000 description 78
- 238000010586 diagram Methods 0.000 description 11
- 238000009826 distribution Methods 0.000 description 9
- 238000005553 drilling Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000006073 displacement reaction Methods 0.000 description 5
- 239000004576 sand Substances 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000007787 solid Substances 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. for interpretation or for event detection
- 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
- G01V2210/632—Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
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
本发明提供一种致密砂岩气水平井压裂改造方案优化方法及装置,包括:根据完钻井的测井和录井数据中选取的物、电性参数曲线确定弹性参数曲线;将物性、电性、弹性参数曲线两两交会,确定最小泊松比与含气饱和度呈负相关,杨氏模量与岩石脆性呈正相关;对AVO模型正演分析确定储层含气性;获取待分析目的层段的三维地震叠前偏移后的CRP道集分角度叠加数据;基于上述理论对地震数据、测井和录井数据进行叠前弹性反演,根据获得的最小泊松比和杨氏模量刻画有效储层空间展布形态和确定井筒周边岩石脆性特征,从而确定改造段数、方向和规模进行压裂改造方案优化。本方案利用最小泊松比和杨氏模量属性优化水平井压裂改造方案,实际应用效果好。
Description
技术领域
本发明涉及水平井压裂技术领域,特别涉及一种致密砂岩气水平井压裂改造方案优化方法及装置。
背景技术
致密砂岩气藏具有“储层薄、低孔、低渗、低压、低丰度、强非均质性”的特征,单井产量低、压力下降快、稳产和开发难度大。为了高效动用地质储量,主要以水平井开发为主。致密砂岩气藏的特性决定了压裂改造是其经济开发的必然选择。
常规水平井压裂改造方案设计主要依托测、录井资料,参考水平段随钻伽玛曲线、气测曲线及录井岩屑资料,通过计算获得岩石的泊松比、杨氏模量等弹性参数,再根据弹性参数与有效储层的相关性,确定压裂改造位置及参数。由于致密气砂岩储层非均质性强,储层纵向相互叠置、横向变化剧烈,空间展布非常复杂,仅凭“一孔之见”(即只有一口水平井的测录井数据,钻孔周边的地质条件没有资料)资料确定压裂改造方案设计,依据不充分,容易造成浪费或者改造不充分等问题。
发明内容
本发明实施例提供了一种致密砂岩气水平井压裂改造方案优化方法及装置,解决了现有技术中仅凭“一孔之见”资料确定压裂改造方案设计,依据不充分,容易造成浪费或者改造不充分的技术问题。
本发明实施例提供了一种致密砂岩气水平井压裂改造方案优化方法,该方法包括:
获取完钻井的测井数据和录井数据,从中选取含气储层段的物性参数曲线和电性参数曲线,所述物性参数曲线和电性参数曲线包括伽玛曲线、纵波时差曲线、横波时差曲线、密度曲线、含气饱和度曲线和岩石脆性曲线;
根据物性参数曲线和电性参数曲线确定含气储层段的弹性参数曲线,弹性参数包括最小泊松比和杨氏模量;
将物性参数曲线和电性参数曲线分别和弹性参数曲线进行两两交会,根据交会结果确定地质情况与地震相关信息规律,所述地质情况与地震相关信息规律为最小泊松比与含气饱和度两者呈负相关,杨氏模量与岩石脆性呈正相关;
根据物性参数曲线和电性参数曲线建立振幅随偏移距变化AVO模型,对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型用于判断储层含气性;
获取待分析目的层段的三维地震叠前偏移后的CRP道集分角度叠加数据;
基于所述正演结果,对待分析目的层段的三维地震叠前偏移后的CRP道集分角度叠加数据、完钻井的测井数据和录井数据进行叠前弹性反演,获得待分析目的层段的弹性参数;
基于地质情况与地震相关信息规律,根据待分析目的层段弹性参数中的最小泊松比刻画有效储层空间展布形态,根据待分析目的层段弹性参数中的杨氏模量确定井筒周边岩石脆性特征;
根据所述有效储层空间展布形态和井筒周边岩石脆性特征确定水平井压裂的改造段数、改造方向和改造规模,基于所述改造段数、改造方向和改造规模对原始水平井的压裂改造方案进行优化;
对于没有实测横波时差的完钻井,采用如下方式确定横波时差曲线:
建立致密砂岩气藏岩石物理模型,根据所述致密砂岩气藏岩石物理模型拟合横波时差曲线。
本发明实施例还提供了一种致密砂岩气水平井压裂改造方案优化装置,该装置包括:
数据获取模块,用于完钻井的测井数据和录井数据,从中选取含气储层段的物性参数曲线和电性参数曲线,所述物性参数曲线和电性参数曲线包括伽玛曲线、纵波时差曲线、横波时差曲线、密度曲线、含气饱和度曲线和岩石脆性曲线;
弹性参数确定模块,用于根据物性参数曲线和电性参数曲线确定含气储层段的弹性参数曲线,弹性参数包括最小泊松比和杨氏模量;
地质情况与地震相关信息规律确定模块,用于将物性参数曲线和电性参数曲线分别和弹性参数曲线进行两两交会,根据交会结果确定地质情况与地震相关信息规律,所述地质情况与地震相关信息规律为最小泊松比与含气饱和度两者呈负相关,杨氏模量与岩石脆性呈正相关;
AVO模型正演模块,用于根据物性参数曲线和电性参数曲线建立振幅随偏移距变化AVO模型,对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型用于判断储层含气性;
所述数据获取模块还用于:获取三维地震叠前偏移后的CRP道集分角度叠加数据;
叠前弹性反演模块,用于基于所述正演结果,对待分析目的层段的三维地震叠前偏移后的CRP道集分角度叠加数据、完钻井的测井数据和录井数据进行叠前弹性反演,获得待分析目的层段的弹性参数;
特性确定模块,用于基于地质情况与地震相关信息规律,根据待分析目的层段弹性参数中的最小泊松比刻画有效储层空间展布形态,根据待分析目的层段弹性参数中的杨氏模量确定井筒周边岩石脆性特征;
水平井压裂改造优化模块,用于根根据所述有效储层空间展布形态和井筒周边岩石脆性特征确定水平井压裂的改造段数、改造方向和改造规模,基于所述改造段数、改造方向和改造规模对原始水平井的压裂改造方案进行优化;
所述数据获取模块还用于:
对于没有实测横波时差的完钻井,采用如下方式确定横波时差曲线:建立致密砂岩气藏岩石物理模型,根据所述致密砂岩气藏岩石物理模型拟合横波时差曲线。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述所述方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述所述方法的计算机程序。
在本发明实施例中,根据完钻井的测井数据和录井数据确定最小泊松比与含气饱和度两者呈负相关,杨氏模量与岩石脆性呈正相关,根据完钻井的测井数据和录井数据建立振幅随偏移距变化AVO模型,对AVO模型正演分析,根据正演结果确定储层含气性,对获取的三维地震叠前偏移后的CRP道集分角度叠加数据、完钻井的测井数据和录井数据进行叠前弹性反演,获得最小泊松比及杨氏模量等岩石物理弹性参数,然后根据上述的结论,精细刻画有效储层空间展布形态、明确井筒周边岩石脆性特征,在其基础上,优化水平井压裂改造方案,优选改造位置、明确改造方向、确定改造规模。在高效动用井间与单井储量的基础上,提高单井产量、实现降本增效。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种致密砂岩气水平井压裂改造方案优化方法流程图;
图2是本发明实施例提供的一种最小泊松比分布图;
图3是本发明实施例提供的一种jing72-60H2井压裂参数设计优化(13段19簇)示意图;
图4是本发明实施例提供的一种致密砂岩气水平井压裂改造方案优化装置结构框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
技术术语解释
测井,也叫地球物理测井或石油测井,简称测井,是利用岩层的电化学特性、导电特性、声学特性、放射性等地球物理特性,测量地球物理参数的方法,属于应用地球物理方法(包括重、磁、电、震、测井)之一。石油钻井时,在钻到设计井深深度后都必须进行测井,又称完井电测,以获得各种石油地质及工程技术资料,作为完井和开发油田的原始资料。
录井是用岩矿分析、地球物理、地球化学等方法,观察、采集、收集、记录、分析随钻过程中的固体、液体、气体等井筒返出物信息,以此建立录井地质剖面、发现油气显示、评价油气层,并为石油工程(投资方、钻井工程、其它工程)提供钻井信息服务的过程。
AVO(Amplitude variation with offset,振幅随偏移距的变化)技术用于研究地震反射振幅随炮点与接收器之间的距离即炮检距(或入射角)的变化特征来探讨反射系数响应随炮检距(或入射角)的变化,进而确定反射界面上覆、下伏介质的岩性特征及物性参数。
CRP(common reflection point)道集是经过叠前时间偏移后的共反射点道集,它基本消除了界面弯曲和倾斜地层对振幅的影响。
在本发明实施例中,提供了一种致密砂岩气水平井压裂改造方案优化方法,如图1所示,该方法包括:
步骤101:获取完钻井的测井数据和录井数据,从中选取含气储层段的物性参数曲线和电性参数曲线,所述物性参数曲线和电性参数曲线包括伽玛曲线、纵波时差曲线、横波时差曲线、密度曲线、含气饱和度曲线和岩石脆性曲线。
对于没有实测横波时差曲线的完钻井,依据Xu-White模型,建立致密砂岩气藏岩石物理模型,根据致密砂岩气藏岩石物理模型拟合横波时差曲线,为叠前反演提供基础井资料。选取的伽马曲线可以判别砂泥岩,从而确定储层和非储层。
步骤102:根据物性参数曲线和电性参数曲线确定含气储层段的弹性参数曲线,弹性参数包括最小泊松比和杨氏模量。
具体的,是根据纵波时差曲线、横波时差曲线、密度曲线确定最小泊松比和杨氏模量,具体公式如下:
E=2×ρ×Vs2×(1+σ);
其中,PI为纵波阻抗,PI=Vp×ρ;SI为横波阻抗,SI=Vs×ρ;Vp为纵波速度,Vs为横波速度,ρ为密度,σ为泊松比,E为杨氏模量。
步骤103:将物性参数曲线和电性参数曲线分别和弹性参数曲线进行两两交会,根据交会结果确定地质情况与地震相关信息规律,所述地质情况与地震相关信息规律为最小泊松比与含气饱和度两者呈负相关,杨氏模量与岩石脆性(岩石受力后变形很小即破裂的性质称为脆性)呈正相关。得出这个结论可以为叠前反演储层预测及压裂方案参数优化提供依据。
步骤104:根据物性参数曲线和电性参数曲线建立振幅随偏移距变化AVO模型,对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型用于判断储层含气性。
具体的,根据纵波时差曲线、横波时差曲线、密度曲线确定含气储层段的纵波阻抗、横波阻抗(所用公式为PI=Vp×ρ,SI=Vs×ρ),根据含气储层段的纵波阻抗、横波阻抗、密度建立振幅随偏移距变化AVO模型,对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型判断储层含气性。
比如,若通过AVO模型正演分析确定致密砂岩气藏AVO表现为第III类异常,即含气储层段顶界地震反射振幅强度随偏移距增大而增强,那么可利用该异常定性判别储层含气性,为叠前弹性反演奠定理论基础。
具体的,根据某一深度含气储层的纵、横波速度及密度参数,建立含气储层段振幅随随偏移距变化的AVO响应模型,即以Zoeppritz方程为基础,求解储层与非储层介质之间的反射界面随偏移距(炮点与接收的点之间的距离)增大而产生的反射系数的变化,再利用Rick子波与反射系数褶积,最终获得含气储层段AVO正演响应特征;进而确定反射界面上覆、下伏介质的岩性特征及物性参数。
步骤105:获取(宽方位、高密度)三维地震叠前偏移后的CRP道集分角度叠加数据。
步骤106:基于所述正演结果,对待分析目的层段的三维地震叠前偏移后的CRP道集分角度叠加数据、完钻井的测井数据和录井数据进行叠前弹性反演,获得待分析目的层段的弹性参数;
具体的,以Zoeppritz方程为理论基础,进行叠前弹性反演。叠前反演是以描述平面波反射和投射的Zoeppritz方程为基础,利用反射系数随入射角变化与地层弹性参数间的关系,求解方程,从叠前地震数据中估算岩石的弹性参数(纵波速度、横波速度、密度、泊松比等),进而利用这些弹性参数进行岩性分析和含油气的预测。实际应用中,把经过叠前精细保幅处理和偏移的共反射点道集,分成多个入射角叠加,形成不同入射角叠加剖面;再利用测井数据分别提取对应不同入射角叠加剖面的子波,保证子波频率、振幅与角道集叠加剖面之间的能量、频率特征相关度达到最佳;再利用解释层位、完钻井测井数据建立低频初始模型;求解Zoeppritz方程,最终获得纵波速度、横波速度及密度参数。
步骤107:基于地质情况与地震相关信息规律,根据待分析目的层段弹性参数中的最小泊松比刻画有效储层空间展布形态,根据待分析目的层段弹性参数中的杨氏模量确定井筒周边岩石脆性特征;
步骤108:根据所述有效储层空间展布形态和井筒周边岩石脆性特征确定水平井压裂的改造段数、改造方向和改造规模,基于所述改造段数、改造方向和改造规模对原始水平井的压裂改造方案进行优化。
具体的,获得了有效储层空间展布形态后,结合随钻测、录井资料,首先优选低伽玛、高气测优质储层段;其次筛选部分仅井筒附近发育3-5米薄泥岩隔层后即为优质储层段的高伽玛、低气测非储层段;最终结合井筒周边岩石脆性特征确定改造段数。
在优选改造段数的基础上,利用有效储层的空间展布三维数据体,以与水平轨迹正交、斜交的方式,截取任意方向的有效储层横截面,落实有效储层发育范围并结合井筒周边岩石脆性特征,最终明确压裂改造方向。
在优选改造段数及方向的基础上,结合有效储层展布范围及井筒周边岩石脆性展布特征,适度增加或减少改造规模,在高效动用井间与单井储量的基础上,提高单井产量、实现降本增效。
下面举例说明本发明方法的效果。
以鄂尔多斯盆地致密砂岩气藏为研究对象,图2是提取的最小泊松比分布图,图中的jing72-60、jing72-60H2、jing72-58H2、jing72-59H2等表示的是水平井号。图2中黑色代表含气饱和度高,白色代表含气饱和度低,因此,利用最小泊松比的分布可预测目的层段含气性分布范围。
图3是水平井水平段纵向与横向杨氏模量属性剖面,图3中的(a)图为jing72-60H2水平井沿水平段方向的杨氏模量属性的纵向剖面,给出了现有压裂方案中的压裂部位及其排量和砂量的数值,其中的h5表示二叠系下石盒子组盒5段底部、h83表示二叠系下石盒子组盒83段底部、s13表示二叠系山西组山13段底部、s23表示二叠系山西组山23段底部等等。基于现有的压裂部位及其排量和砂量的数值,优选杨氏模量较大值并结合水平段含气性确定水平井水平段压裂改造的位置,根据岩石脆性的横向展布并结合有效储层展布特征优选改造方向、确定改造规模,图3中的(b)和(c)图为横向剖面,(b)图为选择(a)图中的第三段裂缝进行改造设计,该段的排量为8m3/min、砂量为60m3,由于该段属于致密带--泥岩段--致密带,对其进行改造的方案是:多簇射孔、大排量压裂,突破泥质砂岩隔层。(c)图为选择(a)图中的第十一段裂缝进行改造设计,该段的排量为7m3/min、砂量为50m3,由于该段属于Ⅰ类储层--致密带--Ⅲ类储层,对其进行改造的方案是:单簇射孔、适度提高排量,实现优质层充分动用。
基于同一发明构思,本发明实施例中还提供了一种致密砂岩气水平井压裂改造方案优化装置,如下面的实施例所述。由于致密砂岩气水平井压裂改造方案优化装置解决问题的原理与致密砂岩气水平井压裂改造方案优化方法相似,因此致密砂岩气水平井压裂改造方案优化装置的实施可以参见致密砂岩气水平井压裂改造方案优化方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图4是本发明实施例的致密砂岩气水平井压裂改造方案优化装置结构框图,如图4所示,包括:
数据获取模块401,用于完钻井的测井数据和录井数据,从中选取含气储层段的物性参数曲线和电性参数曲线,所述物性参数曲线和电性参数曲线包括伽玛曲线、纵波时差曲线、横波时差曲线、密度曲线、含气饱和度曲线和岩石脆性曲线;
弹性参数确定模块402,用于根据物性参数曲线和电性参数曲线确定含气储层段的弹性参数曲线,弹性参数包括最小泊松比和杨氏模量;
地质情况与地震相关信息规律确定模块403,用于将物性参数曲线和电性参数曲线分别和弹性参数曲线进行两两交会,根据交会结果确定地质情况与地震相关信息规律,所述地质情况与地震相关信息规律为最小泊松比与含气饱和度两者呈负相关,杨氏模量与岩石脆性呈正相关;
AVO模型正演模块404,用于根据物性参数曲线和电性参数曲线建立振幅随偏移距变化AVO模型,对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型用于判断储层含气性;
所述数据获取模块401还用于:获取三维地震叠前偏移后的CRP道集分角度叠加数据;
叠前弹性反演模块405,用于基于所述正演结果,对待分析目的层段的三维地震叠前偏移后的CRP道集分角度叠加数据、完钻井的测井数据和录井数据进行叠前弹性反演,获得待分析目的层段的弹性参数;
特性确定模块406,用于基于地质情况与地震相关信息规律,根据待分析目的层段弹性参数中的最小泊松比刻画有效储层空间展布形态,根据待分析目的层段弹性参数中的杨氏模量确定井筒周边岩石脆性特征;
水平井压裂改造优化模块407,用于根据所述有效储层空间展布形态和井筒周边岩石脆性特征确定水平井压裂的改造段数、改造方向和改造规模,基于所述改造段数、改造方向和改造规模对原始水平井的压裂改造方案进行优化。
在本发明实施例中,所述数据获取模块401还用于:
对于没有实测横波时差的完钻井,采用如下方式确定横波时差曲线:建立致密砂岩气藏岩石物理模型,根据所述致密砂岩气藏岩石物理模型拟合横波时差曲线。
在本发明实施例中,按照如下公式确定最小泊松比和杨氏模量:
E=2×ρ×Vs2×(1+σ);
其中,PI为纵波阻抗,PI=Vp×ρ;SI为横波阻抗,SI=Vs×ρ;Vp为纵波速度,Vs为横波速度,ρ为密度,σ为泊松比,E为杨氏模量。
在本发明实施例中,所述AVO模型正演模块404具体用于:
根据纵波时差曲线、横波时差曲线、密度曲线确定含气储层段的纵波阻抗、横波阻抗;
根据含气储层段的纵波阻抗、横波阻抗、密度建立振幅随偏移距变化AVO模型;
对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型判断储层含气性。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述所述方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述所述方法的计算机程序。
综上所述,本发明提出的致密砂岩气水平井压裂改造方案优化方法及装置可以获得如下有益效果:
首次以鄂尔多斯盆地宽方位、高密度三维地震资料为基础,针对致密砂岩气藏储层薄、强非均质性的特征,通过岩石物理分析、AVO正演模拟及叠前弹性反演,获得纵波速度、横波速度、密度、最小泊松比及杨氏模量等多种岩石物理弹性参数,表明最小泊松比与含气饱和度呈负相关,杨氏模量与岩石脆性呈正相关,进而在精细刻画有效储层空间展布形态、明确井筒周边岩石脆性特征的基础上,优化水平井压裂改造方案,优选改造位置、明确改造方向、确定改造规模。在高效动用井间与单井储量的基础上,提高单井产量、实现降本增效。在水平井压裂实际应用中,该方法操作简单,具有较强指导性,效果良好,取得了可观的生产效益和经济效益,充分说明了该方法的可行性和实用性。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种致密砂岩气水平井压裂改造方案优化方法,其特征在于,包括:
获取完钻井的测井数据和录井数据,从中选取含气储层段的物性参数曲线和电性参数曲线,所述物性参数曲线和电性参数曲线包括伽玛曲线、纵波时差曲线、横波时差曲线、密度曲线、含气饱和度曲线和岩石脆性曲线;
根据物性参数曲线和电性参数曲线确定含气储层段的弹性参数曲线,弹性参数包括最小泊松比和杨氏模量;
将物性参数曲线和电性参数曲线分别和弹性参数曲线进行两两交会,根据交会结果确定地质情况与地震相关信息规律,所述地质情况与地震相关信息规律为最小泊松比与含气饱和度两者呈负相关,杨氏模量与岩石脆性呈正相关;
根据物性参数曲线和电性参数曲线建立振幅随偏移距变化AVO模型,对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型判断储层含气性;
获取待分析目的层段的三维地震叠前偏移后的CRP道集分角度叠加数据;
基于所述正演结果,对待分析目的层段的三维地震叠前偏移后的CRP道集分角度叠加数据、完钻井的测井数据和录井数据进行叠前弹性反演,获得待分析目的层段的弹性参数;
基于地质情况与地震相关信息规律,根据待分析目的层段弹性参数中的最小泊松比刻画有效储层空间展布形态,根据待分析目的层段弹性参数中的杨氏模量确定井筒周边岩石脆性特征;
根据所述有效储层空间展布形态和井筒周边岩石脆性特征确定水平井压裂的改造段数、改造方向和改造规模,基于所述改造段数、改造方向和改造规模对原始水平井的压裂改造方案进行优化;
对于没有实测横波时差的完钻井,采用如下方式确定横波时差曲线:
建立致密砂岩气藏岩石物理模型,根据所述致密砂岩气藏岩石物理模型拟合横波时差曲线。
3.如权利要求1所述的致密砂岩气水平井压裂改造方案优化方法,其特征在于,根据物性参数曲线和电性参数曲线建立振幅随偏移距变化AVO模型,对AVO模型正演分析,获得正演结果,包括:
根据纵波时差曲线、横波时差曲线、密度曲线确定含气储层段的纵波阻抗、横波阻抗;
根据含气储层段的纵波阻抗、横波阻抗、密度建立振幅随偏移距变化AVO模型;
对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型判断储层含气性。
4.一种致密砂岩气水平井压裂改造方案优化装置,其特征在于,包括:
数据获取模块,用于完钻井的测井数据和录井数据,从中选取含气储层段的物性参数曲线和电性参数曲线,所述物性参数曲线和电性参数曲线包括伽玛曲线、纵波时差曲线、横波时差曲线、密度曲线、含气饱和度曲线和岩石脆性曲线;
弹性参数确定模块,用于根据物性参数曲线和电性参数曲线确定含气储层段的弹性参数曲线,弹性参数包括最小泊松比和杨氏模量;
地质情况与地震相关信息规律确定模块,用于将物性参数曲线和电性参数曲线分别和弹性参数曲线进行两两交会,根据交会结果确定地质情况与地震相关信息规律,所述地质情况与地震相关信息规律为最小泊松比与含气饱和度两者呈负相关,杨氏模量与岩石脆性呈正相关;
AVO模型正演模块,用于根据物性参数曲线和电性参数曲线建立振幅随偏移距变化AVO模型,对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型判断储层含气性;
所述数据获取模块还用于:获取三维地震叠前偏移后的CRP道集分角度叠加数据;
叠前弹性反演模块,用于基于所述正演结果,对待分析目的层段的三维地震叠前偏移后的CRP道集分角度叠加数据、完钻井的测井数据和录井数据进行叠前弹性反演,获得待分析目的层段的弹性参数;
特性确定模块,用于基于地质情况与地震相关信息规律,根据待分析目的层段弹性参数中的最小泊松比刻画有效储层空间展布形态,根据待分析目的层段弹性参数中的杨氏模量确定井筒周边岩石脆性特征;
水平井压裂改造优化模块,用于根据所述有效储层空间展布形态和井筒周边岩石脆性特征确定水平井压裂的改造段数、改造方向和改造规模,基于所述改造段数、改造方向和改造规模对原始水平井的压裂改造方案进行优化;
所述数据获取模块还用于:
对于没有实测横波时差的完钻井,采用如下方式确定横波时差曲线:建立致密砂岩气藏岩石物理模型,根据所述致密砂岩气藏岩石物理模型拟合横波时差曲线。
6.如权利要求4所述的致密砂岩气水平井压裂改造方案优化装置,其特征在于,所述AVO模型正演模块具体用于:
根据纵波时差曲线、横波时差曲线、密度曲线确定含气储层段的纵波阻抗、横波阻抗;
根据含气储层段的纵波阻抗、横波阻抗、密度建立振幅随偏移距变化AVO模型;
对AVO模型正演分析,根据正演结果确定致密砂岩气藏AVO异常类型,根据所述AVO异常类型判断储层含气性。
7.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至3任一项所述方法。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至3任一项所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010161440.0A CN113376692B (zh) | 2020-03-10 | 2020-03-10 | 致密砂岩气水平井压裂改造方案优化方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010161440.0A CN113376692B (zh) | 2020-03-10 | 2020-03-10 | 致密砂岩气水平井压裂改造方案优化方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113376692A CN113376692A (zh) | 2021-09-10 |
CN113376692B true CN113376692B (zh) | 2023-05-26 |
Family
ID=77568670
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010161440.0A Active CN113376692B (zh) | 2020-03-10 | 2020-03-10 | 致密砂岩气水平井压裂改造方案优化方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113376692B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108398720A (zh) * | 2018-03-07 | 2018-08-14 | 成都理工大学 | 一种基于杨氏模量、泊松比的两项式地震叠前反演方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104500017A (zh) * | 2014-12-12 | 2015-04-08 | 中国石油天然气集团公司 | 一种优化水平井分段压裂位置的方法 |
CN104502970B (zh) * | 2014-12-31 | 2017-09-01 | 中国石油天然气集团公司 | 一种页岩气储层压裂区域选取方法 |
US10310136B2 (en) * | 2015-04-24 | 2019-06-04 | W.D. Von Gonten Laboratories Inc. | Lateral placement and completion design for improved well performance of unconventional reservoirs |
WO2017020267A1 (zh) * | 2015-08-05 | 2017-02-09 | 深圳朝伟达科技有限公司 | 一种水力压裂微地震事件的解释方法 |
-
2020
- 2020-03-10 CN CN202010161440.0A patent/CN113376692B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108398720A (zh) * | 2018-03-07 | 2018-08-14 | 成都理工大学 | 一种基于杨氏模量、泊松比的两项式地震叠前反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113376692A (zh) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
Radwan | Three-dimensional gas property geological modeling and simulation | |
RU2661489C1 (ru) | Способ комплексирования исходных данных для уточнения фильтрационного строения неоднородных карбонатных коллекторов | |
CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
WO2017027433A1 (en) | Method of performing integrated fracture and reservoir operations for multiple wellbores at a wellsite | |
US10359529B2 (en) | Singularity spectrum analysis of microseismic data | |
EP3253946B1 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
Aarre et al. | Seismic detection of subtle faults and fractures | |
Harris | The role of geology in reservoir simulation studies | |
CN105089615A (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
CN112505784B (zh) | 一种煤矿井下盾构机掘进路线选择方法 | |
AlRassas et al. | CO2 storage capacity estimation under geological uncertainty using 3-D geological modeling of unconventional reservoir rocks in Shahejie Formation, block Nv32, China | |
Gu et al. | Shale fracturing characterization and optimization by using anisotropic acoustic interpretation, 3D fracture modeling, and supervised machine learning | |
Trippetta et al. | Carbonate-ramp reservoirs modelling best solutions: Insights from a dense shallow well database in Central Italy | |
CN107831540A (zh) | 储层物性参数直接提取新方法 | |
CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
CN106199712A (zh) | 一种确定压裂套管变形区域的方法及装置 | |
CN115857047B (zh) | 一种地震储层综合预测方法 | |
CN113376692B (zh) | 致密砂岩气水平井压裂改造方案优化方法及装置 | |
CN112282751B (zh) | 一种地质工程三维耦合的致密油气水平井开采检测方法 | |
Soleimani et al. | Permeability variation modeling and reservoir heterogeneity of Bangestan carbonate sequence, Mansouri oilfield, SW Iran | |
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 | |
RU2145101C1 (ru) | Способ оценки эксплуатационных свойств нефтегазовой залежи | |
Xiao et al. | Fine calculation method for fracture morphology based on microseismic events and its application | |
Li et al. | Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China |
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 |