CN109033518A - The water breakthrough time prediction technique and device of bottom water gas condensate reservoir - Google Patents

The water breakthrough time prediction technique and device of bottom water gas condensate reservoir Download PDF

Info

Publication number
CN109033518A
CN109033518A CN201810649707.3A CN201810649707A CN109033518A CN 109033518 A CN109033518 A CN 109033518A CN 201810649707 A CN201810649707 A CN 201810649707A CN 109033518 A CN109033518 A CN 109033518A
Authority
CN
China
Prior art keywords
gas
water
breakthrough time
condensate reservoir
phase
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.)
Pending
Application number
CN201810649707.3A
Other languages
Chinese (zh)
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.)
China Petroleum and Natural Gas Co Ltd
CNPC Jianghan Machinery Research Institute Co Ltd
CNPC Engineering Technology Research Institute Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
CNPC Jianghan Machinery Research Institute Co Ltd
CNPC Engineering Technology Research Institute Co Ltd
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 China Petroleum and Natural Gas Co Ltd, CNPC Jianghan Machinery Research Institute Co Ltd, CNPC Engineering Technology Research Institute Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201810649707.3A priority Critical patent/CN109033518A/en
Publication of CN109033518A publication Critical patent/CN109033518A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Theoretical Computer Science (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The present invention provides the water breakthrough time prediction technique and device of a kind of bottom water gas condensate reservoir, this method comprises: according to skin effect influence factor, gas phase non-Darcy effects influence factor and retrograde condensation function influence factor, the centripetal stream productivity model of gas phase hemispherical of gas phase radial fluid flow productivity model and perforated zone lower part based on perforated zone, establishes the water breakthrough time computation model of bottom water gas condensate reservoir;Obtain the actual production data of bottom water gas condensate reservoir;The water breakthrough time of bottom water gas condensate reservoir is predicted according to the actual production data of bottom water gas condensate reservoir based on the water breakthrough time computation model of the bottom water gas condensate reservoir.Since the program has comprehensively considered the influence of gas phase non-Darcy effects near shaft bottom, skin effect, retrograde condensation effect to the water breakthrough time, foundation more suits the dynamic water breakthrough time computation model of bottom water gas condensate reservoir water cone, and the accuracy of water breakthrough time prediction can be improved.

Description

The water breakthrough time prediction technique and device of bottom water gas condensate reservoir
Technical field
The present invention relates to the water breakthrough risk profile technical field of the natural gas well, in particular to a kind of bottom water gas condensate reservoir is shown in Water time forecasting methods and device.
Background technique
During condensate gas reservoir, when formation pore pressure is reduced to dew-point pressure, liquid will form in the earth formation Condensate flows in stratum for oil gas water three phase, causes flow event increasingly complex.Since there are liquid for bottom water gas condensate reservoir Condensate, therefore its bottom water coning dynamic is different from conventional bottomwater gas field.If calculating mould with the water breakthrough time of conventional Gas Reservoirs out Type is predicted, certainly will will lead to prediction result and error or even mistake occurs.Currently, in the industry to conventional bottom-water reservoir water breakthrough The research of time is more, such as: 1) the water breakthrough time prediction model of bottom water reservoir is established by experiment, with a variety of theory sides Method has studied influence of each factor to oil reservoir bottom water coning.2) the complicated percolation mode based on the dual imperfect well of bottomwater gas field, It derives bottomwater gas field water breakthrough time prediction model, and analyzes influence of the gas reservoir opening degree to the water breakthrough time.3) it considers Hyposmosis bottomwater gas field has starting pressure gradient, derives hypotonic bottomwater gas field water breakthrough time prediction model.4) with artificial mind The water breakthrough time prediction technique of fractured well is established through network technique (ANN) and support vector machines (SVN).5) conformal mapping is used Method derives fractured well new productivity prediction equation, obtains Fractured Gas Wells water breakthrough time model in conjunction with water cone vertex movements equation.Then The influence for considering gas volume factor, each area's seepage state and fracture condudtiviy these three factors, to above-mentioned Fractured Gas Wells Water breakthrough time model is corrected, and improves the precision of hyposmosis bottomwater gas field fractured well water breakthrough time prediction model.6) it pushes away Export band prevention of coning problems in bottom-water reservoirs gas reservoir water breakthrough time prediction model.7) statistical method is used, derives slit formation bottomwater gas field water breakthrough The calculation formula of time, and analyze the influence of gas production and Cracking Thickness to the water breakthrough time.8) fracture hole type bottom water gas is considered The difference of hiding and conventional bottomwater gas field, is managed with homogeneous bottomwater gas field horizontal well water breakthrough time prediction technique and permeability variation By deriving fracture hole type bottomwater gas field water breakthrough time predictor formula.9) it is established with cluster discriminant analysis and Method for Numerical The water breakthrough time prediction model of deep water gas reservoir.And initial stage starting is still in for the research of bottom water gas condensate reservoir water breakthrough time Stage, the influence that some scholars are precipitated by considering condensate are deduced bottom water gas condensate reservoir water breakthrough time prediction model, But its derivation process is the flow model in porous media based on the centripetal stream of spherical surface, does not consider the horizontal radial flowing of perforated zone, does not meet gas The true seepage flow situation of well;And the shaft bottom influence of gas phase non-Darcy effects and skin effect to the water breakthrough time nearby is had ignored, These will lead to the inaccuracy of water breakthrough time prediction.
Summary of the invention
The embodiment of the invention provides a kind of water breakthrough time prediction technique of bottom water gas condensate reservoir and devices, are based on gas water The bottom water coning model that flat radial flow and the centripetal stream of hemispherical combine, and comprehensively considered the shaft bottom non-darcy effect of gas phase nearby It answers, influence of skin effect, the retrograde condensation effect to the water breakthrough time, foundation more suits bottom water gas condensate reservoir water cone and dynamically sees The accuracy of water breakthrough time prediction can be improved in water time computation model.
The water breakthrough time prediction technique of the bottom water gas condensate reservoir includes:
According to skin effect influence factor, gas phase non-Darcy effects influence factor and retrograde condensation function influence factor, it is based on The gas phase radial fluid flow productivity model of perforated zone and the centripetal stream productivity model of the gas phase hemispherical of perforated zone lower part, build The water breakthrough time computation model of vertical bottom water gas condensate reservoir;
Obtain the actual production data of bottom water gas condensate reservoir;
Based on the water breakthrough time computation model of the bottom water gas condensate reservoir, according to the actual production number of bottom water gas condensate reservoir According to the water breakthrough time of prediction bottom water gas condensate reservoir.
The water breakthrough time prediction meanss of the bottom water gas condensate reservoir include:
The water breakthrough time computation model of bottom water gas condensate reservoir establishes module, for according to skin effect influence factor, gas phase Non-Darcy effects influence factor and retrograde condensation function influence factor, the gas phase radial fluid flow productivity model based on perforated zone With the centripetal stream productivity model of gas phase hemispherical of perforated zone lower part, the water breakthrough time computation model of bottom water gas condensate reservoir is established;
Actual production data acquisition module, for obtaining the actual production data of bottom water gas condensate reservoir;
The water breakthrough time prediction module of bottom water gas condensate reservoir is calculated for the water breakthrough time based on the bottom water gas condensate reservoir Model predicts the water breakthrough time of bottom water gas condensate reservoir according to the actual production data of bottom water gas condensate reservoir.
In embodiments of the present invention, the gas phase radial fluid flow productivity model based on perforated zone and perforated zone lower part The centripetal stream productivity model of gas phase hemispherical, and comprehensively considered shaft bottom nearby gas phase non-Darcy effects, skin effect, retrograde condensation The influence to the water breakthrough time is acted on, establishes and more suits the dynamic water breakthrough time computation model of bottom water gas condensate reservoir water cone, it can To improve the accuracy of water breakthrough time prediction.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention without creative efforts, may be used also for those of ordinary skill in the art To obtain other drawings based on these drawings.
Fig. 1 is a kind of water breakthrough time prediction technique flow chart of bottom water gas condensate reservoir provided in an embodiment of the present invention;
Fig. 2 is a kind of bottom water gas condensate reservoir water cone physical model schematic diagram provided in an embodiment of the present invention;
Fig. 3 is a kind of influence signal of the gas condensate reservoir gas pay thickness provided in an embodiment of the present invention to the bottom-water breakthrough time Figure;
Fig. 4 is a kind of influence schematic diagram of the gas condensate reservoir gas production provided in an embodiment of the present invention to the bottom-water breakthrough time;
Fig. 5 is a kind of water breakthrough time prediction meanss structural block diagram of bottom water gas condensate reservoir provided in an embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that the described embodiment is only a part of the embodiment of the present invention, instead of all the embodiments.Based on this Embodiment in invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
In embodiments of the present invention, a kind of water breakthrough time prediction technique of bottom water gas condensate reservoir is provided, as shown in Figure 1, This method comprises:
Step 101: according to skin effect influence factor, gas phase non-Darcy effects influence factor and retrograde condensation function influence Factor, the centripetal miscarriage of gas phase hemispherical of gas phase radial fluid flow productivity model and perforated zone lower part based on perforated zone Energy model, establishes the water breakthrough time computation model of bottom water gas condensate reservoir;
Step 102: obtaining the actual production data of bottom water gas condensate reservoir;
Step 103: the water breakthrough time computation model based on the bottom water gas condensate reservoir, according to the reality of bottom water gas condensate reservoir Border creation data predicts the water breakthrough time of bottom water gas condensate reservoir.
When it is implemented, the water breakthrough time that step 101 specifically establishes bottom water gas condensate reservoir as follows calculates mould Type.
Initially set up bottom water gas condensate reservoir water cone physical model as shown in Figure 2: gas well brill opens part gas condensate reservoir, penetrates Hole interval is the radial fluid flow of gas phase, and perforated zone lower part is the centripetal stream of hemispherical of gas phase;After gas well is gone into operation, water flooding To shaft bottom coning, most short along the gas well axis direction intrusion shaft bottom time, this time is the water breakthrough time of bottom water gas condensate reservoir. In the case where meeting bottom water gas condensate reservoir Liquid Flow rule, make the following assumptions: 1. reservoir uniform thickness and homogeneous;2. air water two In phase flow event, with movement displacement;3. ignoring the influence of capillary force and gravity, viscous force controls air water two-phase on ground Flowing in layer;4. fluid is micro- compressible, gas phase, the density of liquid phase and viscosity are constant;5. gas phase flows full in the earth formation The non-Darcy's law of foot;6. condensate does not flow in the earth formation;7. original gas-water interface approximation can regard a horizontal plane as.
Then following steps are executed:
(1) meet non-Darcy's law according to gas phase and water phase meets Darcy's law, establish air water two phase fluid flow model;
According to two phase fluid flow rule, gas phase meets non-Darcy's law, and water phase meets Darcy's law, then air water two phase fluid flow Model is respectively as follows:
In formula: Pg、PwThe respectively pressure of gas phase and water phase, MPa;μg、μwThe respectively viscosity of gas phase and water phase, mPa·s;Vg、VwThe respectively percolation flow velocity of gas phase and water phase, m/d;Kg、KwThe respectively effective permeability of gas phase and water phase, mD;R is radial radius, m;β is non-Darcy flow coefficient, m-1;ρgFor the density of gas phase, g/cm3
(2) according to air water two phase fluid flow model, liquid and gas percolation flow velocity relational expression is established;
Specifically, it is to be seeped to calculate the time that water cone vertex A reaches shaft bottom according to gas-liquid two-phase that mathematical model, which is established, Stream feature is it is found that the barometric gradient of gas phase and the barometric gradient of water phase are identical at this point, therefore have:
Formula (1), formula (2) are substituted into formula (3) respectively, the relational expression of liquid phase Yu gas phase percolation flow velocity can be obtained:
In formula: MwgFor aqueous vapor mobility ratio, Mwg=(Kww)/(Kgg)。
(3) it establishes the moving distance on the water cone vertex on original gas-water interface and reaches the time history form in shaft bottom;
Specifically, consider the influence of porosity, original water saturation, residual gas saturation and retrograde gas condensate saturation degree, The then relationship of the moving distance of water cone vertex A and time are as follows:
In formula: Φ is porosity, dimensionless;T is time, d;SwiFor original water saturation, dimensionless;SgrFor remnants Gas saturation, dimensionless;SocFor remaining condensate saturation degree, dimensionless.
Formula (5) is deformed, can be obtained:
(4) it is miscarried according to the gas phase planar radial for the perforated zone for considering gas non-Darcy effects and stratum skin effect The centripetal miscarriage energy of gas phase hemispherical of the perforated zone lower part of energy model, consideration gas non-Darcy effects and stratum skin effect Model establishes the gas production of bottom water gas condensate reservoir perforated zone lower part;
Specifically, the gas phase of bottom water gas condensate reservoir perforated zone is to be flowed in the form of radial fluid flow to shaft bottom, it is some Scholar derives the radial fluid flow Productivity Formulae for considering gas non-Darcy effects and stratum skin effect:
In formula (7)~(8): PeFor supply pressure, MPa;PwfFor bottom pressure, MPa;T is reservoir temperature, K;Z is deviation The factor, dimensionless;hpFor bottom water gas condensate reservoir perforated zone thickness, m;Q1For the yield of perforated zone, m3/d;reFor supply half Diameter, m;rwFor wellbore radius, m;S is skin factor, dimensionless;γgFor gas phase relative density, dimensionless;D1For inertia coeffeicent, (m3·d-1)-1
The gas phase of bottom water gas condensate reservoir perforated zone lower part is to be flowed in the form of the centripetal stream of hemispherical to shaft bottom, Yi Xiexue Person derives the centripetal stream Productivity Formulae of hemispherical for considering gas non-Darcy effects and stratum skin effect:
In formula (9)~(10): D2For inertia coeffeicent, (m3·d-1)-1;KsFor effective permeability,mD;Kv For permeability in vertical direction, mD;KhFor permeability in horizontal direction, mD;Q2For bottom water gas condensate reservoir perforated zone lower part Gas production, m3/d。
Formula (7) and formula (9) are subjected to simultaneous, can be obtained:
Assuming that the total gas production of condensate gas well is Q, then have
Q=Q1+Q2 (12)
Formula (11) and formula (12) are subjected to simultaneous, eliminate Q1, obtain related Q2Equation:
Formula (13) is carried out to simplify arrangement:
The then gas production of bottom water gas condensate reservoir perforated zone lower part are as follows:
(5) according to the moving distance on the water cone vertex on liquid and gas percolation flow velocity relational expression, original gas-water interface with Time relationship, the gas production of bottom water gas condensate reservoir perforated zone lower part and the retrograde condensation function influence factor in shaft bottom are reached, is established The water breakthrough time computation model of bottom water gas condensate reservoir.
Specifically, obtaining the expression formula in relation to condensate saturation degree change rate first:
In formula: BgFor the volume factor of gas, dimensionless;H is gas reservoir thickness, m;DC is to coagulate in unit volume condensate gas Condensate oil increment (or reduction amount), m3/m3;rbFor the obstruction radius of condensate saturation degree critical in reservoir gaps, m.
Formula (19) are substituted into formula (6) and are arranged:
Obstruction radius may be expressed as:
In formula: Y is the retrograde condensation factor, m3/(m3·MPa);PRFor original formation pressure, MPa.
If the vertical range in original gas-water interface and shaft bottom is h before gas well liquid loadingb, work as t=tbtWhen, then gas well starts to see Water integrates formula (20), can obtain:
In formula, hbFor the thickness of gas reservoir perforated zone lower part, m.
For bottom water gas condensate reservoir perforated zone lower part, gas makees hemispherical centripetal flow, then the movement speed of gas can It indicates are as follows:
Formula (4), formula (23) are substituted into formula (22), first formula (22) right side of the equal sign is arranged:
It enablesThen formula (24) can be changed by integral are as follows:
Formula (4), formula (21) and formula (23) are substituted into formula (22) again, formula (22) left side of the equal sign is arranged:
Formula (25), formula (26) are substituted into formula (22), can be obtained:
Formula (27) is carried out to simplify arrangement:
Then consider that the bottom water gas condensate reservoir water breakthrough time of skin effect, gas phase non-Darcy effects and retrograde condensation effect calculates mould Type are as follows:
Embodiment
The basic parameter (actual production data) of certain bottom water gas condensate reservoir is as follows: original formation pressure 44.37MPa, gas Hiding temperature is 95.39 DEG C, and the density of gaseous viscosity 0.034mPas, aqueous viscosity 0.23mPas, gas phase are 0.063g/cm3, gas pay thickness 190.31m, perforating depth 88.71m, formation porosity 0.186, gas-bearing formation effective permeability For 175mD, the volume factor of condensate gas is 0.0043, and wellbore radius 0.1m, pressure release radius is 250m, residual gas saturation It is 0.59, remaining condensate saturation degree is 0.44, formation porosity 0.16, skin factor 3, and initial water saturation is 0.56, gas well be averaged gas production be 10.88 × 104m3/d。
The prediction new model derived by the present invention and several common models carry out calculating and precision point to field basic data Analysis can be obtained different bottom water coning models and calculate the resulting water breakthrough time, as shown in table 1.
The comparison of each the model calculation of table 1
Model Calculate time tbt/d Real time/d Relative error/%
Zhang Liehui model 1650 1269 30.02%
King can strong model 2197 1269 73.13%
Yellow complete magnificent model 2004 1269 57.92%
This paper model 1405 1269 10.72%
In table 1, king can the derivation of strong model and Huang Quan China model be based on conventional bottomwater gas field, and Zhang Liehui model is then It is based on bottom water gas condensate reservoir;Wherein, king can strong model have ignored skin effect, gas phase non-Darcy effects and retrograde condensation effect The influence of these three factors, yellow complete magnificent model considers skin effect, gas phase non-Darcy effects, but has ignored retrograde condensation effect Influence, and Zhang Liehui model consider bottom water gas condensate reservoir retrograde condensation effect influence, but the model be based on spherical surface it is centripetal The flow model in porous media of stream, and have ignored the influence of skin effect and gas phase non-Darcy effects.
As seen from Table 1, king can strong model precision of prediction it is minimum, relative error is up to 71.13%, main cause It is to have ignored the influence of skin effect, gas phase non-Darcy effects and retrograde condensation effect, and these three factors can make gas flow rate Increase, the water breakthrough time shortens, therefore king can the bottom water gas condensate reservoir water breakthrough time longest that calculates of strong model.Yellow complete magnificent model is compared For the former, precision of prediction is slightly improved, but relative error is still very big, and nearly 60%, it is primarily due to that takes into account skin effects With gas phase non-Darcy effects, but retrograde condensation effect is had ignored, and retrograde condensation effect can be such that gas flow rate increases, therefore yellow full Hua Mo The water breakthrough time that type calculates is shorter compared to the former.And Zhang Liehui model considers the influence of retrograde condensation effect, but has ignored epidermis Effect and gas phase non-Darcy effects, and the model is the flow model in porous media based on the centripetal stream of spherical surface, when causing to calculate gas well water breakthrough Between the gas production that is taken be gas well total gas production, and confirmed through excessive scholar, it is centripetal based on gas level radial flow and hemispherical The bottom water coning model that stream combines is more reasonable, and gas production should take the centripetal miscarriage amount of hemispherical, therefore in general, Zhang Liehui The calculated result of model is bigger than normal.And with the new model derived of the present invention and meanwhile consider skin effect, gas phase non-Darcy effects and The influence of retrograde condensation effect, the result relative error of calculating is smaller, and about 10%, it is contemplated that the complexity of practical reservoir, and The regular adjustment of yield in production process, this error are acceptables.Compared with above-mentioned common calculation model, precision Improve 19.30%~62.41% and the field working conditions goodness of fit it is preferable, illustrate that computation model of the present invention can Accurate Prediction bottom water The water breakthrough time of gas condensate reservoir effectively instructs live continuous production.
Fig. 3 be different gas pay thicknesses under the conditions of, while ignore skin effect, gas phase non-Darcy effects and retrograde condensation effect, Consider skin effect and gas phase non-Darcy effects, and considers the bottom water gas condensate reservoir water breakthrough time of retrograde condensation effect.It can from Fig. 3 It was found that: while when ignoring skin effect, gas phase non-Darcy effects and retrograde condensation effect, water breakthrough time longest, water breakthrough is the latest;When When gas pay thickness is identical, skin effect, gas phase non-Darcy effects and retrograde condensation, which act on these three factors, will lead to gas condensate reservoir Bottom water accelerates coning, and the water breakthrough time shortens, and the influence of retrograde condensation effect is the most significant;When gas pay thickness is smaller, when water breakthrough Between influenced by skin effect and gas phase non-Darcy effects little, but with the continuous increase of gas pay thickness, the water breakthrough time is by table The influence of skin effect and gas phase non-Darcy effects is gradually increased;Likewise, the water breakthrough time is by condensation when gas pay thickness is smaller The influence of gas reservoir retrograde condensation effect is smaller, but being gradually increased with gas pay thickness, water breakthrough time are influenced gradually by the factor Increase;Regardless of whether consider the influence of these three factors, with the continuous increase of gas pay thickness, the gas well water breakthrough time all can be by It is cumulative to add, and amplitude is gradually increased with the increase of gas pay thickness, and gas well water breakthrough is more late.
Fig. 4 is influence of the different gas production to the bottom water gas condensate reservoir water breakthrough time.As can be seen from Figure 4: with gas production Be gradually increased, the bottom water gas condensate reservoir water breakthrough time constantly shortens;Gas production is bigger, then the producing pressure differential of gas well also increases therewith Greatly, bottom water accelerates coning, and water breakthrough is more early;Meanwhile with the continuous increase of gas production, the water breakthrough time first accelerates to successively decrease, rear slow Successively decrease, in general, lower yield is more advantageous to the coning for inhibiting bottom water, extends the anhydrous gas production phase.
Instance analysis shows that condensate reservoir development process is increasingly complex compared with normal gas pools;Utilize model prediction of the present invention Water breakthrough time is closer practical;Skin effect, gas phase non-Darcy effects and retrograde condensation effect will lead to water breakthrough time shortening; As gas pay thickness increases, gas well water breakthrough is more late;And as gas production increases, the water breakthrough time shows successively decreasing first quick and back slow Trend.
Based on the same inventive concept, a kind of water breakthrough time prediction of bottom water gas condensate reservoir is additionally provided in the embodiment of the present invention Device, as described in the following examples.The principle solved the problems, such as due to the water breakthrough time prediction meanss of bottom water gas condensate reservoir and bottom The water breakthrough time prediction technique of water gas condensate reservoir is similar, therefore the implementation of the water breakthrough time prediction meanss of bottom water gas condensate reservoir can With referring to the implementation of the water breakthrough time prediction technique of bottom water gas condensate reservoir, overlaps will not be repeated.It is used below, art The combination of the software and/or hardware of predetermined function may be implemented in language " unit " or " module ".Although described by following embodiment Device preferably realize that but the combined realization of hardware or software and hardware is also that may and be contemplated with software 's.
Fig. 5 is a kind of structural block diagram of the water breakthrough time prediction meanss of the bottom water gas condensate reservoir of the embodiment of the present invention, is such as schemed Shown in 5, comprising:
The water breakthrough time computation model of bottom water gas condensate reservoir establishes module 501, for according to skin effect influence factor, Gas phase non-Darcy effects influence factor and retrograde condensation function influence factor, the gas phase radial fluid flow production capacity based on perforated zone The centripetal stream productivity model of the gas phase hemispherical of model and perforated zone lower part, the water breakthrough time for establishing bottom water gas condensate reservoir calculate Model;
Actual production data acquisition module 502, for obtaining the actual production data of bottom water gas condensate reservoir;
The water breakthrough time prediction module 503 of bottom water gas condensate reservoir, for the water breakthrough time based on the bottom water gas condensate reservoir Computation model predicts the water breakthrough time of bottom water gas condensate reservoir according to the actual production data of bottom water gas condensate reservoir.
The structure is illustrated below.
It is specifically used for when it is implemented, the water breakthrough time computation model of the bottom water gas condensate reservoir establishes module 501:
Meet non-Darcy's law according to gas phase and water phase meets Darcy's law, establishes air water two phase fluid flow model;
According to air water two phase fluid flow model, liquid and gas percolation flow velocity relational expression is established;
It establishes the moving distance on the water cone vertex on original gas-water interface and reaches the time history form in shaft bottom;
According to the gas phase radial fluid flow production capacity mould for the perforated zone for considering gas non-Darcy effects and stratum skin effect The centripetal stream productivity model of gas phase hemispherical of the perforated zone lower part of type, consideration gas non-Darcy effects and stratum skin effect, Establish the gas production of bottom water gas condensate reservoir perforated zone lower part;
According to the moving distance and arrival on the water cone vertex on liquid and gas percolation flow velocity relational expression, original gas-water interface The time relationship in shaft bottom, the gas production of bottom water gas condensate reservoir perforated zone lower part and retrograde condensation function influence factor, establish bottom water The water breakthrough time computation model of gas condensate reservoir.
It is specifically used for when it is implemented, the water breakthrough time computation model of the bottom water gas condensate reservoir establishes module 501: presses The air water two phase fluid flow model is established according to formula (1)-(2).
It is specifically used for when it is implemented, the water breakthrough time computation model of the bottom water gas condensate reservoir establishes module 501: presses The liquid and gas percolation flow velocity relational expression is established according to formula (3)-(4).
It is specifically used for when it is implemented, the water breakthrough time computation model of the bottom water gas condensate reservoir establishes module 501: presses The moving distance on the water cone vertex on the original gas-water interface is established according to formula (5)-(6) and reaches the time relationship in shaft bottom Formula.
It is specifically used for when it is implemented, the water breakthrough time computation model of the bottom water gas condensate reservoir establishes module 501: presses The gas phase plane diameter of the perforated zone of the consideration gas non-Darcy effects and stratum skin effect is established according to formula (7)-(8) To stream productivity model;
It is established according to formula (9)-(10) under the perforated zone of the consideration gas non-Darcy effects and stratum skin effect The centripetal stream productivity model of the gas phase hemispherical in portion;
The gas production of bottom water gas condensate reservoir perforated zone lower part is established according to formula (11)-(18).
It is specifically used for when it is implemented, the water breakthrough time computation model of the bottom water gas condensate reservoir establishes module 501: presses The water breakthrough time computation model of the bottom water gas condensate reservoir is established according to formula (19)-(31).
In conclusion the present invention is based on the gas of the gas phase radial fluid flow productivity model of perforated zone and perforated zone lower part The centripetal stream productivity model of phase hemispherical, and comprehensively considered shaft bottom nearby gas phase non-Darcy effects, skin effect, retrograde condensation work With the influence to the water breakthrough time, establishes and more suit the dynamic water breakthrough time computation model of bottom water gas condensate reservoir water cone, it can be with Improve the accuracy of water breakthrough time prediction.
It should be understood by those skilled in the art that, the embodiment of the present invention can provide as method, system or computer journey Sequence product.Therefore, complete hardware embodiment, complete software embodiment or combining software and hardware aspects can be used in the present invention The form of embodiment.Moreover, it wherein includes the calculating of computer usable program code that the present invention, which can be used in one or more, The computer program implemented in machine usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) The form of product.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions each in flowchart and/or the block diagram The combination of process and/or box in process and/or box and flowchart and/or the block diagram.It can provide these computers Processor of the program instruction to general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices To generate a machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute For realizing the function of being specified in one or more flows of the flowchart and/or one or more blocks of the block diagram Device.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works, so that instruction stored in the computer readable memory generation includes The manufacture of command device, the command device are realized in one box of one or more flows of the flowchart and/or block diagram Or the function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer Or the instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or box The step of function of being specified in figure one box or multiple boxes.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field For art personnel, the embodiment of the present invention can have various modifications and variations.All within the spirits and principles of the present invention, made Any modification, equivalent substitution, improvement and etc., should all be included in the protection scope of the present invention.

