CN110288227A - A method of fracturing effect Dominated Factors are influenced for evaluating - Google Patents
A method of fracturing effect Dominated Factors are influenced for evaluating Download PDFInfo
- Publication number
- CN110288227A CN110288227A CN201910548367.XA CN201910548367A CN110288227A CN 110288227 A CN110288227 A CN 110288227A CN 201910548367 A CN201910548367 A CN 201910548367A CN 110288227 A CN110288227 A CN 110288227A
- Authority
- CN
- China
- Prior art keywords
- node
- micro
- seismic event
- fracturing effect
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000000694 effects Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000011156 evaluation Methods 0.000 claims abstract description 44
- 238000004458 analytical method Methods 0.000 claims abstract description 42
- 239000007787 solid Substances 0.000 claims abstract description 14
- 238000010586 diagram Methods 0.000 claims abstract description 13
- 239000010755 BS 2869 Class G Substances 0.000 claims description 18
- 238000010276 construction Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 12
- 238000003825 pressing Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000003292 glue Substances 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 6
- 239000007788 liquid Substances 0.000 claims description 5
- 238000005086 pumping Methods 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 2
- 230000029305 taxis Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XQCFHQBGMWUEMY-ZPUQHVIOSA-N Nitrovin Chemical compound C=1C=C([N+]([O-])=O)OC=1\C=C\C(=NNC(=N)N)\C=C\C1=CC=C([N+]([O-])=O)O1 XQCFHQBGMWUEMY-ZPUQHVIOSA-N 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
-
- 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/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Economics (AREA)
- Educational Administration (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geophysics And Detection Of Objects (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
Abstract
The invention discloses a kind of for evaluating the method for influencing fracturing effect Dominated Factors, it is characterized in that, the following steps are included: 1) determine the evaluation index and fracturing effect standard diagrams of the influence fracturing effect factor of analysis well section, wherein, fracturing effect standard diagrams are determined as the SRV volume of pressure break post analysis well section；2) micro-seismic event for obtaining analysis well section carries out clustering to micro-seismic event first with knearest neighbour method, and then the envelope volume of SRV carries out grid dividing, and the SRV volume of pressure break post analysis well section is calculated according to the enveloping solid after grid dividing；3) degree of association of each fracture evaluation index and SRV volume in step 2) in step 1) is calculated using Grey Incidence Analysis, it determines the Dominated Factors for influencing fracturing effect in the fracture evaluation index in step 1), completes the evaluation for influencing fracturing effect Dominated Factors.
Description
Technical field
The present invention relates to a kind of for evaluating the method for influencing fracturing effect Dominated Factors, belongs to oil gas development technology neck
Domain.
Background technique
Horizontal well adds hydraulic fracturing technology to be to speed up the important technique measure of shale gas exploitation, and micro-seismic monitoring is live solution
One of man-made fracture shape (seam length, slit width, seam are high) and effective ways of fracture propagation process are released, microseism validity event
It picks up and crack inversion method is the key point that fracturing fracture form precision is explained in influence microseism.
Research emphasis is that crack and crack extension influence factor are portrayed in research microseism mostly at present, is explained with microseism
Crack result is optimization aim, and the work that evaluation influences fracturing effect composite factor is not yet unfolded, causes in live pressing crack construction
It is unknown to improve the controllable construction parameter guidance method of fracturing effect change.
Summary of the invention
In view of the above-mentioned problems, the object of the present invention is to provide a kind of for evaluating the side for influencing fracturing effect Dominated Factors
Method,.
To achieve the above object, the present invention use following technical scheme, one kind for evaluate influence fracturing effect master control because
The method of element, which comprises the following steps:
1) evaluation index and fracturing effect standard diagrams of the influence fracturing effect factor of analysis well section are determined；
Wherein, fracturing effect standard diagrams are determined as the SRV volume of pressure break post analysis well section；
2) micro-seismic event for obtaining analysis well section carries out grid dividing to the envelope volume of SRV, after grid dividing
Enveloping solid calculate pressure break post analysis well section SRV volume；
3) pass of each fracture evaluation index and SRV volume in step 2) in step 1) is calculated using Grey Incidence Analysis
Connection degree determines the Dominated Factors that fracturing effect is influenced in the fracture evaluation index in step 1), completes to influence fracturing effect master control
The evaluation of factor.
Further, in above-mentioned steps 3) in, in Grey Incidence Analysis, to analyze the post-fracturing SRV body of well section
Product is comparison ordered series of numbers, using the fracture evaluation index in step 1) as reference sequence, calculates fracture evaluation index and post-fracturing SRV
The degree of association of volume, then, degree of being associated sequence obtains the inteerelated order of fracture evaluation index, the fracture evaluation in step 1)
The Dominated Factors that fracturing effect is influenced in index complete the evaluation for influencing fracturing effect Dominated Factors.
Further, in above-mentioned steps 2) in, before the enveloping solid grid dividing for carrying out SRV, utilize knearest neighbour method pair
Micro-seismic event carries out clustering, the noise in the micro-seismic event of removal analysis well section, then to the microseism of every one kind
Case point carries out grid dividing, calculates the SRV volume of the micro-seismic event point of every one kind, and finally superposition obtains total SRV volume.
Further, clustering is carried out to micro-seismic event using knearest neighbour method, detailed process is as follows:
1. classifying to micro-seismic event point, and calculate the distance between class and class matrix D(0),
If any sort is Gx(y), indicate that in xth class include y micro-seismic event point；
2. selecting D(0)In least member, be denoted as Dzk, by least member DzkCorresponding class GzWith class GkIt is merged into one newly
Class is denoted as Gl,
3. calculating new class G according to apart from recurrence formulalWith any sort GxDistance, form new Distance matrix D(1),
It is as follows apart from recurrence formula:
In formula, DlxFor new class GlWith any sort GxDistance, DzxFor class GzWith any sort GxBetween minimum range, DkxFor
Class GkWith any sort GxBetween minimum range；I, j indicates the micro-seismic event point in class；
4. to D(1)It repeats the above steps and 2. and 3. obtains D(2), and so on, the merging of class is carried out,
5. set distance threshold value m and microseism number threshold value n carries out the noise identification in micro-seismic event
When the least member in distance matrix is less than distance threshold m, the final classification of micro-seismic event is completed, when in class
Microseism number when being less than n, point present in such is a small number of noise, removes such.
Further, in above-mentioned steps 2) in, grid dividing is carried out using the enveloping solid that Octree theory carries out SRV, and
The calculating of SRV volume is carried out, detailed process is as follows:
1. according to the X of the micro-seismic event point cloud of analysis well section, it is entire micro- to generate encirclement for the distribution of Y, Z point coordinate
The big bounding box of seismic events point obtains the 0th layer of Octree of root node；
2. setting the target density threshold value of Octree algorithm as t, Octree cube division is carried out to root node, is obtained tree-shaped
Gridding structure,
Point cloud density in each non-empty node is judged according to target density threshold value, if the point Yun Mi in non-empty node
When degree is greater than target density threshold value, then recurrence division is carried out to the node along three change in coordinate axis direction, in each iteration
Cuboid is divided into 8 congruent small cuboids, so that the corresponding child node of a node layer is obtained, until each non-empty
The number of point in child node all meets density threshold, stops loop iteration, and one is at least contained in each non-empty node micro-ly
Shake case point, thus by the tree network that non-empty node forms format structure can completely all micro-seismic event points of envelope,
The empty node in gridding structure, the i.e. non-response area of micro-seismic event in removal root node are rejected,
3. through step, 2. micro-seismic event point is divided into n non-empty child node, the tree that Octree divides
Maximum level be hm；If three-dimensional micro-seismic event is within the scope of work area
M={ mi|mi=(xi, yi, zi), i=1,2 ..., | M | }
In formula, M is three-dimensional micro-seismic event point set, miFor case point subset in i-th after division of non-empty node, xi,
yi, ziIndicate miIn any one occurrence point coordinate；miSize meet setting threshold value, i.e. mi≤t；
The corresponding global number collection of all non-empty child nodes is obtained by Octree multilayer tree to be combined into
B={ bi| i=1,2,3..., n }
In formula, B is the corresponding global index's number of non-empty node, biFor the number of i-th of non-empty node, wherein each
biPositive integer is taken, and meets 1≤bi≤M；
The corresponding relationship of node level and node serial number is obtained by Octree multilayer index structure, obtains the calculating of node level
Function H (x`), to obtain the corresponding level of non-empty child node are as follows:
hi=H (bi)
In formula, hiFor the corresponding division level of i-th of node, and hiRange meet [0, hm], H (x`) is number x`'s
Level calculates function, and H (x`) is single mapping function；
It is the 0th layer since root node is corresponding, if root node is combined into L along the side length collection of three change in coordinate axis direction0,
L0={ l0 j| j=1,2,3 }
In formula, l0 j| j=1,2,3 respectively indicate root node in x, y, the node side length in the direction z,
Then i-th layer of node is expressed as l along the side length in the direction ji j, wherein j=1,2,3 indicate that i-th layer of node exists
The side length in the direction x, y, z；
All be eight equal parts due to dividing each time, be equivalent to the bisection along single dimension, then obtain i-th layer it is single non-
Side length of the gap node along the direction j be
The volume V of corresponding i-th layer single non-empty nodei
Vi=li 1li 2li 3；
4. being added up to obtain the reservoir of total estimation according to different levels to all non-empty child node volumes respectively
Volume, which is transformed, is
Further, above-mentioned steps 2. in, cuboid is used into trisection on stitching high direction in each iteration,
For remaining dimension using halving, 12 congruent small cuboids will be obtained every time by dividing；
Above-mentioned steps 3. in, i-th layer of node is expressed as l along the side length in the direction ji j, wherein j=1,2,3 indicate
I-th layer of node is in x, y, the side length in the direction z, stitches high directional spreding in x, y, either in the direction z upwards；
Due to using trisection on stitching high direction, for remaining dimension using halving, then dividing each time all is ten second-class
Point, then obtaining i-th layer of single non-empty child node along the side length in the direction j is
The side length in the high direction of seam in the direction j is
The side length in the high direction of non-seam in the direction j is,
Further, in above-mentioned steps 3) in, the nothing of data is carried out in Grey Incidence Analysis using initial method
Dimension processing.
Further, in above-mentioned steps 1) in, fracture evaluation index is the controllable parameter in pressing crack construction.
Further, in above-mentioned steps 1) in, fracture evaluation index include perforation length, perforation number of clusters, operational discharge capacity,
Smooth water consumption, pumping liquid measure, 100 mesh haydite dosages, 40/70 haydite dosage, 30/50 haydite dosage, is averaged at linear glue dosage
Sand ratio.
Further, in above-mentioned steps 1) in, analysis well section chooses the fractured well of same fractured interval.
The invention adopts the above technical scheme, and it is a kind of for evaluating influence pressure to have the advantages that the 1, present invention provides
The method for splitting effect Dominated Factors, to the evaluation index (reference sequence) and fracturing effect standard diagrams in pressing crack construction process
(evaluation ordered series of numbers) is determined, and is carried out grid dividing to the enveloping solid of SRV, SRV volume computational accuracy is improved, using grey correlation
Method Calculation Estimation index and the fracturing effect standard diagrams degree of association complete the evaluation for influencing the Dominated Factors of fracturing effect, root
According to the obtained degree of association, the Dominated Factors of the live pressing crack construction of influence are determined, give in live pressing crack construction as raising fracturing effect
Change controllable construction parameter and data support is provided, there is directive significance.
2, the present invention carries out clustering to micro-seismic event using knearest neighbour method, picks before carrying out the calculating of SRV volume
Except the interference due to engineering, monitoring data will appear invalid event point, the computational accuracy of SRV volume is improved, and then improve evaluation
The accuracy of index and fracturing effect standard diagrams calculation of relationship degree, it is accurate to determine the Dominated Factors for influencing live pressing crack construction,
Data support is provided to improve the controllable construction parameter of fracturing effect change.
Detailed description of the invention
Fig. 1 is overall flow schematic diagram of the invention；
Fig. 2 is the 6th section of 48-2 well of micro-seismic event figure, and Fig. 2 a is the micro-seismic event figure before the 6th section of 48-2 well denoising,
Fig. 2 b is the micro-seismic event figure after the 6th section of 48-2 well denoising；
Fig. 3 is the micro-seismic event grid chart after the 6th section of 48-2 well denoising；
Fig. 4 is the micro-seismic event grid chart before the 6th section of 48-2 well denoising；
Fig. 5 is that six sections of well denoising front and back SRV volumes calculate Comparative result histogram.
Specific embodiment
Presently preferred embodiments of the present invention is described in detail below with reference to attached drawing, it is of the invention to be clearer to understand
Objects, features and advantages.It should be understood that embodiment shown in the drawings does not limit the scope of the present invention, and only it is
Illustrate the connotation of technical solution of the present invention.
As shown in Figure 1, the present invention provides a kind of for evaluating the method for influencing fracturing effect Dominated Factors comprising with
Lower step:
1) according to site operation demand, the fracture evaluation index and fracturing effect standard diagrams of analysis well section are determined；
Fracture evaluation index is the controllable parameter in pressing crack construction, such as perforation length, perforation number of clusters, operational discharge capacity, cunning
Slip water consumption, linear glue dosage, pumping liquid measure, 100 mesh haydite dosages, 40/70 haydite dosage, 30/50 haydite dosage, average sand
Than etc. 10 indexs；Fracturing effect standard diagrams are determined as the SRV volume of pressure break post analysis well section.
2) micro-seismic event for obtaining analysis well section carries out grid dividing to the envelope volume of SRV, after grid dividing
Enveloping solid calculate pressure break post analysis well section SRV volume.
In traditional algorithm, often micro-seismic event is wrapped with simple geometry when calculating SRV volume
The size of network, enveloping solid depends on the span of farthest two micro-seismic events point in each dimension in micro-seismic event；This side
Method will lead to the non-response area that will appear blank in enveloping solid, i.e. the region that does not feed through to of hydraulic fracturing, and SRV volume is caused to increase
Greatly, in addition to this, conventional method does not consider the overlapping of micro-seismic event in enveloping solid.
Grid dividing is carried out to the enveloping solid of SRV, the number of grid in the place more than case point quantity is with regard to more, event dot density
Big local grid dividing it is finer, to improve the computational accuracy of SRV volume.
3) each fracture evaluation index and post-fracturing SRV in step 2) in step 1) are calculated using Grey Incidence Analysis
The degree of association of volume determines the Dominated Factors that fracturing effect is influenced in the fracture evaluation index in step 1), completes to influence pressure break
The evaluation of effect Dominated Factors.
In Grey Incidence Analysis, to analyze the post-fracturing SRV volume of well section for comparison ordered series of numbers, in step 1)
Fracture evaluation index is reference sequence, calculates the degree of association of each fracture evaluation index Yu post-fracturing SRV volume, then, is carried out
Relational degree taxis obtains the inteerelated order of fracture evaluation index, and determining in the fracture evaluation index in step 1) influences fracturing effect
Dominated Factors, complete influence fracturing effect Dominated Factors evaluation.
Further, in step 1), analysis well section chooses the fractured well of same fractured interval, in analyzing influence pressure break
The influence of the objective geological conditions such as elasticity modulus, brittleness is excluded during effect Dominated Factors.
Further, when scene is monitored microseism using micro-seismic monitoring means, often due to engineering, prison
The interference of measured data will appear invalid event point, and the presence of invalid event point will affect SRV computational accuracy, therefore, in step 2)
In, before the enveloping solid grid dividing for carrying out SRV, clustering, removal point are carried out to micro-seismic event using knearest neighbour method
The noise (Null Spot case point) in the micro-seismic event of well section is analysed, grid then is carried out to the micro-seismic event point of every one kind and is drawn
Point, the SRV volume of the micro-seismic event point of every one kind is calculated, finally superposition obtains total SRV volume.
Further, clustering is carried out to micro-seismic event using knearest neighbour method, detailed process is as follows:
1. classifying to micro-seismic event point, and calculate the distance between class and class matrix D(0),
If any sort is Gx(y), indicate that in xth class include y micro-seismic event point.
2. selecting D(0)In least member, be denoted as Dzk, by least member DzkCorresponding class GzWith class GkIt is merged into one newly
Class is denoted as Gl,
3. calculating new class G according to apart from recurrence formulalWith any sort GxDistance, form new Distance matrix D(1),
Apart from recurrence formula are as follows:
In formula, DlxFor new class GlWith any sort GxDistance, DzxFor class GzWith any sort GxBetween minimum range, DkxFor
Class GkWith any sort GxBetween minimum range；I, j indicates the micro-seismic event point in class.
4. to D(1)It repeats the above steps and 2. and 3. obtains D(2), and so on, the merging of class is carried out,
5. set distance threshold value m and microseism number threshold value n carries out the noise identification in micro-seismic event,
When the least member in distance matrix is less than distance threshold m, the final classification of micro-seismic event is completed, when in class
Microseism number when being less than n, it is believed that point present in such is a small number of noise, it should remove such.
Further, in step 2), grid dividing is carried out using the enveloping solid that Octree theory carries out SRV, and carry out
SRV volume calculates, and detailed process is as follows:
1. according to the X of the micro-seismic event point cloud of analysis well section, it is entire micro- to generate encirclement for the distribution of Y, Z point coordinate
The big bounding box of seismic events point obtains the 0th layer of Octree of root node；
2. setting the target density threshold value of Octree algorithm as t, recurrence division is carried out to root node, tree network is obtained and formats
Structure, detailed process are as follows:
Point cloud density in each non-empty node is judged according to target density threshold value, if the point Yun Mi in non-empty node
When degree is greater than target density threshold value, then recurrence division is carried out to the node along three change in coordinate axis direction, in each iteration
Cuboid is divided into 8 congruent small cuboids, so that the corresponding child node of a node layer is obtained, until each non-empty
The number of point in child node all meets density threshold, stops loop iteration, and one is at least contained in each non-empty node micro-ly
Shake case point, thus by the tree network that non-empty node forms format structure can completely all micro-seismic event points of envelope,
The empty node in gridding structure, the i.e. non-response area of micro-seismic event in removal root node are rejected,
3. micro-seismic event point is 2. divided into n non-empty child node, the tree-shaped knot that Octree divides through step
The maximum level of structure is hm；If three-dimensional micro-seismic event is within the scope of work area
M={ mi|mi=(xi, yi, zi), i=1,2 ..., | M | }
In formula, M is three-dimensional micro-seismic event point set, miFor case point subset in i-th after division of non-empty node, xi,
yi, ziIndicate miIn any one occurrence point coordinate；miSize meet setting threshold value, i.e. mi≤t；
The corresponding global number collection of all non-empty child nodes is obtained by Octree multilayer tree to be combined into
B={ bi| i=1,2,3..., n }
In formula, B is the corresponding global index's number of non-empty node, biFor the number of i-th of non-empty node, each biIt takes just
Integer, and meet 1≤bi≤M。
The corresponding relationship of node level and node serial number is obtained by Octree multilayer index structure, obtains the calculating of node level
Function H (x`), to obtain the corresponding level of non-empty child node are as follows:
hi=H (bi)
In formula, hiFor the corresponding division level of i-th of node, and hiRange meet [0, hm], H (x`) is number x`'s
Level calculates function, and H (x`) is single mapping function.
It is the 0th layer since root node is corresponding, if root node is combined into L along the side length collection of three change in coordinate axis direction0,
L0={ l0 j| j=1,2,3 } (5)
In formula, l0 j| j=1,2,3 respectively indicate root node in x, y, the node side length in the direction z,
I-th layer of node is expressed as l along the side length in the direction ji j, wherein j=1,2,3, indicate i-th layer of node in x,
The side length in the direction y, z；
All it is eight equal parts due to dividing each time, is equivalent to the bisection along single dimension, then available i-th layer of list
Side length of a non-empty child node along the direction j be
Then obtain the volume V of corresponding i-th layer single non-empty nodei
Vi=li 1li 2li 3
In formula, li 1For the i-th layer of single side length of non-empty node in the direction x, li 2It is i-th layer of single non-empty node in the direction y
Side length, li 3Side length for i-th layer of single non-empty node in the direction z；
4. being added up to obtain the reservoir of total estimation according to different levels to all non-empty child node volumes respectively
Volume, which is transformed, is
Further, cuboid is all divided into eight cuboids by iteration each time in Octree theory, in x, y, the direction z
On be all made of 2iIt is iterated, does not consider the problems of that certain dimension needs to encrypt.
And in fracturing process, to prevent crack from pressing to wear payzone, construction personnel is more concerned about the expansion that fracturing fracture stitches high direction
Exhibition, therefore when calculating SRV volume, the grid for stitching high direction needs to encrypt, and therefore, divides to Octree cube
Method improves, and in each iteration willIt is longCube uses trisection on stitching high direction, remaining dimension, which uses, to be halved,
12 cuboids will be obtained by dividing every time, can be increased the accurate of SRV on the high direction of seam and be portrayed.
Further, detailed process is as follows for improved Octree algorithm:
1. according to the X of the micro-seismic event point cloud of analysis well section, it is entire micro- to generate encirclement for the distribution of Y, Z point coordinate
The big bounding box of seismic events point obtains the 0th layer of Octree of root node；
2. setting the target density threshold value of Octree algorithm as t, grid dividing is carried out to root node, tree network is obtained and formats
Structure,
Point cloud density in each non-empty node is judged according to target density threshold value, if the point Yun Mi in non-empty node
When degree is greater than target density threshold value, then recurrence division is carried out to the node along three change in coordinate axis direction, in each iteration
Cuboid is divided into 12 congruent small cuboids, so that the corresponding child node of a node layer is obtained, until each non-empty
Child node in the number of point all meet density threshold, stop loop iteration, it is at least micro- containing one in each non-empty node
Seismic events point, therefore formatting structure by the tree network that non-empty node forms can completely all micro-seismic event of envelope
Point rejects the empty node in gridding structure, the i.e. non-response area of micro-seismic event in removal root node,
3. through step, 2. micro-seismic event point is divided into n non-empty child node, the tree that Octree divides
Maximum level be hm；If three-dimensional micro-seismic event is within the scope of work area
M={ mi|mi=(xi, yi, zi), i=1,2 ..., | M | }
In formula, M is three-dimensional micro-seismic event point set, miFor case point subset in i-th after division of non-empty node, xi,
yi, ziIndicate miAny one occurrence point coordinate；miSize meet setting threshold value, i.e. mi≤t；
The corresponding global number collection of all non-empty child nodes is obtained by Octree multilayer tree to be combined into
B={ bi| i=1,2,3..., n }
In formula, B is the corresponding global index's number of non-empty node, biFor the number of i-th of non-empty node, wherein each bi
Positive integer is taken, and meets 1≤bi≤M。
The corresponding relationship of node level and node serial number is obtained by Octree multilayer index structure, obtains the calculating of node level
Function H (x`), to obtain the corresponding level of non-empty child node are as follows:
hi=H (bi)
In formula, hiFor the corresponding division level of i-th of node, and hiRange meet [0, hm], H (x`) is number x`'s
Level calculates function, and H (x`) is single mapping function.
It is the 0th layer since root node is corresponding, if root node is combined into L along the side length collection of three change in coordinate axis direction0,
L0={ l0 j| j=1,2,3 }
In formula, l0 j| j=1,2,3 respectively indicate root node in x, y, the node side length in the direction z,
Then i-th layer of node is expressed as l along the side length in the direction ji j, wherein j=1,2,3 indicate that i-th layer of node exists
The side length in the direction x, y, z, stitches high directional spreding in x, y, either in the direction z upwards；
Due to using trisection on stitching high direction, for remaining dimension using halving, dividing each time all is ten bisections,
Then available i-th layer single side length of the non-empty child node along the direction j is
The side length in the direction the j direction Zhong Fenggao is
The side length in the non-high direction of seam is in the direction j
The volume V of corresponding i-th layer single non-empty nodei
Vi=li 1li 2li 3
In formula, li 1For the i-th layer of single side length of non-empty node in the direction x, li 2It is i-th layer of single non-empty node in the direction y
Side length, li 3Side length for i-th layer of single non-empty node in the direction z；
4. being added up to obtain the reservoir of total estimation according to different levels to all non-empty child node volumes respectively
Volume, which is transformed, is
The present invention is explained with specific embodiment below:
By taking reef dam block as an example, illustrate the process of analyzing influence fracturing effect Dominated Factors.
Step (1): the influence in order to exclude the objective geological conditions such as elasticity modulus, brittleness is influencing fracturing effect master control
The fractured well parameter of same fractured interval is chosen in the analysis of factor, is chosen and is walked the pressure break well section that layer position is 3 sections of Longma small stream, it is main
If 1 section of 48-1 well, 22 sections, 9 sections of 48-2 well, 11 sections, 13 sections, 3 sections of 37-3 well, 23 sections of 2-5 well totally 7 sections as analysis well section.
According to site operation needs, the evaluation index for analyzing well section is determined as perforation length, perforation number of clusters, construction row
Amount, smooth water consumption, linear glue dosage, pumping liquid measure, 100 mesh haydite dosages, 40/70 haydite dosage, 30/50 haydite dosage,
10 indexs such as average sand ratio, the results are shown in Table 1 for indicator-specific statistics.Fracturing effect standard diagrams are determined as pressure break post analysis well section
SRV volume.
The evaluation index statistical form of the analysis well section of table 1
Step (2): the SRV volume of analysis well section is calculated
With 13,22 sections of 48-1 well, the 6th, 13 section of 48-2 well, 10,13 sections of 48-3 well etc. 6 sections as the study calculated SRV
Sample first identifies the available point of the micro-seismic event of above-mentioned well section, recycles and carries out to the micro-seismic event after identification
Octree cube divides, and finally calculates SRV volume.
Below for the 6th section of 48-2 well, illustrate entire optimization process,
Micro-seismic event using shortest distance clustering procedure to the 6th section of 48-2 well carries out noise identification, and Fig. 2 a is 48-2 well
Micro-seismic event figure before 6th section of denoising, Fig. 2 b are the micro-seismic event figures after the 6th section of 48-2 well denoising.
After scheming it can be found that being identified by noise, the case point quantity in micro-seismic event figure that the 6th section of 48-2 well
It reduces, the division of Octree cube is carried out to the micro-seismic event available point after denoising, as shown in Figure 3；According to grid dividing knot
It is 177412.3808m that fruit, which calculates the 6th section of 48-2 well of SRV volume after denoising,3。
By it is original do not denoise before micro-seismic event carry out Octree cube division, as shown in figure 4, according to grid dividing knot
It is 1978939.406m that fruit, which calculates and denoises the SRV volume of preceding the 6th section of well of 48-2,3, the SRV of the 6th section of well of 48-2 of comparison denoising front and back
Volume is it is found that the 6th section of 48-2 well of SRV volume reduces by 91.04% after denoising.
Then, the SRV volume of six sections of wells is calculated system by the SRV volume that remaining five groups of well section is calculated according to same thinking
Meter gets up, as shown in table 2,
The contrast table of 2 six sections of well denoising front and back SRV volumes of table
Six sections of well denoising front and back SRV volumes are formed according to 2 data of table and calculate Comparative result histogram, as shown in figure 5, from table
2 can be seen that, the reduction rate of denoising front and back SRV volume is between 45.3%-91.4%.It can be seen that the presence of noise is to SRV's
Calculated result influence is very big, and it is necessary to identify to validity event point before carrying out the calculating of SRV volume.
Therefore, before carrying out the calculating of SRV volume to the analysis well section in step (1), first with shortest distance clustering procedure to 48-
11 section of well, 22 sections, 9 sections of 48-2 well, 11 sections, 13 sections, 3 sections of 37-3 well, 23 sections of 2-5 well totally 7 sections micro-seismic event carry out noise
Then removal calculates the SRV volume of 7 sections of well sections, calculated result is as shown in table 3,
The SRV volume calculation result table of the analysis well section of table 3
Step (3): the SRV calculated result of well section will be analyzed as the factor of evaluation fracturing effect, the i.e. ratio of gray Analysis
Compared with ordered series of numbers, 10 indexs mentioned in step (1) are subjected to correlation analysis using Grey Incidence, analyze result such as table 4
It is shown,
4 Grey Incidence of table analyzes result table
As shown in Table 4, incidence degree sequence is that average sand ratio > operational discharge capacity > linear glue dosage > perforation number of clusters > slippery water is used
Amount > pumping liquid measure > perforation length > 100 haydite dosages > 40/70 haydite dosage > 30/50 haydite dosage, therefore average sand ratio,
Operational discharge capacity and linear glue dosage are to influence the Dominated Factors of live pressing crack construction, and scene can according to actual needs, and emphasis closes
Infuse these three construction parameters.
The present invention is only illustrated with above-described embodiment, and structure, setting position and its connection of each component are all can have
Changed.Based on the technical solution of the present invention, the improvement or equivalent that all principles according to the present invention carry out individual part
Transformation, should not exclude except protection scope of the present invention.
Claims (10)
1. a kind of for evaluating the method for influencing fracturing effect Dominated Factors, which comprises the following steps:
1) evaluation index and fracturing effect standard diagrams of the influence fracturing effect factor of analysis well section are determined,
Wherein, fracturing effect standard diagrams are determined as the SRV volume of pressure break post analysis well section；
2) micro-seismic event for obtaining analysis well section carries out grid dividing to the envelope volume of SRV, according to the packet after grid dividing
The SRV volume of network body calculating pressure break post analysis well section；
3) each fracture evaluation index is calculated in step 1) using Grey Incidence Analysis to be associated with SRV volume in step 2)
Degree, determine in the fracture evaluation index in step 1) influence fracturing effect Dominated Factors, complete influence fracturing effect master control because
The evaluation of element.
2. as described in claim 1 a kind of for evaluating the method for influencing fracturing effect Dominated Factors, it is characterised in that: upper
It states in step 3), in Grey Incidence Analysis, to analyze the post-fracturing SRV volume of well section for comparison ordered series of numbers, with step 1)
In fracture evaluation index be reference sequence, calculate fracture evaluation index and post-fracturing SRV volume the degree of association, then, into
Row relational degree taxis obtains the inteerelated order of fracture evaluation index, and fracturing effect is influenced in the fracture evaluation index in step 1)
Dominated Factors complete the evaluation for influencing fracturing effect Dominated Factors.
3. as described in claim 1 a kind of for evaluating the method for influencing fracturing effect Dominated Factors, it is characterised in that: upper
It states in step 2), before the enveloping solid grid dividing for carrying out SRV, cluster point is carried out to micro-seismic event using knearest neighbour method
The noise in the micro-seismic event of well section is analyzed in analysis, removal, then carries out grid dividing, meter to the micro-seismic event point of every one kind
The SRV volume of the micro-seismic event point of every one kind is calculated, finally superposition obtains total SRV volume.
4. as claimed in claim 3 a kind of for evaluating the method for influencing fracturing effect Dominated Factors, which is characterized in that utilize
Knearest neighbour method carries out clustering to micro-seismic event, and detailed process is as follows:
1. classifying to micro-seismic event point, and calculate the distance between class and class matrix D(0),
If any sort is Gx(y), indicate that in xth class include y micro-seismic event point；
2. selecting D(0)In least member, be denoted as Dzk, by least member DzkCorresponding class GzWith class GkIt is merged into a new class,
It is denoted as Gl,
3. calculating new class G according to apart from recurrence formulalWith any sort GxDistance, form new Distance matrix D(1),
It is as follows apart from recurrence formula:
In formula, DlxFor new class GlWith any sort GxDistance, DzxFor class GzWith any sort GxBetween minimum range, DkxFor class Gk
With any sort GxBetween minimum range；I, j indicates the micro-seismic event point in class；
4. to D(1)It repeats the above steps and 2. and 3. obtains D(2), and so on, the merging of class is carried out,
5. set distance threshold value m and microseism number threshold value n carries out the noise identification in micro-seismic event
When the least member in distance matrix is less than distance threshold m, the final classification of micro-seismic event is completed, when micro- in class
When earthquake number is less than n, point present in such is a small number of noise, removes such.
5. as described in claim 1 a kind of for evaluating the method for influencing fracturing effect Dominated Factors, which is characterized in that upper
It states in step 2), grid dividing is carried out using the enveloping solid that Octree theory carries out SRV, and carry out SRV volume calculating, specific mistake
Journey is as follows:
1. according to the X of the micro-seismic event point cloud of analysis well section, the distribution of Y, Z point coordinate generates and surrounds entire microseism
The big bounding box of case point obtains the 0th layer of Octree of root node；
2. setting the target density threshold value of Octree algorithm as t, Octree cube division is carried out to root node, obtains tree-like mesh
Change structure,
Point cloud density in each non-empty node is judged according to target density threshold value, if the point cloud density in non-empty node is big
When target density threshold value, then recurrence division is carried out to the node along three change in coordinate axis direction, will grown in each iteration
Cube is divided into 8 congruent small cuboids, so that the corresponding child node of a node layer is obtained, until the son section of each non-empty
The number of point in point all meets density threshold, stops loop iteration, at least contains a microseism thing in each non-empty node
Part point, thus by the tree network that non-empty node forms format structure can completely all micro-seismic event points of envelope, reject
Empty node in gridding structure, the i.e. non-response area of micro-seismic event in removal root node,
3. through step, 2. micro-seismic event point is divided into n non-empty child node, and the tree that Octree divides is most
Big level is hm；If three-dimensional micro-seismic event is within the scope of work area
M={ mi|mi=(xi, yi, zi), i=1,2 ..., | M | }
In formula, M is three-dimensional micro-seismic event point set, miFor case point subset in i-th after division of non-empty node, xi, yi, ziTable
Show miIn any one occurrence point coordinate；miSize meet setting threshold value, i.e. mi≤t；
The corresponding global number collection of all non-empty child nodes is obtained by Octree multilayer tree to be combined into
B={ bi| i=1,2,3..., n }
In formula, B is the corresponding global index's number of non-empty node, biFor the number of i-th of non-empty node, wherein each biIt takes
Positive integer, and meet 1≤bi≤M；
The corresponding relationship of node level and node serial number is obtained by Octree multilayer index structure, node level is obtained and calculates function
H (x`), to obtain the corresponding level of non-empty child node are as follows:
hi=H (bi)
In formula, hiFor the corresponding division level of i-th of node, and hiRange meet [0, hm], H (x`) is the level of number x`
Function is calculated, and H (x`) is single mapping function；
It is the 0th layer since root node is corresponding, if root node is combined into L along the side length collection of three change in coordinate axis direction0,
L0={ l0 j| j=1,2,3 }
In formula, l0 j| j=1,2,3 respectively indicate root node in x, y, the node side length in the direction z,
Then i-th layer of node is expressed as l along the side length in the direction ji j, wherein j=1,2,3 indicate i-th layer of node in x, y, z
The side length in direction；
All it is eight equal parts due to dividing each time, is equivalent to the bisection along single dimension, then obtains i-th layer of single non-gap
Side length of the node along the direction j be
The volume V of corresponding i-th layer single non-empty nodei
Vi=li 1li 2li 3；
4. being added up to obtain the reservoir reconstruction of total estimation according to different levels to all non-empty child node volumes respectively
Volume is
6. as claimed in claim 5 a kind of for evaluating the method for influencing fracturing effect Dominated Factors, it is characterised in that:
Above-mentioned steps 2. in, cuboid is used into trisection on stitching high direction in each iteration, remaining dimension uses
It halves, 12 congruent small cuboids will be obtained every time by dividing；
Above-mentioned steps 3. in, i-th layer of node is expressed as l along the side length in the direction ji j, wherein j=1,2,3 indicate i-th layer
Node in x, y, the side length in the direction z stitches high directional spreding in x, y, either in the direction z upwards；
Due to using trisection on stitching high direction, remaining dimension is using halving, then dividing each time all is ten bisections, then
I-th layer of single non-empty child node, which is obtained, along the side length that the side length in the direction j is the high direction of seam in the direction j is
The side length in the high direction of non-seam in the direction j is,
7. as described in claim 1 a kind of for evaluating the method for influencing fracturing effect Dominated Factors, it is characterised in that:
In above-mentioned steps 3) in, the dimensionless processing of data is carried out in Grey Incidence Analysis using initial method.
8. as described in claim 1 a kind of for evaluating the method for influencing fracturing effect Dominated Factors, it is characterised in that: upper
It states in step 1), fracture evaluation index is the controllable parameter in pressing crack construction.
9. as claimed in claim 8 a kind of for evaluating the method for influencing fracturing effect Dominated Factors, it is characterised in that: upper
It states in step 1), fracture evaluation index includes perforation length, perforation number of clusters, operational discharge capacity, smooth water consumption, linear glue use
Amount, pumping liquid measure, 100 mesh haydite dosages, 40/70 haydite dosage, 30/50 haydite dosage, average sand ratio.
10. as described in claim 1 a kind of for evaluating the method for influencing fracturing effect Dominated Factors, it is characterised in that:
Above-mentioned steps 1) in, analysis well section chooses the fractured well of same fractured interval.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910548367.XA CN110288227A (en) | 2019-06-24 | 2019-06-24 | A method of fracturing effect Dominated Factors are influenced for evaluating |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910548367.XA CN110288227A (en) | 2019-06-24 | 2019-06-24 | A method of fracturing effect Dominated Factors are influenced for evaluating |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110288227A true CN110288227A (en) | 2019-09-27 |
Family
ID=68004938
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910548367.XA Pending CN110288227A (en) | 2019-06-24 | 2019-06-24 | A method of fracturing effect Dominated Factors are influenced for evaluating |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110288227A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111859260A (en) * | 2020-07-03 | 2020-10-30 | 中国科学技术大学 | Method for obtaining SRV area by using fracturing construction data |
| CN112855109A (en) * | 2020-12-31 | 2021-05-28 | 西南石油大学 | Fracture acidizing well selecting and layer selecting method based on grey correlation method and analytic hierarchy process |
| CN113687411A (en) * | 2021-09-03 | 2021-11-23 | 西南石油大学 | Earth stress azimuth prediction method based on microseism |
| CN116482752A (en) * | 2023-04-26 | 2023-07-25 | 四川大学 | Rapid three-dimensional positioning method for microseismic event in complex rock mass engineering structure |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104769215A (en) * | 2013-05-29 | 2015-07-08 | 雪佛龙美国公司 | System and method for characterizing uncertainty in subterranean reservoir fracture networks |
| CN105223611A (en) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | Based on the pressure break reservoir leading edge of the three-dimensional point set of micro-seismic monitoring, area and calculation method of physical volume |
| CN107480411A (en) * | 2017-10-20 | 2017-12-15 | 合肥工业大学 | Reservoir fracturing effect evaluation method and evaluation system |
| CN108073786A (en) * | 2017-12-27 | 2018-05-25 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | Shale compressibility evaluation method and device |
| CN108960651A (en) * | 2018-07-11 | 2018-12-07 | 西南石油大学 | A kind of integrated evaluating method of densification oil-gas reservoir multistage fracturing horizontal well completion efficiency |
-
2019
- 2019-06-24 CN CN201910548367.XA patent/CN110288227A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104769215A (en) * | 2013-05-29 | 2015-07-08 | 雪佛龙美国公司 | System and method for characterizing uncertainty in subterranean reservoir fracture networks |
| CN105223611A (en) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | Based on the pressure break reservoir leading edge of the three-dimensional point set of micro-seismic monitoring, area and calculation method of physical volume |
| CN107480411A (en) * | 2017-10-20 | 2017-12-15 | 合肥工业大学 | Reservoir fracturing effect evaluation method and evaluation system |
| CN108073786A (en) * | 2017-12-27 | 2018-05-25 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | Shale compressibility evaluation method and device |
| CN108960651A (en) * | 2018-07-11 | 2018-12-07 | 西南石油大学 | A kind of integrated evaluating method of densification oil-gas reservoir multistage fracturing horizontal well completion efficiency |
Non-Patent Citations (3)
| Title |
|---|
| 刘星 等: "基于八叉树的储层改造体积计算方法及工程应用", 《岩石力学与工程学报》 * |
| 吴爱祥 等: "矿山微地震活动时空分布", 《北京科技大学学报》 * |
| 杨杨: "深层页岩气压裂优化设计及效果评估研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111859260A (en) * | 2020-07-03 | 2020-10-30 | 中国科学技术大学 | Method for obtaining SRV area by using fracturing construction data |
| CN112855109A (en) * | 2020-12-31 | 2021-05-28 | 西南石油大学 | Fracture acidizing well selecting and layer selecting method based on grey correlation method and analytic hierarchy process |
| CN112855109B (en) * | 2020-12-31 | 2022-09-16 | 西南石油大学 | Fracture acidizing well selecting and layer selecting method based on grey correlation method and analytic hierarchy process |
| CN113687411A (en) * | 2021-09-03 | 2021-11-23 | 西南石油大学 | Earth stress azimuth prediction method based on microseism |
| CN116482752A (en) * | 2023-04-26 | 2023-07-25 | 四川大学 | Rapid three-dimensional positioning method for microseismic event in complex rock mass engineering structure |
| CN116482752B (en) * | 2023-04-26 | 2024-01-23 | 四川大学 | Rapid three-dimensional positioning method for microseismic event in complex rock mass engineering structure |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110288227A (en) | A method of fracturing effect Dominated Factors are influenced for evaluating | |
| CN106372402B (en) | The parallel method of fuzzy region convolutional neural networks under a kind of big data environment | |
| CN110674841B (en) | Logging curve identification method based on clustering algorithm | |
| CN104297785B (en) | Petrofacies constraint reservoir physical parameter inversion method and device | |
| CN103775075B (en) | A kind of full well section Lithology Identification Methods | |
| CN104462184B (en) | A kind of large-scale data abnormality recognition method based on two-way sampling combination | |
| CN107894615B (en) | A kind of method of quantitative evaluation 3-D seismics attribute forecast reservoir parameter validity | |
| CN105512366B (en) | The tree-shaped random seam net description method of compact reservoir volume fracturing containing intrinsic fracture | |
| Libohova et al. | Geomorphons: Landform and property predictions in a glacial moraine in Indiana landscapes | |
| CN104101902B (en) | Seismic properties clustering method and device | |
| CN106355011A (en) | Geochemical data element sequence structure analysis method and device | |
| CN105447762B (en) | A kind of calculation method of the low-permeability oil deposit water logging information of fluid replacement | |
| CN106121622B (en) | A kind of Multiple faults diagnosis approach of the Dlagnosis of Sucker Rod Pumping Well based on indicator card | |
| Dikbas et al. | Defining homogeneous regions for streamflow processes in Turkey using a K-means clustering method | |
| CN106934417A (en) | A kind of data stream self-adapting clustering method towards mixed attributes | |
| CN108957587A (en) | Igneous rockss recognition methods and device | |
| Mahjour et al. | Developing a workflow to represent fractured carbonate reservoirs for simulation models under uncertainties based on flow unit concept | |
| CN107301328A (en) | Cancer subtypes based on data stream clustering are precisely found and evolution analysis method | |
| SA521421457B1 (en) | Characterizing Depositional Features by Geologic-Based Seismic Classification | |
| CN108804635A (en) | A kind of method for measuring similarity based on Attributions selection | |
| CN104166806B (en) | A kind of clustering method of inter-well tracer test curve and device | |
| CN109271466A (en) | A kind of weather data analysis method based on hierarchical clustering Yu K mean algorithm | |
| CN112861890A (en) | Reservoir evaluation model construction method and reservoir identification method | |
| Hofmann et al. | Derivation of roof types by cluster analysis in parameter spaces of airborne laserscanner point clouds | |
| Zhao et al. | Data-driven diagenetic facies classification and well-logging identification based on machine learning methods: A case study on Xujiahe tight sandstone in Sichuan Basin |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190927 |