文摘

嵌入离散裂缝模型(EDFM)一直在流行的建模裂缝储层由于其灵活性和效率,同时保持了复杂的裂缝网络的几何。尽管EDFM已经验证了单相流模拟,一些最近的案例表明,EDFM在多相流的情况下可能导致明显的错误。projection-based嵌入离散裂缝模型(pEDFM)和整体嵌入离散裂缝模型(IEDFM)两个最近开发方法,它打算改善EDFM的准确性。在这项研究中,一个projection-based整体嵌入离散裂缝模型(pIEDFM),提出了相结合的优点pEDFM IEDFM。类似于pEDFM, pIEDFM使用一种额外的断裂之间的连接和nonneighboring基质细胞获得更多身体真实的速度场。显著改善,半锥形压力分布,采用IEDFM pIEDFM捕获的断裂表面附近的压力急剧变化。比较基准测试结果和explicit-fracture细网格仿真结果表明,该pIEDFM使用适量的网格提供准确的解决方案。拟议中的pIEDFM也应用于耦合流和地质模拟水库骨折。比较我们的耦合模拟结果与EDFM表明pIEDFM适用于耦合仿真,和不同的方法之间的matrix-fracture遗传性pIEDFM EDFM可能导致应力场地质建模预测的偏差,最终影响石油生产、含水率、含油饱和度资料。

1。介绍

数值模拟方法破碎水库引起了极大关注在过去的几十年。由于显著的不透水岩石渗透率差别矩阵和高导电骨折,建模多相流媒体已经成为具有挑战性的骨折。Barenblatt et al。1)首次提出双孔隙度的概念来描述渗流过程中岩石破裂了。沃伦和根(2]扩展双孔隙度的概念,开发了双孔隙度模型。齐米等。3,4]介绍了双孔隙度模型转化为石油工程和应用裂缝储层建模的方法。之后,一系列的数值方法是双孔隙度模型的扩展,包括双孔隙度双重渗透率模型(5,6),triple-porosity双重渗透率模型(7),和多个相互作用的连续介质模型(8]。所有这些模型可以分为dual-continuum模型。Dual-continuum模型提供了一种有效的方式模拟水库骨折。但是,骨折的几何图形是迷失在双孔隙度模型的假设,并不能代表真正的裂缝网络dual-continuum模型(9]。Panfili et al。10]表明dual-continuum模型中使用的均化可以引入非物质的断开连接之间的裂缝流动区域。Moinfar et al。11]研究了水库的例子复杂裂缝网络,dual-continuum模型不能提供准确的解决方案。

一个离散裂缝模型(DFM) (12- - - - - -15)是发达国家解决dual-continuum模型的局限性。在离散裂缝模型中,每一个裂缝被明确和个人,提供更多的物理和真实的描述骨折,尤其是在复杂的几何图形(16,17]。然而,适应复杂的裂缝网络,非结构化矩阵网格中常用的DFM,在粉磨造成困难。在断裂网络组成的微尺度骨折,生成的大量非结构化网格通常,导致高计算成本,使得DFM不切实际的在实际领域的研究(18]。

嵌入离散裂缝模型(EDFM)提出了李,李19)的一种替代DFM和dual-continuum模型。EDFM,笛卡尔网格用于矩阵离散化方法的效率。明确骨折是离散控制卷,也称为断裂元素。骨折元素嵌入到他们的父母通过matrix-fracture连接矩阵网格。EDFM已成功实现在垂直裂缝情况下(19],nonvertical骨折病例[18),空间的裂缝情况下(20.,21]。一些作者(22,23)结合EDFM dual-continuum方法与多尺度模型水库骨折。Moinfar et al。24)应用的EDFM耦合流和地质模拟水库骨折。

虽然EDFM已经在各种研究验证其准确性的单相流模拟,它可能导致明显的错误在某些情况下对多相模拟。图1显示了一个矩阵块米分割的骨折。两个矩阵的分数米断裂分离的标记为A和b,如图1(一),在现实情况下,传入的水流,走向一个裂缝在注水过程中首先进入分数B,然后分为两个部分,沿着(F₁在(F)和₂分别)骨折。然而,EDFM, A和B两个分数组合作为一个完整的矩阵块而不是被认为是两个单独的网格。因此,跨断裂的水流从分数分数B不能展出。相反,断裂之间的平均流量和矩阵块米评估(F_m-f),如图1 (b)。水从分数A和B可以流同时向裂缝,这就增加了水通量,沿着断裂和结果在一个“非物质的通量分裂”EDFM。因此,水通量沿断裂是高估了,和水通量穿过裂缝被低估了25]。

