买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：中国科学院武汉岩土力学研究所

摘要：本发明提供一种预测利用楔裂压力提高的储层固碳量的方法。利用楔裂压力提高储层固碳量的预测方法包括：建立储层代表单元矿化捕集的三维几何模型；提出CO2矿化封存过程中的化学力学耦合数学模型；基于地层物性参数、应力参数、热化学参数计算单位体积储层岩体在楔裂压力驱动下的固碳量。该发明考虑储层楔裂压力对矿化捕集的促进作用，基于压溶效应原理建立储层代表单元矿化捕集的三维几何模型和储层压溶矿化封存量的计算模型，评估楔裂压力驱动下目标储层的矿物成分、粒度以及孔隙度对矿化封存量的影响，评价目标储层利用楔裂压力加速矿化捕集的潜力。该发明为提高储层矿化固碳量提供了一种可行的评价方法，在碳封存领域有很大的应用前景。

主权项：1.一种预测利用楔裂压力提高的储层固碳量的方法，其特征在于，步骤包括：步骤1，基于目标地层岩样颗粒粒径D50和孔隙率目标地层岩样的n种矿物成分及其百分比αx，建立储层代表单元矿化捕集的三维几何模型，三维几何模型包括立方体及其内嵌的切割球体，其中立方体内部的球体代表矿物颗粒，而立方体内球体外的空间代表孔隙体积，矿物颗粒与相邻的六个立方体分别在三个主应力方向上有六个接触面，初始时刻t＝0时的初始参数包括：切割球体半径r，切割球体初始损失高度切割球体的初始接触面面积切割球体的某一主应力方向上的接触反应面积切割球体沿三个主应力方向的损失体积i＝1，2，3，表示三个主应力方向，x＝1,2,3.....n，n代表矿物种类数量；步骤2，建立CO2矿化封存过程中的化学力学耦合数学模型，基于t时刻内切割球形在某一主应力方向上的损失体积的变化速率等于矿物颗粒在该应力方向上压溶反应体积溶解速率，求出再根据r和公式6更新t时刻沿三个主应力方向的切割球体损失高度 将r和带入公式7更新t时刻沿三个主应力方向的立方体单元面积 将r和带入公式8更新t时刻沿三个主应力方向的接触面积 将和αx带入公式9更新t时刻沿三个主应力方向的矿物反应面积 将t时刻矿物颗粒的接触面沿三个不同主应力方向的楔裂压力与临界平衡压力σc比较，来判断矿物压溶反应启动与否，当时，压溶反应被启动；当时，压溶反应中止，t＝wΔt，w＝1,2…m； 为时刻t某一主应力方向上的楔裂压，表达式如下： σi为目标地层的地应力场的主应力方向的地层压力，p为地下水压力，σc为临界平衡压力，表达式如下： 其中Vm为目标储层岩样的摩尔体积，由查表获得，T为地层温度，从地层物探资料中查询，Em和Tm为熔化热和熔化温度，由查表获得；在三个应力方向上的化学力学耦合数学模型，包括：矿物颗粒与CO2的压溶反应速率：t时刻某矿物颗粒的压溶反应体积溶解速率：其中：k+为溶解速率常数，为矿物摩尔溶解体积，R为气体常数，均由查表获得；步骤3，基于地层物理参数、应力参数、热化学参数计算单位体积储层岩体在楔裂压力驱动下的固碳量，迭代步骤如下：迭代步骤a.通过公式12得到时刻t某矿物颗粒与CO2的反应速率根据公式14即可获得t+Δt时刻矿化捕集的CO2累计摩尔量： 迭代步骤b.获得t+Δt时刻切割球体沿三个主应力方向的损失体积ΔVdi为Δt时步内切割球体某个主应力方向上的溶解体积： 根据公式6更新t+Δt时刻沿三个主应力方向的切割球体损失高度根据公式7更新沿三个主应力方向的立方体单元面积根据公式8更新沿三个主应力方向的接触面积根据公式9更新沿三个主应力方向的矿物反应面积迭代步骤c.以更新后的切割球体损失高度立方体单元面积接触面积作为输入参数，重复迭代步骤a～c，直到矿物压溶反应中止，循环w次，合计历时wΔt，得到最终矿化捕集的CO2累计摩尔量和单位体积储层岩体在压溶效应下的固碳量表达式如下： V0:表示代表单元初始体积，

全文数据：

权利要求：

百度查询： 中国科学院武汉岩土力学研究所 一种预测利用楔裂压力提高的储层固碳量的方法