一种弱凝胶调驱油藏四维试井解释方法及装置

1.本发明属于弱凝胶调驱油藏试井解释领域，涉及一种弱凝胶调驱油藏四维试井解释方法及装置。


背景技术：

2.渤海油田具有储层类型复杂、非均质性严重，主要采用注水提高采收率。但在注入水长期冲刷作用下，注入水沿高渗通道突进，加剧了层间矛盾。常规化学药剂的注入调驱效果较差。弱凝胶深部调驱技术作为一种改善严重非均质油田开发效果的提高采收率技术，调驱效果较好，被渤海油田广泛使用。
3.聚合物溶液交联后形成网状结构，其溶液中的水分子被交联聚合物空间网络“锁住”，失去了一部分液体的流动性，但获得了一部分类似固体的固定特性，最终形成介于液体和固体之间的水凝胶，即弱凝胶。弱凝胶可以顺水流通道进入高渗岩层，成胶后可以在一定程度上封堵优势孔道，从而改变水驱方向，动用低渗岩层所含原油，以此达到调剖目的；弱凝胶在封堵一段时间后，由于驱替压力的增加，可能会在此孔道失去封堵能力，转而随水流向前运移起到驱油效果，直到运移至更小的孔道又能再次起到封堵效果。
4.弱凝胶调驱施工后，其调驱效果通常依据受效井含水率是否下降进行判定，进而确定下一步调驱方案。但制定调驱方案时，由于弱凝胶非均匀突进，其运移规律不清，封堵位置未知，下一步弱凝胶注入量难以确定。因此，针对非均质性较强的弱凝胶调驱油藏，研发一套调驱效果识别方法可以为调驱方案制定提供理论指导。目前调驱效果识别方法通常依据评估对象(注水井、生产井)分为两类。关于生产井调驱见效评价，常用的方法为水驱曲线评价法、递减曲线评价法、净增油降水法；关于注水井调驱见效评价，通常通过吸水剖面、视吸水指数、霍尔曲线、压力变化来评价。
5.现有方法可以高效评价调驱效果，但由于其不能直接评价弱凝胶分布及封堵位置，使得其运用于调驱方案制定时，存在较大局限性。


技术实现要素：