为了解决这个问题,一个projection-based嵌入离散裂缝模型(pEDFM)提出的烦躁和其他人。(26]。江和尤尼斯(25),可以实现更多的物理通量分裂pEDFM,因此修复错误的排水量EDFM预测的过程。Olorode et al。27]扩展pEDFM三维情况下,研究了3 d组分与pEDFM建模。Rao et al。28)修改原始pEDFM和发达micro-translation算法来帮助选择projection-face组合。

EDFM的另一个限制是过于简单化的假设的压力分布矩阵域邻近骨折。EDFM,使用一个近似线性压力分布假设周围骨折。然而,一个锥形分布产生的压力通常是由于大型矩阵和断裂之间的渗透率差异,如图2。EDFM的过于简单化的假设可能导致计算错误的遗传性matrix-fracture连接(29日]。一个整体嵌入离散裂缝模型(IEDFM) (30.)提出了改进EDFM的传递率计算。IEDFM, matrix-fracture连接的遗传性派生的管理系统,得到断裂表面附近的更现实的压力分布,提高建模的准确性流入水库骨折。IEDFM是后来扩展到各向异性裂缝储层的建模(31日]。

一个projection-based整体嵌入离散裂缝模型(pIEDFM)在这项研究中,提出了相结合的优点pEDFM IEDFM。类似于pEDFM,额外matrix-fracture连接添加到pIEDFM骨折nonneighboring矩阵元素和元素之间沿裂缝投影方向。邻近的遗传性和nonneighboring matrix-fracture连接在pIEDFM IEDFM派生管理系统使用方法。pIEDFM验证通过基准测试结果的准确性和细网格仿真结果。拟议中的pIEDFM然后应用于耦合流和地质模拟水库骨折。检查提出的数值方法的适用性。

2。裂缝性储层仿真的控制方程

2.1。质量守恒方程

断裂的媒体给出的质量守恒方程如下: 在哪里代表流体阶段,孔隙度,是饱和的阶段 , 相的密度吗 , 阶段的速度吗 , 通量的阶段吗从井,通量的阶段吗矩阵和断裂之间的元素。

阶段的速度由达西定律定义: 在哪里绝对磁导率,相的相对渗透率吗 , 相的粘度是吗 , 阶段的压力吗 , 重力加速度,是深度,流动相的潜力吗 。

方程(1)是使用控制卷离散的有限差分格式在时间空间和一阶方案,提供以下: 在下标表示元素的值 ,上标 代表当前时间水平,上标代表前面的时间水平,相流项吗元素之间和元素在哪里和相同类型的元素(矩阵或断裂),相流项吗通过matrix-fracture连接元素和元素 ,和通量的阶段吗从井。

流条件和表示如下: 在下标 和 表示适当的平均属性的接口,和和遗传性的连接。

遗传性的连接两个half-transmissibilities的调和平均数: 在哪里和的half-transmissibilities元素1和2,和绝对渗透率,和是距离的中心元素的接口。

2.2。EDFM

水库矩阵使用笛卡尔网格离散,骨折和额外的元素被添加到代表断裂控制卷。因此,有三种元素之间的连接EDFM:矩阵与矩阵,fracture-fracture, matrix-fracture连接。矩阵与矩阵和fracture-fracture连接,传递率可以从两点通量近似导出几何(TPFA)使用(5)。

matrix-fracture连接在2 d储层情况下,矩阵和骨折可以派生的half-transmissibilities如下: 在哪里骨折的表面积,基质渗透率,裂缝渗透率,骨折中心距离接口,孔径等于一半的断裂,然后呢相当于距离矩阵和骨折的元素。

使用近似的线性骨折周围压力分布,给出的等效距离可以如下: 在哪里骨折和距离吗矩阵的体积。

3所示。Projection-Based整体嵌入离散裂缝模型

提出pIEDFM,连接关系建立遵循pEDFM规则,在介绍了额外的连接断裂之间的元素和nonneighboring沿着裂缝投影矩阵元素的方向。如图3,断裂元素有两个投影矩阵元素的边界 ,它的面积吗和 ,分别。骨折的元素连接到相邻矩阵元素和两个nonneighboring矩阵元素和 。选择的标准矩阵的裂缝预测之前,江和尤尼斯的工作25]。矩阵的面孔,更接近断裂中心预计在每个维度面临选择。