Claims (14)

1. a kind of water breakthrough time prediction technique of bottom water gas condensate reservoir characterized by comprising
According to skin effect influence factor, gas phase non-Darcy effects influence factor and retrograde condensation function influence factor, it is based on perforation The gas phase radial fluid flow productivity model of interval and the centripetal stream productivity model of the gas phase hemispherical of perforated zone lower part, establish bottom water The water breakthrough time computation model of gas condensate reservoir;
Obtain the actual production data of bottom water gas condensate reservoir;
Based on the water breakthrough time computation model of the bottom water gas condensate reservoir, according to the actual production data of bottom water gas condensate reservoir, in advance Survey the water breakthrough time of bottom water gas condensate reservoir.
2. the water breakthrough time prediction technique of bottom water gas condensate reservoir as described in claim 1, which is characterized in that according to skin effect Influence factor, gas phase non-Darcy effects influence factor and retrograde condensation function influence factor, the gas phase plane diameter based on perforated zone To the centripetal stream productivity model of gas phase hemispherical of stream productivity model and perforated zone lower part, when establishing the water breakthrough of bottom water gas condensate reservoir Between computation model, comprising:
Meet non-Darcy's law according to gas phase and water phase meets Darcy's law, establishes air water two phase fluid flow model;
According to air water two phase fluid flow model, liquid and gas percolation flow velocity relational expression is established;
It establishes the moving distance on the water cone vertex on original gas-water interface and reaches the time history form in shaft bottom;
According to the gas phase radial fluid flow productivity model for the perforated zone for considering gas non-Darcy effects and stratum skin effect, examine The centripetal stream productivity model of gas phase hemispherical for considering the perforated zone lower part of gas non-Darcy effects and stratum skin effect, establishes bottom The gas production of water gas condensate reservoir perforated zone lower part;
According to the moving distance on the water cone vertex on liquid and gas percolation flow velocity relational expression, original gas-water interface and reach shaft bottom Time relationship, the gas production of bottom water gas condensate reservoir perforated zone lower part and retrograde condensation function influence factor, establish bottom water condensation The water breakthrough time computation model of gas reservoir.
3. the water breakthrough time prediction technique of bottom water gas condensate reservoir as claimed in claim 2, which is characterized in that the air water two-phase Flow model in porous media is established according to following formula:
Wherein, Pg、PwThe respectively pressure of gas phase and water phase, MPa;μg、μwThe respectively viscosity of gas phase and water phase, mPas;Vg、 VwThe respectively percolation flow velocity of gas phase and water phase, m/d;Kg、KwThe respectively effective permeability of gas phase and water phase, mD;R is radial Radius, m;β is non-Darcy flow coefficient, m-1;ρgFor the density of gas phase, g/cm3
4. the water breakthrough time prediction technique of bottom water gas condensate reservoir as claimed in claim 3, which is characterized in that the liquid phase is gentle Phase percolation flow velocity relational expression is established according to following formula:
Wherein, MwgFor aqueous vapor mobility ratio, Mwg=(Kww)/(Kgg)。
5. the water breakthrough time prediction technique of bottom water gas condensate reservoir as claimed in claim 4, which is characterized in that the original air water The moving distance on the water cone vertex on interface and the time history form for reaching shaft bottom are established according to following formula:
Wherein, φ is porosity, dimensionless;T is time, d;SwiFor original water saturation, dimensionless;SgrFor residual gas saturation Degree, dimensionless;SocFor remaining condensate saturation degree, dimensionless.
6. the water breakthrough time prediction technique of bottom water gas condensate reservoir as claimed in claim 5, which is characterized in that the consideration gas The gas phase radial fluid flow productivity model of the perforated zone of non-Darcy effects and stratum skin effect is established according to following formula:
Wherein, PeFor supply pressure, MPa;PwfFor bottom pressure, MPa;T is reservoir temperature, K;Z is deviation factors, dimensionless;hp For gas reservoir perforated zone thickness, m;Q1For the yield of perforated zone, m3/d;reFor drainage radius, m;rwFor wellbore radius, m;S is Skin factor, dimensionless;γgFor gas phase relative density, dimensionless;D1For inertia coeffeicent, (m3·d-1)-1
The centripetal miscarriage energy of gas phase hemispherical of the perforated zone lower part for considering gas non-Darcy effects and stratum skin effect Model is established according to following formula:
Wherein, D2For inertia coeffeicent, (m3·d-1)-1;KsFor effective permeability,mD;KvTo be seeped in vertical direction Saturating rate, mD;KhFor permeability in horizontal direction, mD;Q2For the gas production of gas reservoir perforated zone lower part, m3/d;
The gas production of bottom water gas condensate reservoir perforated zone lower part is established according to following formula:
Wherein, A, B, C are coefficient of dynamics, and Q is the total gas production of bottom water gas condensate reservoir, Q=Q1+Q2
7. the water breakthrough time prediction technique of bottom water gas condensate reservoir as claimed in claim 6, which is characterized in that the bottom water condensation The water breakthrough time computation model of gas reservoir is established according to following formula:
Wherein, tbtFor the water breakthrough time of bottom water gas condensate reservoir, d;F, G, D, E are coefficient of dynamics;hbFor bottom water gas condensate reservoir perforation The thickness of interval lower part, m;DC is condensate increment or reduction amount in unit volume condensate gas, m3/m3;PRTo be primitively laminated Power, MPa;BgFor the volume factor of gas, dimensionless;H is gas reservoir thickness, m;rbFor condensate saturation degree critical in reservoir gaps Obstruction radius, m;Y is the retrograde condensation factor, m3/(m3·MPa)。
8. a kind of water breakthrough time prediction meanss of bottom water gas condensate reservoir characterized by comprising
The water breakthrough time computation model of bottom water gas condensate reservoir establishes module, for according to skin effect influence factor, gas phase is non-reaches Western effects factor and retrograde condensation function influence factor, gas phase radial fluid flow productivity model and perforation based on perforated zone The centripetal stream productivity model of the gas phase hemispherical of interval lower part, establishes the water breakthrough time computation model of bottom water gas condensate reservoir;
Actual production data acquisition module, for obtaining the actual production data of bottom water gas condensate reservoir;
The water breakthrough time prediction module of bottom water gas condensate reservoir calculates mould for the water breakthrough time based on the bottom water gas condensate reservoir Type predicts the water breakthrough time of bottom water gas condensate reservoir according to the actual production data of bottom water gas condensate reservoir.
9. the water breakthrough time prediction meanss of bottom water gas condensate reservoir as claimed in claim 8, which is characterized in that the bottom water condensation The water breakthrough time computation model of gas reservoir is established module and is specifically used for:
Meet non-Darcy's law according to gas phase and water phase meets Darcy's law, establishes air water two phase fluid flow model;
According to air water two phase fluid flow model, liquid and gas percolation flow velocity relational expression is established;
It establishes the moving distance on the water cone vertex on original gas-water interface and reaches the time history form in shaft bottom;
According to the gas phase radial fluid flow productivity model for the perforated zone for considering gas non-Darcy effects and stratum skin effect, examine The centripetal stream productivity model of gas phase hemispherical for considering the perforated zone lower part of gas non-Darcy effects and stratum skin effect, establishes bottom The gas production of water gas condensate reservoir perforated zone lower part;
According to the moving distance on the water cone vertex on liquid and gas percolation flow velocity relational expression, original gas-water interface and reach shaft bottom Time relationship, the gas production of bottom water gas condensate reservoir perforated zone lower part and retrograde condensation function influence factor, establish bottom water condensation The water breakthrough time computation model of gas reservoir.
10. the water breakthrough time prediction meanss of bottom water gas condensate reservoir as claimed in claim 9, which is characterized in that the bottom water is solidifying The water breakthrough time computation model of gassing hiding is established module and is specifically used for:
The air water two phase fluid flow model is established according to following formula:
Wherein, Pg、PwThe respectively pressure of gas phase and water phase, MPa;μg、μwThe respectively viscosity of gas phase and water phase, mPas;Vg、 VwThe respectively percolation flow velocity of gas phase and water phase, m/d;Kg、KwThe respectively effective permeability of gas phase and water phase, mD;R is radial Radius, m;β is non-Darcy flow coefficient, m-1;ρgFor the density of gas phase, g/cm3
11. the water breakthrough time prediction meanss of bottom water gas condensate reservoir as claimed in claim 10, which is characterized in that the bottom water is solidifying The water breakthrough time computation model of gassing hiding is established module and is specifically used for:
The liquid and gas percolation flow velocity relational expression is established according to following formula:
Wherein, MwgFor aqueous vapor mobility ratio, Mwg=(Kww)/(Kgg)。
12. the water breakthrough time prediction meanss of bottom water gas condensate reservoir as claimed in claim 11, which is characterized in that the bottom water is solidifying The water breakthrough time computation model of gassing hiding is established module and is specifically used for:
It is closed according to the time that following formula establishes the moving distance on the water cone vertex on the original gas-water interface and reaches shaft bottom It is formula:
Wherein, φ is porosity, dimensionless;T is time, d;SwiFor original water saturation, dimensionless;SgrFor residual gas saturation Degree, dimensionless;SocFor remaining condensate saturation degree, dimensionless.
13. the water breakthrough time prediction meanss of bottom water gas condensate reservoir as claimed in claim 12, which is characterized in that the bottom water is solidifying The water breakthrough time computation model of gassing hiding is established module and is specifically used for:
The gas phase plane of the perforated zone of the consideration gas non-Darcy effects and stratum skin effect is established according to following formula Radial flow productivity model:
Wherein, PeFor supply pressure, MPa;PwfFor bottom pressure, MPa;T is reservoir temperature, K;Z is deviation factors, dimensionless;hp For gas reservoir perforated zone thickness, m;Q1For the yield of perforated zone, m3/d;reFor drainage radius, m;rwFor wellbore radius, m;S is Skin factor, dimensionless;γgFor gas phase relative density, dimensionless;D1For inertia coeffeicent, (m3·d-1)-1
The gas phase of the perforated zone lower part of the consideration gas non-Darcy effects and stratum skin effect is established according to following formula The centripetal stream productivity model of hemispherical:
Wherein, D2For inertia coeffeicent, (m3·d-1)-1;KsFor effective permeability,mD;KvTo be seeped in vertical direction Saturating rate, mD;KhFor permeability in horizontal direction, mD;Q2For the gas production of gas reservoir perforated zone lower part, m3/d;
The gas production of bottom water gas condensate reservoir perforated zone lower part is established according to following formula:
Wherein, A, B, C are coefficient of dynamics, and Q is the total gas production of bottom water gas condensate reservoir, Q=Q1+Q2
14. the water breakthrough time prediction meanss of bottom water gas condensate reservoir as claimed in claim 13, which is characterized in that the bottom water is solidifying The water breakthrough time computation model of gassing hiding is established module and is specifically used for:
The water breakthrough time computation model of the bottom water gas condensate reservoir is established according to following formula:
Wherein, tbtFor the water breakthrough time of bottom water gas condensate reservoir, d;F, G, D, E are coefficient of dynamics;hbFor bottom water gas condensate reservoir perforation The thickness of interval lower part, m;DC is condensate increment or reduction amount in unit volume condensate gas, m3/m3;PRTo be primitively laminated Power, MPa;BgFor the volume factor of gas, dimensionless;H is gas reservoir thickness, m;rbFor condensate saturation degree critical in reservoir gaps Obstruction radius, m;Y is the retrograde condensation factor, m3/(m3·MPa)。
CN201810649707.3A 2018-06-22 2018-06-22 The water breakthrough time prediction technique and device of bottom water gas condensate reservoir Pending CN109033518A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810649707.3A CN109033518A (en) 2018-06-22 2018-06-22 The water breakthrough time prediction technique and device of bottom water gas condensate reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810649707.3A CN109033518A (en) 2018-06-22 2018-06-22 The water breakthrough time prediction technique and device of bottom water gas condensate reservoir

