CN110288258A - A kind of high water-cut reservoir Tapping Residual Oil method - Google Patents

A kind of high water-cut reservoir Tapping Residual Oil method Download PDF

Info

Publication number
CN110288258A
CN110288258A CN201910589725.1A CN201910589725A CN110288258A CN 110288258 A CN110288258 A CN 110288258A CN 201910589725 A CN201910589725 A CN 201910589725A CN 110288258 A CN110288258 A CN 110288258A
Authority
CN
China
Prior art keywords
field
oil
reservoir
evaluation
flow
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
CN201910589725.1A
Other languages
Chinese (zh)
Other versions
CN110288258B (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.)
China Petroleum and Chemical Corp
Sinopec Jiangsu Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Sinopec Jiangsu Oilfield Co
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 Chemical Corp, Sinopec Jiangsu Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201910589725.1A priority Critical patent/CN110288258B/en
Publication of CN110288258A publication Critical patent/CN110288258A/en
Application granted granted Critical
Publication of CN110288258B publication Critical patent/CN110288258B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK 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/25—Methods for stimulating production
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The present invention relates to a kind of high water-cut reservoir Tapping Residual Oil methods, it include: the flow field geology Dominated Factors by analyzing different oil reservoirs, determine oil reservoir flow field quantitatively characterizing method, and then determine the initial flow-field regularity of distribution and mode, in conjunction with the latter period of oilfield development analysis advantage such as development behavior Flow Field Distribution rule and mode, advantage flow field-initial flow-field method of difference is finally utilized, finds remaining oil, and exploit to remaining oil.High water-cut reservoir Tapping Residual Oil method provided by the invention, on the basis of classification of oil reservoirs, analyse in depth the influence of reservoir, construction feature to initial flow-field, by being reanalysed to development late stage oil reservoir flow field, comparative benefits flow field-initial flow-field, form the understanding to exploitation late residual oil distribution Favorable Areass, targeted Remaining Oil, improve the development degree and ultimate recovery of oil reservoir reserves out of control, effectively reduce Production Decline Prediction of Oilfield amplitude, hence it is evident that improve oil reservoir ultimate recovery.

Description