对于矩阵与矩阵的连接,投影区域的骨折是消除界面区域。修改后的矩阵之间的界面区域和矩阵给出如下: 在哪里是原始矩阵之间的界面面积和矩阵 ,和骨折投射区沿x方向。当裂缝穿过矩阵元素时,矩阵与矩阵的连接将被删除。

图4显示了一个示例连接建立的pIEDFM EDFM。有四个矩阵元素标记为M1, M2, M3, M4。两个骨折离散为若干个骨折段元素的矩阵块边界。骨折段元素由蓝线与红点标记为F1, F2, F3, F4, F5、F6。EDFM,有4矩阵与矩阵的关系,5 fracture-fracture连接,6 matrix-fracture连接。然而,在pIEDFM 3矩阵与矩阵的关系,5 fracture-fracture连接,包括14 matrix-fracture连接。matrix-fracture连接的数量增加pIEDFM由于额外的骨折段之间的连接和nonneighboring矩阵元素。

pIEDFM,矩阵与矩阵和fracture-fracture连接遗传性可以直接导出使用(5)。matrix-fracture连接传递率的计算公式推导的基础上IEDFM,骨折的地方视为一系列点的下沉和传播能力是解决管理系统。内部的压力分布矩阵可以从水槽,骨折,因此繁殖的锥形分布压力提高仿真精度之间的流体交换矩阵和骨折。pIEDFM,传播性骨折和相邻矩阵之间的传播性骨折和nonneighboring矩阵分别计算。

图5的计算提出了一个示意图matrix-fracture pIEDFM遗传性。附近的骨折,点的压力可以叠加的水槽的压力下降造成的问题: 在哪里层的高度,是点下沉,的数量水槽的流量吗 , 点之间的距离吗和水槽 ,和是一个常数与断裂有关的压力。

当骨折被假定为等位,选择几个参考点断裂表面形成一个线性方程系统: 在哪里参考点之间的距离吗和水槽 。

线性方程系统可以改写如下: 在哪里定义如下:

解决线性方程系统给压力点的确切表达式 :

断裂之间的传播性元素和邻近的矩阵元素可以得到如下: 在哪里矩阵的体积吗 。

同样,断裂之间的传播性元素和nonneighboring矩阵元素可以得到如下: 在哪里矩阵的体积吗 , 界面的面积,骨折的面积投影在接口。

4所示。建模的地质力学

4.1。变形的骨折

断裂的变形行为强烈stress-dependent和非线性。基于实验室数据的实证模型用于计算断裂模(32]。在正常应力骨折,杨氏模量如下: 在哪里是初始裂缝正常的刚度,零应力骨折孔径,是有效的法向应力,断裂的最大正常关闭,节理粗糙度系数JRC是,JCS联合抗压强度。

骨折孔径变化在正常压力给出如下:

骨折在剪切应力,剪切模量如下: 在哪里刚度是数字,刚度指数,剪切应力,剪切应力峰值,是失败的比率,残余摩擦角。

变形断裂的剪切应力下的结合水平剪切位移和垂直剪切位移,也称为断裂扩张。横向剪切位移的断裂给出如下: 在哪里是断裂的特征长度。

垂直剪切位移(膨胀)使用一个经验模型给出如下(33]: 在哪里是峰值水平剪切位移,是垂直剪切位移峰值,是一个破坏系数。

骨折在正常和剪切应力,正常的孔径变化和裂缝扩张导致整体断裂孔径变化。因此,骨折孔径在压力下给出如下:

骨折孔径在(23)是两个断裂表面的平均点对点之间的距离,它被定义为“机械”光圈。然而,实际骨折有粗糙的墙壁和可变孔径,机械孔径不合适的渗透系数在计算断裂。“液压”光圈是由流分析和描述裂缝开展能力要好。实证液压孔径和机械孔径之间的关系可以得到如下(34]: 在哪里是动员节理粗糙度系数。

一般来说,液压骨折的孔径小于机械孔径由于粗糙度和曲折断裂的墙壁。

4.2。力学平衡方程

