CN116877039B - 一种燃爆压裂缝数量和长度预测模型的构建方法及设备 - Google Patents
一种燃爆压裂缝数量和长度预测模型的构建方法及设备 Download PDFInfo
- Publication number
- CN116877039B CN116877039B CN202310839336.6A CN202310839336A CN116877039B CN 116877039 B CN116877039 B CN 116877039B CN 202310839336 A CN202310839336 A CN 202310839336A CN 116877039 B CN116877039 B CN 116877039B
- Authority
- CN
- China
- Prior art keywords
- gas
- pressure
- crack
- length
- rock
- 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 58
- 238000005422 blasting Methods 0.000 title claims description 15
- 239000011435 rock Substances 0.000 claims abstract description 44
- 239000011148 porous material Substances 0.000 claims abstract description 23
- 238000004880 explosion Methods 0.000 claims abstract description 21
- 238000002474 experimental method Methods 0.000 claims abstract description 17
- 239000000463 material Substances 0.000 claims abstract description 17
- 230000000977 initiatory effect Effects 0.000 claims abstract description 11
- 208000010392 Bone Fractures Diseases 0.000 claims description 24
- 206010017076 Fracture Diseases 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 23
- 230000035515 penetration Effects 0.000 claims description 13
- 230000015654 memory Effects 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000004891 communication Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 239000000446 fuel Substances 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000012669 compression test Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 230000003746 surface roughness Effects 0.000 claims description 3
- 206010010214 Compression fracture Diseases 0.000 claims description 2
- 230000000644 propagated effect Effects 0.000 claims description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 16
- 238000005516 engineering process Methods 0.000 abstract description 10
- 238000011161 development Methods 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 3
- 239000007789 gas Substances 0.000 description 75
- 238000005336 cracking Methods 0.000 description 6
- 230000003068 static effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000001965 increasing effect Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 238000004200 deflagration Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 230000036962 time dependent Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 239000004568 cement Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000326 densiometry Methods 0.000 description 1
- 238000005474 detonation Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 229910052602 gypsum Inorganic materials 0.000 description 1
- 239000010440 gypsum Substances 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000007800 oxidant agent Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Mathematical Physics (AREA)
- Geology (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Geochemistry & Mineralogy (AREA)
- Algebra (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明属于油气田开发技术领域，并具体公开了一种燃爆压裂缝数量和长度预测模型的构建方法及设备。所述方法包括：获取岩石的力学参数；根据所述力学参数，选择岩石的相似材料作为试样，对试样进行压裂实验，以获取膛压和孔压的升压曲线；根据所述力学参数、膛压和孔压的升压曲线建立裂缝起裂模型，计算裂缝数量和初始长度；基于所述裂缝数量和初始长度构建裂缝扩展模型，计算裂缝扩展长度。本发明建立的燃爆压裂裂缝预测模型力学意义明确，可以准确地预测不同实验条件下产生的裂缝数量和长度，对于甲烷燃爆压裂技术的裂缝控制和参数优化具有重要的理论研究意义和较高的工程应用价值。
Description
技术领域
本发明属于油气田开发技术领域，更具体地，涉及一种燃爆压裂缝数量和长度预测模型的构建方法及设备。
背景技术
尽管水力压裂技术被广泛用于增产技术，但仍存在成本高、环境污染、诱发地震等问题，而且只能产生一条垂直于最小主应力方向的裂缝。由于其强大的瞬态冲击力，爆炸压裂技术可以在钻孔附近快速生成十几条裂缝甚至复杂的裂缝网络。然而，它很容易破坏井筒的完整性并造成安全风险。因而提出利用甲烷原位燃爆来减少水资源消耗、防止地层污染、促进不受应力控制的多条裂缝扩展。
作为一种高效环保的新型无水压裂技术，甲烷原位爆燃压裂技术利用储层中甲烷爆炸反应产生的高温高压气体和氧气等氧化剂冲击页岩储层，形成裂缝网络，为页岩气提供高效的输送通道，最终达到增产的目的。但其作为页岩气储层改造的新方法，在开发模式及参数优化方面仍有很多困难和挑战，需进一步的探究。
在油气开发领域，储层产生裂缝的数量和长度是影响油气井增产效果的关键因素。然而影响裂缝数量和长度的因素较多，包括储层地质静态参数、燃爆压裂施工参数、生产动态参数等方面因素，各参数之间的关系难以用单一表达式表示。通常情况下采用数值模拟的方法(有限元、离散元、扩展有限元、连续非连续单元法等)计算裂缝数量和长度，但存在数值模拟时间长、裂缝参数描述不精确等问题。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种燃爆压裂缝数量和长度预测模型的构建方法及设备，其根据岩石的力学参数进行压裂实验测得膛压曲线和孔压曲线，并据此建立裂缝起裂模型计算裂缝数量和初始裂缝长度，最后建立裂缝扩展模型计算裂缝长度，因此本发明建立的燃爆压裂裂缝预测模型力学意义明确，可以准确地预测不同实验条件下产生的裂缝数量和长度，对于甲烷燃爆压裂技术的裂缝控制和参数优化具有重要的理论研究意义和较高的工程应用价值。
为实现上述目的，按照本发明的一个方面，提出了一种燃爆压裂缝数量和长度预测模型的构建方法及设备，包括以下步骤：
S100获取岩石的力学参数；
S200根据所述力学参数，选择岩石的相似材料作为试样，对试样进行压裂实验，以获取膛压和孔压的升压曲线；
S300根据所述力学参数、膛压和孔压的升压曲线建立裂缝起裂模型，计算裂缝数量和初始长度；
S400基于所述裂缝数量和初始长度构建裂缝扩展模型，计算裂缝扩展长度。
作为进一步优选的，步骤S100包括以下步骤：
进行多种室内实验获取岩石的力学参数；
具体的，首先测得岩石密度；然后进行单轴压缩试验测得岩石的抗压强度、弹性模量和泊松比；最后进行巴西劈裂试验测得岩石的抗拉强度。
作为进一步优选的，步骤S200包括以下步骤：
S201根据岩石的力学参数选择岩石的相似材料，并浇筑试样；
S202对试样进行燃爆压裂实验并测量致裂管内的膛压曲线和作用在孔壁上的孔压曲线；
S203根据膛压和孔压的升压曲线分别得到膛压和孔压的峰值压力和升压速率。
作为进一步优选的，步骤S300包括以下步骤：
S301根据所述力学参数、膛压和孔压的升压曲线计算燃爆生成的气体密度，并基于裂缝起裂模型计算裂缝的初始长度和宽度；
S302计算流经裂缝入口的气体质量；
S303根据气体质量守恒计算裂缝数量。
作为进一步优选的，所述气体密度为：
ρ＝PM/RT
式中，T是峰值温度，P是膛压峰值压力，M是气体的分子质量，R是理想气体常数；
所述裂缝的初始长度为：
式中，σtd是动态拉伸强度；Pp是孔压峰值压力，rw是炮孔半径，α是应力波衰减指数，b＝1-β2，μd＝0.8μ，μ是静态泊松比，β是岩石的横波速度与纵波速度的比值；
裂缝的初始宽度为：
w0＝πrs 2/h，或者
式中，rs是射孔通道的半径，h是试样高度，E为弹性模量；
所述气体质量为：
式中，Cp材料纵波速度，αs是比例系数，γ是孔压升压速率，σ是地应
所述裂缝数量为：
作为进一步优选的，步骤S400包括以下步骤：
S401将得到孔压曲线划分为N个时间步，得到的裂缝数量N0、初始长度L0以及岩石力学参数均为输入数据；
S402假设裂缝尖端与气体尖端不重合，设定气体贯入深度和气体尖端与裂缝入口压力比的初值满足：
0＜θ″(t)＜1,p″(θ”,t)＝0
式中，θ”(t)是气体贯入深度；p”(θ”,t)是气体尖端与裂缝入口压力比；
S403根据应力强度因子判据计算得到气体压力分布指数m；
S404计算得到第i时间步气体密度、裂缝张开位移、气体流动速度、裂缝扩展速度以及气体滤失速度，N≥i＞0；
S405判断气体在裂缝中的流动过程满足质量守恒，若是，进入步骤S406，否则返回至步骤S403，改变气体贯入深度的值；
S406重复步骤S403至步骤S405，进行下一时间步的裂缝扩展计算，直至N个时间步计算完成。
作为进一步优选的，所述步骤S405还包括：若气体在裂缝中的流动过程不满足质量守恒：
且当θ″(t)≤0时仍不满足守恒条件，则设定气体贯入深度和气体尖端与裂缝入口压力比的初值满足：
0＜p″(θ”,t)＜1,θ″(t)＝0
并执行步骤重复步骤S403至步骤S404，判断气体在裂缝中的流动过程满足质量守恒，若是，进入步骤S406，否则，改变p”(θ”,t)的值，返回步骤S403；
若改变p”(θ”,t)的值，且p″(θ,t)≥1时仍不满足守恒条件，则第i时间步裂缝不会扩展，进入步骤S406，以进行下一时间步的裂缝扩展计算。
作为进一步优选的，所述气体压力分布指数m为：
式中，K为应力强度因子，Lf为裂缝长度，p(θ,t)是压力分布函数,σ是地应力,θ是计算点的相对位置，即横坐标值与裂缝长度之比；
步骤S404中，所述气体密度为：
式中，Mg是分子气体质量，pk(θ)和Tk(θ)分别是相对位置θ处的气体压力和温度，R是理想气体常数；
所述裂缝张开位移为：
式中，μ是气体粘度系数，为第k时间步的裂缝长度，θ1和θ2分别是裂缝扩展过程中的瞬时长度和微段长度；
所述气体流动速度为：
式中，ε是表面粗糙度，a，b均为实验测量值，是缝内气体压力梯度；
所述裂缝扩展速度为：
式中，Cr是瑞利波速度，KIC是岩石的断裂韧性，p为缝内气体压力；
所述气体滤失速度为：
式中，β是气体压缩系数，k0是渗透率，Φ是孔隙率，μ是气体粘度系数，L1是气体穿透深度；
步骤S405中判断气体在裂缝中的流动过程质量守恒的模型为：
式中，是气体流入质量流量，/>是压差作用下的气体过滤质量流量，/>和分别是第k和第k-1时间步计算得到的裂缝中的气体质量。
按照本发明的另一个方面，还提供了一种燃爆压裂缝数量和长度预测模型的构建系统，包括：
第一主控模块，用于获取岩石的力学参数；
第二主控模块，用于根据所述力学参数，选择岩石的相似材料作为试样，对试样进行压裂实验，以获取膛压和孔压的升压曲线；
第三主控模块，用于根据所述力学参数、膛压和孔压的升压曲线建立裂缝起裂模型，计算裂缝数量和初始长度；
第四主控模块，用于基于所述裂缝数量和初始长度构建裂缝扩展模型，计算裂缝扩展长度。
按照本发明的另一个方面，还提供了一种电子设备，包括：
至少一个处理器、至少一个存储器和通信接口；其中，
所述处理器、存储器和通信接口相互间进行通信；
所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令，以执行上述任意一种燃爆压裂缝数量和长度预测模型的构建方法。
按照本发明的另一个方面，还提供了一种非暂态计算机可读存储介质，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行上述任意一种燃爆压裂缝数量和长度预测模型的构建方法。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，主要具备以下的技术优点：
1.本发明根据岩石的力学参数进行压裂实验测得膛压曲线和孔压曲线，并据此建立裂缝起裂模型计算裂缝数量和初始裂缝长度，最后建立裂缝扩展模型计算裂缝长度，因此本发明建立的燃爆压裂裂缝预测模型力学意义明确，可以准确地预测不同实验条件下产生的裂缝数量和长度，对于甲烷燃爆压裂技术的裂缝控制和参数优化具有重要的理论研究意义和较高的工程应用价值。
2.本发明的模型参数物理意义明确，可以快速准确地计算出燃爆压裂产生的裂缝数量和长度，对于甲烷燃爆压裂技术的裂缝控制和参数优化具有重要的理论研究意义和较高的工程应用价值。
3.本发明的模型精度和理论严谨性，与试验数据和与以往模型的比较，强于拟合方法，与数学模型相比使用更少的训练样本获得更高的预测精度，与数值模拟方法相比计算时间更少，体现了模型的优越性。
4.本发明方法构建的模型，不但可以应用于甲烷燃爆压裂，还可以应用到其他压裂技术领域(水力裂缝中层内的爆燃技术和高能气体压裂技术等)，体现了模型的适用性。
附图说明
图1是本发明实施例涉及的一种燃爆压裂缝数量和长度预测模型的构建方法的流程图；
图2是本发明实施例涉及的基于裂缝起裂模型的裂缝数量和初始长度计算流程图；
图3是本发明实施例涉及的基于裂缝扩展模型的裂缝长度计算流程图；
图4中的(a)为本发明实施例计算得到的裂缝长度随时间变化的曲线，图4中的(b)为本发明实施例计算得到的压力分布随时间变化曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
下面对本发明第一实施方式的燃爆压裂裂缝数量和长度预测模型的构建方法的实现细节进行说明，以下内容仅为方便理解而提供的实现细节，并非实施本方案的必须。
本实施方式的具体流程图如图1所示，具体包括：
S100进行多种室内实验获取岩石的力学参数。
具体地，多种室内试验包括岩石密度测定试验、单轴压缩试验、巴西劈裂试验。相应地，岩石的力学参数包括密度、抗压强度、抗拉强度、弹性模量、泊松比和剪切模量。
S200选择力学参数相似材料，进行压裂实验获取膛压和孔压的升压曲线。具体地，步骤S200包括以下过程：
S201根据岩石的力学参数选择岩石相似材料，并浇筑试样。
更具体地，根据测得的岩石密度、抗压强度和抗拉强度选择相似材料，浇筑的试样尺寸最小为直径和高度均为1m的圆柱，并在其中心预留贯穿孔模拟井筒。
S202进行燃爆压裂实验并测量致裂管内的膛压曲线和作用在孔壁上的孔压曲线。
更具体地，进行燃爆压裂实验地步骤包括：将确保密封性合格的致裂管装置放入模拟井筒内，检查致裂管顶盖部件的安装，用粗砂和速凝石膏填充模拟井筒，连接高压脉冲点火系统和传感器采集系统，测试和调整传感器和瞬态信号采集仪的采集参数设置，使用区分的管路依次充装甲烷和氧气至指定压力，充装完成后及时撤出充装系统，启动电源进行起爆，保存测得的压力曲线并记录裂缝参数。
更具体地，采集系统中两个压力传感器分别放置在致裂管端部和致裂管与井壁之间测量致裂管内的膛压时程曲线和作用在井壁上的孔压时程曲线。
S203根据膛压和孔压的升压曲线分别得到膛压和孔压的峰值压力和升压速率。
S300结合得到的力学参数和压力参数建立裂缝起裂模型，计算裂缝数量和初始长度。具体地，如图2所示，步骤S300包括以下过程：
S301根据现有参数计算燃爆生成的气体密度，裂缝的初始长度和宽度：
ρ＝PM/RT (1)
其中，T是峰值温度；P是膛压峰值压力；M是气体的分子质量；R是理想气体常数。
其中，σtd是动态拉伸强度；Pp是孔压峰值压力；rw是炮孔半径；α是应力波衰减指数，α＝2-μd/(1-μd)；b＝1-β2，μd是动态泊松比，β是岩石的横波速度与纵波速度的比值。
w0＝πrs 2/h (3)
其中，rs是射孔通道的半径；h是试样高度。
更具体地，动态抗拉强度和动态泊松比可由下式计算：
σtd＝Ktσt，μd＝0.8μ
其中，Kt是岩石在动态载荷下的抗拉强度增加系数；σt是静态抗拉强度；μ是静态泊松比。
可选地，如果钻孔压力不易测量，可以用冲击载荷幅值和地应计算初始裂缝长度。
可选地，裂缝初始宽度也可按照下式计算：
其中，E为弹性模量。
S302计算流经裂缝入口的气体质量：
其中，Cp材料纵波速度；αs是比例系数；γ是孔压升压速率；σ是地应力。
S303根据气体质量守恒计算裂缝数量：
S400结合得到的裂缝数量和初始长度建立裂缝扩展模型，计算裂缝扩展长度。具体地，如图3所示，步骤S400包括以下过程：
S401将得到孔压曲线划分为N个时间步，得到的裂缝数量N0、初始长度L0以及岩石力学参数均为输入数据；
S402先假设裂缝尖端与气体尖端不重合，指定θ”(t)和p”(θ”,t)的初值满足：
0＜θ″(t)＜1,p″(θ”,t)＝0 (6)
其中，θ”(t)是气体贯入深度；p”(θ”,t)是气体尖端与裂缝入口压力比。
S403根据应力强度因子判据计算得到气体压力分布指数m：
其中，p(θ,t)是压力分布函数。
其中，p(0,t)是该时刻裂缝入口处的压力；p(θ”,t)是该时刻气体流动尖端处的压力。
更具体地，求解m时可假设燃爆裂缝内爆生气体压力分布指数m的初值为1，运用牛顿迭代法寻找满足裂缝起裂扩展判据的缝内爆生气体压力分布指数m。
S404计算得到该时刻气体密度，裂缝张开位移，气体流动速度，裂缝扩展速度以及气体滤失速度：
其中，Mg是分子气体质量，pk(θ)和Tk(θ)是相对位置θ处的气体压力和温度。
其中，θ1和θ2是裂缝扩展过程中的瞬时长度和微段长度；G是剪切模量。
其中，ε是表面粗糙度；a，b是Huit实验测量值，a＝0.1，b＝0.5。
其中，Cr是瑞利波速度；KIC是岩石的断裂韧性。
其中，β是气体压缩系数，k0是渗透率，Φ是孔隙率，μ是气体粘度系数，L1是气体穿透深度。
S405判断气体在裂缝中的流动过程满足质量守恒：
其中，是气体流入质量流量，/>是压差作用下的气体过滤质量流量，/>是裂缝中的气体质量：
若满足守恒方程，令进行下一时间步的裂缝扩展计算；若不满足，改变θ”(t)的值(0<θ”(t)<1)，返回步骤S403。
更具体地，令p”(θ”,t)＝0，并从θ”(t)＝0时开始迭代，求取压力分布以及S404中各种参数，寻找满足守恒方程的气体贯入深度θ”(t)。
S406如果当θ″(t)≥1时仍不满足守恒条件，则假设裂缝尖端与气体尖端重合指定θ”(t)和p”(θ”,t)的初值满足：
0＜p″(θ”,t)＜1,θ″(t)＝0 (18)
执行步骤S403-S404，判断气体在裂缝中的流动过程满足质量守恒。若满足守恒方程，令进行下一时间步的裂缝扩展计算；若不满足，改变p”(θ”,t)的值(0<p”(θ”,t)<1)，返回步骤S403。
更具体地，令θ”(t)＝1，并从p”(θ”,t)＝0时开始迭代，求取压力分布以及S404中各种参数，寻找满足守恒方程的气体尖端与裂缝入口处的压力之比p”(θ”,t)。
S407如果当p″(θ,t)≥1时仍不满足守恒条件，则说明该时间步裂缝不会扩展，则进行下一时间步的裂缝扩展计算；
S408重复步骤S402-S407，直至全部时间步计算完成。
为了进一步的说明本申请的技术方案，以下以举例的形式进行说明(应当理解其仅为多种实施例中的一种)：
先进行多种室内试验测定岩石力学参数，见表1。
表1岩石力学参数
| 力学参数 | 值 | 力学参数 | 值 |
| 材料密度(kg/m3) | 2280 | 弹性模量(GPa) | 22.52 |
| 抗压强度(MPa) | 84.08 | 泊松比 | 0.24 |
| 抗拉强度(MPa) | 6.84 | 剪切模量(GPa) | 9.08 |
再根据力学参数选取水泥灌浆料作为相似材料，浇筑直径和高度均为1米的圆柱试样进行燃爆压裂试验。实验后观察到试样产生了9条裂缝，裂缝长度范围为1.25～1.5m。
根据冲击波传感器测得的压力曲线提取到的压力参数如表2所示。
表2压力参数
将力学参数和压力参数带入起裂模型中，计算得到其裂缝数量为9，初始长度为0.0218m。最后将上述参数带入裂缝扩展模型中，计算得到的裂缝长度和气体分布随时间变化如图4所示。计算得到的裂缝长度为1.277m，且在整个裂隙扩展过程中，裂隙尖端和气体尖端呈现出重合、非重合、最终重合的变化过程。
进一步的，基于上述实施例的内容，作为一种可选的实施例，本发明实施例中提供的一种燃爆压裂缝数量和长度预测模型的构建系统，包括：
第一主控模块，用于获取岩石的力学参数；第二主控模块，用于根据所述力学参数，选择岩石的相似材料作为试样，对试样进行压裂实验，以获取膛压和孔压的升压曲线；第三主控模块，用于根据所述力学参数、膛压和孔压的升压曲线建立裂缝起裂模型，计算裂缝数量和初始长度；第四主控模块，用于基于所述裂缝数量和初始长度构建裂缝扩展模型，计算裂缝扩展长度。该系统用于实现上述任意实施例涉及的方法。
本发明实施例的方法是依托电子设备实现的，因此对相关的电子设备有必要做一下介绍。基于此目的，如本发明的实施例提供了一种电子设备，该电子设备包括：至少一个处理器(processor)、通信接口(Communications Interface)、至少一个存储器(memory)和通信总线，其中，至少一个处理器，通信接口，至少一个存储器通过通信总线完成相互间的通信。至少一个处理器可以调用至少一个存储器中的逻辑指令，以执行前述各个方法实施例提供的方法的全部或部分步骤。
此外，上述的至少一个存储器中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个方法实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (9)
1.一种燃爆压裂缝数量和长度预测模型的构建方法，其特征在于，包括以下步骤：
S100获取岩石的力学参数；
S200根据所述力学参数，选择岩石的相似材料作为试样，对试样进行压裂实验，以获取膛压和孔压的升压曲线；
S300根据所述力学参数、膛压和孔压的升压曲线建立裂缝起裂模型，计算裂缝数量和初始长度；
S301根据现有参数计算燃爆生成的气体密度，裂缝的初始长度和宽度：
ρ＝PM/RT (1)
其中，T是峰值温度；P是膛压峰值压力；M是气体的分子质量；R是理想气体常数；
S302计算流经裂缝入口的气体质量：
其中，Cp材料纵波速度；αs是比例系数；γ是孔压升压速率；σ是地应力；
S303根据气体质量守恒计算裂缝数量：
其中：
σtd是动态拉伸强度；Pp是孔压峰值压力；rw是炮孔半径；α是应力波衰减指数，α＝2-μd/(1-μd)；b＝1-β2，μd是动态泊松比，β是岩石的横波速度与纵波速度的比值；
w0＝πrs 2/h (5)
其中，rs是射孔通道的半径；h是试样高度；S400基于所述裂缝数量和初始长度构建裂缝扩展模型，计算裂缝扩展长度。
2.根据权利要求1所述的一种燃爆压裂缝数量和长度预测模型的构建方法，其特征在于，步骤S100包括以下步骤：
进行多种室内实验获取岩石的力学参数；
具体的，首先测得岩石密度；然后进行单轴压缩试验测得岩石的抗压强度、弹性模量和泊松比；最后进行巴西劈裂试验测得岩石的抗拉强度。
3.根据权利要求1所述的一种燃爆压裂缝数量和长度预测模型的构建方法，其特征在于，步骤S200包括以下步骤：
S201根据岩石的力学参数选择岩石的相似材料，并浇筑试样；
S202对试样进行燃爆压裂实验并测量致裂管内的膛压曲线和作用在孔壁上的孔压曲线；
S203根据膛压和孔压的升压曲线分别得到膛压和孔压的峰值压力和升压速率。
4.根据权利要求1所述的一种燃爆压裂缝数量和长度预测模型的构建方法，其特征在于，步骤S400包括以下步骤：
S401将得到孔压曲线划分为N个时间步，得到的裂缝数量N0、初始长度L0以及岩石力学参数均为输入数据；
S402假设裂缝尖端与气体尖端不重合，设定气体贯入深度和气体尖端与裂缝入口压力比的初值满足：
0＜θ″(t)＜1,p″(θ”,t)＝0
式中，θ”(t)是气体贯入深度；p”(θ”,t)是气体尖端与裂缝入口压力比；
S403根据应力强度因子判据计算得到气体压力分布指数m；
S404计算得到第i时间步气体密度、裂缝张开位移、气体流动速度、裂缝扩展速度以及气体滤失速度，N≥i＞0；
S405判断气体在裂缝中的流动过程满足质量守恒，若是，进入步骤S406，否则返回至步骤S403，改变气体贯入深度的值；
S406重复步骤S403至步骤S405，进行下一时间步的裂缝扩展计算，直至N个时间步计算完成。
5.根据权利要求4所述的一种燃爆压裂缝数量和长度预测模型的构建方法，其特征在于，所述步骤S405还包括：若气体在裂缝中的流动过程不满足质量守恒：
且当θ″(t)≤0时仍不满足守恒条件，则设定气体贯入深度和气体尖端与裂缝入口压力比的初值满足：
0＜p″(θ”,t)＜1,θ″(t)＝0
并执行步骤重复步骤S403至步骤S404，判断气体在裂缝中的流动过程满足质量守恒，若是，进入步骤S406，否则，改变p”(θ”,t)的值，返回步骤S403；
若改变p”(θ”,t)的值，且p″(θ,t)≥1时仍不满足守恒条件，则第i时间步裂缝不会扩展，进入步骤S406，以进行下一时间步的裂缝扩展计算。
6.根据权利要求5所述的一种燃爆压裂缝数量和长度预测模型的构建方法，其特征在于，所述气体压力分布指数m为：
式中，K为应力强度因子；Lf为裂缝长度；p(θ,t)是压力分布函数；σ是地应力；θ是计算点的相对位置，即横坐标值与裂缝长度之比；
步骤S404中，所述气体密度为：
式中，Mg是分子气体质量，pk(θ)和Tk(θ)分别是相对位置θ处的气体压力和温度，R是理想气体常数；
所述裂缝张开位移为：
式中，μ是气体粘度系数，为第k时间步的裂缝长度，θ1和θ2分别是裂缝扩展过程中的瞬时长度和微段长度；
所述气体流动速度为：
式中，ε是表面粗糙度，a，b均为实验测量值，是缝内气体压力梯度；
所述裂缝扩展速度为：
式中，Cr是瑞利波速度，KIC是岩石的断裂韧性，p为缝内气体压力；
所述气体滤失速度为：
式中，β是气体压缩系数，k0是渗透率，Φ是孔隙率，μ是气体粘度系数，L1是气体穿透深度；
步骤S405中判断气体在裂缝中的流动过程质量守恒的模型为：
式中，是气体流入质量流量，/>是压差作用下的气体过滤质量流量，/>和/>分别是第k和第k-1时间步计算得到的裂缝中的气体质量。
7.一种燃爆压裂缝数量和长度预测模型的构建系统，其特征在于，应用如权利要求1-6中任一项所述的燃爆压裂缝数量和长度预测模型的构建方法实现，包括：
第一主控模块，用于获取岩石的力学参数；
第二主控模块，用于根据所述力学参数，选择岩石的相似材料作为试样，对试样进行压裂实验，以获取膛压和孔压的升压曲线；
第三主控模块，用于根据所述力学参数、膛压和孔压的升压曲线建立裂缝起裂模型，计算裂缝数量和初始长度；
第四主控模块，用于基于所述裂缝数量和初始长度构建裂缝扩展模型，计算裂缝扩展长度。
8.一种电子设备，其特征在于，包括：
至少一个处理器、至少一个存储器和通信接口；其中，
所述处理器、存储器和通信接口相互间进行通信；
所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令，以执行权利要求1至6任一项权利要求所述的方法。
9.一种非暂态计算机可读存储介质，其特征在于，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行权利要求1至6中任一项权利要求所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310839336.6A CN116877039B (zh) | 2023-07-10 | 2023-07-10 | 一种燃爆压裂缝数量和长度预测模型的构建方法及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310839336.6A CN116877039B (zh) | 2023-07-10 | 2023-07-10 | 一种燃爆压裂缝数量和长度预测模型的构建方法及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116877039A CN116877039A (zh) | 2023-10-13 |
| CN116877039B true CN116877039B (zh) | 2024-03-22 |
Family
ID=88259847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310839336.6A Active CN116877039B (zh) | 2023-07-10 | 2023-07-10 | 一种燃爆压裂缝数量和长度预测模型的构建方法及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116877039B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101173603A (zh) * | 2006-04-13 | 2008-05-07 | 施蓝姆伯格技术公司 | 通过快速氧化增强地下煤层甲烷井 |
| CN107545113A (zh) * | 2017-09-08 | 2018-01-05 | 西南石油大学 | 非常规油气藏水力压裂复杂缝网形成过程模拟方法 |
| CN109598029A (zh) * | 2018-11-13 | 2019-04-09 | 山西潞安环保能源开发股份有限公司常村煤矿 | 一种应用于高瓦斯煤巷co2高压气体致裂的设计方法 |
| CN111579579A (zh) * | 2020-05-26 | 2020-08-25 | 中国地质大学（武汉） | 多点全方位co2相变致裂压力时程曲线测试系统与方法 |
| CN112036098A (zh) * | 2020-09-15 | 2020-12-04 | 中国石油大学(华东) | 一种深层油气藏水力裂缝扩展数值模拟的方法 |
| CN112364518A (zh) * | 2020-11-20 | 2021-02-12 | 中国石油大学(华东) | 非常规油气地质工程一体化开发运筹学优化定量决策方法 |
| CN113107363A (zh) * | 2021-05-31 | 2021-07-13 | 中国矿业大学 | 一种促进裂隙错动提升自支撑能力的水平井施工方法 |
| CN113281182A (zh) * | 2021-05-25 | 2021-08-20 | 中国科学院武汉岩土力学研究所 | 一种多手段集成的压裂缝定量评价方法 |
| CN115577645A (zh) * | 2022-12-08 | 2023-01-06 | 中国石油大学(华东) | 一种燃爆压裂破裂范围预测模型的构建方法及预测方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7730951B2 (en) * | 2008-05-15 | 2010-06-08 | Halliburton Energy Services, Inc. | Methods of initiating intersecting fractures using explosive and cryogenic means |
| WO2012003027A1 (en) * | 2010-06-28 | 2012-01-05 | Exxonmobil Upstream Research Company | Method and system for modeling fractures in ductile rock |
| WO2016061179A1 (en) * | 2014-10-16 | 2016-04-21 | Reme, L.L.C. | Smart lower end |
-
2023
- 2023-07-10 CN CN202310839336.6A patent/CN116877039B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101173603A (zh) * | 2006-04-13 | 2008-05-07 | 施蓝姆伯格技术公司 | 通过快速氧化增强地下煤层甲烷井 |
| CN107545113A (zh) * | 2017-09-08 | 2018-01-05 | 西南石油大学 | 非常规油气藏水力压裂复杂缝网形成过程模拟方法 |
| CN109598029A (zh) * | 2018-11-13 | 2019-04-09 | 山西潞安环保能源开发股份有限公司常村煤矿 | 一种应用于高瓦斯煤巷co2高压气体致裂的设计方法 |
| CN111579579A (zh) * | 2020-05-26 | 2020-08-25 | 中国地质大学（武汉） | 多点全方位co2相变致裂压力时程曲线测试系统与方法 |
| CN112036098A (zh) * | 2020-09-15 | 2020-12-04 | 中国石油大学(华东) | 一种深层油气藏水力裂缝扩展数值模拟的方法 |
| CN112364518A (zh) * | 2020-11-20 | 2021-02-12 | 中国石油大学(华东) | 非常规油气地质工程一体化开发运筹学优化定量决策方法 |
| CN113281182A (zh) * | 2021-05-25 | 2021-08-20 | 中国科学院武汉岩土力学研究所 | 一种多手段集成的压裂缝定量评价方法 |
| CN113107363A (zh) * | 2021-05-31 | 2021-07-13 | 中国矿业大学 | 一种促进裂隙错动提升自支撑能力的水平井施工方法 |
| CN115577645A (zh) * | 2022-12-08 | 2023-01-06 | 中国石油大学(华东) | 一种燃爆压裂破裂范围预测模型的构建方法及预测方法 |
Non-Patent Citations (2)
| Title |
|---|
| 我院高能气体压裂技术十年发展综述;秦发动, 吴晋军;西安石油学院学报(自然科学版)(第03期);14-17+52 * |
| 油气层爆燃压裂裂缝动态延伸模型;陈德春;吴晓东;李海波;吴飞鹏;;中国石油大学学报(自然科学版)(第04期);103-107 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116877039A (zh) | 2023-10-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Stress-dependent fluid flow and permeability in fractured media: from lab experiments to engineering applications | |
| CN109374867B (zh) | 一种基于离散元的砂砾岩水力压裂模拟方法 | |
| Younessi et al. | Sand production simulation under true-triaxial stress conditions | |
| CN103206206A (zh) | 用于计算超深井地层破裂压力的方法 | |
| CN104655820A (zh) | 隧道硬质岩岩爆的判断、等级划分及处理方法 | |
| Wang et al. | Case study of modified H–B strength criterion in discrimination of surrounding rock loose circle | |
| Wang et al. | Three-dimensional experimental and numerical investigations on fracture initiation and propagation for oriented limited-entry perforation and helical perforation | |
| Xie et al. | Dynamic response of rock mass subjected to blasting disturbance during tunnel shaft excavation: a field study | |
| Mohamadi et al. | Comparison of the classical and fracture mechanics approaches to determine in situ stress/hydrofracturing method | |
| Liu et al. | Nondestructive testing on cumulative damage of watery fractured rock mass under multiple cycle blasting | |
| Guo et al. | Study on fracture morphological characteristics of refracturing for longmaxi shale formation | |
| Jiang et al. | Calculation of fracture number and length formed by methane deflagration fracturing technology | |
| CN116877039B (zh) | 一种燃爆压裂缝数量和长度预测模型的构建方法及设备 | |
| Park et al. | Investigation of the development of an excavation damaged zone and its influence on the mechanical behaviors of a blasted tunnel | |
| Fan et al. | Influence of water depth on the range of crushed zones and cracked zones for underwater rock drilling and blasting | |
| Hou et al. | Experimental study on deformation and acoustic emission characteristics of arch roadway under different unloading rates | |
| CN114737925B (zh) | 一种水压致裂煤岩体瓦斯渗流模拟装置及抽采量预测方法 | |
| Wen et al. | The use of Hoek Brown failure criterion on determination of the geo-mechanical parameters of a grouting consolidation body | |
| Lv et al. | Empty-hole effect on fracture propagation under blasting load | |
| Yang et al. | Response characteristics of coal measure strata subjected to hydraulic fracturing: insights from a field test | |
| Xu et al. | A Study on Triaxial Unloading Test of Columnar‐Jointed‐Rock‐Mass‐Like Material with AW Velocity Analysis | |
| Liu et al. | A new method for wellbore stability evaluation based on fractured carbonate reservoir rock breaking degree | |
| Deng et al. | A Model for Estimating Penetration Length under Different Conditions | |
| Jiang et al. | Impact analysis of multiple parameters on fracture formation during volume fracturing in coalbed methane reservoirs | |
| Yang et al. | Analysis of explosion stress field of air radial decoupling charge |
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 |