A kind of high water-cut reservoir Tapping Residual Oil method
Technical field
The invention belongs to oil-gas field development technical fields, and in particular to a kind of high water-cut reservoir Tapping Residual Oil method.
Background technique
(1) connotation in oil reservoir flow field
The static distribution feature and its dynamic rule of fluid and fluid reservoir are collectively referred to as oil reservoir flow field in oil reservoir. Advantage flow field generally refers to the flow field intensity in porous media body and is apparently higher than the intracorporal stream of the comparable neighbouring medium of volume therewith Field intensity.The apparent feature in advantage flow field drain amount between injection-production well is larger, and drain amount and water injection well swept volume between well Ratio is also larger.Furthermore, it is also possible to be accompanied by the features such as fluid flow rate is fast, amount of water injected is big, moisture content is high.What initial flow-field referred to It is oil reservoir flow field situation when oil field is untapped.
(2) high water cut development stage oil reservoir Study on Flow Field status
Currently, the research to oil reservoir initial flow-field more lacks.For grinding for high water cut development stage oil field advantage flow field The method of studying carefully specifically includes that oil-water well dynamic-analysis method, logging interpretation method, well-logging method, tracer method and Later development Fuzzy comprehensive evoluation mathematical method etc..
Interwell tracer numerical simulation technology: being in the water phase one or more of tracers being injected into water injection well, Then in a period of simulation trial, tracer trace routine is tracked in three-dimensional grid space, and then determines water filling well water The moving situation of phase component.By the output of the 3-D graphic to water injection well tracer each in each layer institute affected area, can analyze Out the connection situation of injection-production well, involve range, employ in plane and non-producing status, and the size and range of adjustment potentiality. Then the tracer for drawing each producing well in each well group adds up quantum of output and the relation curve between the production time.Due to tracer Producing well is always charged into first along high permeability zone or macropore, while being observed producing well tracer and being charged into sequencing, and Tracer changes in the cumulative production amount of each producing well, and the tracer cumulative amount curve rate of climb is fast, is advantage flow field direction, most Shown afterwards by calculating water injection well to each producing well tracer partition coefficient quantitatively characterizing advantage flow field size.
Streamline simulation method: basic ideas are the pressure field and velocity field for first finding out fluid in porous media, then seek stream Line (flow trace of fluid and time), finally finds out the remaining oil saturation value of any point in any streamline.The streamline model is not There are numerical dispersions and unstability, faster than traditional numerical simulation speed, allow node more, can fully demonstrate fine geology mould Influence of the type to computational fluid dynamics, flow characteristics reaches fine description waterflooding reservoir high water-cut stage remaining oil distributed in three dimensions Purpose.
Fuzzy comprehensive evaluation method: current high water cut development stage numerical model reservoir porosity, permeability are extracted first, is contained Then water saturation etc. handles parameter normalization by the method for normalized to obtain the blurring mapping square of single grid Battle array, and the weight matrix that contribution determines each factor is formed for advantage flow field according to the above Selected Factors, finally by square Blurring mapping matrix × weight matrix is obtained the advantage flow field fuzzy evaluation value of single grid by battle array multiplication.By to the block The fuzzy evaluation value that each grid is obtained after all grids progress fuzzy evaluations generates research block advantage flow field evaluation of estimate distribution Figure, and range restriction is carried out by the judge value to the whole district, finally obtain advantage Flow Field Distribution.
(3) importance of high water-cut reservoir Study on Flow Field
China's most oilfields entered in High water cut mining phase, waterflooding extraction ratio is high, meanwhile, main force's oil at this stage Field has carried out the subdivision and adjustment of excessive round and note adopts Perfect the work, continues the potentiality very little for segmenting infilling adjustment from now on.But Remaining oil enrichment region is still had in high water-cut reservoir, finds the emphasis that these remaining oil enriched area are lower part developments. Due to lack be directed to undeveloped reserves initial flow-field characterization scheme, cause to exploitation the late residual oil regularity of distribution understanding compared with For difficulty, cause development late stage to Tapping Residual Oil difficulty.The characterizing method of undeveloped reserves initial flow-field is urgently to solve at present Technological difficulties certainly.
The prior art lacks the research for being directed to undeveloped reserves initial flow-field, for remaining oil distribution area flow field change feelings Condition lacks understanding, this seriously affects the adjustment of oil-gas field development tail period technical policy and Tapping Residual Oil.From current From the point of view of domestic and international research conditions, more lack for the evaluation method of untapped period initial flow-field, when evaluation influences initial flow The many factors of field are difficult to evaluate simultaneously；Meanwhile for development late stage oil reservoir, lack in conjunction with behavioral characteristics and other means of production The method for carrying out comprehensive analysis.
Summary of the invention
For above-mentioned problems of the prior art, it can avoid above-mentioned skill occur the purpose of the present invention is to provide one kind The high water-cut reservoir Tapping Residual Oil method of art defect.
In order to achieve the above-mentioned object of the invention, technical solution provided by the invention is as follows:
A kind of high water-cut reservoir Tapping Residual Oil method, comprising: oil reservoir initial flow-field evaluation；The evaluation of oil reservoir advantage flow field.
Further, which comprises
By analyzing the flow field geology Dominated Factors of different oil reservoirs, oil reservoir flow field quantitatively characterizing method is determined, and then determine The initial flow-field regularity of distribution and mode, in conjunction with the latter period of oilfield development analysis advantage such as development behavior Flow Field Distribution rule and Mode finally utilizes advantage flow field-initial flow-field method of difference, finds remaining oil, and exploit to remaining oil；It specifically includes:
(1) binding area geologic feature, preferably reservoir heterogeneity, coefficient of permeability variation, crack, glue at permeability The evaluation index of native content, viscosity of crude, stratigraphic dip as initial flow-field；Using the fuzzy synthesis for solving the problems, such as Fuzzy Engineering Law theory is judged to establish the overall evaluation system of initial flow-field, and determines each evaluation criterion weight using analytic hierarchy process (AHP) Coefficient；
(2) high water-cut reservoir is divided into four classes, respectively with respect to integrity reservoir, narrow strip oil reservoir, lithological pool and extremely multiple Miscellaneous fault block oil reservoir；Utilize the method for dynamically associating, fuzzy comprehensive evaluation method, streamline simulation fado method comprehensive distinguishing development late stage advantage Flow field；
(3) it is directed to mid-later development phase high water-cut reservoir, using advantage flow field-initial flow-field method of difference, finds remaining oil.
Further, oil reservoir initial flow-field evaluation includes: by evaluation index to final evaluation goal parameter influence degree Ambiguity is indicated with fuzzy set form, is formed the evaluations matrix that can directly carry out judging operation, is obtained one by blurring mapping The evaluation result of a fuzzy set.
Further, the evaluation of oil reservoir initial flow-field includes: to choose the factor compositing factor for influencing effect of reservoir development first Collection, then determine the Comment gathers for evaluating its result, then determine each factor for the weight sets of evaluation result；It is chosen by tentative calculation Suitable membership function then calculates the synthetic evaluation matrix of evaluation object；Mould is carried out to weight sets and synthetic evaluation matrix Paste transformation, obtains comprehensive evaluation result.
Further, oil reservoir initial flow-field, which is evaluated, includes:
(1) set of object is determined；The set of the determination each layer of each well to be calculated in block to be evaluated, is denoted as:
A=(A1, A2, A3... An)；
(2) set of index factor, the i.e. set of each attribute of object are determined；Attribute set refers to influence initial flow-field Various heterogeneity parameters, be denoted as:
P=(P1, P2, P3... Pm)；
(3) determine that index matrix, matrix element are property parameters；
(4) index value, element a are normalizedijIt can be normalized as the following formula as bij
After providing each evaluation index, nondimensionalization need to be carried out to each parameter, to reject the influence between different dimensions；
Parameter is normalized；Matrix after normalization are as follows:
(5) determination of weight
The corresponding weight of each index is provided, weight vectors are obtained:
W=(w1, w2, w3, L, wm)；
(6) initial flow-field evaluation number I is calculated
I=BW；
Acquire initial flow-field evaluation number vector I=(i1, i2, i3, L in), wherein (0,1) i ∈, I is bigger, and flow field is got over By force；I is smaller, and flow field is weaker.
Further, evaluation index is comprised determining that in the evaluation of oil reservoir initial flow-field；Standardize initial data；Determine weight； Calculate initial flow-field intensity index；Determine initial flow-field discriminant classification standard.
Further, the evaluation of oil reservoir advantage flow field includes: and utilizes to dynamically associate method, fuzzy comprehensive evaluation method, streamline simulation Three kinds of methods of method evaluate research area.
Further, the method for dynamically associating includes:
(1) analytical sequence is determined；
On the basis of qualitative analysis, a dependent variable factor and multiple independent variable factors are determined；If dependent variable data structure At reference sequences X0', each argument data composition ratio is compared with sequence Xi' (i=1,2 ... m), and m+1 data sequence forms following square Battle array:
Wherein Xi '=(xi ' (1), xi ' (2) ..., xi ' (n))T, the m of i=0,1,2 ...；N is the length of Variables Sequence.
(2) nondimensionalization is carried out to Variables Sequence；
Respectively because prime sequences form following matrix after nondimensionalization:
Nondimensionalization is carried out using first value method:
(3) difference sequence, maximum difference and lowest difference are asked
Form absolute difference matrix:
Wherein: △0i(k)=| x0(k)-xi(k)|
I=0,1,2 ..., m；K=1,2 ..., n；
Maximum number and minimum number are maximum difference and lowest difference in absolute difference battle array:
(4) calculate correlation coefficient
Such as down conversion is made to data in absolute difference battle array:
Obtain incidence coefficient matrix:
Resolution ratio ρ is in (0,1) interior value in formula；
Incidence coefficient ξ0i(k) 1 positive number, Δ are no more than0i(k) smaller, ξ0i(k) bigger；
(5) calculating correlation
Averaging obtains XiWith X0The degree of association:
The degree of association is adopted by calculating note and describes the power of relationship between water injection well and producing well, judges whether to generate between well excellent Potential-flow field.
Further, fuzzy comprehensive evaluation method includes:
Permeability and clay content are chosen as static evaluation index；
For water injection well, three daily water-injection rate, water filling oil pressure and unit thickness cumulative water injected indexs are chosen as dynamic Evaluation index；For producing well, daily fluid production rate, moisture content, unit thickness cumulative liquid production are chosen as dynamic assessment index；
Nondimensionalization processing is carried out to parameter；
Using Weight of Coefficient through Analytic Hierarchy Process:
Firstly, determining permeability and shale content to the relative importance for forming advantage flow field according to each influence factor Weight constitutes influence factor weight vectors [ω1,ω2]T, guarantee that the sum of weight is equal to 1；
Secondly, according to oil well, water injection well dynamic performance characteristic to the relative importance in advantage flow field, determine each index Weight；
Both finally, according to the importance that the influence factor in advantage flow field and performance characteristic identify existence, determine Weight [ωa,ωb], guarantee that the sum of weight is equal to 1.
Further, initial flow-field discriminant classification standard are as follows:
If composite index I < 0.4, for weak current field；If composite index 0.4≤I < 0.6, for middle flow field；If composite index 0.6≤I < 0.8 is then high current field；If composite index 0.8≤I < 1, for advantage flow field.
High water-cut reservoir Tapping Residual Oil method provided by the invention analyses in depth reservoir on the basis of classification of oil reservoirs Influence of the heterogeneous, construction feature to initial flow-field, by being reanalysed to development late stage oil reservoir flow field, comparative benefits flow field- Initial flow-field, forms the understanding to exploitation late residual oil distribution Favorable Areass, and targeted Remaining Oil improves oil The development degree and ultimate recovery for hiding reserves out of control, effectively reduce Production Decline Prediction of Oilfield amplitude, extend the stable production period in oil field Limit, hence it is evident that improve oil reservoir ultimate recovery, achieve obvious development effectiveness and benefit, practical application can be met well It needs.
Detailed description of the invention
Fig. 1 is each factor stream field percentage contribution histogram；
Fig. 2 is G6 fault block substratum initial flow-field distribution map；
Fig. 3 is G6 fault block substratum remaining oil saturation figure；
Fig. 4 is Z67 well intake profile figure；
Fig. 5 is G6-39 well, the flat 3 well production curve figures of G6；
Fig. 6 is that G6 fault block substratum dynamically associates method predominant current field intensity distribution map；
Fig. 7 is G6 fault block substratum fuzzy comprehensive evaluation method predominant current field intensity distribution map；
Fig. 8 is the phase percolation curve of G6-104 well difference displacement multiple；
Fig. 9 is middle and high infiltration core permeability and pressue-graph；
Figure 10 is the matched curve figure of G6 fault block comprehensive water cut；
Figure 11 is the matched curve figure of high 6 fault block day oil-producing；
Figure 12 is high 6 fault block substratum streamline simulation method predominant current field intensity distribution map.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawing and specific implementation The present invention will be further described for example.It should be appreciated that described herein, specific examples are only used to explain the present invention, and does not have to It is of the invention in limiting.Based on the embodiments of the present invention, those of ordinary skill in the art are not making creative work premise Under every other embodiment obtained, shall fall within the protection scope of the present invention.
Sandstone oil reservoir classification is completed by rock core information, chemical examination data, drilling data etc. for oil reservoir untapped period Flow field geology Dominated Factors are analyzed in evaluation, are summarized oil reservoir flow field quantitatively characterizing method, are completed the initial flow-field research of oil reservoir；Again Advantage flow field is identified by a variety of methods, is summarized advantage Flow Field Distribution rule and mode, is completed the table in development late stage advantage flow field Sign；Eventually by comparison initial flow-field and development late stage flow field, remaining oil enriched area is found.
Conducted a research by the evaluation of initial flow-field to high water-cut reservoir and the evaluation of advantage flow field, using initial flow-field with it is excellent Existing difference between potential-flow field, excavation initial flow-field intensity is higher, but the remaining oil of predominant current field intensity lower region.
A kind of high water-cut reservoir Tapping Residual Oil method, comprising the following steps:
Step 1: oil reservoir initial flow-field is evaluated
1. oil reservoir initial flow-field evaluation method principle
Fuzzy comprehensive evaluation method is to use ambiguity of the evaluation index to final evaluation goal parameter influence degree to obscure Aggregate form indicates, forms the evaluations matrix that can directly carry out judging operation, obtains a fuzzy set by blurring mapping Evaluation result.The essence of fuzzy comprehensive evaluation method is built upon a kind of fuzzy linear transformation in the fuzzy foundation of mathematics.
The factor compositing factor collection for influencing effect of reservoir development is chosen first, then determines the Comment gathers for evaluating its result, so Determine each factor for the weight sets of evaluation result afterwards.Suitable membership function is chosen by tentative calculation, then calculates evaluation The synthetic evaluation matrix of object.Blurring mapping is carried out to weight sets and synthetic evaluation matrix, obtains comprehensive evaluation result.
(1) set of object is determined.The set of the determination each layer of each well to be calculated in block to be evaluated, is denoted as:
A=(A1, A2, A3... An) (1)
(2) set of index factor, the i.e. set of each attribute of object are determined.Here attribute set refers to that influence is initial The various heterogeneity parameters in flow field, such as permeability coefficient of variation parameter set, be denoted as:
P=(P1, P2, P3... Pm) (2)
(3) determine that index matrix, matrix element are property parameters.
A is also referred to as AT=[A1, A2, A3, L, An], wherein
(4) index value, element a are normalizedijIt can be normalized as the following formula as bij
After providing each evaluation index, nondimensionalization need to be carried out to each parameter, to reject the influence between different dimensions.This Secondary selection is maximum, minimum value method carries out nondimensionalization processing to each parameter.
It is bigger for the value of parameter, the stronger parameter in characterization flow field, such as permeability parameter, it is normalized using formula 4； Bigger for the value of parameter, the weaker parameter in characterization flow field is normalized using formula 5.Matrix after normalization are as follows:
(5) determination of weight
The corresponding weight of each index is provided, weight vectors are obtained:
W=(w1, w2, w3, L, wm) (7)
(6) initial flow-field evaluation number I is calculated
I=BW (8)
Acquire initial flow-field evaluation number vector I=(i1, i2, i3, L in), wherein (0,1) i ∈, I is bigger, and flow field is got over By force；I is smaller, and flow field is weaker.
2. the selection of evaluation index
The screening of evaluation index should follow following principle:
(1) importance
The index filtered out has a conclusive effect to the differentiation in oil reservoir flow field, rather than those dispensable fingers Mark.
(2) independence
There is opposite independence, guarantee will not filter out repetition and lengthy and jumbled index between the index filtered out.
(3) practicability
The index filtered out has certain operability, be as far as possible can directly survey index (such as log well hole infiltration) or By being simply calculated.
(4) comprehensive
Selected index should describe the influence factor of initial flow-field comprehensively as far as possible, cannot omit.
Above four principles are followed in oil reservoir initial flow-field evaluation index screening process, can be obtained both has representative Property, independence, while again more practical, succinct index system.
According to the above principle, determine on initial flow-field influence apparent reservoir heterogeneity, permeability, sedimentary facies, crack, Clay content, viscosity of crude, stratigraphic dip are evaluation index (table 1, Fig. 1).
1. reservoir heterogeneity (H)
Reservoir heterogeneity reflects the variation of sand body permeability, and reservoir heterogeneity will lead to waterline and promote uneven, cause Keep remaining oil distribution uneven, or injecting water and advancing by leaps and bounds makes oil well sudden water flooding influence development effectiveness, directly affects water flood field Water flooding effectiveness.
2. permeability (K)
Permeability is the principal element of grease plane motion, and the principal element of control Flow Field Distribution.Permeability is higher Well, the flow for injecting water is bigger, more serious to washing away for rock, also advantage flow field more easy to form.Middle and high infiltration development block Permeability contrast at different well points is larger, is easy to form advantage flow field between the high well group of permeability.
3. sedimentary facies (S)
Depositional environment is the foundation factors for influencing reservoir storage and collection performance.Sedimentary micro is advantageous, and sand body is formed in hydrodynamic force item Part strength, is transformed by water flow, and sandstone maturity is high, and shale grey matter content is low, and reservoir properties are good.And the equal deposition in water channel lateral margin Facies tract, hydrodynamic condition is poor, and flow action is weak, and especially at shallow water calmness, deposited particles are thin, and shale content is high, reservoir properties Difference.
The high 6 fault block initial flow-field parameter value table of table 1
4. microcrack/pressure-break (F)
Due to the requirement of laboratory sample, crack is seeped etc. in hole and can not be embodied in physical property tests, but these cracks not only can be with Improve the connectivity and permeability of reservoir, controls oil well initial productivity and well water rises rule, rise and determine in seepage direction Sexual factor.The effect of pressure-break in the earth formation, be equivalent in homogeneous formation entrainment of a high permeable strip, fracture guide Ability is higher, and the heterogeneity of oil reservoir is stronger.According to the development degree in crack, difference value 0.8,0.6,0.4,0.2,0.Pressure break The fractured layer position value 1.0 of well.Numerical value is bigger, and influence of the crack to initial flow-field is bigger.
5. clay mineral (C)
Influence of the variation of clay mineral to waterflooding extraction mainly generates reservoir properties due to the variation of clay mineral It influences, to influence the efficiency of extraction and the water filling of oil gas indirectly, what clay mineral was filled the water washes away, easy water-swellable or fortune It moves, to change the permeability of reservoir and inject the direction of water.
6. viscosity of crude (V)
Viscosity of crude is bigger, and viscosity ratio of oil and water is bigger, with the increase of injection rate, the change of divergence of high and low permeable formation It is more obvious, injection water is easier in high permeability zone fingering, to exacerbate the inhomogeneities of oil reservoir Flow Field Distribution.Conversely, if crude oil Viscosity is smaller, and viscosity ratio of oil and water is smaller, then the effect of piston displacement is more obvious, and the uniformity of oil reservoir Flow Field Distribution can then obtain Improve.
7. stratigraphic dip (D)
There are the oil fields of stratigraphic dip, and the oil well below water injection well is early compared with the oil well water breakthrough being located above water injection well, Rate of water cut increase is very fast after water breakthrough；With the reduction of stratigraphic dip, gravity is gradually become smaller, when stratigraphic dip is less than 5 ° When, influence of the stratigraphic dip to waterflooding extraction can be ignored.
3. initial data standardizes
Assuming that there is each layer of each well of N mouthfuls of wells to participate in classification, each well has K with reference to oil deposit parameter, constitutes matrix:
For " becoming big " (value is the bigger the better) index, data normalization method:
For " becoming small " (value is the smaller the better) index, using standardized method:
4. weight determines
Determine that the common method of weight sets X has: statistic law, Solution of Fuzzy Relation Equation method, similar sample determine method, level Analytic approach etc. is several.Think that analytic hierarchy process (AHP) computational accuracy is higher, simple and easy to do by Optimization Analysis, it is theoretical mature, adaptive Property is good (when required index lacks a certain item or a few, as long as will remove a certain item or a few in the weight coefficient table acquired The weight coefficient of item, then normalizes again, can acquire new weight coefficient table), practicability is preferable.Therefore, using level Analytic approach determines weight sets.
Hierarchy Analysis Method is traditional weighing computation method, and this method is the beginning of the seventies by Univ. of Pittsburgh Professor A.L.Satty foundes.Complicated problem is divided into several levels by analytic hierarchy process (AHP), then according to certain extension Judgement indicates, i.e. described Judgement Matricies with regard to the relative importance quantitative of each level each element.It should by solving The Maximum characteristic root of matrix and corresponding feature vector determine the weight of each level each element relative importance.
For scheme evaluation, problem is divided into several levels, each level includes a quantitative target.According to visitor The judgement of sight compares and quantifies the relative importance between each level, determines weighted value by judgment matrix.
Analytic hierarchy process (AHP) (The Analytic Hierarchy Process, abbreviation AHP) is a kind of quantifies and qualitative phase In conjunction with the method that the subjective judgement of people is expressed and handled with Number Forms.
When carrying out decision with AHP, it can substantially be divided into 4 steps and carry out:
(1) hierarchical model
The relationship of each factor in analysis system, and the general objective of defining issue requirement and the factor for influencing general objective, and According to one recursive hierarchy structure of dominance relation constructed in groups between factor.
(2) judgment matrix is established
The lower layer's factors dominated to each factor are compared two-by-two, compare their important journeys to upper layer factor Degree, and is come out with digital representation, thus to the judgment matrix two-by-two of each upper layer factors composition one.
Provide giving a mark two-by-two according to the result of the actual development experience in oil field and expert discussion as a result, construction compares two-by-two Matrix.
(3) it solves judgment matrix and carries out consistency desired result, the underproof judgment matrix of consistency desired result is carried out again Assignment or amendment finally acquire lowest level factor to the synthetic weight of general objective.
1. the consistency of judgment matrix
A, judge consistency
If each element for judgment matrix meets aik×akj=aij(wherein i, j, k=1,2 ... n), then claim to sentence Disconnected matrix, which has, judges consistency, when constructing discrimination matrix, since the factor compared two-by-two is too many, be difficult formula this formula all at It is vertical.If judgment matrix A, if its each column vector has identical ranking results, A is referred to as to judge consistency matrix, in reality In operation, when judgment matrix, which has, judges consistency, it is believed that policymaker's thinking is clear orderly.Therefore, judge one Cause property is reasonable as a reference standard of acceptable matrix.
B, localized delivery matrix
Positive reciprocal matrix A=(aij) in n × n, if having a to any kik≥ajkOr to any k, there is aik≤ajk, then claim There is local sequence transitivity between i row element and j row element.Since single criterion decision is mainly used to choose the best alternatives, because If this largest component for acquiring right main characteristic vector is wt, and there is local sequence transitivity in t row and remaining each row element, then still Receive the judgment matrix.
2. the consistency desired result method of judgment matrix
A, column vector normalizes
Have following several method: row vector summation column vector normalizing, row vector even multiply out Nth power column vector normalizing.
B, Maximum characteristic root is sought
Approximate algorithm:
Calculate coincident indicator and relative uniformity index
CRA=CIA/RIACompared with 0.1 or 0.01, consistency is judged.
3. common judgment matrix error correcting method
Common method has following several:
A, Nth column coefficient returns 1, observes modification method.
B, judgment matrix is reconfigured, deviation maximal term modifies method.
C, judgment matrix is reconfigured, sum of square of deviations maximal term modifies method.
D, multi-expert weighted mean method.
(4) it is calculated by comparison element by judgment matrix for the relative weighting of the criterion.
There are many ways to calculating weight by judgment matrix, common handled with method, consistent for one judgement Matrix is exactly corresponding weight vectors after its each row normalization.It is similar to weight after each row normalization when A is inconsistent Vector and method are exactly the arithmetic mean of instantaneous value using this n column vector as weight vectors.Therefore have
Steps are as follows for its calculating:
Step 1: the element of A presses row normalization；Step 2: each column after normalization are added；Step 3: after will add up Vector divided by n up to weight vectors.
2 initial flow-field factor important ratio of table is compared with quantitative data table
Factor of evaluation Reservoir heterogeneity Permeability Sedimentary facies Crack Clay mineral Viscosity of crude Stratigraphic dip
Reservoir heterogeneity 1.00 2.00 3.00 4.00 5.00 6.00 7.00
Permeability 0.50 1.00 2.00 3.00 4.00 5.00 6.00
Sedimentary facies 0.33 0.50 1.00 2.00 3.00 4.00 5.00
Crack 0.25 0.33 0.50 1.00 2.00 3.00 4.00
Clay mineral 0.20 0.25 0.33 0.50 1.00 2.00 3.00
Viscosity of crude 0.17 0.20 0.25 0.33 0.50 1.00 2.00
Stratigraphic dip 0.14 0.17 0.20 0.25 0.33 0.50 1.00
3 influence factor index of table and its weight
Factor Weight
Reservoir heterogeneity H 0.35
Permeability K 0.24
Sedimentary facies S 0.16
Crack F 0.10
Clay mineral C 0.07
Viscosity of crude V 0.04
Stratigraphic dip D 0.03
Above-mentioned matrix is handled according to formula (12) respectively, and then calculates each level index for a upper level Weight.
Determine the relative importance that different factors influence evaluation goal parameter.For initial flow-field evaluation, seven Major class geologic(al) factor importance definitive result table 2.
According to the relative importance that different factors influence evaluation goal parameter, different factor normalized weights can be acquired Coefficient results are as shown in table 3.Reservoir heterogeneity and permeability account for principal element, remaining parameter influences relatively small.
5. initial flow-field intensity index (I) calculates
It is obtained just by the normalization of 7 quantitative parameters and the weight result of study of parameters selected in the studies above The calculation formula of beginning flow field intensity index I:
I=0.35H+0.24K+0.16S+0.1F+0.07C+0.07V+0.03D
By the calculating to initial flow-field intensity index R, the initial flow-field plane point of each substratum of each oil reservoir may be implemented The production of Butut, it will thus be seen that the size of initial flow-field intensity, and the advantageous of advantage Flow Field Distribution can be predicted accordingly Zone.
According to random consistency calculation method and test stone (random consistency ratio CR < 0.1), according to the weight of analysis The random consistency ratio CR=0.023 of coefficient illustrates that the weight coefficient value for the initial flow-field geologic(al) factor that this is determined is to compare Correctly.
The initial flow-field value of each well point series of strata is acquired according to comprehensive weight.High 6 fault block E1f2 3-1For non-active layers system, initially Flow field is low, and remaining oil saturation is high；E1f2 3-5For main force's series of strata, initial flow-field is high, and remaining oil saturation is low (Fig. 2, Fig. 3).Initially Flow field low value area, remaining oil are relatively enriched with, and initial flow-field can preferably reflect the flowing of water filling and the distribution of remaining oil.
6. initial flow-field discriminant classification standard
Advantage flow field discrimination standard is determined in conjunction with pertinent literature according to geology characteristic:
If composite index I < 0.4, for weak current field；
If composite index 0.4≤I < 0.6, for middle flow field；
If composite index 0.6≤I < 0.8, for high current field；
If composite index 0.8≤I < 1, for advantage flow field.
Step 2: oil reservoir advantage flow field is evaluated
Oil reservoir flow field constantly changes with the progress of oil field development.Mid-late oilfield development, reservoir is in injection water Long-term souring under, along with the gravitational differentiation of grease, the flow field intensity in local reservoir is apparently higher than neighbouring storage Layer forms the inefficient or inefficient cycle of injection water, influences Development Response of Oilfield.Therefore the stream field that needs to be association of activity and inertia is commented Valence identifies advantage flow field, to instruct the adjustment in flow field.From G6-2 well, Z67 well different times intake profile from the point of view of, with Waterflooding extraction gos deep into, and the water absorption of main water absorption layer is increasing, and advantage flow field gradually shows (Fig. 4).
By the investigation to domestic and foreign literature etc., the common recognition methods in advantage flow field mainly dynamically associates method, obscures The methods of Comprehensive Evaluation, streamline simulation method, well testing method, injection profile logging data, the monitoring of well head falloff curve, due to various Method is each advantageous, also there is its limitation, needs to consider that a variety of methods are mutually authenticated, overall merit.This research and utilization dynamic Three kinds of correlation method, fuzzy comprehensive evoluation, streamline simulation method methods evaluate research area.
1. dynamically associating method
After advantage flow field is formed, the injection rate of water injection well has stronger correlation with the Liquid output of producing well, uses Grey correlation theory adopts degree of association parameter with note to describe in same flooding unit the strong of relationship between water injection well and different producing wells It is weak.The degree of association is bigger, and flow field is stronger.
The Liquid output of the injection rate of water injection well and producing well has certain correlation, and note is adopted between the apparent well group of response Flow field is strong, and it is weak that note adopts flow field intensity between the weak well group of response.If G6-39 well is after 2017 stop note, corresponding oil well G6 puts down 3 wells Liquid measure decline is obvious, and after G6-39 well restores water filling within 2018, the Liquid output that G6 puts down 3 wells is also risen with it.Illustrate G6-39 well with It is strong (Fig. 5) that G6 equals dynamically associating property between 3 wells.
Grey correlation analysis (GRA) method is a kind of multifactor statistical analysis technique, between factor each in solving system Prevailing relationship, find out an important factor for influencing target value, to grasp the main feature of things, become to a System Development Change the quantitative description and comparative approach of situation.It using the sample data of each factor as foundation, described with grey relational grade it is each because The power of relationship, size between element.If sample data series reflects the situation (direction, size, speed etc.) of two factors variation Almost the same, then the degree of association between them is larger；Conversely, the degree of association is smaller.It mainly calculates step are as follows:
(1) analytical sequence is determined
On the basis of to studied a question qualitative analysis, a dependent variable factor and multiple independent variable factors are determined.If Dependent variable data constitute reference sequences X0', each argument data composition ratio is compared with sequence Xi' (i=1,2 ... m), m+1 data sequence Column form following matrix:
Wherein Xi '=(xi ' (1), xi ' (2) ..., xi ' (n))T, the m of i=0,1,2 ...；N is the length of Variables Sequence.
(2) nondimensionalization is carried out to Variables Sequence
Under normal circumstances, original variable sequence has different dimension or the order of magnitude, in order to guarantee to analyze the reliable of result Property, need to carry out Variables Sequence nondimensionalization, respectively because prime sequences form following matrix after nondimensionalization:
Nondimensionalization is carried out using first value method:
(3) difference sequence, maximum difference and lowest difference are asked
The absolute difference of first row (reference sequences) and remaining column (compared with sequence) in calculating formula (16) is formed following absolute Matrix of differences:
Wherein:
△0i(k)=| x0(k)-xi(k)|
I=0,1,2 ..., m；K=1,2 ..., n
Maximum number and minimum number are maximum difference and lowest difference in absolute difference battle array:
(4) calculate correlation coefficient
Such as down conversion is made to data in absolute difference battle array:
Obtain incidence coefficient matrix:
Resolution ratio ρ is in (0,1) interior value in formula, and under normal circumstances mostly in 0.1 to 0.5 value, ρ is smaller more to mention the pass G Contact several difference.
Incidence coefficient ξ0i(k) 1 positive number, Δ are no more than0i(k) smaller, ξ0i(k) bigger, it, which reflects i-th, compares sequence Arrange XiWith reference sequences X0In the correlation degree of kth phase.
(5) calculating correlation
Compare sequence XiWith reference sequences X0Correlation degree be to be come instead by N number of incidence coefficient (i.e. in (22) formula i-th column) It reflects, averaging can be obtained by XiWith X0The degree of association:
The degree of association is adopted by calculating note and describes the power of relationship between water injection well and producing well, judges whether to generate between well excellent Potential-flow field.Note by calculating G6-41 well group, which adopts the degree of association and the fltting speed of tracer, has good correlation (table 4), Illustrate that it is reasonable for adopting the degree of association for judging flow field power using note.
Table 4G 6-41 well group note adopts the degree of association
Producing well G6-39 G6-40 G6-73
The degree of association 1 0.78 0.56
Tracer speed, m/d 1.57 1.25 1.01
Creation data is split by the monitoring materials of two sections and assigns to single sand body, it is strong to calculate each sand body flow field of G6 fault block It spends (Fig. 6).Advantage flow field is mainly distributed between G6-19 well and G6-75 well, between G6-6 well and G6-59, G6-13 and G6-99 well Between, between G6-97 and G6-18, side G6-14 well, between G6-29 well and G6-108 well.
The identification for dynamically associating method progress advantage flow field does not account for geologic(al) factor to its shadow although parameter easily obtains It rings, creation data splits that point there are subjectivities, influences the evaluation in single sand body flow field.
2. fuzzy comprehensive evaluation method
The power in flow field is to influence each other caused by many factors, and between each factor, mutually restrict.In addition, each Factor is difficult to be stated with accurate concept to whole influence and correlation, that is to say, that entire evaluation procedure is non-linear 's.Due to advantage flow field forming process particular/special requirement and fuzzy comprehensive evoluation the characteristics of, here use fuzzy comprehensive evoluation The evaluation in method progress flow field.
The index for influencing flow field has Static implicit method and dynamic factor, and Static implicit method describes the geology to form advantage flow field Basis is the genetic prerequisite that advantage flow field is formed, and can show that the potentiality that advantage flow field is formed.The direct characterization in oil reservoir flow field refers to Mark is exactly dynamic indicator, and value reflects the size of oil reservoir flow field actual strength.
(1) Static implicit method
The heterogeneity on stratum itself leads to inject water during waterflooding extraction preferentially along high permeability zone position or structure Make low area flowing.This long-term unbalanced flowing causes the washing degree at Thief zone position obviously than hyposmosis position Washing degree it is high, and this species diversity is gradually expanding with the increase of injected slurry volume multiple.Water is injected also just along these Low-resistance position gradually forms advantage flow field.Permeability size, formation heterogeneity, sedimentary rhythm, porosity size, formation cementation Degree, core intersection, grease viscosity differences and clay content are the major reasons to form advantage flow field.
Permeability is to restrict the main factor of injectability, and the higher well of permeability, displacement flow quantity is bigger, to rock Wash away more serious, also advantage flow field more easy to form.Permeability contrast at middle and high infiltration development block difference well point is larger, more It is easy to form advantage flow field between the high well group of permeability.
What clay was filled the water washes away, easy water-swellable or migration, to change the permeability of reservoir and inject the side of water To.Permeability and clay content are thus chosen as static evaluation index.
(2) dynamic factor
After advantage flow field is formed between oil-water well, injection water just along predominant current place layer position and band flow into recover the oil Well.Due to the layer position of other hyposmosis and band seldom involving even without injection water, injection water flows in reservoir During the entire process of filtrational resistance all very littles, this has different degrees of embodiment in the dynamic indicator of oil-water well.
For water injection well, three daily water-injection rate, water filling oil pressure and unit thickness cumulative water injected indexs are chosen as dynamic Evaluation index.For producing well, three indexs such as daily fluid production rate, moisture content, unit thickness cumulative liquid production are chosen as dynamic Evaluation index.
(3) parameter nondimensionalization is handled
Index is divided into using " opposite optimization " by " more bigger more excellent type " and " smaller more excellent type ".According to the spy of practical problem Point, these parameters such as permeability, daily water-injection rate, moisture content are bigger, and reflection flow field is stronger, and trapezoidal profile has been selected to be analyzed. The water filling oil pressure of water injection well is smaller, and reflection flow field is stronger, is distributed using lower semi-trapezoid.
(4) weight determines
Index weights use analytic hierarchy process (AHP).Analytic hierarchy process (AHP) is a kind of qualitative analysis and quantitatively calculates the system combined Analysis method.Complicated problem is divided into several levels by it, and then basis is to certain extension judgement, according to each level The relative importance of each element provides the quantificational expression of its weight.
Firstly, according to each influence factor to the relative importance for forming advantage flow field, determine permeability and shale content this The weight of two indexs constitutes influence factor weight vectors [ω1,ω2]T, guarantee that the sum of weight is equal to 1, concrete outcome such as table 5.
The static influence factor index of table 5 and its weight
Static implicit method Weight
Permeability 0.55
Shale content 0.45
Secondly, according to oil well, water injection well dynamic performance characteristic to the relative importance in advantage flow field, determine each index Weight, concrete outcome such as table 6, table 7.
Each index of 6 oil well of table and its weight
Anomalous character parameter factor Weight
Daily fluid production rate 0.25
Moisture content 0.30
Unit thickness cumulative liquid production 0.45
Each index of 7 well of table and its weight
Water injection well dynamic parameter factor Weight
Daily water-injection rate 0.35
Fill the water oil pressure 0.35
Unit thickness cumulative water injected 0.30
Both finally, according to the importance that the influence factor in advantage flow field and performance characteristic identify existence, determine Weight [ωa,ωb], guarantee that the sum of weight is equal to 1, concrete outcome such as table 8.
The weight of 8 Static implicit method of table, dynamic factor
Project Weight
Static implicit method 0.45
Dynamic factor 0.55
By selected index and weight, according to Comprehensive Evaluation principle, (figure is identified to each well group advantage of G6 fault block flow field 7), advantage flow field be mainly distributed between G6-70 and CG6-14 well, between G6-12 and G6-111 well, between G6-64 and G6-2 well, G6-9 The regions such as between G6-20, G6-16 well.
3. streamline simulation method
Streamline numerical simulation can intuitively describe the direction and region that fluid flows in three-dimensional space, pass through dredging for streamline The variation of close and streamline direction reflection advantage flow field direction and the power in flow field, streamline is more intensive, and flow field is stronger.
Streamline in streamline simulation indicates the throughflow of liquid, and dilute power to judge flow field can be dredged by streamline, flows Line is closeer, and fluid throughput is bigger, and flow field is stronger；Streamline is diluter, and liquid measure throughput is smaller, and flow field is weaker.
With the increase of displacement in flooding multiple, reservoir characteristics change high water-cut reservoir, and convenient value analogue technique is difficult Accurately to describe this variation.Targetedly method is taken in this research, improves the accuracy of numerical simulation, is the fine of flow field Description is laid a solid foundation.
(1) processing of high Injection volume phase percolation curve
Show that strong water drive can achieve the higher displacement of reservoir oil from the oil displacement efficiency calculated result of sealed coring well ZJ4 well sample product Efficiency, for sample segment oil displacement efficiency up to 73%, residual oil saturation only 19.9% is below Initial Experimental Stage result of study: remaining Oily saturation degree 24%~26% or so.After illustrating displacement multiple increase, residual oil saturation value is reduced.It is different from G6-104 well From the point of view of the phase percolation curve of displacement multiple, amount of water injected is increased to 50 times by conventional 20 times, and residual oil is dropped to from 18.7% 11.6%, the water phase permeability under residual oil rises to 0.42 (Fig. 8) from 0.29.
The positive rhythm theoretical model adopted by establishing the very poor note one for being 5, when well water 98% stops operation, utilization The water phase percolation ability at conventional phase percolation curve model oil reservoir top bottom differs 15 times, oil saturation difference 6.7%；And utilize high note Enter multiple phase percolation curve model oil reservoir top bottom water phase percolation ability and differs 100 times, oil saturation difference 25%.From 4 wells of true inspection, From the point of view of the sealing core drilling of 1 well of old inspection, at the top of positive rhythmic unit, oil saturation is higher, and part is in a state of nature, it is seen that high The numerical simulation of Injection volume phase percolation curve more can really reflect ultra-high water cut stage oil reservoir remaining oil distribution situation.
(2) time-varying of physical parameter
In waterflooding reservoir, since injection water constantly washes away, the clay mineral fragment in reservoir is made to be broken up, be moved It moves, and has part as Produced Liquid is taken out of, so as to cause reservoir permeability time-varying, influence oil-water movement and remaining oil distribution rule Rule.
It is that physical parameter approximation is modified by point water-cut stage in fact that the numerical simulation of traditional physical parameter time-varying, which describes method, It is existing.Since the region of High water cut is just existing near waterflooding development initial stage water injection well, and constantly expand, using dividing water-cut stage Simulation has its limitation.This research obtains hypotonic, middle height by establishing saturation degree, effective relationship for crossing water multiple, permeability It seeps area's difference and crosses the permeability variation multiple under water multiple, be truly reflected physical property time-varying.
Hypotonic area: K=0.9215 × e-0.004×PV (24)
Middle and high infiltration area: K=0.0231 × ln (PV) (25)
The advantages of this time-varying numerical simulation, is: establishing saturation degree, effective relationship for crossing water multiple, permeability, truly Dynamic physical property time-varying is realized, is greatly improved and perfect traditional time-varying numerical simulation technology of physical property stage by stage.
(3) processing of quick effect is answered
In the process of development, as strata pressure is gradually reduced, formation rock is compressed, and the small venturi in rock closes It closes, so as to cause the reduction of permeability, influences underground movement disciplinarian of oil and water (Fig. 9), it is therefore desirable to consider Ying Minxiao in digital-to-analogue It answers.Voltage-sensitive effect is considered by introducing related keyword (ROCKCOMP, OVERBURD, ROCKOPTS and ROCKTAB etc.) and realizing Reservoir numerical simulation.
By fine history matching, the whole district, individual well index have obtained preferable fitting (Figure 10, Figure 11), true reflection oil Geology, exploitation, the underground fluid seepage characteristics of hiding.The characterization in advantage flow field is carried out on the basis of fine analog.
Such as G6 fault block main force sand body E1f2 3-5Most of region flow field intensity is big, and stagnant area flow field is weaker between local well, by force Flow field is mainly distributed between G6-70 and CG6-14 well, between G6-12 and G6-111 well, between G6-64 and G6-2 well, G6-9 and G6- 25, regions (Fig. 7) between G6-16 well etc..Non- main force's sand body E1f1 1-2Whole flow field is weaker (Figure 12).
It all include following three wellblocks: the wellblock G6-2, G6-57 well using the advantage flow field area that above-mentioned three kinds of sides do not identify Area, the wellblock CG6-14.These three wellblocks will be adjusted, preferentially in lower successive step to reduce inefficient water injection, suggestions for improvement.
The differentiation in advantage flow field is carried out by three kinds of methods for the same substratum, it can be seen that three kinds of weak current field areas and The advantage flow field regions goodness of fit with higher.Advantage flow field can be preferably differentiated by multi-method synthesis.
Step 3: " the advantage flow field-initial flow-field method of difference " Tapping Residual Oil is analyzed
Research area's oil reservoir advantage flow field is realized using the method for dynamically associating, fuzzy comprehensive evaluation method and streamline simulation method to know Not, on the basis of the quantitative assessment of sandstone oil reservoir initial flow-field and advantage flow field identify, proposition is by initial flow-field intensity and at present The comparison of predominant current field intensity result of study, and then the method for finding potentiality block, to propose sandstone oil reservoir " advantage flow field- The initial flow-field method of difference " Technology of Tapping Residual Oil.
According to the initial flow-field value sought.In most cases initial flow-field is high, and remaining oil saturation is low, and initial flow-field is low It is worth area, remaining oil is relatively enriched with, and initial flow-field can preferably reflect the flowing of water filling and the distribution of remaining oil.But different small It often will appear initial flow-field height in layer, remaining oil saturation also high situation illustrates to study area's later period note in this case Occur the water filling of the area Liao Gai during water and does not involve or involve the lower situation of degree, thus the area belongs to the master to tap the latent power in next step Want target area.
Under the guidance of above-mentioned thinking, this, which is researched and proposed, utilizes the preferred method in the Study of recognition of advantage flow field, Using dynamically associating, the methods of method, fuzzy comprehensive evaluation method, streamline simulation method carry out the identification of advantage flow field and distribution is ground in work area Study carefully, by numerical simulation, the advantage Flow Field Distribution signature analysis in realization area combines initial flow-field intensity point on this basis Analysis is as a result, can be very good discovery REGION OF WATER INJECTION OILFIELD.Such as the wellblock G6-119 primary flow field intensity is higher, but at present flow field intensity compared with Low, G6-119 well metaideophone, weak current field is enhanced, and effect is good (Fig. 2, Fig. 7, Figure 12).
The working principle of the invention is on the basis of classification of oil reservoirs evaluation, by the flow field geology master for analyzing different oil reservoirs Control factor determines oil reservoir flow field quantitatively characterizing method, and then determines the initial flow-field regularity of distribution and mode, in conjunction with development behavior Equal latter period of oilfield development analysis advantage Flow Field Distribution rule and mode.Detailed process:
(1) binding area geologic feature, preferably reservoir heterogeneity, coefficient of permeability variation, crack, glue at permeability The evaluation index of native content, viscosity of crude, stratigraphic dip as initial flow-field.Using the fuzzy synthesis for solving the problems, such as Fuzzy Engineering Law theory is judged to establish the overall evaluation system of initial flow-field, and determines each evaluation criterion weight using analytic hierarchy process (AHP) Coefficient.
(2) high water-cut reservoir is divided into four classes, is respectively opposite integrity reservoir, narrow strip oil reservoir, lithological pool and extremely multiple Miscellaneous fault block oil reservoir.Utilize the method for dynamically associating, fuzzy comprehensive evaluation method, streamline simulation fado method comprehensive distinguishing development late stage advantage Flow field.
(3) it is directed to mid-later development phase high water-cut reservoir, using advantage flow field-initial flow-field method of difference, finds remaining oil, and Remaining oil is exploited.
The invention proposes the methods of the flow field of the untapped state of oil reservoir characterization；Propose that " advantage flow field-initial flow-field is poor Different method " Technology of Tapping Residual Oil.The present invention can needle it is untapped when oil reservoir, carry out the understanding of initial flow-field；The present invention proposes Sandstone oil reservoir " the advantage flow field-initial flow-field method of difference " Technology of Tapping Residual Oil scheme, improve development late stage it is superpower it is non- Remaining oil cognitive ability in matter reservoir；On the basis of classification of oil reservoirs, reservoir, construction feature are analysed in depth to initial The influence in flow field, by being reanalysed to development late stage oil reservoir flow field, comparative benefits flow field-initial flow-field, after being formed to exploitation The understanding of phase remaining oil distribution Favorable Areass, targeted Remaining Oil improve the development degree of oil reservoir reserves out of control With ultimate recovery, Production Decline Prediction of Oilfield amplitude is effectively reduced, extends the stable yields time limit in oil field, hence it is evident that it is final to improve oil reservoir Recovery ratio achieves obvious development effectiveness and benefit.
Embodiments of the present invention above described embodiment only expresses, the description thereof is more specific and detailed, but can not Therefore limitations on the scope of the patent of the present invention are interpreted as.It should be pointed out that for those of ordinary skill in the art, Without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to protection model of the invention It encloses.Therefore, the scope of protection of the patent of the invention shall be subject to the appended claims.