6.本发明目的在于克服现有技术的不足，提供一种弱凝胶调驱油藏四维试井解释方法及装置，将沿径向非均匀分布的流体粘度考虑到渗流数学模型中，建立了一种弱凝胶调驱油藏四维试井解释方法。解决了通过现有方法不能直接识别注水井波及范围内弱凝胶分布导致下一步调驱方案难以制定等问题。
7.实现本发明目的的技术方案为：
8.本发明的第一方面是提供了一种弱凝胶调驱油藏四维试井解释方法，使用弱凝胶调驱后的试井实测压力数据反演井周探测半径范围内流体参数分布，进而求取调驱后弱凝胶粘度分布；通过弱凝胶粘度与浓度关系的实验数据，用弱凝胶粘度分布反求井周探测半径范围内弱凝胶浓度分布；并通过分布图中弱凝胶浓度峰值，判定弱凝胶调驱封堵的位置。包括以下步骤：
9.步骤一、依据弱凝胶调驱后的试井实测压力数据，由(61)式计算tn时刻探测半径rn对应的从而获得流体参数分布；
10.计算流体参数分布的表达式为：
[0011][0012]
其中tn为第n个时间步，d；rn为探测半径，m；k为储层渗透率，md；h为储层厚度，m；μ为流体粘度，mpa
·
s；n为时间步；δ为微变量。
[0013]
为第n个时间步的流体参数；
[0014]
步骤二、通过可求取调驱后tn时刻探测半径rn对应的弱凝胶粘度μn分布；
[0015]
步骤三、依据步骤二中弱凝胶粘度μn分布曲线，绘制注入井周储层中弱凝胶粘度μn随探测半径rn分布的平面图；
[0016]
步骤四、通过实验测试不同浓度弱凝胶稳定后的粘度，绘制弱凝胶稳定粘度与浓度关系曲线，并回归粘度与浓度间的关系式；
[0017]
步骤五、结合步骤三中求取的调驱后tn时刻探测半径rn对应的弱凝胶粘度μn分布、步骤四中回归的粘度与浓度间的关系式，反算出调驱后地层中弱凝胶的浓度随探测半径rn分布的曲线；
[0018]
步骤六、通过调驱后地层中弱凝胶的浓度分布曲线，绘制出注入井周储层中弱凝胶浓度随探测半径rn分布的平面图，依据该图可以确定储层中弱凝胶浓度大小，明确弱凝胶浓度最大的区域。
[0019]
进一步地，所述探测半径rn的计算公式为：
[0020][0021]
式中各参数的含义：为第n个时间步对应的探测前缘渗透率，md；tn为第n个时间步，d；φ为孔隙度，小数；μ为流体粘度，mpa
·
s；c
t
为综合压缩系数，小数。
[0022]
本发明的第二方面是提供了实现上述弱凝胶调驱油藏四维试井解释方法的装置，包括：
[0023]
流体参数计算模块，依据弱凝胶调驱后的试井实测压力数据，计算tn时刻探测半径rn对应的从而获得流体参数分布；
[0024]
粘度分布计算模块，求取调驱后tn时刻探测半径rn对应的弱凝胶粘度μn分布，绘制注入井周储层中弱凝胶粘度μn随探测半径rn分布的平面图；
[0025]
粘度与浓度换算模块，通过实验测试不同浓度弱凝胶稳定后的粘度，绘制弱凝胶稳定粘度与浓度关系曲线，并回归粘度与浓度间的关系式；
[0026]
浓度分布计算模块，结合求取的调驱后tn时刻探测半径rn对应的弱凝胶粘度μn分布、步骤四中回归的粘度与浓度间的关系式，反算出调驱后地层中弱凝胶的浓度随探测半径rn分布的曲线；通过调驱后地层中弱凝胶的浓度分布曲线，绘制出注入井周储层中弱凝胶浓度随探测半径rn分布的平面图；
[0027]
封堵位置计算模块，依据弱凝胶浓度随探测半径rn分布的平面图确定储层中弱凝胶浓度大小，明确弱凝胶浓度最大的区域。
[0028]
本发明的第三方面是提供了一种电子设备，包括存储器和处理器；其中：
[0029]
存储器：用于存储处理器可执行的指令；
[0030]
处理器：所述处理器被配置执行：使用弱凝胶调驱后的试井实测压力数据反演井周探测半径范围内流体参数分布，进而求取调驱后弱凝胶粘度分布；通过弱凝胶粘度与浓度关系的实验数据，用弱凝胶粘度分布反求井周探测半径范围内弱凝胶浓度分布；并通过分布图中弱凝胶浓度峰值，判定弱凝胶调驱封堵的位置。
[0031]
本发明的第四方面是提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序用于使所述计算机执行所述的弱凝胶调驱油藏四维试井解释方法。
[0032]
本发明的优点和有益效果：
[0033]
依据本发明弱凝胶调驱油藏四维试井解释方法，结合弱凝胶粘度与浓度实验数据及试井资料，可以解释注入井波及范围内弱凝胶分布，该方法可以为下一步调驱方案制定提供理论指导。
附图说明
[0034]
图1为注入井压力及压力导数曲线图；
[0035]
图2为注入井附近储层弱凝胶粘度与半径关系反演曲线图；
[0036]
图3为注入井附近储层弱凝胶粘度分布场图；
[0037]
图4为实验测试不同浓度弱凝胶粘度随时间变化曲线图；
[0038]
图5为实验测试弱凝胶稳定粘度与浓度回归关系图；
[0039]
图6为注入井附近储层弱凝胶浓度与半径关系反演曲线图；
[0040]
图7为注入井附近储层弱凝胶粘度分布场图。
具体实施方式
[0041]
下面通过具体实施例对本发明作进一步详述，以下实施例只是描述性的，不是限定性的，不能以此限定本发明的保护范围。
[0042]
本发明提供了一种弱凝胶调驱油藏四维试井解释方法，该方法包括如下步骤：
[0043]
步骤一、设定在无限大径向非均质油藏中一口井定产量生产，上、下边界封闭，储层流体流度是径向距离的函数，地层孔隙度、厚度及岩石压缩系数、流体压缩系数为常数，忽略重力和毛管力影响。
[0044]
步骤二、建立无因次渗流数学模型。依据步骤一中假设条件，建立无限大径向非均质油藏中一口井定产量生产的无因次渗流数学模型如下：
[0045][0046]
初始条件如下：
[0047][0048]
内边界条件如下：
[0049][0050]
外边界条件如下：
[0051][0052]
式中：rd为无因次径向距离；kd(rd)为无因次地层渗透率，是无因次径向距离的函数；pd为无因次地层压力；td为无因次时间；pd(rd,td)为无因次地层压力，是无因次径向距离与无因次时间的函数。
[0053]
步骤三、定义无因次渗流数学模型中的无因次变量。步骤二中无因次渗流数学模型的无因次变量定义如下：
[0054]
无因次时间：
[0055][0056]
无因次半径：
[0057][0058]
无因次压力：
[0059][0060]
无因次井底压力：
[0061][0062]
无因次井底压力导数：
[0063][0064]
无因次径向渗透率
[0065][0066]
其中，为参考渗透率。
[0067]
式中：为参考渗透率，md；t为时间，d；φ为孔隙度，小数；c
t
为综合压缩系数，小数；μ为天然气粘度，mpa
·
s；rw为井筒半径，m；r为径向距离，m；h为储层厚度，m；pi为初始地
层压力，mpa；p(r,t)为地层压力，是径向距离与时间的函数，mpa；q为天然气产量，m3/d；b为天然气地层体积系数，小数；p
wd
为无因次井底流压；pw为井底流压，mpa；p’wd
为无因次井底流压导数；k(r)为地层渗透率，是径向距离的函数，md。
[0068]
步骤四、求解渗流数学模型。利用扰动理论和拉普拉斯理论求解径向非均质渗流扩散方程组，得到任意渗透率分布下的近似压力和压力导数的解析解。
[0069]
首先定义如下渗透率：
[0070][0071]
式中ε是一小量，代表了地层渗透率在参考渗透率附近作较小的波动。
[0072]
将无因次扩散方程写成算子的形式：
[0073][0074]
根据拉普拉斯理论，式(12)的拉普拉斯变换为：
[0075][0076]
根据扰动理论，拉氏空间中的井底压力解可以表示为：
[0077][0078]
该方程表示渗透率和角度不发生变化的均质油藏的扩散方程，为无因次均质储层压力在拉氏空间中的零阶解。利用van everdingen和hurst的研究，零阶解可以表示为：
[0079][0080]
式中：s—拉氏变量，顶部标有“～”符号的字符表示拉氏空间变量
[0081]
k0是第二类零阶修正贝塞尔函数；
[0082]
k1是第二类一阶修正贝塞尔函数。
[0083]
根据平面非均质数学模型的分析，可以方便地推导出无因次第一次扰动压力在拉氏空间中的一阶解
[0084][0085]
当rd＝1时，得到井底压力：
[0086][0087]
利用下列关系简化ψ
d1
的表达式：
[0088][0089][0090]
令t
→
∞(s
→
0)，即当时间足够长时，通过拉氏数值反演，可以得到真实空间中的无因次井底压力和无因次井底压力导数的近似解：
[0091][0092][0093]
式中：s为拉氏变量，顶部标有“～”符号的字符表示拉氏空间变量；k0是第二类零阶修正贝塞尔函数；k1是第二类一阶修正贝塞尔函数。
[0094]
步骤五、建立出弱凝胶调驱前后储层粘度等参数场的反演算法。
[0095]
定义k(r)为位置r的k(r，θ)的调和平均：
[0096][0097]
无因次变量定义为：
[0098][0099][0100][0101][0102]
oliver的解可以表示为：
[0103][0104]
变形后的核函数：
[0105][0106]
尽管oliver在利用扰动理论推导时假设渗透率在参考渗透率k附近微小扰动，但研究表明，当渗透率变化很大时，这种扰动解仍可能给出合理的近似解。
[0107]
当时，即在均质油藏条件下，由(27)式可以得到：
[0108][0109]
另一种情况是渗透率为常数k1的均质油藏时，(29)式变为：
[0110]
依据试井曲线特征，径向流阶段的导数值大小与渗透率呈正相关关系，因此，压力导数曲线特征满足如下关系式
[0111][0112]
由(27)和(30)式可以导出：
[0113][0114]
为下面推导的方便，令k
1*
(rd，td)＝2k
*
(rd，td)，则：
[0115][0116]
将(32)式代入(27)式，则得到：
[0117][0118]
由(25)、(26)式无因次压力导数又可以表示为：
[0119][0120]
结合(26)、(33)和(34)式，得到
[0121][0122]
从(35)式可以看出，参考渗透率仅隐含在无因次时间中，结合(35)式的特点重新
定义参考流体参数，记为
[0123][0124]
由于与某一特定时刻的压力导数相关，称为瞬时流体参数。
[0125]
由(35)式可以推导出
[0126][0127]
(37)式表明，瞬时流体参数是带权重函数的调和平均，因而对于流体参数随径向距离变化的地层平均流体参数(上面的推导过程中，认为流体参数是径向距离r的函数，不考虑随角度θ的变化)的正确定义应该运用调和平均。
[0128]
由于包含了更多的油藏信息，研究表明，如果用代替无因次时间td中的由(37)式得出的结果会更加准确。
[0129]
定义无因次拟时间：
[0130][0131]
为无因次拟时间，为探测前缘渗透率，t为时间，rw为井筒半径。
[0132]
则(37)式变为：
[0133][0134][0135]
为了从中求得渗透率k(r)，假设k(r)可以表示为离散的分段函数，即当r
i-1
＜r≤ri(i＝1，2，
…
，n)，k(r)＝ki。这样，这个反演算法的实质就是求解在tn时刻探测半径rn对应下的探测半径rn定义为：
[0136]
[0137]
为第n个时间步对应的探测前缘渗透率，tn为第n个时间步，φ为孔隙度，μ为流体粘度，c
t
为综合压缩系数。
[0138]
定义：
[0139][0140][0141]
在t＝tn时刻，(39)式变形为：
[0142][0143]
由(32)式，可以得到：
[0144][0145]
从而(44)式可以改写为：
[0146][0147]
由(42)式可以得到：
[0148][0149]
将(47)式代入(46)式，得到：
[0150][0151]
可以首先求出
[0152][0153]
(49)式中的t1表示存在可靠的压力导数值的第一时间点，再利用(3-40)式可求得(48)式可变形为：
[0154][0155]
因为
[0156][0157]
(50)式中的积分式可表达为：
[0158]
[0159]
根据无因次boltzmann变量定义，相应的瞬时无因次boltzmann变量定义为：
[0160][0161]
由(52)式，得：
[0162][0163]
将(53)式代入(51)式，得到：
[0164][0165]
oliver研究表明，在区域以外的地方核函数k
*
(rd,td)≈0。这里定义如下的oliver无因次探测内半径：
[0166][0167]
相应的外探测半径为：
[0168][0169]
对于任意函数f(rd)，存在如下关系式：
[0170][0171]
(57)式的精度取决于函数的表达式。研究证实，当函数f(rd)是常数时(57)式有足够高的精度。
[0172]
综合(52)～(57)式的分析，(51)式可以进一步表示为：
[0173][0174]
若令：
[0175][0176]
则(58)式可以简写为：
[0177][0178]
根据(42)式的定义，计算流体参数分布的表达式为：
[0179][0180]
tn为第n个时间步，rn为探测半径，k为储层渗透率，h为储层厚度，μ为流体粘度，n为时间步，δ为微变量；
[0181]
步骤六、利用弱凝胶调驱后的试井实测压力数据，由(34)式计算瞬时流体参数
[0182]
步骤七、依据(41)式，计算探测半径rn。
[0183]
步骤八、结合图1(注入井压力及压力导数曲线)横坐标对应的时间，由(38)式计算无因次拟时间。
[0184]
步骤九、由(61)式计算流体参数便可以得到tn时刻探测半径rn对应的从而获得流体参数分布。
[0185]
步骤十、通过可求取调驱后tn时刻探测半径rn对应的弱凝胶粘度μn分布(图2)。
[0186]
步骤十一、依据步骤十中弱凝胶粘度μn分布曲线，绘制注入井周储层中弱凝胶粘度μn随探测半径rn分布的平面图(图3)。
[0187]
步骤十二、通过实验测试不同浓度弱凝胶稳定后的粘度，本次实验测试了弱凝胶浓度分别为1000～4000mg/l对应的弱凝胶稳定后的粘度(图4)。
[0188]
步骤十三、依据步骤十二中实验测试的不同浓度弱凝胶稳定后的粘度，绘制了弱凝胶稳定粘度与浓度关系曲线，并回归了粘度与浓度间的关系式(图5)。
[0189]
步骤十四、结合步骤十中求取的调驱后tn时刻探测半径rn对应的弱凝胶粘度μn分布(图2)、步骤十三中回归的粘度与浓度间的关系式(图5)，即可反算出调驱后地层中弱凝胶的浓度随探测半径rn分布的曲线(图6)。
[0190]
步骤十五、通过调驱后地层中弱凝胶的浓度分布曲线(图6)，可绘制出注入井周储层中弱凝胶浓度随探测半径rn分布的平面图(图7)。依据该图中图例，可以确定储层中弱凝胶浓度大小，即可以明确红色区域是弱凝胶封堵效果最好的区域，即为距离井16～30m左右
的区域。
[0191]
本发明进一步提供了一种实现上述弱凝胶调驱油藏四维试井解释方法的装置，包括：
[0192]
流体参数计算模块，依据弱凝胶调驱后的试井实测压力数据，计算tn时刻探测半径rn对应的从而获得流体参数分布；
[0193]
粘度分布计算模块，求取调驱后tn时刻探测半径rn对应的弱凝胶粘度μn分布，绘制注入井周储层中弱凝胶粘度μn随探测半径rn分布的平面图；
[0194]
粘度与浓度换算模块，通过实验测试不同浓度弱凝胶稳定后的粘度，绘制弱凝胶稳定粘度与浓度关系曲线，并回归粘度与浓度间的关系式；
[0195]
浓度分布计算模块，结合求取的调驱后tn时刻探测半径rn对应的弱凝胶粘度μn分布、步骤四中回归的粘度与浓度间的关系式，反算出调驱后地层中弱凝胶的浓度随探测半径rn分布的曲线；通过调驱后地层中弱凝胶的浓度分布曲线，绘制出注入井周储层中弱凝胶浓度随探测半径rn分布的平面图；
[0196]
封堵位置计算模块，依据弱凝胶浓度随探测半径rn分布的平面图确定储层中弱凝胶浓度大小，明确弱凝胶浓度最大的区域。
[0197]
本发明进一步提供了一种电子设备，包括存储器和处理器；其中：
[0198]
存储器：用于存储处理器可执行的指令；
[0199]
处理器：所述处理器被配置执行：使用弱凝胶调驱后的试井实测压力数据反演井周探测半径范围内流体参数分布，进而求取调驱后弱凝胶粘度分布；通过弱凝胶粘度与浓度关系的实验数据，用弱凝胶粘度分布反求井周探测半径范围内弱凝胶浓度分布；并通过分布图中弱凝胶浓度峰值，判定弱凝胶调驱封堵的位置。
[0200]
本发明进一步提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序用于使所述计算机执行上述弱凝胶调驱油藏四维试井解释方法。
[0201]
以上所述的仅是本发明的优选实施方式，应当指出，对于本领域的普通技术人员来说，在不脱离发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。