地质力学的控制方程,也称为力学平衡方程,得到从动量守恒定律: 在哪里 是流体混合物的密度,岩石骨架的密度,是岩石骨架的位移矢量,总应力张量,是身体力量,重力在这项工作。在静态分析中,动态项左边的(25可以省略)。

在多孔介质骨折、矩阵和骨折是两个单独的多孔包含流体的空间,和双孔隙度有效应力法给出如下(35]: 在哪里有效应力张量,矩阵块的平均压力,骨折的平均压力,是毕奥系数矩阵,骨折的毕奥系数,是单位矩阵。

方程(25)是离散使用有限元方法(FEM),提供以下: 在哪里 节点刚度矩阵,是strain-nodal位移矩阵,是弹性刚度矩阵,是节点位移增量向量,是体积刚度向量,和矩阵的平均压力增加和骨折,然后呢是外部荷载增量向量。

弹性刚度矩阵的建立骨折的多孔介质使用等效连续体的方法。局部坐标系的骨折,合规矩阵可以写成: 在哪里是合规矩阵的等效裂隙岩体,和岩石的合规矩阵是矩阵和骨折,是一个系数矩阵,是岩石的体积分数矩阵。

合规矩阵在全球坐标系统从坐标变换可以得到: 在哪里坐标变换矩阵相关断裂方向。

弹性刚度矩阵合规矩阵的反转吗 :

4.3。耦合策略

采用迭代耦合策略。流和地质模块调用顺序在每一个迭代循环。在每个迭代中,首先执行流仿真。压力和饱和度信息流动模拟结果随后被传递到地质模块计算应力、应变和位移的结果。水力特性的变化,如孔隙度和渗透率、更新流仿真之前,下一个迭代的步骤。液压的属性矩阵的更新和骨折分别执行。对矩阵,使用实证关系(36]: 在哪里残余孔隙度,是零压力孔隙度,是平均有效压力,是零应力渗透率,和是更新参数。

除了孔隙度和渗透率的更新功能,毛细管压力的岩石莱弗里特矩阵更新的功能:

更新裂缝的孔隙度和渗透率的孔径变化: 在哪里是零应力裂缝孔隙度和零应力裂缝渗透率。

骨折刚度也更新使用(16)和(19下一个迭代循环之前)。时停止迭代过程收敛流和地质模块都达到标准。然后,下一个时间步的仿真过程的开始。

5。pIEDFM的验证

5.1。与基准测试结果验证

pIEDFM的仿真结果比较的基准测试结果在工作Karimi-Fard et al。37]。如图6,一个简单的裂缝性储层具有水平和垂直裂缝。

矩阵的孔隙度是0.20,渗透率是1。骨折孔径是0.1毫米。块最初饱和油。左下角的水注入。油和水的混合物产生的右上角。油和水的粘度0.45 cP和1.0 cP。岩石矩阵和骨折的相对渗透率曲线是线性的。毛细管压力是被忽视的矩阵和骨折。模拟执行,直到2 PV的水注入。

图7显示了良好的协议之间的累积产油量pIEDFM结果和细网格的结果Karimi-Fard et al .图的工作8显示了在0.1 PV后储层含水饱和度资料,0.3 PV和0.5 PV的水注入。相比之下,EDFM模拟的结果了。指出一个更好的协议达成pIEDFM结果和基准测试结果。水流在骨折EDFM被低估了,而pIEDFM精确再现了含水饱和度剖面和细微的差别。

5.2。验证与细网格模拟

拟议中的pIEDFM进一步验证储层情况下与一个断裂,pIEDFM的表演和传统EDFM进行了比较。储层为0.15米0.15米1.0年的大小和离散3030 pIEDFM和EDFM笛卡尔网格,如图9。矩阵的孔隙度是0.20。矩阵的渗透率是1。骨折孔径是0.5毫米。水库的初始应力是x 10.0 MPay方向和25.0 MPaz方向。

水库最初饱和油的密度为800.0公斤/米³0.45 cP。水的粘度和密度为1000.0公斤/米³和1.0的粘度cP左下角的注入。注射速率为0.1 m³/天。油和水的混合物产生的右上角固定井底压力为10.0 MPa。线性相对渗透率曲线用于矩阵和骨折,毛细管压力是被忽视的。储层的初始压力为10.0 MPa。仿真运行为1800秒。

