【技术实现步骤摘要】
本申请涉及地下氢气存储，尤其涉及一种岩心氢气扩散系数模拟计算方法以及相关设备。
技术介绍
1、中国天然气产业迅速发展，天然气消费持续增长，在国家能源体系中的重要性不断提高，同时，天然气等能源需求季节性差异强，充足的储备是保证稳定供应的关键，因此，储气库建设迫在眉睫。另外，近年来，随着新能源的发展和国家能源转型战略的实施，氢气等清洁能源消费和供给需求快速上涨，气体长期、大规模储存是能源产业中至关重要的一环。
2、目前建造和运营的储气库主要为枯竭油气藏储气库、含水层储气库及盐穴储气库。对于枯竭油气藏储气库和含水层储气库，主要储存空间为储层多孔介质岩石。而对于盐穴储气库，其主要储存空间为岩盐地层溶腔，溶腔围岩主要为岩盐或夹层泥砂岩，同样属于多孔介质。因此在储氢库运行过程中，均涉及多孔介质气体运移。
3、气体流通性和渗透性强，在岩石中的扩散作用不容忽视。扩散系数是表征岩石介质中气体扩散的主要参数。而扩散系数通常与气体的特性及岩石多孔介质的孔隙结构有关。气体扩散计算主要有欧拉法和拉格朗日法。相比于欧拉法，拉格朗日法由于其计算方法简便、运算要求低等特点，更适合于强非均质介质单元的流体流动计算。现有的岩石气体扩散系数数值测试模型主要岩石孔隙结构或孔隙迂曲度计算，仅关注气体在孔隙空间中的扩散和运移，忽视了气体在岩石矿物中的扩散，同时，现有模型无法计算岩石中气体扩散各向异性。例如，专利号为cn202111284159.7公布的一种计算岩石迂曲度和气体扩散系数的方法仅考虑岩石中孔隙扩散，并未考虑矿物内扩散，因此无法获得准确的
技术实现思路
1、为了解决上述问题，本申请实施例提供了一种岩心氢气扩散系数模拟计算方法，在ct扫描微观岩石样品重构模型基础上，充分考虑岩石孔隙和矿物内气体扩散特征，有效解决复杂强非均质和各向异性岩石的扩散系数计算问题，所述技术方案如下：
2、本申请第一方面提供一种岩心氢气扩散系数模拟计算方法，包括以下步骤：
3、s1获取岩石样品ct扫描灰度数据，按照矿物及孔隙灰度阈值进行阈值分割；
4、s2基于时间域扩散气体扩散系数反演方法，模拟气体岩心单元体中扩散运移，获取某时刻在某方向上距离与气体浓度比数据；
5、s3根据解析模型，拟合某方向上气体扩散距离与浓度比的关系方程，计算岩石样品有效扩散系数。
6、例如，在一个实施例提供的岩心氢气扩散系数模拟计算方法中，所述s1中，获取岩石样品ct扫描灰度数据，按照矿物及孔隙灰度阈值进行阈值分割，划分像素单元类型i，确定样品孔隙结构和矿物分布，并导出划分结果数据，其中，像素单元类型i＝0，1，2…n。
7、例如，在一个实施例提供的岩心氢气扩散系数模拟计算方法中，所述基于时间域扩散气体扩散系数反演方法，模拟气体岩心单元体中扩散运移，获取某时刻在某方向上距离与气体浓度比数据，包括以下步骤：
8、s2.1读取岩石样品ct扫描划分结果数据，根据需要截取预设尺寸长度的样品数据，获取样品坐标，并根据所述像素单元类型为其赋予扩散系数di和孔隙度φi；
9、s2.2在样品数据外边界添加密封边界单元，并修正样品坐标；
10、s2.3设置样品数据外边界某一面为气体运移初始位置，记初始边界处气体扩散单元组数为c0，并按照随机分布方式布置气体扩散初始位置点；
11、s2.4开始模拟每个气体扩散单元组扩散运移情况，并保持初始位置处气体扩散单元组数不变，计算迭代次数为a，计算气体从像素单元i扩散像素单元j的概率pij和时间tij；
12、s2.5计算每个气体扩散单元组的累计扩散时间；
13、s2.6统计某一累计时刻t时气体扩散单元组所在位置，计算该时刻测试方向的每个法向面上的气体扩散单元组数ci并计算其与初始边界处气体扩散单元组数c0之比，即浓度比ci/c0。
14、例如，在一个实施例提供的岩心氢气扩散系数模拟计算方法中，所述s2.1中根据需要截取预设尺寸长度的样品数据，样品坐标为(x，y，z)，样品坐标范围为([0:s],[0:p],[0:q])，所述s2.2中，在样品数据外边界添加密封边界单元，密封边界单元扩散系数d0＝0，孔隙度φ0＝0，样品坐标范围修正为([0:s+2a],[0:p+2a],[0:q+2a])，其中，a为密封边界单元的宽度值。
15、例如，在一个实施例提供的岩心氢气扩散系数模拟计算方法中，所述s2.4中气体从像素单元i扩散像素单元j的概率pij和时间tij满足以下关系式：
16、
17、
18、其中，sij为样品像素单元i和像素单元j的交界面面积；(φd)ij为像素单元i和像素单元j的扩散系数d和孔隙度φ乘积的调和平均数；lij为像素单元i和像素单元j中心之间的距离；vi为像素单元i的体积；u01为[0,1]的均匀分布随机数。
19、例如，在一个实施例提供的岩心氢气扩散系数模拟计算方法中，所述s2.6中测试方向的每个法向面上的气体扩散单元组数ci与初始边界处气体扩散单元组数c0之比满足以下关系式：
20、
21、其中，li为气体沿某向的扩散距离，de为岩石样品有效扩散系数。
22、本申请第二方面提供一种岩心氢气扩散系数模拟计算装置，包括：ct扫描单元、阈值分割单元、模拟单元和计算单元，所述ct扫描单元用于对岩石样品进行ct扫描，获取岩石样品ct扫描灰度数据；所述阈值分割单元用于按照矿物及孔隙灰度阈值进行阈值分割，划分像素单元类型，确定样品孔隙结构和矿物分布，并导出划分结果数据；所述模拟单元基于时间域扩散气体扩散系数反演方法，模拟气体岩心单元体中扩散运移，获取某时刻在某方向上距离与气体浓度比数据；所述计算单元根据解析模型，拟合某方向上气体扩散距离与浓度比的关系方程，计算岩石样品有效扩散系数。
23、本申请第三方面提供一种电子设备，包括：存储器、处理器，所述处理器用于执行存储器中存储的计算机管理类程序时实现如上述所述的岩心氢气扩散系数模拟计算方法。
24、本申请第四方面提供一种计算机可读存储介质，其上存储有计算机管理类程序，所述计算机管理类程序被处理器执行时实现上述所述的岩心氢气扩散系数模拟计算方法。
25、本申请一些实施例提供的一种岩心氢气扩散系数模拟计算方法带来的有益效果为：本申请在ct扫描微观岩石样品重构模型基础上，充分考虑岩石孔隙和矿物内气体扩散特征，构建气体扩散计算模拟方法可以便捷有效地确定岩石扩散系数，有效解决复杂强非均质和各向异性岩石的扩散系数计算问题，并通过拟合反演，方便、快捷、准确的计算岩样气体扩散系数，并且充分考虑气体真实扩散特性，有效完善了岩石气体扩散系数计算理论。本申请可为岩石气体扩散系数计算提供新的思路，对于储气库建设运行过程中注采能力和密封性评价具有积极的指导意义。
本文档来自技高网...【技术保护点】
1.一种岩心氢气扩散系数模拟计算方法，其特征在于，包括以下步骤：
2.根据权利要求1所述岩心氢气扩散系数模拟计算方法，其特征在于，所述S1中，获取岩石样品CT扫描灰度数据，按照矿物及孔隙灰度阈值进行阈值分割，划分像素单元类型i，确定样品孔隙结构和矿物分布，并导出划分结果数据，其中，像素单元类型i＝0，1，2…n。
3.根据权利要求1所述岩心氢气扩散系数模拟计算方法，其特征在于，所述基于时间域扩散气体扩散系数反演方法，模拟气体岩心单元体中扩散运移，获取某时刻在某方向上距离与气体浓度比数据，包括以下步骤：
4.根据权利要求3所述岩心氢气扩散系数模拟计算方法，其特征在于，所述S2.1中根据需要截取预设尺寸长度的样品数据，样品坐标为(X，Y，Z)，样品坐标范围为([0:s],[0:p],[0:q])，所述S2.2中，在样品数据外边界添加密封边界单元，密封边界单元扩散系数D0＝0，孔隙度φ0＝0，样品坐标范围修正为([0:s+2a],[0:p+2a],[0:q+2a])，其中，a为密封边界单元的宽度值。
5.根据权利要求3所述岩心氢气扩散系数
6.根据权利要求3所述岩心氢气扩散系数模拟计算方法，其特征在于，所述S2.6中测试方向的每个法向面上的气体扩散单元组数Ci与初始边界处气体扩散单元组数C0之比满足以下关系式：
7.一种岩心氢气扩散系数模拟计算装置，其特征在于，包括：
8.一种电子设备，其特征在于，包括：
9.一种计算机可读存储介质，其上存储有计算机管理类程序，其特征在于：所述计算机管理类程序被处理器执行时实现如权利要求1至6中任意一项所述的岩心氢气扩散系数模拟计算方法。
...【技术特征摘要】
1.一种岩心氢气扩散系数模拟计算方法，其特征在于，包括以下步骤：
2.根据权利要求1所述岩心氢气扩散系数模拟计算方法，其特征在于，所述s1中，获取岩石样品ct扫描灰度数据，按照矿物及孔隙灰度阈值进行阈值分割，划分像素单元类型i，确定样品孔隙结构和矿物分布，并导出划分结果数据，其中，像素单元类型i＝0，1，2…n。
3.根据权利要求1所述岩心氢气扩散系数模拟计算方法，其特征在于，所述基于时间域扩散气体扩散系数反演方法，模拟气体岩心单元体中扩散运移，获取某时刻在某方向上距离与气体浓度比数据，包括以下步骤：
4.根据权利要求3所述岩心氢气扩散系数模拟计算方法，其特征在于，所述s2.1中根据需要截取预设尺寸长度的样品数据，样品坐标为(x，y，z)，样品坐标范围为([0:s],[0:p],[0:q])，所述s2.2中，在样品数据外边界添加密封边界单元，密封边界单...
【专利技术属性】
技术研发人员：宋睿，宋宇家，刘建军，
申请(专利权)人：中国科学院武汉岩土力学研究所，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。