CN104949899A - Method for measuring effective viscosity of oil displacing polymer in porous medium - Google Patents

Method for measuring effective viscosity of oil displacing polymer in porous medium Download PDF

Info

Publication number
CN104949899A
CN104949899A CN201510255465.6A CN201510255465A CN104949899A CN 104949899 A CN104949899 A CN 104949899A CN 201510255465 A CN201510255465 A CN 201510255465A CN 104949899 A CN104949899 A CN 104949899A
Authority
CN
China
Prior art keywords
permeability
formula
water
polymer
sand
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.)
Granted
Application number
CN201510255465.6A
Other languages
Chinese (zh)
Other versions
CN104949899B (en
Inventor
黄文章
张鹏
周成裕
陈世兰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing yukeyuan Environmental Protection Technology Co.,Ltd.
Original Assignee
Chongqing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN201510255465.6A priority Critical patent/CN104949899B/en
Publication of CN104949899A publication Critical patent/CN104949899A/en
Application granted granted Critical
Publication of CN104949899B publication Critical patent/CN104949899B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses a method for measuring the effective viscosity of an oil displacing polymer in a porous medium. The method is characterized by comprising the following steps: a, preparing an oil displacing polymer solution; b, through a relation curve of shear stress and shear rate, performing regression to obtain a consistency coefficient H and a flow behavior index n of the polymer solution; c, filling a sand packed model, and measuring the permeability k and porosity phi of the sand packed model; d, flooding the sand packed model with the polymer solution; e, performing subsequent water flooding, namely, measuring the permeability to water k<f> and a residual resistance factor F<RR> after polymer flooding; f, correcting a Blake-kozeny equation to obtain an effective viscosity calculation formula (5) or (6); and g, substituting the consistency coefficient H, the flow behavior index n, the permeability k, the permeability to water k<f> after polymer flooding and the porosity phi or the residual resistance factor F<RR> into the formula to calculate the effective viscosity.

Description

The assay method of a kind of Polymer Used For Oil Displacement virtual viscosity in porous medium
Technical field
The invention belongs to Polymer Used For Oil Displacement viscosimetric analysis technical field in tertiary oil production in oil field, specifically, relate to the measure and calculation method of a kind of Polymer Used For Oil Displacement virtual viscosity in porous medium.
Background technology
Current polymer flooding has become a kind of important means of Old Eastern oilfield stimulation stable yields, the displacing agent overwhelming majority that polymer flooding uses is acrylamide copolymer, such as partial hydrolysis type polyacrylamide, cationic-type polyacrylamide, hydrophobically associated polyacrylamide etc.It is generally acknowledged, the dominant mechanism of polymer flooding reduces aqueous phase mobility, improves water-oil mobility ratio, and then improve sweep efficiency, reaches the object improving recovery ratio.The reduction of aqueous phase mobility has two approach usually: increase aqueous viscosity and reduce water phase permeability (Zhang Peng, Wang Yefei, Zhang Jian etc., several influence factors of hydrophobic association polyacrylamide oil-displacement capacity, oilfield chemistry, 72 volume 4 phase in 2010 462 ~ 468 pages).Therefore the viscosity of polymer solution weighs one of important index of of displacing agent performance.Need it is noted that said aqueous viscosity is not apparent viscosity here, but the virtual viscosity of polymer solution in pore media during seepage flow, this viscosity not only contains the viscosity factor of polymer solution in porous medium, further comprises elastic reaction.Therefore, the virtual viscosity of polymkeric substance reflects viscosity when polymkeric substance flows in pore media more really, design and implementation for mining site scheme more has reference value (Cheng Jiecheng, the flow characteristics of BP16 solution in porous medium, grand celebration oil geology and exploitation, 8 volume 2 phase in 1989 47 ~ 52 pages).
At present, the computing formula of virtual viscosity mainly contains two kinds, and formula (7) is wherein a kind of, μ effor virtual viscosity; Δ P is polymer solution displacement pressure reduction; L is fill out sand tube length; v dfor the percolation flow velocity of polymer solution; K ffor polymer solution rinses permeability by rear water.
&mu; ef = K f &Delta;P v D L - - - ( 7 )
Use formula (7) calculates virtual viscosity and take into account the delay institute of polymkeric substance in porous medium to viscosity role, but does not embody the non-newtonian flow volume property of polymer solution in this formula.
Another computing formula calculating virtual viscosity is Blake-kozeny equation, and see formula (8), in this formula, H is consistency index, n flow stance index, k and φ is respectively permeability and the factor of porosity of rock.Although this computing method take into account the non-newtonian flow volume property of polymkeric substance, do not consider the impact of the delay of polymkeric substance in hole on permeability and factor of porosity.
In the calculating of virtual viscosity, should consider the non-Newtonian fluid mechanic character of polymer solution, also will embody polymer flow impact on permeability and factor of porosity after pore media, the virtual viscosity calculated like this could closing to reality situation more.
Summary of the invention
For solving the problems of the technologies described above, the object of the present invention is to provide the measure and calculation method of a kind of Polymer Used For Oil Displacement virtual viscosity in porous medium, what consider the impact that polymkeric substance non-Newtonian fluid mechanic character and polymkeric substance be detained in pore media sets up virtual viscosity measure and calculation method.
The present invention seeks to the measure and calculation method of a kind of Polymer Used For Oil Displacement virtual viscosity in porous medium realized like this, it is characterized in that being made up of following steps:
The preparation of a, Polymer Used For Oil Displacement solution: the total mineralization injecting water according to Oil Field, uses distilled water and chlorate preparation simulated injection water, re-uses the water-soluble solution Polymer Used For Oil Displacement of simulated injection and obtain Polymer Used For Oil Displacement solution;
B, mensuration consistency index and flow stance index: the rheological curve using rheometer measurement polymer solution, obtains the relation curve of shear stress and shear rate, return out consistency index H and the flow stance index n of polymer solution;
C, filling sand-filling tube model: use silica sand filling fill out sand tube, fill out sand tube diameter is 2.5cm, measures permeability k and the factor of porosity of sand-filling tube model , and draw permeability k and factor of porosity relational expression formula (1), wherein r is pore radius, and τ is tortuosity:
D, polymer solution displacement fill out sand tube: setting constant current pump capacity is 0.23ml/min, and injection of polymer solution is until pressure is steady;
E, sequent water flooding: setting constant current pump capacity is 0.23ml/min, injects simulated injection water, the water measured after poly-driving surveys permeability k fand residual resistance factor F rR, drive rear hole radius r with tortuosity τ does not change assuming that poly-, so can draw and gather the factor of porosity after driving drive rear water survey permeability k with poly- fpass be formula (2); Residual resistance factor F rRwith permeability k and poly-drive rear water and survey permeability k fpass be formula (3):
F RR = k k f - - - ( 3 ) ;
F, calculating virtual viscosity:
Update equation formula: cause permeability and factor of porosity to reduce because polymkeric substance exists delay character in pore media, therefore uses to gather and drives rear water survey permeability replace the permeability of pore media , use the factor of porosity k after poly-driving freplace original porosity k, so, revised Blake-kozeny equation is formula (3):
Formula (2), formula (3) are substituted into formula (4), and the computing formula finally obtaining virtual viscosity is formula (5) or formula (6):
Or
Substitute into and calculate: then by consistency index H, flow stance index n, permeability k, poly-drive rear water and survey permeability k f, factor of porosity or residual resistance factor F rRsubstitution formula (5) or formula (6) result of calculation are virtual viscosity.
In the preparation of above-mentioned steps a Polymer Used For Oil Displacement solution: described chlorate is made up of sodium chloride, lime chloride, magnesium chloride.
Above-mentioned steps c loads in sand-filling tube model: the granularity of described silica sand is 80-100 order, and fill out sand tube length is 30-50cm.
Beneficial effect:
The measure and calculation method of a kind of Polymer Used For Oil Displacement of the present invention virtual viscosity in porous medium, consider the impact of reduction on virtual viscosity that polymkeric substance non-Newtonian fluid mechanic character and polymkeric substance are detained caused permeability and factor of porosity in pore media, solve the problem that prior art accurately cannot obtain Polymer Used For Oil Displacement virtual viscosity and computing method in pore media.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of the inventive method;
Fig. 2 is that in the present invention, massfraction is the shear stress-rate profile of the polymer solution of 0.1%;
Fig. 3 is that in the present invention, massfraction is the shear stress-rate profile of the polymer solution of 0.2%.
Embodiment
Below in conjunction with embodiment and accompanying drawing, the invention will be further described.
Embodiment
A measure and calculation method for Polymer Used For Oil Displacement virtual viscosity in porous medium, as shown in Figure 1:
1, the simulated injection water that distilled water preparation total mineralization is 8625mg/L (calcium ions and magnesium ions 602mg/L) is used.1097.71g sodium chloride, 68.13g lime chloride, 297.28g Magnesium dichloride hexahydrate crystal is dissolved successively in 100mL distilled water.Use simulated injection water-soluble solution displacement of reservoir oil acrylamide copolymer, preparation massfraction is respectively the polymer solution of 0.1% and 0.2%.
2, HAAKE MARS type III flow graph (rotor: PZ38) is used to measure the rheological curve of polymkeric substance.Arranging and measuring temperature is 80 DEG C, and shear rate is 0-800s -1, measure shear stress and shear rate change relation, the shear stress-rate curve of the polymer solution of 0.1% and 0.2% is shown in Fig. 2, Fig. 3 respectively.The consistency index H of polymer solution and flow stance index n is gone out respectively in table 1 by curvilinear regression.
The consistency index H of polymer solution and flow stance index n under table 1 different quality mark
The massfraction of polymer solution Consistency index H Flow stance index n
0.1% 0.2115 0.5150
0.2% 0.3154 0.4794
3, use 80-100 object silica sand to load 6 fill out sand tube, fill out sand tube diameter is 2.5cm, and length is 30cm.People's Republic of China's oil and gas industry standard (SY/T6576-2003) is used to measure permeability and the factor of porosity of sand-filling tube model for the method improved in the recommendation way of the evaluation of polymer of oil recovery.The permeability of six roots of sensation sand-filling tube model and factor of porosity are in table 2.
The permeability k of table 2 sand-filling tube model and factor of porosity
4, service property (quality) mark is polymer solution displacement 6 groups of sand-filling tube models of 0.1% and 0.2% respectively.Experimental temperature is 80 DEG C, and setting constant current pump capacity is 0.23ml/min.After inlet pressure is steady, then use simulated injection water drive for sand-filling tube model by identical discharge capacity, record entry pressure when inlet pressure no longer reduces, calculate to gather by Darcy formula and drive rear water survey permeability k f, through type (3) calculates residual resistance factor F rR.
6, respectively by consistency index H, flow stance index n, permeability k, poly-drive rear water and survey permeability k f, factor of porosity or residual resistance factor F rRsubstitution formula (5) or formula (6) calculate virtual viscosity, and what obtain the results are shown in Table 3.
Or
The various data that table 3 measures in testing

Claims (3)

1. a measure and calculation method for Polymer Used For Oil Displacement virtual viscosity in porous medium, is characterized in that being made up of following steps:
The preparation of a, Polymer Used For Oil Displacement solution: the total mineralization injecting water according to Oil Field, uses distilled water and chlorate preparation simulated injection water, re-uses the water-soluble solution Polymer Used For Oil Displacement of simulated injection and obtain Polymer Used For Oil Displacement solution;
B, mensuration consistency index and flow stance index: the rheological curve using rheometer measurement polymer solution, obtains the relation curve of shear stress and shear rate, return out consistency index H and the flow stance index n of polymer solution;
C, filling sand-filling tube model: use silica sand filling fill out sand tube, fill out sand tube diameter is 2.5cm, measures permeability k and the factor of porosity of sand-filling tube model
D, polymer solution displacement fill out sand tube: setting constant current pump capacity is 0.23ml/min, and injection of polymer solution is until pressure is steady;
E, sequent water flooding: setting constant current pump capacity is 0.23ml/min, injects simulated injection water, the water measured after poly-driving surveys permeability k fand residual resistance factor F rR, through type (2) draws the factor of porosity after poly-driving through type (3) draws residual resistance factor F rR, formula (2) and formula (3) as follows:
F RR = k k f - - - ( 3 ) ;
F, calculating virtual viscosity: by consistency index H, flow stance index n, permeability k, gather and drive rear water survey permeability k f, factor of porosity or residual resistance factor F rRsubstitution formula (5) or formula (6) result of calculation are virtual viscosity, formula (5) or formula (6) as follows:
Or
2. the measure and calculation method of a kind of Polymer Used For Oil Displacement according to claim 1 virtual viscosity in porous medium, is characterized in that in the preparation of step a Polymer Used For Oil Displacement solution: described chlorate is made up of sodium chloride, lime chloride, magnesium chloride.
3. the measure and calculation method of a kind of Polymer Used For Oil Displacement according to claim 1 virtual viscosity in porous medium, is characterized in that step c loads in sand-filling tube model: the granularity of described silica sand is 80-100 order, and fill out sand tube length is 30-50cm.
CN201510255465.6A 2015-05-19 2015-05-19 A kind of assay method of Polymer Used For Oil Displacement effective viscosity in porous media Active CN104949899B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510255465.6A CN104949899B (en) 2015-05-19 2015-05-19 A kind of assay method of Polymer Used For Oil Displacement effective viscosity in porous media

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510255465.6A CN104949899B (en) 2015-05-19 2015-05-19 A kind of assay method of Polymer Used For Oil Displacement effective viscosity in porous media

Publications (2)

Publication Number Publication Date
CN104949899A true CN104949899A (en) 2015-09-30
CN104949899B CN104949899B (en) 2017-09-08

Family

ID=54164724

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510255465.6A Active CN104949899B (en) 2015-05-19 2015-05-19 A kind of assay method of Polymer Used For Oil Displacement effective viscosity in porous media

Country Status (1)

Country Link
CN (1) CN104949899B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105571988A (en) * 2015-12-14 2016-05-11 中国石油天然气股份有限公司 Polymer thermal stability detecting equipment and detecting method
CN107271328A (en) * 2017-06-13 2017-10-20 大港油田集团有限责任公司 A kind of polymer solution anti-shear performance evaluating apparatus and method
CN107869338A (en) * 2017-11-03 2018-04-03 中国石油化工股份有限公司 Chemical flooding common heavy oil dynamic viscosity reduction evaluation method
CN110059438A (en) * 2019-04-28 2019-07-26 东北石油大学 The poly- calculation method for driving the deformation of microchannel residual oil film
CN110457857A (en) * 2019-08-21 2019-11-15 中国海洋石油集团有限公司 A method of polymer flooding effect under the conditions of prediction different crude oils viscosity
CN114324066A (en) * 2021-12-30 2022-04-12 中海油田服务股份有限公司 Evaluation method of modified nano material for rock core
CN115306370A (en) * 2022-08-10 2022-11-08 重庆科技学院 Experimental device and method for on-line monitoring of residual resistance coefficient of polymer

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10197436A (en) * 1996-12-28 1998-07-31 Soubu Namakon Kk Method for testing fluidity characteristic of highly fluid concrete and apparatus therefor
CN101058723A (en) * 2007-04-26 2007-10-24 中国石油大学(华东) Re-use agent for polymer after driving and injecting method thereof
CN101545368A (en) * 2009-04-27 2009-09-30 中国石油大学(华东) On-line depth profile control method for improving oil displacement effect of polymer of offshore production well
CN102278090A (en) * 2011-08-16 2011-12-14 中国石油化工股份有限公司 Method and control system for stopping leakage of reservoir fracture in drilling
CN102455274A (en) * 2010-10-18 2012-05-16 袁俊海 Method for measuring intrinsic viscosity of polymer
CN102980835A (en) * 2012-12-13 2013-03-20 中国石油化工股份有限公司 Method for measuring viscosity-reducing effect of oil-soluble viscosity reducer suitable for super heavy oil and application of method
CN103196796A (en) * 2013-04-15 2013-07-10 中国石油大学(华东) Experimental device and method for researching viscosity property of supercritical carbon dioxide in stratum
CN103267710A (en) * 2013-05-30 2013-08-28 中国石油大学(北京) Measurement and calculation method for effective viscosity of VES (viscoelastic surfactant) variable-viscosity acid in porous medium
CN103344525A (en) * 2013-08-02 2013-10-09 成都理工大学 Method and device for testing effective viscosity of foams in pore medium
CN103760065A (en) * 2013-06-28 2014-04-30 中国石油化工股份有限公司 Test method and test system for polymer solution flow effective apparent viscosity
CN103775078A (en) * 2014-02-14 2014-05-07 武汉科技大学 Method for forecasting effective viscosity of power law fluid in porous medium

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10197436A (en) * 1996-12-28 1998-07-31 Soubu Namakon Kk Method for testing fluidity characteristic of highly fluid concrete and apparatus therefor
CN101058723A (en) * 2007-04-26 2007-10-24 中国石油大学(华东) Re-use agent for polymer after driving and injecting method thereof
CN101545368A (en) * 2009-04-27 2009-09-30 中国石油大学(华东) On-line depth profile control method for improving oil displacement effect of polymer of offshore production well
CN102455274A (en) * 2010-10-18 2012-05-16 袁俊海 Method for measuring intrinsic viscosity of polymer
CN102278090A (en) * 2011-08-16 2011-12-14 中国石油化工股份有限公司 Method and control system for stopping leakage of reservoir fracture in drilling
CN102980835A (en) * 2012-12-13 2013-03-20 中国石油化工股份有限公司 Method for measuring viscosity-reducing effect of oil-soluble viscosity reducer suitable for super heavy oil and application of method
CN103196796A (en) * 2013-04-15 2013-07-10 中国石油大学(华东) Experimental device and method for researching viscosity property of supercritical carbon dioxide in stratum
CN103267710A (en) * 2013-05-30 2013-08-28 中国石油大学(北京) Measurement and calculation method for effective viscosity of VES (viscoelastic surfactant) variable-viscosity acid in porous medium
CN103760065A (en) * 2013-06-28 2014-04-30 中国石油化工股份有限公司 Test method and test system for polymer solution flow effective apparent viscosity
CN103344525A (en) * 2013-08-02 2013-10-09 成都理工大学 Method and device for testing effective viscosity of foams in pore medium
CN103775078A (en) * 2014-02-14 2014-05-07 武汉科技大学 Method for forecasting effective viscosity of power law fluid in porous medium

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
佟曼丽 等: "聚合物稀溶液流经孔隙介质时的粘弹效应及其表征", 《油田化学》 *
张鹏 等: "疏水缔合聚丙烯酰胺驱油能力的几种影响因素", 《油田化学》 *
杨清彦 等: "聚合物驱相对渗透率曲线计算方法研究", 《石油化学》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105571988A (en) * 2015-12-14 2016-05-11 中国石油天然气股份有限公司 Polymer thermal stability detecting equipment and detecting method
CN107271328A (en) * 2017-06-13 2017-10-20 大港油田集团有限责任公司 A kind of polymer solution anti-shear performance evaluating apparatus and method
CN107271328B (en) * 2017-06-13 2019-11-22 大港油田集团有限责任公司 A kind of polymer solution anti-shear performance evaluating apparatus and method
CN107869338A (en) * 2017-11-03 2018-04-03 中国石油化工股份有限公司 Chemical flooding common heavy oil dynamic viscosity reduction evaluation method
CN110059438A (en) * 2019-04-28 2019-07-26 东北石油大学 The poly- calculation method for driving the deformation of microchannel residual oil film
CN110059438B (en) * 2019-04-28 2022-10-11 东北石油大学 Calculation method for residual oil film deformation of polymer-driven microchannel
CN110457857A (en) * 2019-08-21 2019-11-15 中国海洋石油集团有限公司 A method of polymer flooding effect under the conditions of prediction different crude oils viscosity
CN110457857B (en) * 2019-08-21 2022-11-15 中国海洋石油集团有限公司 Method for predicting polymer flooding effect under different crude oil viscosity conditions
CN114324066A (en) * 2021-12-30 2022-04-12 中海油田服务股份有限公司 Evaluation method of modified nano material for rock core
CN114324066B (en) * 2021-12-30 2023-09-12 中海油田服务股份有限公司 Evaluation method of modified nanomaterial for core
CN115306370A (en) * 2022-08-10 2022-11-08 重庆科技学院 Experimental device and method for on-line monitoring of residual resistance coefficient of polymer
CN115306370B (en) * 2022-08-10 2023-11-24 重庆科技学院 Experimental device and method for online monitoring of residual resistance coefficient of polymer

Also Published As

Publication number Publication date
CN104949899B (en) 2017-09-08

Similar Documents

Publication Publication Date Title
CN104949899A (en) Method for measuring effective viscosity of oil displacing polymer in porous medium
CN102720486B (en) Device for testing shale gas fracture network diversion capacity and working method of device
Patton et al. Prediction of polymer flood performance
CN105715241B (en) A kind of measuring method of Curve of Polymer Displacement Relative Permeability
WO2012083068A3 (en) Method of determining reservoir pressure
CN105670584B (en) A kind of temporary stifled type kill-job fluid being suitable for hypotonic low-pressure gas well well workover and preparation method
CN204402469U (en) Hydraulic fracturing geostress survey die system
CN109444339B (en) A kind of shale gas horizontal well refracturing fiber temporarily blocks up experimental test procedures
RU2006138042A (en) METHODS AND DEVICES FOR EVALUATING THE PHYSICAL PARAMETERS OF TANKS USING THE CURVE PRESSURE RESTORATION CURVES METHOD WHEN TESTING BY A PRESSURE / RESET TEST
CN109033504B (en) Oil-water well casing damage prediction method
US20190368349A1 (en) Method for determining gelation time in a core plug
CN106520096A (en) A well killing fluid, a preparing method thereof and applications of the well killing fluid
US4848461A (en) Method of evaluating fracturing fluid performance in subsurface fracturing operations
CN103867178A (en) Method for effectively evaluating damage of fracturing fluid to support crack flow diversion capability
CN106204304A (en) A kind of Conglomerate Reservoir gathers the determination method driving permeability saturation curve
CN102425386B (en) Drilling fluid rheological parameter control method meeting power law mode
CN106404634A (en) Method of evaluating polymer solution lowering the permeability of oil field storage layer
CN105257278A (en) Method for acquiring insertion depth of propping agent
CN103334725A (en) Method and device for evaluating low permeability reservoir displacement effectiveness
CN107679338B (en) Reservoir fracturing effect evaluation method and evaluation system based on flowback data
CN108952668B (en) Fracturing method of normal-pressure shale gas reservoir
CN111963149B (en) Post-fracturing stratum pressure solving method taking earth stagnation amount pressurization into consideration
CN103217360A (en) Method for detecting molecular weight of hydrolyzed polyacrylamide
CN106295235A (en) A kind of computational methods of shale formation intrinsic fracture shearing slip amount
CN109918769A (en) Utilize the method for instantaneous equation calculation fracture-pore reservoir unstable state water enchroachment (invasion) water influx

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220321

Address after: 400050 502-8, 5th floor, Julong steel market complex building, hub Park, international logistics park, Shapingba District, Chongqing

Patentee after: Chongqing yukeyuan Environmental Protection Technology Co.,Ltd.

Address before: No. 20, East Road, University City, Chongqing, Shapingba District, Chongqing

Patentee before: Chongqing University of Science & Technology