技术特征：
1.一种弱凝胶调驱油藏四维试井解释方法，其特征在于，使用弱凝胶调驱后的试井实测压力数据反演井周探测半径范围内流体参数分布，进而求取调驱后弱凝胶粘度分布；通过弱凝胶粘度与浓度关系的实验数据，用弱凝胶粘度分布反求井周探测半径范围内弱凝胶浓度分布；并通过分布图中弱凝胶浓度峰值，判定弱凝胶调驱封堵的位置。2.根据权利要求1所述的弱凝胶调驱油藏四维试井解释方法，其特征在于，包括以下步骤：步骤一、依据弱凝胶调驱后的试井实测压力数据，由(61)式计算t
n
时刻探测半径r
n
对应的从而获得流体参数分布；计算流体参数分布的表达式为：其中：t
n
为第n个时间步，r
n
为探测半径，k为储层渗透率，h为储层厚度，μ为流体粘度，n为时间步，δ为微变量；为第n个时间步的流体参数；步骤二、通过可求取调驱后t
n
时刻探测半径r
n
对应的弱凝胶粘度μ
n
分布；步骤三、依据步骤二中弱凝胶粘度μ
n
分布曲线，绘制注入井周储层中弱凝胶粘度μ
n
随探测半径r
n
分布的平面图；步骤四、通过实验测试不同浓度弱凝胶稳定后的粘度，绘制弱凝胶稳定粘度与浓度关系曲线，并回归粘度与浓度间的关系式；步骤五、结合步骤三中求取的调驱后t
n
时刻探测半径r
n
对应的弱凝胶粘度μ
n
分布、步骤四中回归的粘度与浓度间的关系式，反算出调驱后地层中弱凝胶的浓度随探测半径r
n
分布的曲线；步骤六、通过调驱后地层中弱凝胶的浓度分布曲线，绘制出注入井周储层中弱凝胶浓度随探测半径r
n
分布的平面图，依据该图可以确定储层中弱凝胶浓度大小，明确弱凝胶浓度最大的区域。3.根据权利要求2所述的弱凝胶调驱油藏四维试井解释方法，其特征在于，所述探测半径r
n
的计算公式为：其中：为第n个时间步对应的探测前缘渗透率，t
n
为第n个时间步，φ为孔隙度，μ为流体粘度，c
t
为综合压缩系数。4.根据权利要求1～3任一权利要求所述的弱凝胶调驱油藏四维试井解释方法的装置，其特征在于，包括：流体参数计算模块，依据弱凝胶调驱后的试井实测压力数据，计算t
n
时刻探测半径r
n
对
应的从而获得流体参数分布；粘度分布计算模块，求取调驱后t
n
时刻探测半径r
n
对应的弱凝胶粘度μ
n
分布，绘制注入井周储层中弱凝胶粘度μ
n
随探测半径r
n
分布的平面图；粘度与浓度换算模块，通过实验测试不同浓度弱凝胶稳定后的粘度，绘制弱凝胶稳定粘度与浓度关系曲线，并回归粘度与浓度间的关系式；浓度分布计算模块，结合求取的调驱后t
n
时刻探测半径r
n
对应的弱凝胶粘度μ
n
分布、步骤四中回归的粘度与浓度间的关系式，反算出调驱后地层中弱凝胶的浓度随探测半径r
n
分布的曲线；通过调驱后地层中弱凝胶的浓度分布曲线，绘制出注入井周储层中弱凝胶浓度随探测半径r
n
分布的平面图；封堵位置计算模块，依据弱凝胶浓度随探测半径r
n
分布的平面图确定储层中弱凝胶浓度大小，明确弱凝胶浓度最大的区域。5.一种电子设备，其特征在于，包括存储器和处理器；其中：存储器：用于存储处理器可执行的指令；处理器：所述处理器被配置执行：使用弱凝胶调驱后的试井实测压力数据反演井周探测半径范围内流体参数分布，进而求取调驱后弱凝胶粘度分布；通过弱凝胶粘度与浓度关系的实验数据，用弱凝胶粘度分布反求井周探测半径范围内弱凝胶浓度分布；并通过分布图中弱凝胶浓度峰值，判定弱凝胶调驱封堵的位置。6.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序用于使所述计算机执行权利要求1-3的任一项所述的弱凝胶调驱油藏四维试井解释方法。

技术总结
本发明公开了一种弱凝胶调驱油藏四维试井解释方法及装置，该方法是基于多孔介质中流体渗流理论，依据数学物理方法，建立并求解的一种弱凝胶调驱油藏四维试井正演及反演算法。依据该方法，可以使用弱凝胶调驱后的试井实测压力数据反演井周探测半径范围内流体参数分布，进而求取调驱后弱凝胶粘度分布。通过弱凝胶粘度与浓度关系的实验数据，可以用弱凝胶粘度分布反求井周探测半径范围内弱凝胶浓度分布。并通过分布图中弱凝胶浓度峰值，可以判定弱凝胶调驱封堵的位置。本发明建立的弱凝胶调驱油藏四维试井正演及反演算法，解决了弱凝胶调驱油藏封堵位置的弱凝胶浓度分布难以确定的问题，为下一步调驱方案制定提供理论依据。为下一步调驱方案制定提供理论依据。为下一步调驱方案制定提供理论依据。


技术研发人员：冯青 李成勇 李胜胜 胡鹏 刘蜀东 孙艳妮 张小蓉
受保护的技术使用者：成都理工大学
技术研发日：2023.06.20
技术公布日：2023/10/15