Claims (10)

1. a kind of high water-cut reservoir Tapping Residual Oil method characterized by comprising the evaluation of oil reservoir initial flow-field；Oil reservoir advantage Flow field evaluation.
2. high water-cut reservoir Tapping Residual Oil method according to claim 1, which is characterized in that the described method includes: logical The flow field geology Dominated Factors for analyzing different oil reservoirs are crossed, determine oil reservoir flow field quantitatively characterizing method, and then determine initial flow-field point Cloth rule and mode, in conjunction with the latter period of oilfield development analysis advantage such as development behavior Flow Field Distribution rule and mode, finally Using advantage flow field-initial flow-field method of difference, remaining oil is found, and remaining oil is exploited；It specifically includes:
(1) binding area geologic feature, preferably reservoir heterogeneity, permeability, coefficient of permeability variation, crack, clay contain The evaluation index of amount, viscosity of crude, stratigraphic dip as initial flow-field；Using the fuzzy comprehensive evoluation for solving the problems, such as Fuzzy Engineering Law theory establishes the overall evaluation system of initial flow-field, and each evaluation criterion weight system is determined using analytic hierarchy process (AHP) Number；
(2) high water-cut reservoir is divided into four classes, respectively with respect to integrity reservoir, narrow strip oil reservoir, lithological pool and extremely complicated disconnected Block oil reservoir；Utilize the method for dynamically associating, fuzzy comprehensive evaluation method, streamline simulation fado method comprehensive distinguishing development late stage predominant current ?；
(3) it is directed to mid-later development phase high water-cut reservoir, using advantage flow field-initial flow-field method of difference, finds remaining oil.
3. high water-cut reservoir Tapping Residual Oil method according to claim 1, which is characterized in that the evaluation of oil reservoir initial flow-field Include: to indicate ambiguity of the evaluation index to final evaluation goal parameter influence degree with fuzzy set form, is formed direct The evaluations matrix that can carry out judging operation, the evaluation result of a fuzzy set is obtained by blurring mapping.
4. high water-cut reservoir Tapping Residual Oil method according to claim 1, which is characterized in that the evaluation of oil reservoir initial flow-field It include: to choose the factor compositing factor collection for influencing effect of reservoir development first, then determine the Comment gathers for evaluating its result, then really Weight sets of the fixed each factor for evaluation result；Suitable membership function is chosen by tentative calculation, then calculates evaluation object Synthetic evaluation matrix；Blurring mapping is carried out to weight sets and synthetic evaluation matrix, obtains comprehensive evaluation result.
5. high water-cut reservoir Tapping Residual Oil method according to claim 4, which is characterized in that the evaluation of oil reservoir initial flow-field Include:
(1) set of object is determined；The set of the determination each layer of each well to be calculated in block to be evaluated, is denoted as:
A=(A1, A2, A3... An)；
(2) set of index factor, the i.e. set of each attribute of object are determined；Attribute set, which refers to, influences each of initial flow-field Kind heterogeneity parameter, is denoted as:
P=(P1, P2, P3... Pm)；
(3) determine that index matrix, matrix element are property parameters；
(4) index value, element a are normalizedijIt can be normalized as the following formula as bij
After providing each evaluation index, nondimensionalization need to be carried out to each parameter, to reject the influence between different dimensions；
Parameter is normalized；Matrix after normalization are as follows:
(5) determination of weight
The corresponding weight of each index is provided, weight vectors are obtained:
W=(w1, w2, w3, L, wm)；
(6) initial flow-field evaluation number I is calculated
I=BW；
Acquire initial flow-field evaluation number vector I=(i1, i2, i3, L in), wherein (0,1) i ∈, I is bigger, and flow field is stronger；I is got over Small, flow field is weaker.
6. high water-cut reservoir Tapping Residual Oil method according to claim 1, which is characterized in that commented in oil reservoir initial flow-field Valence comprises determining that evaluation index；Standardize initial data；Determine weight；Calculate initial flow-field intensity index；Determine initial flow-field Discriminant classification standard.
7. high water-cut reservoir Tapping Residual Oil method according to claim 1, which is characterized in that the evaluation of oil reservoir advantage flow field It include: to be evaluated using three kinds of method of dynamically associating, fuzzy comprehensive evaluation method, streamline simulation method methods research area.
8. high water-cut reservoir Tapping Residual Oil method according to claim 1, which is characterized in that the method for dynamically associating includes:
(1) analytical sequence is determined；
On the basis of qualitative analysis, a dependent variable factor and multiple independent variable factors are determined；If dependent variable data constitute ginseng Examine sequence X0', each argument data composition ratio is compared with sequence Xi' (i=1,2 ... m), and m+1 data sequence forms following matrix:
Wherein Xi '=(xi ' (1), xi ' (2) ..., xi ' (n))T, the m of i=0,1,2 ...；N is the length of Variables Sequence.
(2) nondimensionalization is carried out to Variables Sequence；
Respectively because prime sequences form following matrix after nondimensionalization:
Nondimensionalization is carried out using first value method:
(3) difference sequence, maximum difference and lowest difference are asked
Form absolute difference matrix:
Wherein: △0i(k)=| x0(k)-xi(k)|
I=0,1,2 ..., m；K=1,2 ..., n；
Maximum number and minimum number are maximum difference and lowest difference in absolute difference battle array:
(4) calculate correlation coefficient
Such as down conversion is made to data in absolute difference battle array:
Obtain incidence coefficient matrix:
Resolution ratio ρ is in (0,1) interior value in formula；
Incidence coefficient ξ0i(k) 1 positive number, Δ are no more than0i(k) smaller, ξ0i(k) bigger；
(5) calculating correlation
Averaging obtains XiWith X0The degree of association:
The degree of association is adopted to describe the power of relationship between water injection well and producing well by calculating note, judges whether generate predominant current between well ?.
9. high water-cut reservoir Tapping Residual Oil method according to claim 1, which is characterized in that fuzzy comprehensive evaluation method packet It includes:
Permeability and clay content are chosen as static evaluation index；
For water injection well, daily water-injection rate, water filling three indexs of oil pressure and unit thickness cumulative water injected are chosen as dynamic evaluation Index；For producing well, daily fluid production rate, moisture content, unit thickness cumulative liquid production are chosen as dynamic assessment index；
Nondimensionalization processing is carried out to parameter；
Using Weight of Coefficient through Analytic Hierarchy Process:
Firstly, the weight of permeability and shale content is determined to the relative importance for forming advantage flow field according to each influence factor, Constitute influence factor weight vectors [ω1,ω2]T, guarantee that the sum of weight is equal to 1；
Secondly, according to oil well, water injection well dynamic performance characteristic to the relative importance in advantage flow field, determine the power of each index Weight；
Finally, determining the weight of the two according to the importance that the influence factor in advantage flow field and performance characteristic identify existence [ωa,ωb], guarantee that the sum of weight is equal to 1.
10. high water-cut reservoir Tapping Residual Oil method according to claim 6, which is characterized in that initial flow-field classification is sentenced Other standard are as follows:
If composite index I < 0.4, for weak current field；If composite index 0.4≤I < 0.6, for middle flow field；If composite index 0.6 ≤ I < 0.8 is then high current field；If composite index 0.8≤I < 1, for advantage flow field.
CN201910589725.1A 2019-07-02 2019-07-02 High-water-content oil reservoir residual oil submergence excavating method Active CN110288258B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910589725.1A CN110288258B (en) 2019-07-02 2019-07-02 High-water-content oil reservoir residual oil submergence excavating method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910589725.1A CN110288258B (en) 2019-07-02 2019-07-02 High-water-content oil reservoir residual oil submergence excavating method