Publications (1)

Publication Number Publication Date
CN109033518A true CN109033518A (en) 2018-12-18

Family

ID=64609960

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810649707.3A Pending CN109033518A (en) 2018-06-22 2018-06-22 The water breakthrough time prediction technique and device of bottom water gas condensate reservoir

Country Status (1)

Country Link
CN (1) CN109033518A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109884726A (en) * 2019-03-07 2019-06-14 中国石油大学(北京) The gas breakthrough time prediction technique and device of gas-drive pool
CN111364979A (en) * 2020-03-23 2020-07-03 中国石油大学(华东) Underground gas invasion monitoring system based on ultrasonic waves
CN112446118A (en) * 2019-08-13 2021-03-05 中国石油化工股份有限公司 Prediction method for water breakthrough time of water gas reservoir at bottom of complex biological reef
CN113356817A (en) * 2021-07-19 2021-09-07 中海油田服务股份有限公司 Perforation optimization design method for bottom water reservoir horizontal well based on prolonging anhydrous oil production period

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2293214C2 (en) * 2005-01-25 2007-02-10 Иван Яковлевич Клюшин Method of action on pre-bottom zone of well in hydrocarbon field with bottom water and recovery of oil and water by pumps-compressors with separate intake for coneless operation of well
CN106499371A (en) * 2015-09-06 2017-03-15 中国石油天然气股份有限公司 A kind of standing column well gas injection improves the method and device of Condensate Gas Reservoir recovery ratio
CN107045671A (en) * 2017-03-22 2017-08-15 重庆科技学院 Water-producing gas well hydrops Risk Forecast Method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2293214C2 (en) * 2005-01-25 2007-02-10 Иван Яковлевич Клюшин Method of action on pre-bottom zone of well in hydrocarbon field with bottom water and recovery of oil and water by pumps-compressors with separate intake for coneless operation of well
CN106499371A (en) * 2015-09-06 2017-03-15 中国石油天然气股份有限公司 A kind of standing column well gas injection improves the method and device of Condensate Gas Reservoir recovery ratio
CN107045671A (en) * 2017-03-22 2017-08-15 重庆科技学院 Water-producing gas well hydrops Risk Forecast Method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
明瑞卿 等: "边水凝析气藏高产井见水时间预测新模型" *
黄全华 等: "带隔板底水气藏见水时间预测方法" *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109884726A (en) * 2019-03-07 2019-06-14 中国石油大学(北京) The gas breakthrough time prediction technique and device of gas-drive pool
CN109884726B (en) * 2019-03-07 2020-08-04 中国石油大学(北京) Gas-drive reservoir gas-visible time prediction method and device
CN112446118A (en) * 2019-08-13 2021-03-05 中国石油化工股份有限公司 Prediction method for water breakthrough time of water gas reservoir at bottom of complex biological reef
CN112446118B (en) * 2019-08-13 2022-08-12 中国石油化工股份有限公司 Prediction method for water breakthrough time of water gas reservoir at bottom of complex biological reef
CN111364979A (en) * 2020-03-23 2020-07-03 中国石油大学(华东) Underground gas invasion monitoring system based on ultrasonic waves
CN111364979B (en) * 2020-03-23 2023-05-23 中国石油大学(华东) Underground gas invasion monitoring system based on ultrasonic waves
CN113356817A (en) * 2021-07-19 2021-09-07 中海油田服务股份有限公司 Perforation optimization design method for bottom water reservoir horizontal well based on prolonging anhydrous oil production period

