本文在Langmuir等温吸附理论的基础上,建立了考虑应力敏感作用的裂缝性页岩气藏在不同边界条件下的双重介质渗流数学模型。考虑到模型存在很强的非线性,在求解过程中将解吸作用考虑成基质等温压缩系数的增加,降低了模型的非线性程度。通过Laplace变换、摄动变换最终得到了数学模型在Laplace空间的解析解,最后通过编程绘制了压力动态曲线,分析了压力动态曲线的特征,并对影响压力动态曲线的因素进行了参数敏感性分析,研究结果可为页岩气藏的试井分析提供科学依据。 Based on the theory of Langmuir isothermal adsorption, the mathematical model of double media seepage in fractured shale gas reservoirs which considers stress sensitivity under different boundary conditions is established in this paper. Taking the strong nonlinearity of the model into account, desorption effect can be considered as the increase of isothermal compressibility of matrix during the resolution to decrease the degree of nonlinearity of the model. The analytical solution of the mathematical model in Laplace space is obtained by means of Laplace transformation and perturbation transformation. Finally, pressure transient curves are drawn with the aid of computer programs, and its characteristics have been analyzed, and parameter sensitivity analysis for the contributing factors affecting the dynamic curve of pressure also has been done. The results can provide scientific basis for well test analysis of shale gas reservoirs.
赵梓瑜,赵二猛
东北石油大学石油工程学院,黑龙江 大庆
收稿日期:2016年8月1日;录用日期:2016年8月19日;发布日期:2016年8月22日
本文在Langmuir等温吸附理论的基础上,建立了考虑应力敏感作用的裂缝性页岩气藏在不同边界条件下的双重介质渗流数学模型。考虑到模型存在很强的非线性,在求解过程中将解吸作用考虑成基质等温压缩系数的增加,降低了模型的非线性程度。通过Laplace变换、摄动变换最终得到了数学模型在Laplace空间的解析解,最后通过编程绘制了压力动态曲线,分析了压力动态曲线的特征,并对影响压力动态曲线的因素进行了参数敏感性分析,研究结果可为页岩气藏的试井分析提供科学依据。
关键词 :页岩气,压力动态,应力敏感
页岩气是新兴的能源资源,对于页岩气的渗流机理及压力动态分析尤为重要。国内外学者围绕页岩气的渗流机理做了很多研究工作 [
根据Langmuir方程:
式中VE为等温吸附量,m3/m3;VL为Langmuir体积,m3/m3;pL为Langmuir压力,MPa;p为气体压力,MPa。
由Pedrosa [
式中K为地层渗透率,m2;Ki为原始地层渗透率,m2;ψ为气藏拟压力,Pa2/(Pa.s);ψi为原始气藏拟压力,Pa2/(Pa.s);γ为渗透率模数,(Pa.s)/Pa2。
并定义如下无因次变量:无因次拟压力
可得无因次的渗流数学模型:
式中(a)、(b)、(c)分别对应无限大、圆形封闭以及圆形定压边界三种情况。
该数学模型是一个强非线性方程组,无法直接求解,为了简化方程组,引入如下变换:
关于γD的正则扰动:
根据摄动技术可得到其近似解。根据同登科等人 [
其中:
上述方程的形式与常规气藏类似,可以得到Laplace空间中不同边界条件下的解。
1) 无限大外边界
因此通解为
2) 圆形封闭边界
因此其通解为
3) 圆形定压外边界:
因此其通解为
令rD = 1便可得到井底压力解η0wD。如果要考虑井筒存储系数和表皮效应对压力以及产能的影响,可以通过以下关系式对井底压力解进行修正:
由于上面求得的井底拟压力是拉普拉斯空间的,所以要根据拉普拉斯空间的反演公式进行反演,得到真实空间中
如图1,考虑应力敏感后,压力导数曲线从井筒储集阶段之后的过渡阶段便开始出现上移。尤其在后期阶段,无因次渗透率模数γD越大,曲线上翘现象越明显,表现出与封闭边界类似的特征。因此可以得出结论,在页岩气藏中,应力敏感作用的存在会使得裂缝性气藏压力下降速度变快,即使油藏范围很大,单井的控制范围也有限。
由图2可知,无限大地层时,边界的影响非常小,因此压力曲线接近于一条直线,压力导数曲线为一条水平线。外边界为封闭边界时,当封闭边界影响到井筒以后,压力曲线末端压力趋近45˚线而压力导数为45˚线。定压边界时,后期流动将达到稳定流动阶段,因此曲线末端为压力水平线而压力导数急剧下滑。
通过(图3)曲线对比可以得到如下结论:由于吸附气主要存在于基质系统中,在压力导数曲线上表现为过渡段的“凹子”变深,随着Langmuir体积的不断增大,气藏中的吸附气的解吸量增大,补偿井底压力损失明显,井底压力降低缓慢。
图1. 不同渗透率模数下压力动态曲线
图2. 不同边界类型页岩气藏压力动态曲线
图3. 不同Langmuir体积下压力动态曲线
1) 建立并求解了考虑应力敏感作用的裂缝性页岩气藏在不同边界条件下的双重介质渗流数学模型。
2) 页岩气藏应力敏感作用的存在使得气藏压力降落速度变快,并使得单井控制程度有限,所以对于应力敏感性较强的储层,要合理控制压差,以防对储层造成损害。
3) 不同的外边界条件前期压力动态曲线是一样的,只是在流动阶段后期有不同的表现。供给半径的大小只影响外边界反映阶段,供给半径越大,探测到外边界的时间越晚。
4) 随着Langmuir体积的不断增大,在压力导数曲线上表现为过渡段的“凹子”变深,补偿井底压力损失明显,井底压力降低缓慢。
赵梓瑜,赵二猛. 应力敏感页岩气藏压力动态分析Pressure Dynamic Analysis of Stress-Sensitive Shale Gas Reservoir[J]. 自然科学, 2016, 04(03): 327-332. http://dx.doi.org/10.12677/OJNS.2016.43038