Publications (2)

Publication Number Publication Date
CN110288258A true CN110288258A (en) 2019-09-27
CN110288258B CN110288258B (en) 2021-05-18

Family

ID=68021785

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910589725.1A Active CN110288258B (en) 2019-07-02 2019-07-02 High-water-content oil reservoir residual oil submergence excavating method

Country Status (1)

Country Link
CN (1) CN110288258B (en)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110778314A (en) * 2019-10-08 2020-02-11 中国石油化工股份有限公司 Reasonable mechanical recovery system efficiency measuring and calculating method based on oil reservoir conditions
CN110851982A (en) * 2019-11-12 2020-02-28 中国海洋石油集团有限公司 Method for analyzing feasibility of oxygen reduction air flooding of medium-low permeability oil reservoir
CN111101924A (en) * 2019-11-15 2020-05-05 中国石油天然气股份有限公司大港油田分公司 Lithologic reservoir dominant facies band prediction method and device
CN111119875A (en) * 2019-12-24 2020-05-08 中海石油(中国)有限公司 Water flooding front edge monitoring method based on crude oil alkylphenol compound content change
CN111236901A (en) * 2020-01-09 2020-06-05 西南石油大学 Water drive reservoir seepage field heterogeneity evaluation method based on flow field diagnosis
CN111400972A (en) * 2020-03-24 2020-07-10 西南石油大学 Semi-closed fault block oil reservoir productivity analysis method
CN112001055A (en) * 2019-11-07 2020-11-27 中海石油(中国)有限公司 Low-amplitude thin oil reservoir water content prediction method based on microstructure
CN112561279A (en) * 2020-12-09 2021-03-26 中国石油化工股份有限公司胜利油田分公司孤岛采油厂 Identification method and system for high-water-consumption zone
CN112576238A (en) * 2020-12-02 2021-03-30 中国石油大学(华东) System, method and application for determining position and content of residual oil in low-permeability reservoir
CN113029892A (en) * 2020-03-17 2021-06-25 中国海洋石油集团有限公司 Method for evaluating reasonability of oil-water relative permeability curve based on regional statistical law
CN113434812A (en) * 2021-05-20 2021-09-24 青岛海洋地质研究所 Hydrate trial production target optimization method based on fuzzy comprehensive evaluation
CN113761697A (en) * 2020-06-01 2021-12-07 中国石油化工股份有限公司 Method for rapidly evaluating water-drive reservoir transfer field development effect
CN113818849A (en) * 2020-06-18 2021-12-21 中国石油化工股份有限公司 Method for improving oil extraction speed of oil-water reverse migration development with low cost
CN114066666A (en) * 2021-11-17 2022-02-18 西南石油大学 Method for analyzing connectivity among wells through injection-production profile monitoring data
CN114810006A (en) * 2021-01-27 2022-07-29 中国石油化工股份有限公司 Potential evaluation method for regulating and controlling high-water-consumption zone by separate-layer water injection after heterogeneous flooding
CN115018225A (en) * 2021-03-04 2022-09-06 中国石油化工股份有限公司 Oil reservoir bottom water volume multiple evaluation method and device, computer equipment and medium
CN118194204A (en) * 2024-05-14 2024-06-14 合肥城市云数据中心股份有限公司 Industrial data feature selection and outlier detection method, system, equipment and medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102346806A (en) * 2010-08-06 2012-02-08 中国石油天然气股份有限公司 Selection method and device for high-water-cut oil field layer system recombination scheme
CN103628850A (en) * 2012-08-23 2014-03-12 中国石油化工股份有限公司 Overall profile control and water plugging decision-making method for water-flooding oil field
CN104778378A (en) * 2015-05-05 2015-07-15 中国石油大学(华东) Method for analyzing influence factors of production decline of oil and gas field
RU2631451C1 (en) * 2016-07-29 2017-09-22 федеральное государственное автономное образовательное учреждение высшего образования "Национальный исследовательский ядерный университет МИФИ" (НИЯУ МИФИ) Method to increase oil recovery of formation with high viscosity oil
CN107844653A (en) * 2017-11-10 2018-03-27 中国石油天然气股份有限公司 Comprehensive evaluation method and device for reservoir water flooding development potential
CN108733632A (en) * 2018-05-17 2018-11-02 成都北方石油勘探开发技术有限公司 Hypotonic high water-cut reservoir refracturing selects well evaluation method in one kind

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102346806A (en) * 2010-08-06 2012-02-08 中国石油天然气股份有限公司 Selection method and device for high-water-cut oil field layer system recombination scheme
CN103628850A (en) * 2012-08-23 2014-03-12 中国石油化工股份有限公司 Overall profile control and water plugging decision-making method for water-flooding oil field
CN104778378A (en) * 2015-05-05 2015-07-15 中国石油大学(华东) Method for analyzing influence factors of production decline of oil and gas field
RU2631451C1 (en) * 2016-07-29 2017-09-22 федеральное государственное автономное образовательное учреждение высшего образования "Национальный исследовательский ядерный университет МИФИ" (НИЯУ МИФИ) Method to increase oil recovery of formation with high viscosity oil
CN107844653A (en) * 2017-11-10 2018-03-27 中国石油天然气股份有限公司 Comprehensive evaluation method and device for reservoir water flooding development potential
CN108733632A (en) * 2018-05-17 2018-11-02 成都北方石油勘探开发技术有限公司 Hypotonic high water-cut reservoir refracturing selects well evaluation method in one kind

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110778314A (en) * 2019-10-08 2020-02-11 中国石油化工股份有限公司 Reasonable mechanical recovery system efficiency measuring and calculating method based on oil reservoir conditions
CN112001055A (en) * 2019-11-07 2020-11-27 中海石油(中国)有限公司 Low-amplitude thin oil reservoir water content prediction method based on microstructure
CN112001055B (en) * 2019-11-07 2024-04-09 中海石油(中国)有限公司 Low-amplitude thin oil reservoir water content prediction method based on microstructure
CN110851982A (en) * 2019-11-12 2020-02-28 中国海洋石油集团有限公司 Method for analyzing feasibility of oxygen reduction air flooding of medium-low permeability oil reservoir
CN111101924A (en) * 2019-11-15 2020-05-05 中国石油天然气股份有限公司大港油田分公司 Lithologic reservoir dominant facies band prediction method and device
CN111119875A (en) * 2019-12-24 2020-05-08 中海石油(中国)有限公司 Water flooding front edge monitoring method based on crude oil alkylphenol compound content change
CN111119875B (en) * 2019-12-24 2023-04-07 中海石油(中国)有限公司 Water flooding front edge monitoring method based on crude oil alkylphenol compound content change
CN111236901A (en) * 2020-01-09 2020-06-05 西南石油大学 Water drive reservoir seepage field heterogeneity evaluation method based on flow field diagnosis
CN113029892A (en) * 2020-03-17 2021-06-25 中国海洋石油集团有限公司 Method for evaluating reasonability of oil-water relative permeability curve based on regional statistical law
CN113029892B (en) * 2020-03-17 2022-12-13 中国海洋石油集团有限公司 Method for evaluating reasonability of oil-water relative permeability curve based on regional statistical rule
CN111400972A (en) * 2020-03-24 2020-07-10 西南石油大学 Semi-closed fault block oil reservoir productivity analysis method
CN111400972B (en) * 2020-03-24 2022-02-15 西南石油大学 Semi-closed fault block oil reservoir productivity analysis method
CN113761697A (en) * 2020-06-01 2021-12-07 中国石油化工股份有限公司 Method for rapidly evaluating water-drive reservoir transfer field development effect
CN113818849A (en) * 2020-06-18 2021-12-21 中国石油化工股份有限公司 Method for improving oil extraction speed of oil-water reverse migration development with low cost
CN113818849B (en) * 2020-06-18 2023-08-01 中国石油化工股份有限公司 Method for improving oil extraction speed by reverse migration of oil and water at low cost
CN112576238A (en) * 2020-12-02 2021-03-30 中国石油大学(华东) System, method and application for determining position and content of residual oil in low-permeability reservoir
CN112561279B (en) * 2020-12-09 2023-02-24 中国石油化工股份有限公司胜利油田分公司孤岛采油厂 Identification method and system for high-water-consumption zone
CN112561279A (en) * 2020-12-09 2021-03-26 中国石油化工股份有限公司胜利油田分公司孤岛采油厂 Identification method and system for high-water-consumption zone
CN114810006A (en) * 2021-01-27 2022-07-29 中国石油化工股份有限公司 Potential evaluation method for regulating and controlling high-water-consumption zone by separate-layer water injection after heterogeneous flooding
CN114810006B (en) * 2021-01-27 2024-05-31 中国石油化工股份有限公司 Heterogeneous flooding post-treatment separate layer water injection regulation and control high water consumption layer potential evaluation method
CN115018225A (en) * 2021-03-04 2022-09-06 中国石油化工股份有限公司 Oil reservoir bottom water volume multiple evaluation method and device, computer equipment and medium
CN113434812A (en) * 2021-05-20 2021-09-24 青岛海洋地质研究所 Hydrate trial production target optimization method based on fuzzy comprehensive evaluation
CN114066666A (en) * 2021-11-17 2022-02-18 西南石油大学 Method for analyzing connectivity among wells through injection-production profile monitoring data
CN118194204A (en) * 2024-05-14 2024-06-14 合肥城市云数据中心股份有限公司 Industrial data feature selection and outlier detection method, system, equipment and medium

Also Published As

Publication number Publication date
CN110288258B (en) 2021-05-18

Similar Documents

Publication Publication Date Title
CN110288258A (en) A kind of high water-cut reservoir Tapping Residual Oil method
CN112343587B (en) Identification characterization method for dominant seepage channel of ultralow permeability oil reservoir
CN112561144B (en) Method for evaluating main control factor of productivity of tight oil fracturing horizontal well and predicting productivity
CN106837297B (en) Method for identifying connectivity among wells and predicting oil-water dynamic state
CN110276827B (en) Effectiveness evaluation method based on shale reservoir
CN102041995B (en) System for monitoring complicated oil deposit flooding conditions
Wang et al. Quantitative prediction of oil and gas prospects of the Sinian-Lower Paleozoic in the Sichuan Basin in central China
CN103352693B (en) A kind of output prediction system based on oily content and method thereof
CN110717301B (en) Flow unit information classification and identification method based on support vector machine algorithm
CN109815516A (en) Method and device for predicting productivity of shale gas well
CN110608023B (en) Adaptability boundary analysis and evaluation method for stratified steam injection of thickened oil
CN105740563A (en) Preferential channel identification technology for secondary development of mature oil field
CN113743023B (en) Hierarchical characterization method for carbon dioxide flooding gas channeling channel
CN110162922A (en) A kind of integrated recognition method of water-drive pool dominant flowing path
CN108343430A (en) Fault sealing property fuzzy mathematics comprehensive quantitative evaluation method
CN109403962A (en) Oil reservoir block Monitoring Indexes association analysis method
CN106651060A (en) Method for quantitatively evaluating degree of enrichment of contactor-relay control natural gas reservoirs
CN111444621A (en) High-water-content oil reservoir flow unit dividing method based on dynamic and static seepage interface
CN108561129A (en) Craven fault block reservoirs water to oil area remaining oil fast appraisement method
Zhang et al. Experimental study on natural gas migration and accumulation mechanism in sweet spots of tight sandstones
CN114139242A (en) Water flooded layer well logging evaluation method based on lithofacies
CN111287739B (en) Residual oil distribution prediction method based on stratum crude oil viscosity
CN113971528A (en) Method for identifying dominant seepage channel
CN108490484B (en) Method for dissection of blank space of mature exploration area of fractured basin
CN114076992A (en) Classification evaluation method for micro pore throats of meandering stream reservoir based on lithofacies

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
GR01 Patent grant
GR01 Patent grant