Similar Documents

Publication Publication Date Title
CN109033518A (en) The water breakthrough time prediction technique and device of bottom water gas condensate reservoir
CN111236931B (en) Method and system for generating unsteady well testing chart of gas reservoir highly-deviated well
CN105354639B (en) The complete period PRODUCTION FORECASTING METHODS and its device of fine and close oil multi-dielectric seepage coupling
CN105089595B (en) Numerical reservoir simulation method and device under the effect of horizontal fractures fracture guide
CN105089582B (en) Numerical reservoir simulation method and device based on downhole flow control device
CN104879104B (en) Oil reservoir water injection method
CN108763657A (en) Horizontal gas well continuously takes liquid critical flow and determines method and device
CN112434426B (en) Shale gas multistage fracturing horizontal well step gradient pressure drop development method and device
CN106156439B (en) A kind of heavy crude reservoir dissolved gas drive method for numerical simulation with foam oil phenomenon
CN109933860A (en) Shale gas pressure break horizontal well cluster spacing optimization method
CN104834807A (en) Stress sensitive reservoir stratum relative permeability calculation method based on fractal theory
CN108131122A (en) Improve the CO2 amounts of sealing up for safekeeping and the method for oil recovery factor
CN105899756A (en) Optimizing flow control device properties on both producer and injector wells in coupled injector-producer liquid flooding systems
CN108280312B (en) A kind of horizontal well in segments design method for taking into account control water and sand control
CN109958404B (en) Bottom water coning and water plugging profile control agent parameter calculation method for bottom water reservoir
CN114427432B (en) Method for determining development potential of residual gas in gas reservoir
CN108729912B (en) Yield splitting method suitable for numerical reservoir simulation
CN106529184A (en) Method for calculating productivity of water-producing gas well of tilted water-bearing gas reservoir
CN108830410A (en) The water breakthrough time prediction technique and device of hyposmosis bottomwater gas field
CN107869309B (en) Track control method and device in the case of small hole deviation
CN109611073A (en) A kind of densification grease horizontal well mud stone can the preferred method and device in fractured layer position
CN105765161A (en) Optimizing flow control device properties on injector wells in liquid flooding systems
CN112836438B (en) Method and device for predicting water breakthrough time of gas well of heterogeneous edge water-gas reservoir
AU2013405168B2 (en) Optimizing flow control device properties on a producer well in coupled injector-producer liquid flooding systems
CN112377182A (en) Method and device for determining parameters of large-size hole crack type carbonate rock reservoir body

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20181218