CN105715241B - 一种聚合物驱相对渗透率曲线的测量方法 - Google Patents
一种聚合物驱相对渗透率曲线的测量方法 Download PDFInfo
- Publication number
- CN105715241B CN105715241B CN201610024682.9A CN201610024682A CN105715241B CN 105715241 B CN105715241 B CN 105715241B CN 201610024682 A CN201610024682 A CN 201610024682A CN 105715241 B CN105715241 B CN 105715241B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- data
- polymer
- curve
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
一种聚合物驱相对渗透率曲线的测量方法，包括：首先进行聚合物岩心驱替实验，得到不同时刻的生产实验数据，然后以得到的生产实验数据作为拟合数据，建立拟合目标函数和聚合物驱相对渗透率曲线模型，最后通过数值模拟器，利用自动历史拟合算法不断调整相对渗透率曲线模型参数，直至拟合目标函数值达到误差允许范围内，得到最优的聚合物驱相对渗透率曲线，该方法简单、易操作，为研究聚合物驱的渗流规律提供帮助。
Description
技术领域
本发明涉及一种聚合物驱相对渗透率曲线的测量方法，属于石油开发的技术领域。
背景技术
聚合物驱是一种重要的三次采油技术，该技术将水溶性的高分子聚合物加入注入水中作为驱油剂，增加注入水相的粘度，达到降低水相渗透率，改善油水流度比，提高原油采收率的目的。目前，聚合物驱油技术在矿场广泛应用。其中，聚合物驱相对渗透率曲线是聚合物驱油藏分析以及数值模拟研究的基础，准确获取聚合物驱相对渗透率曲线尤为重要。
目前，主要是通过传统非稳态法或直接降低水相相对渗透率的方法获得聚合物驱的油水相对渗透率曲线。但在实际应用过程中，由于实验条件无法满足相对渗透率曲线数学模型的假设条件，按照传统非稳态法测得的实验数据难以真实反映聚合物驱油的渗流规律，因此计算出的聚合物驱相对渗透率曲线可靠性差。此外，通过直接降低水相相对渗透率获取聚合物驱相对渗透率曲线的方法，忽略了聚合物驱过程中的影响因素，因而误差很大，实用性差。综上可知目前尚没有一种可靠的聚合物驱相对渗透率曲线的获取方法，因此有必要提出一种获取更加真实、准确的聚合物驱相对渗透率曲线的方法，为研究聚合物驱油提供技术支持。
发明内容
针对现有技术的不足，本发明提供一种聚合物驱相对渗透率曲线的测量方法，其在聚合物驱替实验的基础上，结合油藏数值模拟器对油水相对渗透率进行准确计算。
本发明的技术方案如下：
一种聚合物驱相对渗透率曲线的测量方法，过程具体如下：
a.测量聚合物驱岩心模型的基本参数，包括岩心长度、岩心直径、孔隙度、渗透率、束缚水饱和度、残余油饱和度；测量所用聚合物的基本物化参数，包括聚合物溶液粘度、吸附量、残余阻力系数、不可及孔隙体积。
b.通过现有的聚合物驱替实验，记录生产试验数据包括注入端压力数据、累积产油数据和产出聚合物浓度数据，并将实验数据进行归一化处理：
在式(I)中，i＝123...N，为实验数据序号，N为实验数据个数；l代表数据类型，l＝p时代表压力，l＝o时代表累积产油，l＝cp代表产出聚合物浓度；Dl为数据l经归一化处理后的值；De,l为数据l的实验测量值；De,max,l为数据l的实验测量的最大值；De,min,l为数据l的实验测量的最小值。
c.以聚合物驱替实验得到的不同时刻对应的注入端压力数据、累积产油数据和产出聚合物浓度数据作为动态拟合数据，建立目标函数，具体为：
在式(II)中，P为建立的目标函数；为拟合参数向量；Dpred,l为由数值模拟器计算得到的数据l的模型预测值；为l的数据权重。
其中，依据主观赋权法，根据实际问题中数据l的重要程度确定，这里给出各数据权重的取值范围为：且满足
d.利用幂律模型作为聚合物驱相对渗透率曲线的表征模型，建立相对渗透率曲线数学模型，具体为：
在式(III)中，krw为水相相对渗透率；krow为油相相对渗透率；Sw为含水饱和度；Swc为束缚水饱和度；Sorw为残余油饱和度；krwro为残余油饱和度下的水相相对渗透率；krocw为束缚水饱和度下的油相相对渗透率；nw为水相相对渗透率曲线的指数参数；no油相相对渗透率曲线的指数参数。
其中，Swc和Sorw均可由聚合物驱替实验获得，krwro，krocw，nw，no均为需要拟合的参数，即拟合的参数向量为拟合参数的取值范围为：krwro＝[0.1,1]，krocw＝[0.1,1]，nw＝[1.5,4.5]，no＝[1.5,4.5]。
e.利用聚合物驱替实验测定的模型基本参数和初始相对渗透率曲线建立岩心数学模型，调用数值模拟器计算生产动态数据：包括注入端压力、累积产油量和产出聚合物浓度，并进行归一化处理。其中，初始相对渗透率曲线通过给出步骤d中拟合的参数初始值，利用相对渗透率曲线模型得到。
f.计算目标函数值，通过自动历史拟合算法调整拟合参数向量得到新的相对渗透率曲线，循环迭代，直至目标函数值达到预设的允许范围，拟合结束，得到最优的聚合物驱相对渗透率曲线。选用的自动历史拟合为阻尼最小二乘法，属现有算法，具体实现过程可参见文献：S.Qadeer，K.Dehghani，D.O.gbe.Correcting Oil/Water RelativePermeability Data for Capillary End Effect in Displacement Experiment[J].SPE17423,1988。
根据本发明优选的，步骤e中所述预设的允许范围取为0～1×10-5。
本发明的优点在于：
本发明所述一种聚合物驱相对渗透率曲线的测量方法，在聚合物驱替实验的基础上，结合数值模拟器对相对渗透率曲线进行修正，准确性强，实用性好，为研究聚合物驱油提供技术支持。
附图说明
图1：注入端压力和累积产油量拟合效果；
图2：产出聚合物浓度拟合效果；
图3：计算得到的相对渗透率曲线。
具体实施方式
为进一步公开本发明的技术方案，下面结合附图1、2、3对本发明的实施例作详细说明，但不限于此。
实施例1：测量聚合物驱岩心模型的基本参数，包括岩心长度、岩心直径、孔隙度、渗透率、束缚水饱和度、残余油饱和度；测量所用聚合物的基本物化参数，包括聚合物溶液粘度、吸附量、残余阻力系数、不可及孔隙体积。
将岩心装入岩心夹持器中，利用恒温箱将其恒温(45℃)8小时，设定围压为3.5Mpa, 回压为300Kpa,用分子量为1500╳104的聚合物配成浓度为1000mg/L的实验用聚合物溶液，进行实验操作，测量岩心以及聚合物的基本实验参数，测量结果见表1。
表1基本参数
| 参数名称 | 参数值 | 参数名称 | 参数值 |
| 岩心长度，cm | 30 | 聚合物溶液粘度，mPa.s | 15.3 |
| 岩心直径，cm | 3.8 | 聚合物吸附量，μg*g-1 | 86.45 |
| 孔隙度 | 0.315 | 聚合物不可及孔隙体积 | 0.18 |
| 渗透率，md | 1064 | 聚合物残余阻力系数 | 2.51 |
实施例2：进行聚合物驱替实验，记录生产试验数据。具体过程为：
a.将实施例1中实验结束后的岩心进行冲洗，饱和地层水，
b.以驱替速度由小及大的方式进行油驱水，直至产水量不再增加，计算出束缚水饱和度Swc＝0.285；
c.以1ml/min的恒定驱替速度进行聚合物驱油过程，记录不同时刻的注入端压力以及对应时刻的累积产油量和产出聚合物浓度，产油量不再增加时驱替过程结束，计算出聚合物驱残余油饱和度Sorw＝0.215。
实施例3：使用实施例2的生产动态数据进行自动历史拟合，得出最优的聚合物驱相对渗透率曲线。具体步骤如下：
a.首先将生产动态数据进行归一化处理，以注入端压力为例进行归一化。注入端压力最大值De,maxp＝303.62KPa，De,min,p＝296.33KPa，第10个压力数据De,p(10)＝300.96KPa，则归一化处理为：
利用同样的方法对每个生产动态进行归一化处理后用于后续计算；
b.利用归一化后的生产动态数据建立目标函数，具体为：
其中，
c.以幂律模型作为聚合物驱油对渗透率曲线的表征模型，建立相对渗透率曲线数学模型，具体为：
其中，Swc＝0.285，Sorw＝0.215，krwro，krocw，nw，no为待拟合参数。
d.利用聚合物驱替实验测定的基本参数和初始相对渗透率曲线模型参数值建立岩心数学模型，调用数值模拟器计算动态数据：包括注入端压力和累积产油量，并进行归一化处理。其中，给出的参数初始值为：krwro＝0.5，krocw＝1.0，nw＝2，no＝2。
e.利用步骤a中的测得的生产动态数据和步骤d中得到的生产动态数据，计算步骤b中的目标函数值，通过自动历史拟合算法调整拟合参数向量得到新的相对渗透率曲线，循环迭代，直至目标函数值达到预设的允许范围，拟合结束，得到最优的聚合物驱相对渗透率曲线。选用的自动历史拟合为阻尼最小二乘法，属现有算法，具体实现过程可参见文献：S.Qadeer，K.Dehghani，D.O.gbe.Correcting Oil/Water Relative Permeability Datafor Capillary End Effect in Displacement Experiment[J].SPE17423,1988。
f.本例中的允许误差范围为小于1╳10-5，经过64次迭代，目标函数值达到0.98╳10-5，在允许的误差范围之内，迭代过程结束。最终的相对渗透率曲线模型参数值为：
krwro＝0.28，krocw＝0.88，nw＝2.36，no＝3.42
注入端压力和累积产油量拟合效果如图1所示，产出聚合物浓度的拟合效果如图2所示，最终得到的聚合物驱油、水相相对渗透率曲线如图3所示。
Claims (4)
1.一种聚合物驱相对渗透率曲线的测量方法，其特征在于，该方法过程具体为：
a.测量聚合物驱岩心模型的基本参数，包括岩心长度、岩心直径、孔隙度、渗透率、束缚水饱和度、残余油饱和度；测量所用聚合物的基本物化参数，包括聚合物溶液粘度、吸附量、残余阻力系数、不可及孔隙体积；
b.通过现有的聚合物驱替实验，记录生产试验数据包括注入端压力数据、累积产油数据和产出聚合物浓度数据，并将实验数据进行归一化处理：
在式(I)中，i＝1，2，3...N，为实验数据序号，N为实验数据个数；l代表数据类型，l＝p时代表压力，l＝o时代表累积产油，l＝cp代表产出聚合物浓度；Dl为数据l经归一化处理后的值；De,l为数据l的实验测量值；De,max,l为数据l的实验测量的最大值；De,min,l为数据l的实验测量的最小值；
c.以聚合物驱替实验得到的不同时刻对应的注入端压力数据、累积产油数据和产出聚合物浓度数据作为动态拟合数据，建立目标函数，具体为：
在式(II)中，P为建立的目标函数；为拟合参数向量；Dpred,l为由数值模拟器计算得到的数据l的模型预测值；为l的数据权重；
d.利用幂律模型作为聚合物驱相对渗透率曲线的表征模型，建立相对渗透率曲线数学模型，具体为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>r</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>r</mi>
<mi>w</mi>
<mi>r</mi>
<mi>o</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>S</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>w</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>o</mi>
<mi>r</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>w</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>n</mi>
<mi>w</mi>
</msub>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>c</mi>
<mi>w</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>o</mi>
<mi>r</mi>
<mi>w</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>w</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>o</mi>
<mi>r</mi>
<mi>w</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>n</mi>
<mi>o</mi>
</msub>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>I</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
在式(III)中，krw为水相相对渗透率；krow为油相相对渗透率；Sw为含水饱和度；Swc为束缚水饱和度；Sorw为残余油饱和度；krwro为残余油饱和度下的水相相对渗透率；krocw为束缚水饱和度下的油相相对渗透率；nw为水相相对渗透率曲线的指数参数；no油相相对渗透率曲线的指数参数，其中，Swc和Sorw均可由聚合物驱替实验获得，krwro，krocw，nw，no均 为需要拟合的参数，即拟合的参数向量为
e.利用聚合物驱替实验测定的模型基本参数和初始相对渗透率曲线建立岩心聚合物驱数学模型，调用数值模拟器计算生产动态数据，所述生产动态数据包括注入端压力数据、累积产油量数据和产出聚合物浓度数据，所述初始相对渗透率曲线通过给出步骤c中拟合的参数初始值，利用相对渗透率曲线模型得到；
f.计算目标函数值，通过自动历史拟合算法调整拟合参数向量得到新的相对渗透率曲线，循环迭代，直至目标函数值达到预设的允许范围，拟合结束，得到最优的聚合物驱相对渗透率曲线。
2.根据权利要求1所述的一种聚合物驱相对渗透率曲线的测量方法，其特征在于，所述步骤c中依据主观赋权法，根据实际问题中数据l的重要程度确定，这里给出各数据权重的取值范围为：且满足
3.根据权利要求1所述的一种聚合物驱相对渗透率曲线的测量方法，其特征在于，所述步骤d中拟合的参数的取值范围为；krwro＝[0.1,1]，krocw＝[0.1,1]，nw＝[1.5,4.5]，no＝[1.5,4.5] 。
4.根据权利要求1所述的一种聚合物驱相对渗透率曲线的测量方法，其特征在于，所述步骤f中，目标函数值预设的允许范围为0～10-5。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610024682.9A CN105715241B (zh) | 2016-01-15 | 2016-01-15 | 一种聚合物驱相对渗透率曲线的测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610024682.9A CN105715241B (zh) | 2016-01-15 | 2016-01-15 | 一种聚合物驱相对渗透率曲线的测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105715241A CN105715241A (zh) | 2016-06-29 |
| CN105715241B true CN105715241B (zh) | 2017-11-17 |
Family
ID=56147107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610024682.9A Active CN105715241B (zh) | 2016-01-15 | 2016-01-15 | 一种聚合物驱相对渗透率曲线的测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105715241B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106055843B (zh) * | 2016-07-01 | 2019-05-14 | 北京金士力源科技有限公司 | 一种聚驱过程中产聚曲线的量化分析方法及装置 |
| CN106932324B (zh) * | 2017-03-10 | 2019-09-17 | 中国石油化工股份有限公司 | 一种确定高含水砂岩油藏储层渗透率变化规律的方法 |
| CN109357967B (zh) * | 2018-10-24 | 2019-05-28 | 中国石油大学(华东) | 一种油田污水配制的聚合物溶液表观粘度损失量化评价方法 |
| CN109632604B (zh) * | 2019-01-04 | 2021-06-15 | 中国海洋石油集团有限公司 | 一种孔隙尺度到岩心尺度聚合物驱相对渗透率粗化方法 |
| CN110439518B (zh) * | 2019-08-14 | 2020-05-12 | 中国石油大学(华东) | 三元复合驱替效果定量化评价方法 |
| CN114167030A (zh) * | 2020-09-10 | 2022-03-11 | 中国石油天然气股份有限公司 | 泥质砂岩储层的含水率的确定方法、装置和设备 |
| CN112145139B (zh) * | 2020-10-26 | 2022-02-01 | 西南石油大学 | 一种原油黏度合理映射聚合物驱使用浓度的方法 |
| CN112081588B (zh) * | 2020-10-26 | 2022-02-15 | 西南石油大学 | 一种判断相对渗透率曲线中化学体系黏度变化的方法 |
| CN114508331B (zh) * | 2022-01-20 | 2022-08-12 | 东北石油大学 | 一种基于霍尔曲线导数确定聚合物注入能力新方法 |
| CN117592387A (zh) * | 2023-05-25 | 2024-02-23 | 中国石油大学(北京) | 低渗致密油藏浸润调控渗流规律表征方法、装置及设备 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103573234B (zh) * | 2012-08-06 | 2016-06-08 | 中国石油化工股份有限公司 | 一种确定完整油水相对渗透率曲线的方法 |
| US9341557B2 (en) * | 2012-11-14 | 2016-05-17 | Kuwait Oil Company (K.S.C.) | Method and system for permeability calculation using production logs for horizontal wells, using a downhole tool |
| CN104834807B (zh) * | 2015-03-18 | 2017-11-03 | 成都北方石油勘探开发技术有限公司 | 一种基于分形理论的应力敏感储层相对渗透率计算方法 |
-
2016
- 2016-01-15 CN CN201610024682.9A patent/CN105715241B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105715241A (zh) | 2016-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105715241B (zh) | 一种聚合物驱相对渗透率曲线的测量方法 | |
| CN104568694B (zh) | 一种致密岩心气‑水相对渗透率的测试方法 | |
| CN104834807B (zh) | 一种基于分形理论的应力敏感储层相对渗透率计算方法 | |
| CN103674799B (zh) | 一种测定气体在多孔隙介质中轴向扩散系数的装置及方法 | |
| CN102809526B (zh) | 一种测量二氧化碳在饱和油岩心中扩散系数的方法 | |
| CN107038268B (zh) | 一种确定非均质储层五点井网水驱波及系数的方法 | |
| CN107526891B (zh) | 一种聚合物驱大孔道油藏试井分析方法 | |
| CN105893679B (zh) | 低产水平井续流修正试井解释方法 | |
| CN106501144A (zh) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 | |
| CN104297126B (zh) | 低渗透储层气体渗流启动压力梯度测量装置及测量方法 | |
| CN105784567B (zh) | 一种测试岩心相对渗透率的设备和方法 | |
| CN112727424B (zh) | 模拟压裂液注入的井筒-裂缝实验系统以及实验方法 | |
| CN102915406A (zh) | 径向流条件下油水相对渗透率曲线的计算方法 | |
| CN113484216B (zh) | 一种评估致密砂岩气藏水相返排率及合理返排压差的方法 | |
| CN106204304A (zh) | 一种砾岩油藏聚驱相对渗透率曲线的确定方法 | |
| CN109900614A (zh) | 测定超低渗岩心渗透率的方法 | |
| CN105443093A (zh) | 用于注聚井的井口组合测试装置及其方法 | |
| Liu et al. | An inversion method of relative permeability curves in polymer flooding considering physical properties of polymer | |
| Zhang et al. | Tight rock permeability measurement by pressure pulse decay and modeling | |
| CN105804713B (zh) | 一种快速确定注水井各小层井口注水启动压力的方法 | |
| CN106404600B (zh) | 判别粘弹性颗粒驱油剂在多孔介质中渗流行为的方法 | |
| CN111963149B (zh) | 一种考虑滞地液量增压的压裂后地层压力求取方法 | |
| CN104729970B (zh) | 泡沫驱气液相对渗透率曲线的测量方法 | |
| CN113008752B (zh) | 一种油藏型储气库库容计算的有效孔隙体积确定方法 | |
| CN110331966A (zh) | 点坝砂体双向水驱油模拟实验的剩余油饱和度测定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Hou Jian Inventor after: Liu Yongge Inventor after: Liu Lingling Inventor after: Zhou Kang Inventor after: Fu Hongfei Inventor after: Du Qingjun Inventor after: Li Shuxia Inventor after: Yao Chuanjin Inventor after: Xia Zhizeng Inventor before: Hou Jian Inventor before: Fu Hongfei Inventor before: Du Qingjun Inventor before: Li Shuxia Inventor before: Yao Chuanjin Inventor before: Xia Zhizeng Inventor before: Liu Yongge Inventor before: Yang Yong Inventor before: Cao Xulong Inventor before: Liu Lingling Inventor before: Zhou Kang Inventor before: Dai Tao Inventor before: Guo Lanlei Inventor before: Cao Weidong |