pIEDFM和EDFM仿真结果与细网格explicit-fracture结果相比,认为是准确的。在细网格explicit-fracture仿真,水库300年离散300年笛卡儿网格。网格的大小等于孔径骨折,骨折是显式地处理一系列高度渗透的网格。1800秒后的含油饱和度剖面注入如图10和用作参考的结果在下面讨论。

在图11的含油饱和度资料pIEDFM和EDFM经过1800秒的注入与细网格含油饱和度剖面。支持的讨论,在细网格上执行一个升级含油饱和度资料。它可以发现的含油饱和度剖面pEDFM显示更好的协议比EDFM细网格的结果。图12礼物的资料绝对pIEDFM含油饱和度误差和EDFM仿真结果与细网格的结果,这也表明pIEDFM优于EDFM再现现实的含油饱和度分布。EDFM结果,含油饱和度的中心周围骨折高于细网格的结果,这表明水断裂位移是低估了。此外,断裂的上边缘周围的含油饱和度低于细网格的结果,这表明水沿裂缝位移是高估了。

仿真场景中重复使用细网格网格(6060)。图13和图14显示出含油饱和度资料和含油饱和度的错误配置文件pIEDFM EDFM仿真后1800秒的注入。网格细化后,pIEDFM饱和度剖面显示一个更好的与细网格的结果。然而,EDFM的改善结果是有限的,这表明网格细化不能有效减少错误造成的非物质的通量分裂EDFM。通过添加附加骨折和nonneighboring基质细胞之间的连接,pIEDFM成功抓住了现实的通量分裂和减少了错误在预测含油饱和度分布。图15显示了流场中的油相细网格,pIEDFM,和EDFM结果,这表明pIEDFM可以获得更现实的速度场比周围的EDFM断裂表面。

执行网格敏感性分析检查pIEDFM的性能在不同的网格密度。的规范代表介绍了含油饱和度的总体误差: 在哪里采样点的数量和吗是高档细网格含油饱和度。图16比较了规范pIEDFM和EDFM在不同网格大小。结果表明,误差pIEDFM随矩阵网格的细化,演示了该方法的收敛性。与细网格pIEDFM有更好的协议结果比EDFM网格大小,显示出该方法的优越性。

证明的改善pIEDFM字段在捕捉现实的压力,压力分布在三片的水库,通过降低边缘骨折,骨折的中心,和骨折的上边缘,分别调查了,如图17。图18显示了压力分布的细网格的三片,结果是明显的非线性。

图19显示错误的分布的油压pIEDFM和EDFM结果与细网格的结果。它可以发现pIEDFM具有显著的优势比EDFM预测压力分布。EDFM,骨折不能被认为是附近非线性压力变化,导致断裂表面附近的明显错误,而pIEDFM更好地捕捉压力急剧变化附近的骨折。

6。应用程序流和地质模拟pIEDFM的耦合

拟议中的pIEDFM应用于耦合流和地质油藏模拟。仿真场景中重复使用的原始EDFM比较水库的地质pIEDFM和EDFM预测的行为。

2 d储层天然裂缝穿透整个储层研究的深度,和图所示的配置20.。原则假定水库的应力方向与轴方向一致,和零内摩擦角被认为是。因此,平均攻击角骨折从轴的方向。骨折的平均长度是10米。水库是由30离散30笛卡尔网格。矩阵的孔隙度是0.15。矩阵的零应力渗透率是5。零应力骨折孔径0.3毫米。初始压力是6.0 MPa在x方向上,y方向和18.0 MPaz方向。在初始应力条件下,初始基质渗透率是2.023 mD,初始裂缝孔径为0.222毫米。

水库最初饱和油的密度为800.0公斤/米³和5.0 cP。水的粘度从注入井注入位于左下角的恒定井底压力20.0 MPa。水的密度是1000.0公斤/ m³,水粘度是0.9 cP。油和水的混合物产生的右上角固定井底压力为5.0 MPa。骨折的相对渗透率是假定为线性的,和Brooks-Corey模型用于矩阵的相对渗透率,如图21。岩石的毛管压力曲线给出了矩阵在图22骨折,而毛细管压力是被忽视的。储层的初始压力为20.0 MPa。列出了力学参数表1。模拟生产运行20000天。

图23比较的含油饱和度资料pIEDFM和EDFM模拟5000年之后,10000年,20000天的生产。含油饱和度资料表明,水流在骨折是pIEDFM更加突出的结果。由于捕捉适当的限制通过骨折通量分裂,EDFM不能代表准确多相断裂位移。

