CN106703776A - Method for optimizing fracturing parameters - Google Patents
Method for optimizing fracturing parameters Download PDFInfo
- Publication number
- CN106703776A CN106703776A CN201611202818.7A CN201611202818A CN106703776A CN 106703776 A CN106703776 A CN 106703776A CN 201611202818 A CN201611202818 A CN 201611202818A CN 106703776 A CN106703776 A CN 106703776A
- Authority
- CN
- China
- Prior art keywords
- parameter
- fracturing
- effect
- fractured well
- reservoir
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000001965 increasing effect Effects 0.000 claims abstract description 75
- 238000011156 evaluation Methods 0.000 claims abstract description 48
- 238000010276 construction Methods 0.000 claims abstract description 11
- 238000013461 design Methods 0.000 claims abstract description 7
- 230000000694 effects Effects 0.000 claims description 79
- 238000004519 manufacturing process Methods 0.000 claims description 72
- 239000011159 matrix material Substances 0.000 claims description 50
- 238000012360 testing method Methods 0.000 claims description 24
- 239000004576 sand Substances 0.000 claims description 20
- 238000002474 experimental method Methods 0.000 claims description 17
- 239000007788 liquid Substances 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000003786 synthesis reaction Methods 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 6
- 241000208340 Araliaceae Species 0.000 claims description 4
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 4
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000010219 correlation analysis Methods 0.000 claims description 4
- 235000008434 ginseng Nutrition 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 3
- 238000012417 linear regression Methods 0.000 claims description 3
- 229940050561 matrix product Drugs 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 3
- 238000005457 optimization Methods 0.000 abstract 1
- 238000003860 storage Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000001537 neural effect Effects 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010141 design making Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- E21B43/26—Methods for stimulating production by forming crevices or fractures
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention discloses a method for optimizing fracturing parameters, and belongs to the field of oil and gas reservoir exploration and development. The method includes steps of collecting a plurality of groups of reservoir parameters, fracturing parameters and yield increasing effect data of fractured wells and reservoir parameters of to-be-fractured wells; analyzing correlation between the collected reservoir parameters and the collected yield increasing effect data of the fractured wells and correlation between the collected fracturing parameters and the collected yield increasing effect data of the fractured wells so as to obtain mathematical statistics results; applying orthogonal experimental scheme designs and combining the orthogonal experimental scheme designs with the fracturing parameters of the fractured wells and the reservoir parameters of the to-be-fractured wells to form a plurality of groups of fracturing parameter alternative schemes for the to-be-fractured wells; converting fuzzy comprehensive evaluation matrixes of the alternative schemes into fuzzy comprehensive evaluation scores, sorting the fuzzy comprehensive evaluation scores from high to low according to maximum membership degree principles, sorting construction risks from low to high and determining ultimate optimal schemes. The method has the advantage that the reasonableness and the effectiveness of fracturing parameter optimization and yield increasing effects can be improved.
Description
Technical field
The present invention relates to reservoir exploration and development field, more particularly to a kind of fracturing parameter method for optimizing.
Background technology
Hydraulic fracturing is hyposmosis, the key technology of fine and close reservoir exploration and development.Effect of increasing production after well fracturing
In addition to being influenceed by reservoir parameter (effective thickness, porosity, gas saturation, compensated neutron, compensation density), pressure break ginseng
Number (operational discharge capacity, average sand ratio, preposition liquid proportional, plus sand intensity) be preferably also determine effect of increasing production quality it is crucial because
Element.
Conventional fracturing parameter method for optimizing is using a large amount of reservoir parameters of fractured well of oil field block, pressure at present
Parameter and effect of increasing production are split, set up reservoir parameter, fracturing parameter by regression analysis or neural net method imitates with volume increase
Mathematic(al) representation between fruit.Then for wanting the well of preferred fracturing parameter, because the reservoir parameter of the well is it has been determined that logical
Cross and change different fracturing parameter combinations, using the mathematic(al) representation set up, the fractured well calculated under different parameters combination increases
Effect is produced, and selects the best parameter combination of effect of increasing production as preferred fracturing parameter.
Realize it is of the invention during, inventor find prior art at least there is problems with:
The numerical value of each reservoir parameter is taken from the deterministic numerical value such as the limiting value of each parameter, average value, does not examine fully
Each reservoir parameter is considered due to uncertainty or by laboratory apparatus, experimental technique, calculating analysis equal error, and causes these parameters
Actual value of the parameter under reservoir conditions is not necessarily reflected, have ignored these parameters has the characteristic of certain fuzziness.Its
It is secondary, using the reservoir parameter of fractured well, fracturing parameter and effect of increasing production, set up by regression analysis or neural net method
During the method for mathematical relationship, when reservoir parameter, fracturing parameter and effect of increasing production do not exist certain mathematical relationship, can cause excellent
The fracturing parameter deviation increase of choosing.Meanwhile, alternative fracturing parameter preferred scheme very little, finally preferably out actually
Local optimum fracturing parameter may be just corresponded to ask rather than global optimum's fracturing parameter truly.
The content of the invention
In order to solve problem of the prior art, a kind of fracturing parameter method for optimizing is the embodiment of the invention provides.The skill
Art scheme is as follows:
A kind of fracturing parameter method for optimizing, comprises the following steps:
S100 collects the reservoir parameter of multigroup fractured well, fracturing parameter, effect of increasing production data and treats the storage of fractured well
Layer parameter;
S200 will be collected into described in the reservoir parameter of fractured well and effect of increasing production data and fracturing parameter and volume increase
Effect data carries out correlation analysis, draws mathematical statistics result;
S300 is based on the mathematical statistics result, sets up reservoir parameter and effect of increasing production data and fracturing parameter and volume increase
The membership function relation and single factor test fuzzy evaluating matrix of effect data;
S400 be based on the mathematical statistics result, using grey correlation theory set up reservoir parameter account for effect of increasing production data with
And fracturing parameter accounts for the weight matrix of effect of increasing production data;
S500 application orthogonal experiment conceptual designs, in combination with fracturing parameter and the storage for treating fractured well of the fractured well
Layer parameter, forms multigroup fracturing parameter alternative for treating fractured well;
S600 is based on the fracturing parameter alternative for treating fractured well, and the membership function set up using fractured well is closed
System and the single factor test fuzzy evaluating matrix, the single factor test fuzzy evaluating matrix and the weight matrix are combined, and are obtained
To the fuzzy overall evaluation matrix of the alternative;
The fuzzy overall evaluation matrix of the alternative is converted into fuzzy overall evaluation score by S700, is subordinate to according to maximum
Category degree principle, final preferred scheme is determined by fuzzy synthesis scoring from low to high from high to low and with reference to construction risk.
Alternatively, in the step S100 reservoir parameter, fracturing parameter, effect of increasing production data include that reservoir is effectively thick
Degree, porosity, gas saturation, compensated neutron, compensation density, operational discharge capacity, preposition liquid proportional, average sand ratio, plus sand intensity
And effect of increasing production data.
Alternatively, the step S200 is specially:
The reservoir parameter of fractured well and effect of increasing production data and fracturing parameter and effect of increasing production described in being collected into
Data carry out single factor test linear regression analysis, it is determined that the reservoir parameter of the fractured well, fracturing parameter and effect of increasing production are category
In positive correlation or negative correlativing relation, and line parameter minimum value is entered to each parameter according to positive correlation and negative correlativing relation
The opinion rating of 4 class deciles and maximum between.
Alternatively, the step S300 is specially:
The mathematical statistics of the reservoir parameter according to the fractured well, the single factor test numerical value between fracturing parameter and effect of increasing production
As a result, the membership function of parameters is set up, approximate normal distribution membership function is selected,
In formula:uv- membership function, zero dimension;The occurrence of d-critical parameter, zero dimension;V-opinion rating, altogether
" I ", " II ", " III ", " IV " 4 class opinion rating, zero dimension;The index parameter of a, b-each opinion rating, zero dimension;
The membership function for determining same opinion rating is,
In formula:uvj(di)—diIn opinion rating vjOn degree of membership, zero dimension;di- i-th evaluating it is specific
Value, zero dimension;
From formula (2):Work as diDuring=a, uvjA ()=1, represents that the evaluating belongs to being subordinate to for some opinion rating
Degree is maximum, so a is the mathematic expectaion of certain grade, i.e.,
A=(d1+d2)/2 (3)
The equal point of adjacent opinion rating degree of membership is referred to as transition point, and the degree of membership for taking transition point is 0.5, i.e.
In formula:d1、d2The interval upper limit value and lower limit value of-j-th opinion rating, the d of each opinion rating of each evaluating1、d2
Value is determined by mathematical statistics result, after a, b value determine, you can it is determined that known membership function relation.
The single factor test fuzzy evaluating matrix corresponding to each parameter is can be calculated by formula (1)-(5), if total m is commented
Valency grade, n evaluating, then can obtain a single factor test fuzzy evaluating matrix R for n × m ranks,
Alternatively, the step S500 is specially:
According to fractured well fracturing parameter, fractured well fracturing parameter scope is counted, minimum value in each fracturing parameter,
Determine three 3 horizontal factors between maximum, form the operational discharge capacity of fractured well fracturing parameter, preposition liquid proportional, average
Sand than, plus sand intensity scheme, this scheme is one and has four factors, three experimental programs of 3 horizontal factors, sets up preferred
The preferred form of orthogonal experiment;
What cooperation was collected into treats fractured well reservoir effective thickness, porosity, gas saturation, compensated neutron, compensation density
Parameter, with reference to orthogonal experiment form, the alternative of fractured well fracturing parameter is treated in formation.
Alternatively, the step S600 is specially:
According to effective thickness, porosity, gas saturation, compensated neutron, compensation density parameter and the pressure break for the treatment of fractured well
The preferred form combination of the preferred orthogonal experiment of parameter, applying equation (6) determines the single factor evaluation matrix R of each parameter;
According to the single factor evaluation matrix R, using matrix product summation algorithm, by fuzzy weight vector A and single factor test mould
Paste Evaluations matrix R is combined, and considers all parameter contributions and draws fuzzy overall evaluation matrix B,
Alternatively, the step S700 is specially:
According to maximum membership grade principle, the fuzzy overall evaluation matrix of each alternative is converted into fuzzy overall evaluation and is obtained
Point, fraction is higher, and effect of increasing production is better after representing alternative implementation;When the fuzzy synthesis scoring of different alternatives is identical
When, consider pressure break risk and operating cost determines final preferred scheme;
Computing formula is as follows,
Wherein:vI、vII、vIII、vIVThe corresponding fraction of opinion rating is respectively 100,75,50,25.
The beneficial effect that technical scheme provided in an embodiment of the present invention is brought is:
Fracturing parameter method for optimizing has taken into full account the parameter of many factors of influence effect of increasing production, acquisition with fuzzy
Property, the problem of effect of increasing production and many factors without clear and definite functional relation expression formula, in the prior art may so as to overcome
Fracturing parameter local optimum is only realized rather than all optimal problems, the preferred reasonability of fracturing parameter, effectively is improve
Property and effect of increasing production.
Brief description of the drawings
Technical scheme in order to illustrate more clearly the embodiments of the present invention, below will be to that will make needed for embodiment description
Accompanying drawing is briefly described, it should be apparent that, drawings in the following description are only some embodiments of the present invention, for
For those of ordinary skill in the art, on the premise of not paying creative work, other can also be obtained according to these accompanying drawings
Accompanying drawing.
Fig. 1 is fracturing parameter method for optimizing flow chart provided in an embodiment of the present invention;
Fig. 2 a are the graphs of a relation of fractured well effect of increasing production provided in an embodiment of the present invention and effective thickness;
Fig. 2 b are the graphs of a relation of fractured well effect of increasing production provided in an embodiment of the present invention and porosity;
Fig. 2 c are the graphs of a relation of fractured well effect of increasing production provided in an embodiment of the present invention and gas saturation;
Fig. 2 d are the graphs of a relation of fractured well effect of increasing production provided in an embodiment of the present invention and compensated neutron;
Fig. 2 e are the graphs of a relation of fractured well effect of increasing production provided in an embodiment of the present invention and compensation density;
Fig. 2 f are the graphs of a relation of fractured well effect of increasing production provided in an embodiment of the present invention and operational discharge capacity;
Fig. 2 g be fractured well effect of increasing production provided in an embodiment of the present invention with average sand than graph of a relation;
Fig. 2 h are the graphs of a relation of fractured well effect of increasing production provided in an embodiment of the present invention and preposition liquid proportional.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with accompanying drawing to embodiment party of the present invention
Formula is described in further detail.
A kind of fracturing parameter method for optimizing is the embodiment of the invention provides, referring to Fig. 1, is comprised the following steps:
S100 collects the reservoir parameter of multigroup fractured well, fracturing parameter, effect of increasing production data and treats the storage of fractured well
Layer parameter.
S200 will be collected into described in the reservoir parameter of fractured well and effect of increasing production data and fracturing parameter and volume increase
Effect data carries out correlation analysis, draws mathematical statistics result.
S300 is based on the mathematical statistics result, sets up reservoir parameter and effect of increasing production data and fracturing parameter and volume increase
The membership function relation and single factor test fuzzy evaluating matrix of effect data.
S400 be based on the mathematical statistics result, using grey correlation theory set up reservoir parameter account for effect of increasing production data with
And fracturing parameter accounts for the weight matrix of effect of increasing production data.
S500 application orthogonal experiment conceptual designs, in combination with fracturing parameter and the storage for treating fractured well of the fractured well
Layer parameter, forms multigroup fracturing parameter alternative for treating fractured well.
S600 is based on the fracturing parameter alternative for treating fractured well, and the membership function set up using fractured well is closed
System and the single factor test fuzzy evaluating matrix, the single factor test fuzzy evaluating matrix and the weight matrix are combined, and are obtained
To the fuzzy overall evaluation matrix of the alternative.
The fuzzy overall evaluation matrix of the alternative is converted into fuzzy overall evaluation score by S700, is subordinate to according to maximum
Category degree principle, final preferred scheme is determined by fuzzy synthesis scoring from low to high from high to low and with reference to construction risk.
In the present embodiment, in step S100 reservoir parameter, fracturing parameter, effect of increasing production data specifically include reservoir has
Effect thickness, porosity, gas saturation, compensated neutron, compensation density, operational discharge capacity, preposition liquid proportional, average sand ratio, plus sand
Intensity and effect of increasing production data.
In the present embodiment, step S200 is specially:
The reservoir parameter of fractured well, fracturing parameter and effect of increasing production data are carried out into single factor test linear regression analysis, really
The numerical value for determining effect of increasing production is to increase with the increase of the reservoir parameter, fracturing parameter of fractured well or increase and reduce.
Further determine that the reservoir parameter of fractured well, fracturing parameter and effect of increasing production are belonging to positive correlation and still bear
Dependency relation, and according to positive correlation and negative correlativing relation each parameter is entered between line parameter minimum value and maximum " I ",
The opinion rating of " II ", " III ", " IV " totally 4 class deciles, parameter " I " level opinion rating of wherein positive correlation represents parameter
Greatly, on the contrary the parameter of negative correlativing relation " I " level evaluation represents that parameter is small.
In the present embodiment, step S300 is specially:
The mathematical statistics knot of the reservoir parameter according to fractured well, the single factor test numerical value between fracturing parameter and effect of increasing production
Really, the membership function of parameters is set up, approximate normal distribution membership function is selected,
In formula:uv- membership function, zero dimension;The occurrence of d-critical parameter, zero dimension;V-opinion rating, altogether
" I ", " II ", " III ", " IV " 4 class opinion rating, zero dimension;The index parameter of a, b-each opinion rating, zero dimension.
The membership function for determining same opinion rating is,
In formula:uvj(di)—diIn opinion rating vjOn degree of membership, zero dimension;di- i-th evaluating it is specific
Value, zero dimension.
From formula (2):Work as diDuring=a, uvjA ()=1, represents that the evaluating belongs to being subordinate to for some opinion rating
Degree is maximum, so a is the mathematic expectaion of certain grade, i.e.,
A=(d1+d2)/2 (3)
The equal point of adjacent opinion rating degree of membership is referred to as transition point, and it is 0.5 to take the degree of membership of transition point, i.e.,:
In formula:d1、d2The interval upper limit value and lower limit value of-j-th opinion rating, the d of each opinion rating of each evaluating1、d2
Value is determined by mathematical statistics result, after a, b value determine, you can it is determined that known membership function relation.
The single factor test fuzzy evaluating matrix corresponding to each parameter is can be calculated by formula (1)-(5), if total m is commented
Valency grade, n evaluating, then can obtain a single factor test fuzzy evaluating matrix R for n × m ranks,
In the present embodiment, step S500 is specially:
According to fractured well fracturing parameter, fractured well fracturing parameter scope is counted, minimum value in each fracturing parameter,
Determine three 3 horizontal factors between maximum, form the operational discharge capacity of fractured well fracturing parameter, preposition liquid proportional, average
Sand ratio, plus sand intensity, this scheme is one and has four factors, three experimental programs of 3 horizontal factors, sets up preferably orthogonal reality
Preferred form is tested, as shown in table 1, the data in each space are obtained according to the pressure break data statistics homogenous segmentations of fractured well.
Table 1 treats the fractured well fracturing parameter preferred form of preferred orthogonal experiment
What cooperation was collected into treats fractured well reservoir effective thickness, porosity, gas saturation, compensated neutron, compensation density
Parameter, with reference to the preferred form of above-mentioned orthogonal experiment, the alternative of fractured well fracturing parameter is treated in formation.
In the present embodiment, step S600 is specially:
According to effective thickness, porosity, gas saturation, compensated neutron, compensation density parameter and the pressure break for the treatment of fractured well
The preferred form combination of the preferred orthogonal experiment of parameter, applying equation (6) determines the single factor evaluation matrix R of each parameter.
According to single factor evaluation matrix R, using matrix product summation algorithm, fuzzy weight vector A is obscured with single factor test and is commented
Valency matrix R is combined, and considers all parameter contributions and draws fuzzy overall evaluation matrix B,
In the present embodiment, step S700 is specially:
According to maximum membership grade principle, the fuzzy overall evaluation matrix of each alternative is converted into fuzzy overall evaluation and is obtained
Point, fraction is higher, and effect of increasing production is better after representing alternative implementation;When the fuzzy synthesis scoring of different alternatives is identical
When, consider pressure break risk and operating cost determines final preferred scheme.
Computing formula is as follows
Wherein:vI、vII、vIII、vIVThe corresponding fraction of opinion rating is respectively 100,75,50,25.
Using above method step, the present embodiment implementation process is as follows:
For a specific oil-gas reservoir, it is considered to which fractured well is in fracturing technology, equipment, fracturing fluid and the constant feelings of proppant
Under condition, entered to subsequently treating the fracturing parameter of fractured well using the reservoir basic parameter of fractured well, fracturing parameter and effect of increasing production
Row is preferred.In the factor of numerous influence effect of increasing production, effect of increasing production is sensitive to some Parameters variations, and unwise to some parameters
Sense.By fractured well data, using correlation analysis and the live expertise of combination, final choice determines 15 mouthfuls of wells to be had
Effect thickness, porosity, gas saturation, compensated neutron, compensation density, operational discharge capacity, preposition liquid proportional, average sand ratio, plus sand
Intensity and effect of increasing production data are shown in Table 2 as Research foundation.Wherein numbering 1-15 is that fractured well basic parameter is imitated with volume increase
Really, numbering 16 be need preferred construction parameter treat fractured well basic parameter.
Table 2 fractured well reservoir parameter, construction parameter and effect of increasing production statistics
Referring to Fig. 2, fractured well effect of increasing production and effective thickness, porosity, gas saturation, compensated neutron, plus sand intensity
Into positive correlation, in other words as the numerical value of these parameters increases, fractured well effect of increasing production improves;And with compensation density,
Into negative correlativing relation, i.e., as the numerical value of these parameters increases, fractured well increases production for operational discharge capacity, average sand ratio, preposition liquid proportional
Effect is deteriorated.
According to each parameter and the dependency relation statistics of effect of increasing production, each parameter is divided into 4 opinion ratings, such as table 3
It is shown.
The by stages opinion rating of each evaluating of table 3
Using fractured well reservoir parameter, fracturing parameter and effect of increasing production data, parameter is determined using grey correlation theory
To the weight shared by effect of increasing production, the weight matrix A shared by each parameter is asked for, as shown in table 4.
Each factor of evaluation of table 4 is to the weight shared by effect of increasing production
A=[0.1329 0.1284 0.1274 0.1248 0.0906 0.0879 0.0881 0.0850 0.1344] (9)
Using orthogonal experiment conceptual design, with reference to the fracturing parameter scope of fractured well, selection orthogonal experiment form, design
Orthogonal experiment parameter and scheme;It is simultaneously comprehensive to treat that the reservoir parameter and orthogonal experiment schematic design making of fractured well treat the pressure break of fractured well
Parameter alternative, as shown in table 5.
Table 5 fractured well reservoir parameter, construction parameter and effect of increasing production statistics
According to each parameter of fracturing parameter alternative for treating fractured well, using each parameter membership function that fractured well is set up,
The single factor evaluation function of each parameter of fractured well alternative is treated in evaluation, further by each single factor evaluation function and weight matrix
It is combined, is obtained the fuzzy overall evaluation matrix of each alternative, by taking each basic parameter in Protocol Numbers 1 as an example, according to
Formula (1)-(5), with reference to table 3, the data of table 4, the single factor evaluation matrix of numerical procedure 1 is,
Above-mentioned single factor evaluation matrix only reflects each factor in each interval degree of membership, it is impossible to embody each parameter
To the combined influence result of effect of increasing production.By taking effective thickness as an example, effective thickness is subordinate to " I " class evaluation knot in reflecting scheme 1
The degree of membership of fruit is 0.5000, and it is 0 to be subordinate to " II ", " III ", the degree of membership of " IV " class evaluation result, according to maximum membership degree
Principle, effective thickness is under the jurisdiction of " I " class evaluation result.
Weight and single factor evaluation matrix are combined, Multifactor Comprehensive Evaluation matrix is obtained.
According to maximum membership grade principle, the comprehensive evaluation result of arrangement and method for construction I belongs to " I " class evaluation result, according to formula (8)
Fuzzy evaluating matrix is converted into 83 points of evaluation score.
Repeat the above steps, the comprehensive evaluation result of scheme 2-9 in fractured well alternative is treated in completion, and calculates fuzzy commenting
Point, as shown in table 6.
The comprehensive evaluation result of the different schemes of table 6
According to 9 kinds of comprehensive grading results of scheme in table 6, the fuzzy overall evaluation fraction of scheme 3,4,5,8 is all 84 points,
Identical best result can be all obtained after showing this 4 kinds of scheme constructions.Further consider pressure break risk and operating cost is true
Fixed final preferred scheme, because used in scheme 5 plus sand intensity is minimum, only 1.2m3/ m, advantageously reduces cost and carries
Construction safety high, final preferred scheme 5 is used as arrangement and method for construction.
In sum, the embodiment of the present invention has taken into full account that many factors of influence effect of increasing production, the parameter of acquisition have
Ambiguity, effect of increasing production and problem of many factors without clear and definite functional relation expression formula, so as to overcome in the prior art
Fracturing parameter local optimum may only be realized rather than all optimal problems, improve the preferred reasonability of fracturing parameter,
Validity and effect of increasing production.
The foregoing is only presently preferred embodiments of the present invention, be not intended to limit the invention, it is all it is of the invention spirit and
Within principle, any modification, equivalent substitution and improvements made etc. should be included within the scope of the present invention.
Claims (7)
1. a kind of fracturing parameter method for optimizing, it is characterised in that comprise the following steps:
S100 collects the reservoir parameter of multigroup fractured well, fracturing parameter, effect of increasing production data and treats the reservoir ginseng of fractured well
Number;
S200 is by the reservoir parameter of collected fractured well and effect of increasing production data and fracturing parameter and effect of increasing production number
According to correlation analysis are carried out, mathematical statistics result is drawn;
S300 is based on the mathematical statistics result, sets up reservoir parameter and effect of increasing production data and fracturing parameter and effect of increasing production
The membership function relation and single factor test fuzzy evaluating matrix of data;
S400 is based on the mathematical statistics result, sets up reservoir parameter using grey correlation theory and accounts for effect of increasing production data and pressure
Split the weight matrix that parameter accounts for effect of increasing production data;
S500 application orthogonal experiment conceptual designs, in combination with fracturing parameter and the reservoir ginseng for treating fractured well of the fractured well
Number, forms multigroup fracturing parameter alternative for treating fractured well;
S600 is based on the fracturing parameter alternative for treating fractured well, the membership function relation set up using fractured well and
The single factor test fuzzy evaluating matrix, the single factor test fuzzy evaluating matrix and the weight matrix are combined, and obtain institute
State the fuzzy overall evaluation matrix of alternative;
The fuzzy overall evaluation matrix of the alternative is converted into fuzzy overall evaluation score by S700, according to maximum membership degree
Principle, final preferred scheme is determined by fuzzy synthesis scoring from low to high from high to low and with reference to construction risk.
2. fracturing parameter method for optimizing according to claim 1, it is characterised in that the reservoir ginseng in the step S100
Number, fracturing parameter, effect of increasing production data include reservoir effective thickness, porosity, gas saturation, compensated neutron, compensation density,
Operational discharge capacity, preposition liquid proportional, average sand ratio, plus sand intensity and effect of increasing production data.
3. fracturing parameter method for optimizing according to claim 2, it is characterised in that the step S200 is specially:
The reservoir parameter of fractured well and effect of increasing production data and fracturing parameter and effect of increasing production data described in being collected into
Single factor test linear regression analysis is carried out, it is determined that the reservoir parameter of the fractured well, fracturing parameter and effect of increasing production are belonging to just
Dependency relation or negative correlativing relation, and line parameter minimum value and most is entered to each parameter according to positive correlation and negative correlativing relation
The opinion rating of 4 class deciles between big value.
4. fracturing parameter method for optimizing according to claim 3, it is characterised in that the step S300 is specially:
The mathematical statistics knot of the reservoir parameter according to the fractured well, the single factor test numerical value between fracturing parameter and effect of increasing production
Really, the membership function of parameters is set up, approximate normal distribution membership function is selected,
In formula:uv- membership function, zero dimension;The occurrence of d-critical parameter, zero dimension;V-opinion rating, common " I ",
" II ", " III ", " IV " 4 class opinion rating, zero dimension;The index parameter of a, b-each opinion rating, zero dimension;
The membership function for determining same opinion rating is,
In formula:uvj(di)—diIn opinion rating vjOn degree of membership, zero dimension;diThe occurrence of-i-th evaluating, nothing
Dimension;
From formula (2):Work as diDuring=a, uvjA ()=1, represents that the evaluating belongs to the degree of membership of some opinion rating most
Greatly, so a is the mathematic expectaion of certain grade, i.e.,
A=(d1+d2)/2 (3)
The equal point of adjacent opinion rating degree of membership is referred to as transition point, and it is 0.5 to take the degree of membership of transition point, i.e.,
In formula:d1、d2The interval upper limit value and lower limit value of-j-th opinion rating, the d of each opinion rating of each evaluating1、d2Value by
Mathematical statistics result determines, after a, b value determine, you can it is determined that known membership function relation;
The single factor test fuzzy evaluating matrix corresponding to each parameter is can be calculated by formula (1)-(5), if total m evaluation etc.
Level, n evaluating, then can obtain a single factor test fuzzy evaluating matrix R for n × m ranks,
5. fracturing parameter method for optimizing according to claim 4, it is characterised in that the step S500 is specially:
According to fractured well fracturing parameter, fractured well fracturing parameter scope is counted, minimum value, maximum in each fracturing parameter
Determine three 3 horizontal factors between value, formed the operational discharge capacity of fractured well fracturing parameter, preposition liquid proportional, average sand ratio,
Plus the scheme of sand intensity, this scheme is one and has four factors, three experimental programs of 3 horizontal factors, sets up L9(34) preferably
The preferred form of orthogonal experiment;
What cooperation was collected into treats fractured well reservoir effective thickness, porosity, gas saturation, compensated neutron, compensation density parameter,
With reference to the preferred form of preferred orthogonal experiment, the alternative of fractured well fracturing parameter is treated in formation.
6. fracturing parameter method for optimizing according to claim 5, it is characterised in that the step S600 is specially:
According to effective thickness, porosity, gas saturation, compensated neutron, compensation density parameter and the fracturing parameter for the treatment of fractured well
It is preferred that the preferred form of orthogonal experiment is combined, applying equation (6) determines the single factor evaluation matrix R of each parameter;
According to the single factor evaluation matrix R, using matrix product summation algorithm, fuzzy weight vector A is obscured with single factor test and is commented
Valency matrix R is combined, and considers all parameter contributions and draws fuzzy overall evaluation matrix B,
7. fracturing parameter method for optimizing according to claim 6, it is characterised in that the step S700 is specially:
According to maximum membership grade principle, the fuzzy overall evaluation matrix of each alternative is converted into fuzzy overall evaluation score,
Fraction is higher, and effect of increasing production is better after representing alternative implementation;When the fuzzy synthesis scoring of different alternatives is identical,
Consider pressure break risk and operating cost determines final preferred scheme;
Computing formula is as follows
Wherein:vI、vII、vIII、vIVThe corresponding fraction of opinion rating is respectively 100,75,50,25.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611202818.7A CN106703776B (en) | 2016-12-23 | 2016-12-23 | Fracturing parameter optimization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611202818.7A CN106703776B (en) | 2016-12-23 | 2016-12-23 | Fracturing parameter optimization method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106703776A true CN106703776A (en) | 2017-05-24 |
CN106703776B CN106703776B (en) | 2020-04-21 |
Family
ID=58903043
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611202818.7A Active CN106703776B (en) | 2016-12-23 | 2016-12-23 | Fracturing parameter optimization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106703776B (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108222910A (en) * | 2018-01-22 | 2018-06-29 | 成都西油联合石油天然气工程技术有限公司 | A kind of fixed displacement slug formula fracturing technology |
CN108446797A (en) * | 2018-03-06 | 2018-08-24 | 西南石油大学 | A kind of compact oil reservoir horizontal well volume fracturing initial productivity prediction technique |
CN108442911A (en) * | 2018-02-28 | 2018-08-24 | 西南石油大学 | A kind of shale gas horizontal well refracturing hydraulic fracture method for optimally designing parameters |
CN108843296A (en) * | 2018-06-25 | 2018-11-20 | 成都北方石油勘探开发技术有限公司 | A kind of individual well refracturing effect prediction method based under multifactor impact |
CN109083641A (en) * | 2018-10-15 | 2018-12-25 | 西南石油大学 | Measure job execution effect evaluation method in a kind of underground workover treatment |
CN109209321A (en) * | 2018-09-03 | 2019-01-15 | 中国石油大学(华东) | It is a kind of based on fracturing potential to pressure break horizontal well fracturing design method and device |
CN109815542A (en) * | 2018-12-21 | 2019-05-28 | 东软集团股份有限公司 | Recognition methods, device, storage medium and the electronic equipment of sand plugging after fracturing |
CN110955985A (en) * | 2019-12-19 | 2020-04-03 | 长江大学 | Method and device for optimizing fracturing construction parameters and readable storage medium |
CN111894545A (en) * | 2019-05-06 | 2020-11-06 | 中国石油化工股份有限公司 | Method for determining proppant pumping scheme |
CN112529443A (en) * | 2020-12-18 | 2021-03-19 | 河南理工大学 | Coal bed gas development target area optimization method based on fuzzy pattern recognition |
CN113033114A (en) * | 2021-03-04 | 2021-06-25 | 中国石油大学(北京) | Optimization method for mineral kinetic parameters in reservoir water rock reaction simulation |
CN113685162A (en) * | 2021-07-22 | 2021-11-23 | 中国石油大学(北京) | Fracturing parameter determination method, device, equipment and storage medium |
CN113833446A (en) * | 2020-06-23 | 2021-12-24 | 中国石油天然气股份有限公司 | Fracturing effect evaluation method and device |
CN113833445A (en) * | 2020-06-08 | 2021-12-24 | 中国石油天然气股份有限公司 | Fracture oil reservoir transformation process |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104018831A (en) * | 2014-06-24 | 2014-09-03 | 西南石油大学 | Method for evaluating reservoir of fractured well |
CN104376420A (en) * | 2014-11-20 | 2015-02-25 | 中国石油天然气股份有限公司 | Water breakthrough risk evaluation method and evaluation device for gas well with water gas reservoir |
CA2945467A1 (en) * | 2014-05-23 | 2015-11-26 | Halliburton Energy Services, Inc. | Enhancing reservoir characterization using real-time srv and fracture evolution parameters |
CN105626009A (en) * | 2015-12-17 | 2016-06-01 | 中国石油天然气股份有限公司 | Quantitative evaluation method for single-well water injection oil replacement effect of fracture-cave carbonate reservoir |
-
2016
- 2016-12-23 CN CN201611202818.7A patent/CN106703776B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2945467A1 (en) * | 2014-05-23 | 2015-11-26 | Halliburton Energy Services, Inc. | Enhancing reservoir characterization using real-time srv and fracture evolution parameters |
CN104018831A (en) * | 2014-06-24 | 2014-09-03 | 西南石油大学 | Method for evaluating reservoir of fractured well |
CN104376420A (en) * | 2014-11-20 | 2015-02-25 | 中国石油天然气股份有限公司 | Water breakthrough risk evaluation method and evaluation device for gas well with water gas reservoir |
CN105626009A (en) * | 2015-12-17 | 2016-06-01 | 中国石油天然气股份有限公司 | Quantitative evaluation method for single-well water injection oil replacement effect of fracture-cave carbonate reservoir |
Non-Patent Citations (1)
Title |
---|
肖勇: ""合川气田须家河组压裂选井与施工参数优选研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108222910A (en) * | 2018-01-22 | 2018-06-29 | 成都西油联合石油天然气工程技术有限公司 | A kind of fixed displacement slug formula fracturing technology |
CN108442911A (en) * | 2018-02-28 | 2018-08-24 | 西南石油大学 | A kind of shale gas horizontal well refracturing hydraulic fracture method for optimally designing parameters |
CN108446797A (en) * | 2018-03-06 | 2018-08-24 | 西南石油大学 | A kind of compact oil reservoir horizontal well volume fracturing initial productivity prediction technique |
CN108446797B (en) * | 2018-03-06 | 2020-02-07 | 西南石油大学 | Method for predicting capacity of tight oil reservoir horizontal well at initial volume fracturing stage |
CN108843296A (en) * | 2018-06-25 | 2018-11-20 | 成都北方石油勘探开发技术有限公司 | A kind of individual well refracturing effect prediction method based under multifactor impact |
CN109209321A (en) * | 2018-09-03 | 2019-01-15 | 中国石油大学(华东) | It is a kind of based on fracturing potential to pressure break horizontal well fracturing design method and device |
CN109209321B (en) * | 2018-09-03 | 2020-01-21 | 中国石油大学(华东) | Fracturing potential-based horizontal well fracturing design method and device to be fractured |
CN109083641A (en) * | 2018-10-15 | 2018-12-25 | 西南石油大学 | Measure job execution effect evaluation method in a kind of underground workover treatment |
CN109815542A (en) * | 2018-12-21 | 2019-05-28 | 东软集团股份有限公司 | Recognition methods, device, storage medium and the electronic equipment of sand plugging after fracturing |
CN111894545A (en) * | 2019-05-06 | 2020-11-06 | 中国石油化工股份有限公司 | Method for determining proppant pumping scheme |
CN111894545B (en) * | 2019-05-06 | 2022-07-22 | 中国石油化工股份有限公司 | Method for determining proppant pumping scheme |
CN110955985A (en) * | 2019-12-19 | 2020-04-03 | 长江大学 | Method and device for optimizing fracturing construction parameters and readable storage medium |
CN113833445A (en) * | 2020-06-08 | 2021-12-24 | 中国石油天然气股份有限公司 | Fracture oil reservoir transformation process |
CN113833445B (en) * | 2020-06-08 | 2023-10-27 | 中国石油天然气股份有限公司 | Fracture reservoir reconstruction process |
CN113833446A (en) * | 2020-06-23 | 2021-12-24 | 中国石油天然气股份有限公司 | Fracturing effect evaluation method and device |
CN113833446B (en) * | 2020-06-23 | 2023-08-22 | 中国石油天然气股份有限公司 | Fracturing effect evaluation method and device |
CN112529443A (en) * | 2020-12-18 | 2021-03-19 | 河南理工大学 | Coal bed gas development target area optimization method based on fuzzy pattern recognition |
CN112529443B (en) * | 2020-12-18 | 2024-05-28 | 河南理工大学 | Coalbed methane development target area optimization method based on fuzzy pattern recognition |
CN113033114A (en) * | 2021-03-04 | 2021-06-25 | 中国石油大学(北京) | Optimization method for mineral kinetic parameters in reservoir water rock reaction simulation |
CN113033114B (en) * | 2021-03-04 | 2022-02-11 | 中国石油大学(北京) | Optimization method for mineral kinetic parameters in reservoir water rock reaction simulation |
CN113685162A (en) * | 2021-07-22 | 2021-11-23 | 中国石油大学(北京) | Fracturing parameter determination method, device, equipment and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN106703776B (en) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106703776A (en) | Method for optimizing fracturing parameters | |
CN111625922B (en) | Large-scale oil reservoir injection-production optimization method based on machine learning agent model | |
CN107358021A (en) | DO prediction model establishment method based on BP neural network optimization | |
CN104018831B (en) | A kind of fractured well reservoir evaluation methods | |
CN105930931A (en) | Electric power engineering cost management method | |
Enke et al. | A hybrid neuro-fuzzy model to forecast inflation | |
CN103440525B (en) | Lake and reservoir water bloom emergency treatment decision-making method based on Vague value similarity measurement improved algorithm | |
CN105931116A (en) | Automated credit scoring system and method based on depth learning mechanism | |
Xiao et al. | The heterogeneous impacts of human capital on green total factor productivity: Regional diversity perspective | |
CN106355540A (en) | Small- and medium-sized reservoir dam safety evaluating method based on GRA-BP (grey relational analysis and back propagation) neural network | |
Liu et al. | Comprehensive evaluation of water resources security in the Yellow River basin based on a fuzzy multi-attribute decision analysis approach | |
CN112070362A (en) | Seasonal river ecological corridor function evaluation method suitable for plain area | |
CN112801687A (en) | Overhead line engineering-based cost prediction model construction method | |
CN106777987A (en) | River water cluster and quality evaluation method based on hydrodynamic force, water quality and SOM | |
Shan et al. | Study on river health assessment weight calculation | |
CN105550804A (en) | Machine tool product manufacturing system energy efficiency evaluation method based on gray fuzzy algorithm | |
CN103942604A (en) | Prediction method and system based on forest discrimination model | |
CN108711100A (en) | A kind of system of the P2P platform operation risk assessment based on neural network | |
CN102789546B (en) | Reference lake quantitative determination method based on human disturbance intensity | |
Liu et al. | Improved design of risk assessment model for PPP project under the development of marine architecture | |
Du et al. | Operating mechanism and set pair analysis model of a sustainable water resources system | |
CN115906669A (en) | Dense residual error network landslide susceptibility evaluation method considering negative sample selection strategy | |
CN106447112A (en) | Construction land scale prediction method of multiple cities based on stack limited Boltzmann machine | |
CN114595628A (en) | Differential transformation method for volume fracturing of horizontal well | |
CN114066166A (en) | Cltobalite low-permeability reservoir CO2Miscible flooding screening method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |