CN112302640A - 水平井压裂工程风险预测方法及装置 - Google Patents
水平井压裂工程风险预测方法及装置 Download PDFInfo
- Publication number
- CN112302640A CN112302640A CN202011172344.2A CN202011172344A CN112302640A CN 112302640 A CN112302640 A CN 112302640A CN 202011172344 A CN202011172344 A CN 202011172344A CN 112302640 A CN112302640 A CN 112302640A
- Authority
- CN
- China
- Prior art keywords
- fracture
- index
- azimuth
- dimensional seismic
- crack
- 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 36
- 230000000694 effects Effects 0.000 claims abstract description 46
- 238000005553 drilling Methods 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims description 39
- 101000934489 Homo sapiens Nucleosome-remodeling factor subunit BPTF Proteins 0.000 claims description 20
- 101100326803 Neurospora crassa (strain ATCC 24698 / 74-OR23-1A / CBS 708.71 / DSM 1257 / FGSC 987) fac-2 gene Proteins 0.000 claims description 20
- 102100025062 Nucleosome-remodeling factor subunit BPTF Human genes 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 12
- 101100066315 Trypanosoma brucei brucei (strain 927/4 GUTat10.1) FAZ1 gene Proteins 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 18
- 230000035882 stress Effects 0.000 description 17
- 238000004458 analytical method Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism 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
- 230000003334 potential effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
本发明公开了一种水平井压裂工程风险预测方法及装置,其中该方法包括:获取观测区的钻井轨迹方位、区域应力方位和三维地震数据体;根据三维地震数据体计算裂缝指数和裂缝方位;根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数;根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数;根据裂缝影响指数和裂缝活动指数计算工程风险指数,所述工程风险指数用于预测水平井压裂工程风险。本发明应用区域应力评价井旁裂缝活动指数,应用裂缝方位评价裂缝对套管的影响指数,应用裂缝活动指数和影响指数综合评价压裂工程风险,以降低或减少压裂工程风险,降低压裂工程成本,提高单井产量。
Description
技术领域
本发明涉及井中地球物理解释技术领域,尤其涉及水平井压裂工程风险预测方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
水平井在钻井、压裂等工程实施过程中,常常会发生套管变形等风险,严重影响工程实施的进度,增加工程成本。压裂是一种有效的储层改造手段,但常常伴随工程风险,对于工程风险的防治,首先需要预测工程风险位置,评价风险程度。因此,有效的进行水平井压裂工程风险的预测意义重大。但是,目前对压裂工程风险的研究主要集中在机理研究和风险的防治方法研究,尚未见有关应用地球物理技术预测压裂工程风险的公开文献或报道,也未公开有关应用地球物理技术预测压裂工程风险的详细描述和具体细节。
发明内容
本发明实施例提供一种水平井压裂工程风险预测方法,用以,该方法包括:
获取观测区的钻井轨迹方位、区域应力方位和三维地震数据体;
根据三维地震数据体计算裂缝指数和裂缝方位;
根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数;
根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数;
根据裂缝影响指数和裂缝活动指数计算工程风险指数,所述工程风险指数用于评价水平井压裂工程风险。
本发明实施例还提供一种水平井压裂工程风险预测装置,用以,该装置包括:
数据采集模块,用于获取观测区的钻井轨迹方位、区域应力方位和三维地震数据体;
裂缝方位与裂缝指数计算模块,用于根据三维地震数据体计算裂缝指数和裂缝方位;
裂缝影响指数计算模块,用于根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数;
裂缝活动指数计算模块,用于根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数;
工程风险指数计算模块,用于根据裂缝影响指数和裂缝活动指数计算工程风险指数,所述工程风险指数用于评价水平井压裂工程风险。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述水平井压裂工程风险预测方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述水平井压裂工程风险预测方法的计算机程序。
本发明实施例中,应用区域应力评价井旁裂缝活动指数,应用裂缝方位评价裂缝对套管的影响指数,应用裂缝活动指数和影响指数综合评价压裂工程风险,以降低或减少压裂工程风险,降低压裂工程成本,提高单井产量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中水平井压裂工程风险预测方法流程图;
图2为本发明实施例中水平井压裂工程风险预测方法具体流程图;
图3为本发明实施例中裂缝强度示意图;
图4为本发明实施例中裂缝强度与裂缝方位示意图;
图5为本发明实施例中方位权重系数和权重指数示意图;
图6为本发明实施例中裂缝计算的邻域和相对方位示意图;
图7为本发明实施例中方位加权前后裂缝强度示意图;
图8为本发明实施例中裂缝强度和裂缝指数;
图9为本发明实施例中压裂工程风险指数评价图
图10为本发明实施例中水平井压裂工程风险预测装置结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
本发明提供了一种水平井压裂工程风险预测方法,采用双重影响分析,包括应力对裂缝的影响分析:即裂缝活动指数FAI分析,以及裂缝对钻井的影响分析:即裂缝影响指数FII分析,评价水平压裂工程潜在的风险,降低了风险评价的多解性,使风险评价更精准。
图1为本发明实施例中水平井压裂工程风险预测方法流程图,如图1所示,该方法包括:
步骤101:获取观测区的钻井轨迹方位WAZ、区域应力方位SAZ和三维地震数据体;
步骤102:根据三维地震数据体计算裂缝指数FAC和裂缝方位FAZ;
步骤103:根据钻井轨迹方位、裂缝指数FAC和裂缝方位FAZ,计算裂缝影响指数FII,评价裂缝对钻井轨迹的影响;
步骤104:根据区域应力方位、裂缝指数FAC和裂缝方位FAZ,计算裂缝活动指数FAI,评价裂缝的潜在活动性;
步骤105:根据裂缝影响指数FII和裂缝活动指数FAI计算工程风险指数FRI,所述工程风险指数用于评价水平井压裂工程风险。
在本发明实施例中,所述三维地震数据体包括三维地震叠加数据体和三维地震横波数据体,还包括三维地震速度比或泊松比数据体;
如图2所示,根据三维地震数据体计算裂缝指数和裂缝方位,包括:
根据三维地震叠加数据体计算第一个体裂缝指数FAC1和第一裂缝方位FAZ1;
根据三维地震横波数据体计算第二个体裂缝指数FAC2和第二裂缝方位FAZ2;
根据三维地震速度比或泊松比数据体计算第三个体裂缝指数FAC3和第三裂缝方位FAZ3;
将第一个体裂缝指数、第二个体裂缝指数和第三个体裂缝指数融合,获得全局裂缝指数FAC;
将第一裂缝方位、第二裂缝方位和第三裂缝方位融合,获得全局裂缝方位FAZ。
步骤103和104利用的都是全局裂缝指数FAC和全局裂缝方位FAZ。
在本发明实施例中,步骤102:根据三维地震数据体计算裂缝指数FACi和裂缝方位FAZi(i=1,2,3),包括:
S1:根据三维地震数据体计算曲率CURi;即在某个待计算点的邻域内,应用三维地震数据体的幅值,拟合三维曲面,计算三维曲面的曲度即为该点的曲率CURi。
S2:根据曲率CURi计算裂缝方位FAZi。即在某个待计算点的邻域内,应用三维空间的曲率CURi,拟合三维曲面,并计算三维曲面的方位梯度,梯度最小的方向即为裂缝方位FAZi。
S3:根据裂缝方位FAZi对曲率CURi进行加权处理,计算裂缝强度FSTi;采用如下关系式计算:
其中,Wj为加权系数;Q为加权指数,Q∈(0,4];加权系数和加权指数的应用使得裂缝强度在裂缝方位上分辨率更高,信噪比更高。加权指数越小,加权系数随着方位变化越小;加权指数越大,加权系数越偏重于裂缝方位,即沿裂缝方位加权系数较大,其它方位加权系数逐渐减小,垂直于裂缝方位的加权系数为0。当Q=4时,裂缝方位10度范围以外的加权系数小于0.5。
FSTip表示第p点的裂缝强度,表示第j点曲率,k为指数,值越小越彰显细节,建议取值0.25到3.00之间;FAZip表示第p点的裂缝方位,其中,i=1,2,3,分别表示根据三维地震叠加数据体和三维地震横波数据体、三维地震速度比或泊松比数据体计算得到的裂缝强度、曲率、裂缝方位;DAZjp表示第j点到第p点的空间方位;s是第p点的邻域,第j点位于s内;angj表示观测点到计算点j方位相对于裂缝方位的方位角。
如图3展示了裂缝强度,颜色越深代表值越大,即裂缝强度越大;图4展示了裂缝方位,图中直线的朝向即为裂缝方位。图5展示了不同加权指数下的方位加权系数。图中W为加权系数,ANG为观测点到计算点方位相对于裂缝方位的方位角,Q为加权指数。图6展示了空间一点的邻域及邻域内某点相对于该点的方位示意图。图7展示了采用不同加权指数,加权处理前后的裂缝强度示意图,其中,第一幅为加权处理前的裂缝强度原始图,第2、3和4是采用加权指数Q=0.2、0.8、0.4处理后的裂缝强度示意图,裂缝信噪比获得提高。
S4:根据裂缝强度计算裂缝指数。采用如下关系式:
其中,FACi表示裂缝指数,λ表示阈值,可根据实际设置。
所述的裂缝指数取值范围为0和1,当数值为0是,指示无裂缝;当数值为1是,表示裂缝发育。
所述的裂缝方位取值范围为0和360度。
图8展示了裂缝强度与裂缝指数,表示裂缝指数的图中灰度深黑色的为0,灰度浅接近于白色为1。
在本发明实施例中,采用如下关系式根据三维地震数据体计算裂缝指数和裂缝方位:
FAC=FAC1|FAC2|FAC3;
其中,|竖线表示或运算,当FAC1、FAC2、FAC3中任意一个数值为1时,FAC为1,否则为0;FAC表示全局裂缝指数;
当FAC=0时,FAZ=0;当FAC≠0&&FAC1≠0时,FAZ=FAZ1;当FAC≠0&&FAC1=0时,FAZ=(FAZ2×FAC2+FAZ3×FAC3)/(FAC2+FAC3);其中,&&表示并且的意思;FAZ表示全局裂缝方位。
在本发明实施例中,采用如下关系式根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数:
FII=FAC×{sin(asin(|sin(WAZ-FAZ)|))};
其中,FII表示裂缝影响指数;FAZ表示全局裂缝方位;FAC表示全局裂缝指数;WAZ表示钻井轨迹方位,asin()是sin()的反函数。
FII值域从0到1,当FII=0时,表示裂缝对水平井的压裂无影响。FII值越大,表示裂缝对水平井的压裂影响程度越大。WAZ为水平井井轨迹方位。
在本发明实施例中,采用如下关系式根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数:
FAI=FAC×{cos(|asin(|sin(SAZ-FAZ)|)-(π/4-φ/2)|)-cos(π/4-φ/2)};
其中,FAI表示裂缝活动指数;FAZ表示全局裂缝方位;FAC表示全局裂缝指数;SAZ表示区域应力方位。asin()是sin()的反函数。
FAI值域从-1到+1,当FAI<=0时,表示裂缝不活动,当FAI>0时,表示裂缝活动,FAI值越大,裂缝活动强度越大。
SAZ为区域应力方位,可通过(但不限于)区域地质分析等方式获得;φ为裂缝内摩擦角,可通过(但不限于)岩芯测试等方式获得。
在本发明实施例中,采用如下关系式根据裂缝影响指数和裂缝活动指数计算工程风险指数:
其中,FRI表示工程风险指数;FAI表示裂缝活动指数;FII表示裂缝影响指数。
FRI的取值范围为[0,1]和[10,11],个位和小数位表征风险程度,十位数表征裂缝活动程度。当十位数为0时,表征裂缝活动程度低于影响程度;当十位数为1时,表征裂缝活动程度高于影响程度。
假设区域应力方位为110度,水平井轨迹方位为20度,储层的内摩擦角为30度,如图9展示了风险指数示意图。
本发明实施例中还提供了一种水平井压裂工程风险预测装置,如下面的实施例所述。由于该装置解决问题的原理与水平井压裂工程风险预测方法相似,因此该装置的实施可以参见水平井压裂工程风险预测方法的实施,重复之处不再赘述。
图10为本发明实施例中水平井压裂工程风险预测装置结构框图,如图10所示,该装置包括:
数据采集模块02,用于获取观测区的钻井轨迹方位、区域应力方位和三维地震数据体;
裂缝方位与裂缝指数计算模块04,用于根据三维地震数据体计算裂缝指数和裂缝方位;
裂缝影响指数计算模块06,用于根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数;
裂缝活动指数计算模块08,用于根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数;
工程风险指数计算模块10,用于根据裂缝影响指数和裂缝活动指数计算工程风险指数,所述工程风险指数用于评价水平井压裂工程风险。
在本发明实施例中,所述三维地震数据体包括三维地震叠加数据体和三维地震横波数据体,还包括三维地震速度比或泊松比数据体;
裂缝方位与裂缝指数计算模块具体用于:
根据三维地震叠加数据体计算第一个体裂缝指数和第一裂缝方位;
根据三维地震横波数据体计算第二个体裂缝指数和第二裂缝方位;
根据三维地震速度比或泊松比数据体计算第三个体裂缝指数和第三裂缝方位;
将第一个体裂缝指数、第二个体裂缝指数和第三个体裂缝指数融合,获得全局裂缝指数;
将第一裂缝方位、第二裂缝方位和第三裂缝方位融合,获得全局裂缝方位。
在本发明实施例中,裂缝方位与裂缝指数计算模块具体用于:
按照如下方式根据三维地震数据体计算裂缝指数和裂缝方位:
根据三维地震数据体计算曲率;
根据曲率计算裂缝方位;
根据裂缝方位对曲率进行加权处理,计算裂缝强度;
根据裂缝强度计算裂缝指数。
在本发明实施例中,裂缝方位与裂缝指数计算模块具体用于:
采用如下关系式根据裂缝方位对曲率进行加权处理,计算裂缝强度:
其中,Wj为加权系数;Q为加权指数,Q∈(0,4];FSTip表示第p点的裂缝强度,表示第j点曲率,k为指数,值越小越彰显细节,建议取值0.25到3.00之间;FAZip表示第p点的裂缝方位,其中,i=1,2,3,分别表示根据三维地震叠加数据体和三维地震横波数据体、三维地震速度比或泊松比数据体计算得到的裂缝强度、曲率、裂缝方位;DAZjp表示第j点到第p点的空间方位;s是第p点的邻域,第j点位于s内;angj表示观测点到计算点j方位相对于裂缝方位的方位角。
在本发明实施例中,裂缝方位与裂缝指数计算模块具体用于:
采用如下关系式根据裂缝强度计算裂缝指数:
其中,FACi表示裂缝指数,λ表示阈值;其中,i=1,2,3,分别表示根据三维地震叠加数据体和三维地震横波数据体、三维地震速度比或泊松比数据体计算得到的裂缝指数。
在本发明实施例中,裂缝方位与裂缝指数计算模块具体用于:
采用如下关系式根据三维地震数据体计算裂缝指数和裂缝方位:
FAC=FAC1|FAC2|FAC3;
其中,|竖线表示或运算,当FAC1、FAC2、FAC3中任意一个数值为1时,FAC为1,否则为0;FAC表示全局裂缝指数;FAC1、FAC2、FAC3分别表示根据三维地震叠加数据体和三维地震横波数据体、三维地震速度比或泊松比数据体计算得到的裂缝指数;
当FAC=0时,FAZ=0;当FAC≠0&&FAC1≠0时,FAZ=FAZ1;当FAC≠0&&FAC1=0时,FAZ=(FAZ2×FAC2+FAZ3×FAC3)/(FAC2+FAC3);其中,&&表示并且的意思;FAZ表示全局裂缝方位;FAZ1、FAZ2、FAZ3分别表示根据三维地震叠加数据体和三维地震横波数据体、三维地震速度比或泊松比数据体计算得到的裂缝方位。
在本发明实施例中,裂缝影响指数计算模块具体用于:
采用如下关系式根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数:
FII=FAC×{sin(asin(|sin(WAZ-FAZ)|))};
其中,FII表示裂缝影响指数;FAZ表示全局裂缝方位;FAC表示全局裂缝指数;WAZ表示钻井轨迹方位,asin()是sin()的反函数。
在本发明实施例中,裂缝活动指数计算模块具体用于:
采用如下关系式根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数:
FAI=FAC×{cos(|asin(|sin(SAZ-FAZ)|)-(π/4-φ/2)|)-cos(π/4-φ/2)};
其中,FAI表示裂缝活动指数;FAZ表示全局裂缝方位;FAC表示全局裂缝指数;SAZ表示区域应力方位,asin()是sin()的反函数。
在本发明实施例中,工程风险指数计算模块具体用于:
采用如下关系式根据裂缝影响指数和裂缝活动指数计算工程风险指数:
其中,FRI表示工程风险指数;FAI表示裂缝活动指数;FII表示裂缝影响指数。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述水平井压裂工程风险预测方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述水平井压裂工程风险预测方法的计算机程序。
本发明实施例中,提供的一种水平井压裂工程风险预测方法,采用双重影响分析,即应力对裂缝的影响(裂缝活动指数FAI)和裂缝对钻井的影响(裂缝影响指数FII),评价水平井压裂工程潜在的风险,降低了多解性,使风险评价更精确。以期降低或减少压裂工程风险,提高单井压裂工程时效,降低成本,提高单井产量。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (20)
1.一种水平井压裂工程风险预测方法,其特征在于,包括:
获取观测区的钻井轨迹方位、区域应力方位和三维地震数据体;
根据三维地震数据体计算裂缝指数和裂缝方位;
根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数;
根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数;
根据裂缝影响指数和裂缝活动指数计算工程风险指数,所述工程风险指数用于评价水平井压裂工程风险。
2.如权利要求1所述的水平井压裂工程风险预测方法,其特征在于,所述三维地震数据体包括三维地震叠加数据体和三维地震横波数据体,还包括三维地震速度比或泊松比数据体;
根据三维地震数据体计算裂缝指数和裂缝方位,包括:
根据三维地震叠加数据体计算第一个体裂缝指数和第一裂缝方位;
根据三维地震横波数据体计算第二个体裂缝指数和第二裂缝方位;
根据三维地震速度比或泊松比数据体计算第三个体裂缝指数和第三裂缝方位;
将第一个体裂缝指数、第二个体裂缝指数和第三个体裂缝指数融合,获得全局裂缝指数;
将第一裂缝方位、第二裂缝方位和第三裂缝方位融合,获得全局裂缝方位。
3.如权利要求2所述的水平井压裂工程风险预测方法,其特征在于,根据三维地震数据体计算裂缝指数和裂缝方位,包括:
根据三维地震数据体计算曲率;
根据曲率计算裂缝方位;
根据裂缝方位对曲率进行加权处理,计算裂缝强度;
根据裂缝强度计算裂缝指数。
6.如权利要求2所述的水平井压裂工程风险预测方法,其特征在于,采用如下关系式根据三维地震数据体计算裂缝指数和裂缝方位:
FAC=FAC1|FAC2|FAC3;
其中,|竖线表示或运算,当FAC1、FAC2、FAC3中任意一个数值为1时,FAC为1,否则为0;FAC表示全局裂缝指数;FAC1、FAC2、FAC3分别表示根据三维地震叠加数据体和三维地震横波数据体、三维地震速度比或泊松比数据体计算得到的裂缝指数;
当FAC=0时,FAZ=0;当FAC≠0&&FAC1≠0时,FAZ=FAZ1;当FAC≠0&&FAC1=0时,FAZ=(FAZ2×FAC2+FAZ3×FAC3)/(FAC2+FAC3);其中,&&表示并且的意思;FAZ表示全局裂缝方位;FAZ1、FAZ2、FAZ3分别表示根据三维地震叠加数据体和三维地震横波数据体、三维地震速度比或泊松比数据体计算得到的裂缝方位。
7.如权利要求1所述的水平井压裂工程风险预测方法,其特征在于,采用如下关系式根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数:
FII=FAC×{sin(asin(|sin(WAZ-FAZ)|))};
其中,FII表示裂缝影响指数;FAZ表示全局裂缝方位;FAC表示全局裂缝指数;WAZ表示钻井轨迹方位,asin()是sin()的反函数。
8.如权利要求1所述的水平井压裂工程风险预测方法,其特征在于,采用如下关系式根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数:
FAI=FAC×{cos(|asin(|sin(SAZ-FAZ)|)-(π/4-φ/2)|)-cos(π/4-φ/2)};
其中,FAI表示裂缝活动指数;FAZ表示全局裂缝方位;FAC表示全局裂缝指数;SAZ表示区域应力方位;φ为裂缝内摩擦角,asin()是sin()的反函数。
10.一种水平井压裂工程风险预测装置,其特征在于,包括:
数据采集模块,用于获取观测区的钻井轨迹方位、区域应力方位和三维地震数据体;
裂缝方位与裂缝指数计算模块,用于根据三维地震数据体计算裂缝指数和裂缝方位;
裂缝影响指数计算模块,用于根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数;
裂缝活动指数计算模块,用于根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数;
工程风险指数计算模块,用于根据裂缝影响指数和裂缝活动指数计算工程风险指数,所述工程风险指数用于评价水平井压裂工程风险。
11.如权利要求10所述的水平井压裂工程风险预测装置,其特征在于,所述三维地震数据体包括三维地震叠加数据体和三维地震横波数据体,还包括三维地震速度比或泊松比数据体;
裂缝方位与裂缝指数计算模块具体用于:
根据三维地震叠加数据体计算第一个体裂缝指数和第一裂缝方位;
根据三维地震横波数据体计算第二个体裂缝指数和第二裂缝方位;
根据三维地震速度比或泊松比数据体计算第三个体裂缝指数和第三裂缝方位;
将第一个体裂缝指数、第二个体裂缝指数和第三个体裂缝指数融合,获得全局裂缝指数;
将第一裂缝方位、第二裂缝方位和第三裂缝方位融合,获得全局裂缝方位。
12.如权利要求11所述的水平井压裂工程风险预测装置,其特征在于,裂缝方位与裂缝指数计算模块具体用于:
按照如下方式根据三维地震数据体计算裂缝指数和裂缝方位:
根据三维地震数据体计算曲率;
根据曲率计算裂缝方位;
根据裂缝方位对曲率进行加权处理,计算裂缝强度;
根据裂缝强度计算裂缝指数。
15.如权利要求11所述的水平井压裂工程风险预测装置,其特征在于,裂缝方位与裂缝指数计算模块具体用于:
采用如下关系式根据三维地震数据体计算裂缝指数和裂缝方位:
FAC=FAC1|FAC2|FAC3;
其中,|竖线表示或运算,当FAC1、FAC2、FAC3中任意一个数值为1时,FAC为1,否则为0;FAC表示全局裂缝指数;FAC1、FAC2、FAC3分别表示根据三维地震叠加数据体和三维地震横波数据体、三维地震速度比或泊松比数据体计算得到的裂缝指数;
当FAC=0时,FAZ=0;当FAC≠0&&FAC1≠0时,FAZ=FAZ1;当FAC≠0&&FAC1=0时,FAZ=(FAZ2×FAC2+FAZ3×FAC3)/(FAC2+FAC3);其中,&&表示并且的意思;FAZ表示全局裂缝方位;FAZ1、FAZ2、FAZ3分别表示根据三维地震叠加数据体和三维地震横波数据体、三维地震速度比或泊松比数据体计算得到的裂缝方位。
16.如权利要求10所述的水平井压裂工程风险预测装置,其特征在于,裂缝影响指数计算模块具体用于:
采用如下关系式根据钻井轨迹方位、裂缝指数和裂缝方位,计算裂缝影响指数:
FII=FAC×{sin(asin(|sin(WAZ-FAZ)|))};
其中,FII表示裂缝影响指数;FAZ表示全局裂缝方位;FAC表示全局裂缝指数;WAZ表示钻井轨迹方位,asin()是sin()的反函数。
17.如权利要求10所述的水平井压裂工程风险预测装置,其特征在于,裂缝活动指数计算模块具体用于:
采用如下关系式根据区域应力方位、裂缝指数和裂缝方位,计算裂缝活动指数:
FAI=FAC×{cos(|asin(|sin(SAZ-FAZ)|)-(π/4-φ/2)|)-cos(π/4-φ/2)};
其中,FAI表示裂缝活动指数;FAZ表示全局裂缝方位;FAC表示全局裂缝指数;SAZ表示区域应力方位;φ为裂缝内摩擦角,asin()是sin()的反函数。
19.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至9任一所述方法。
20.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至9任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011172344.2A CN112302640A (zh) | 2020-10-28 | 2020-10-28 | 水平井压裂工程风险预测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011172344.2A CN112302640A (zh) | 2020-10-28 | 2020-10-28 | 水平井压裂工程风险预测方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112302640A true CN112302640A (zh) | 2021-02-02 |
Family
ID=74331457
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011172344.2A Pending CN112302640A (zh) | 2020-10-28 | 2020-10-28 | 水平井压裂工程风险预测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112302640A (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080183393A1 (en) * | 2007-01-25 | 2008-07-31 | Rundle John B | Method and system for forecasting earthquakes |
US20160069182A1 (en) * | 2014-09-10 | 2016-03-10 | Fracture ID, Inc. | Apparatus and method using measurements taken while drilling to map mechanical boundaries and mechanical rock properties along a borehole |
CN105952427A (zh) * | 2016-05-04 | 2016-09-21 | 中国石油大学(北京) | 一种低渗透油藏注水诱导裂缝的预测与评价方法 |
WO2017020267A1 (zh) * | 2015-08-05 | 2017-02-09 | 深圳朝伟达科技有限公司 | 一种水力压裂微地震事件的解释方法 |
CN107942400A (zh) * | 2017-11-10 | 2018-04-20 | 中国石油天然气股份有限公司 | 一种预测钻井风险的方法及装置 |
CN109543338A (zh) * | 2019-01-04 | 2019-03-29 | 陈建永 | 一种地震波作用下基坑稳定性安全系数的计算方法 |
WO2020048028A1 (zh) * | 2018-09-03 | 2020-03-12 | 中国石油大学(华东) | 一种基于压裂潜力的待压裂水平井压裂设计方法和装置 |
CN111257946A (zh) * | 2018-11-30 | 2020-06-09 | 中国石油化工股份有限公司 | 一种地球物理指导钻井方法及更新地层地震速度的方法 |
-
2020
- 2020-10-28 CN CN202011172344.2A patent/CN112302640A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080183393A1 (en) * | 2007-01-25 | 2008-07-31 | Rundle John B | Method and system for forecasting earthquakes |
US20160069182A1 (en) * | 2014-09-10 | 2016-03-10 | Fracture ID, Inc. | Apparatus and method using measurements taken while drilling to map mechanical boundaries and mechanical rock properties along a borehole |
WO2017020267A1 (zh) * | 2015-08-05 | 2017-02-09 | 深圳朝伟达科技有限公司 | 一种水力压裂微地震事件的解释方法 |
CN105952427A (zh) * | 2016-05-04 | 2016-09-21 | 中国石油大学(北京) | 一种低渗透油藏注水诱导裂缝的预测与评价方法 |
CN107942400A (zh) * | 2017-11-10 | 2018-04-20 | 中国石油天然气股份有限公司 | 一种预测钻井风险的方法及装置 |
WO2020048028A1 (zh) * | 2018-09-03 | 2020-03-12 | 中国石油大学(华东) | 一种基于压裂潜力的待压裂水平井压裂设计方法和装置 |
CN111257946A (zh) * | 2018-11-30 | 2020-06-09 | 中国石油化工股份有限公司 | 一种地球物理指导钻井方法及更新地层地震速度的方法 |
CN109543338A (zh) * | 2019-01-04 | 2019-03-29 | 陈建永 | 一种地震波作用下基坑稳定性安全系数的计算方法 |
Non-Patent Citations (1)
Title |
---|
储仿东 等: "《利用维地震识别水力压裂井旁天然裂缝》", 《石油地球物理勘探》, vol. 53, no. 2, 30 December 2018 (2018-12-30), pages 143 - 147 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Stork et al. | Application of machine learning to microseismic event detection in distributed acoustic sensing data | |
Sinha et al. | Spectral decomposition of seismic data with continuous-wavelet transform | |
Rebetsky et al. | Development of inversion methods on fault slip data: Stress state in orogenes of the Central Asia | |
Jin et al. | An accurate and adaptative cutterhead torque prediction method for shield tunneling machines via adaptative residual long-short term memory network | |
US9261613B2 (en) | Passive monitoring method for seismic events | |
US20150316673A1 (en) | Systems and Methods for 3D Seismic Data Depth Conversion Utilizing Artificial Neural Networks | |
US20170108604A1 (en) | Denoising seismic data | |
Riedel et al. | Underground vertical seismic profiling with conventional and fiber-optic systems for exploration in the Kylylahti Polymetallic Mine, Eastern Finland | |
Cui et al. | VP and VS prediction from digital rock images using a combination of U-Net and convolutional neural networks | |
Emery et al. | Truncated Gaussian simulation of discrete-valued, ordinal coregionalized variables | |
Cortina et al. | Neural networks applied to the wave-induced fatigue analysis of steel risers | |
Fathi et al. | High-quality fracture network mapping using high frequency logging while drilling (LWD) data: MSEEL case study | |
Wojtecki et al. | An attempt to use machine learning algorithms to estimate the rockburst hazard in underground excavations of hard coal mine | |
Liu et al. | Convolutional neural network‐based classification of microseismic events originating in a stimulated reservoir from distributed acoustic sensing data | |
Riedel et al. | A hybrid exploration approach for the prediction of geological changes ahead of mechanized tunnel excavation | |
CN112302640A (zh) | 水平井压裂工程风险预测方法及装置 | |
Li et al. | Reconstruction of critical coalbed methane logs with principal component regression model: A case study | |
Maleki et al. | Comparison of several different methods of in situ stress determination | |
Wang et al. | Attribute-guided target data separation network for DAS VSP data | |
CA2780539A1 (en) | System and method for analyzing and transforming geophysical and petrophysical data | |
Ortega et al. | Automatic selection of dispersion curves based on a weighted probability scheme | |
Yang et al. | SLKNet: An attention-based deep-learning framework for downhole distributed acoustic sensing data denoising | |
Panetta et al. | 3D visualization of integrated geologic and geophysical subsurface data using open-source programming: A case study using data from the MSEEL project | |
Datta et al. | Application of multivariate change detection in automated lithofacies classification from well-log data in a nonstationary subsurface | |
Mahmoudi et al. | A pattern recognition approach to evaluate the uncertain soil stratum in front of the tunnel face |
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 |