图24显示了生产速度历史生产井的油和水。耦合与非耦合模拟,水生产速率预测的pIEDFM大于EDFM,这可以归因于低估EDFM水跨断裂位移过程的预测。由于更加突出pIEDFM模拟水流在骨折,更多的石油从水库流离失所,导致稳态油率就越高。图25介绍了累积油和水生产历史,表明水驱替过程之间的差异pIEDFM和EDFM不仅影响流动模拟的结果也是地质储层的性能。在15000天的生产,石油产量累计在耦合情况下比非耦合情况下是92.7% pIEDFM EDFM模拟、仿真和93.9%,累计生产水的比例是20.8% pIEDFM EDFM仿真模拟和19.0%。

图26比较了含水pIEDFM和EDFM模拟的历史。在非耦合模拟,水突破时间是6112天pIEDFM EDFM和5904天。前面的水突破的主要原因EDFM模拟可以归因于流沿着骨折的高估超过一半的骨折以来东方向生产井。当地质力学认为,水突破时间将推迟由于岩石的渗透率降低矩阵和关闭骨折。在耦合模拟,水突破时间是10926天pIEDFM EDFM和11157天。相比之下,水突破时间的预测是早些时候pIEDFM模拟,这表明EDFM仿真更受到地质力学的耦合效应的影响。

平均有效压力资料pIEDFM和EDFM模拟后20000天的生产比较图27。平均有效压力资料表明,储层更巩固在EDFM模拟中的损耗,这也是支持的等高线图呈现在图28。

在耦合模拟,有效应力场中的错误导致pIEDFM渗透率分布差异和EDFM结果。图29日礼物的渗透率资料pIEDFM和EDFM模拟生产的20000天,这表明储层更容易渗透减少EDFM模拟。图30.显示了不同渗透率的分布场pIEDFM和EDFM之间的结果,和等高线地图比较图31日。在大多数地区的储层,渗透率pIEDFM结果高于EDFM结果,和主要的差异渗透在生产井附近出现。

明显的延迟水突破时间,有效应力的差异,和渗透率分布表明EDFM可能导致结果的偏差耦合流和地质模拟,和pIEDFM可以有效降低这些偏差。

7所示。结论

提出了一个projection-based整体嵌入离散裂缝模型。pIEDFM,添加额外的连接断裂nonneighboring矩阵元素和元素之间提供一个更现实的表示的通量分裂水流入跨断裂的过程。骨折之间的连接元素的传播性和邻近的矩阵元素之间的连接和传播性骨折元素和附加nonneighboring矩阵元素分别计算使用半锥形压力分布在断裂表面。pIEDFM的准确性验证基准测试结果和explicit-fracture细网格的结果。pIEDFM的表演和EDFM比较和讨论通过单一裂缝注水情况。pIEDFM应用于耦合流和地质模拟,结果与EDFM相比。结论如下:(1)良好的协议达成的饱和概要文件pIEDFM模拟和基准测试结果,这地址的准确性pIEDFM建模多相流过程中断裂的媒体。(2)比较EDFM和pIEDFM结果与explicit-fracture细网格的模拟结果表明,pIEDFM身体可以获得一个更真实的速度场和更好的预测多相流动过程中断裂的水库。pIEDFM的含油饱和度与细网格的结果显示了良好的协议使用适量的网格。(3)pIEDFM具有显著优势在预测比EDFM压力分布。与EDFM结果相比,周围的压力错误pIEDFM断裂表面明显减少,表明非线性压力改变附近的骨折可以被捕捉到。(4)应用pIEDFM耦合流和地质模拟表明,预测水驱替过程之间的差异不仅pIEDFM和EDFM影响生产和饱和度资料的结果也是地质储层的性能。偏见水驱替过程EDFM导致水突破时间的预测偏差,有效应力场,stress-dependent渗透率分布。因此,它承诺将pIEDFM耦合流和地质模拟裂缝储层。

数据可用性

所有的数据都可以发现在我们的手稿。

的利益冲突

作者宣称没有利益冲突有关这项研究的出版物。

确认

这个项目是由国家自然科学基金委的共同基金为企业创新和发展(批准号U19B6003-02-06)和中国国家自然科学基金(批准号51674010)。