CN114594526A - Quantitative evaluation method for inter-well plane connectivity of river facies sand body - Google Patents

Quantitative evaluation method for inter-well plane connectivity of river facies sand body Download PDF

Info

Publication number
CN114594526A
CN114594526A CN202011397661.4A CN202011397661A CN114594526A CN 114594526 A CN114594526 A CN 114594526A CN 202011397661 A CN202011397661 A CN 202011397661A CN 114594526 A CN114594526 A CN 114594526A
Authority
CN
China
Prior art keywords
well
sand
parameter
connectivity
inter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011397661.4A
Other languages
Chinese (zh)
Other versions
CN114594526B (en
Inventor
张海娜
杜玉山
李健
翟亮
徐鹏晔
刘丽
姜书荣
王志伟
胡晨彬
崔晴晴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN202011397661.4A priority Critical patent/CN114594526B/en
Publication of CN114594526A publication Critical patent/CN114594526A/en
Application granted granted Critical
Publication of CN114594526B publication Critical patent/CN114594526B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/40—Protecting water resources
    • Y02A20/402—River restoration

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Geophysics (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention provides a quantitative evaluation method for inter-well plane connectivity of river facies sand bodies, which comprises the following steps: step 1, identifying and dividing sand bodies of single river channels, and establishing a well-to-sample detailed data table; step 2, setting 5 key parameters closely related to connectivity, and calculating 5 key parameters of a well pair; step 3, carrying out data processing and weight calculation on the 5 key parameters of the well pair; and 4, calculating a quantitative evaluation parameter M of the inter-well connectivity, and carrying out quantitative evaluation on the inter-well plane connectivity of the river facies sand body. Compared with the prior art, the method for quantitatively evaluating the inter-well plane connectivity of the river-phase sand body has better operability, innovativeness and practicability, is beneficial to popularization, and provides a simple, quick and feasible method for quantitatively evaluating the inter-well plane connectivity.

Description

Quantitative evaluation method for inter-well plane connectivity of river facies sand body
Technical Field
The invention relates to the technical field of oil and gas field development, in particular to a quantitative evaluation method for inter-well plane connectivity of river facies sand bodies.
Background
In the middle and later stages of development of oil fields, research on connectivity of reservoirs among wells plays an important decision-making guidance role in prediction of residual oil among wells, adjustment of well pattern series, well position arrangement of a newly drilled and adjusted well, formulation of production and perforation schemes, production allocation and injection allocation of injection and production wells and the like. The existing plane connectivity evaluation method of the interwell reservoir stratum has more qualitative evaluation methods, has a method for dividing connectivity levels according to a certain parameter or a plurality of parameters in the aspect of geology, and has a method for analyzing according to a sediment type and a sand body superposition mode, and also has a pressure drop analysis method, an injection and production dynamic data correlation analysis method, an interference well testing method, a tracer agent and other methods in the aspect of oil reservoir engineering. The quantitative evaluation method is few, and mainly comprises machine learning, numerical simulation and the like. The machine learning method needs to establish a prediction model, while the numerical simulation method needs to establish a geological model, the accuracy requirement on the established model is high, the prediction result of the method depends on the established three-dimensional geological model and the well logging interpretation result of the pore permeability parameter, and the method needs long research time and high implementation cost.
In the application No.: 201711257300.8, to a method and apparatus for determining the connectivity of sand bodies. The method comprises the following steps: respectively determining membership degree relations between each standard index and each specified connectivity grade; setting a weight matrix corresponding to the standard indexes according to the standard indexes, and respectively determining the target weight value of each standard index in the weight matrix; and determining a target index and a target index parameter value according to the first geological parameter information and the second geological parameter information, and determining sand body connectivity between the first single sand body and the second single sand body according to the target index parameter value, the target weight value and the membership degree relation of the standard index. According to the method, the sand-to-ground ratio, the interlayer density, the sand permeability and the porosity of two sand bodies are considered, three levels (level I, level II and level III) of connectivity are divided according to the product of the parameters, a standard index is determined for each level, and the connectivity level evaluation between the two single sand bodies in a work area is realized. There is a lack of quantitative assessment of the connectivity between well pairs within a sand body from well to well.
In the application No.: 201811002234.4, to a method, apparatus and system for determining sand body connectivity. The method comprises the steps of scanning a reservoir outcrop analogue body of a target work area by laser to obtain a digital outcrop section of the reservoir outcrop analogue body; extracting sand body thickness data from the digital outcrop section; and determining the sand body connectivity of the target work area according to the extracted sand body thickness data and the logging sand body thickness data. The method utilizes the outcrop sand body communication relation to restrict the underground sand body hooking mode to determine the sand body connectivity of the target work area, and is a qualitative evaluation method of the connectivity.
In the application No.: 201910658689.X relates to a method and a device for comparing underground single-river sand bodies among wells. The method comprises the following steps: acquiring detection information of a tracer injected from an injection well, which is detected in a monitoring well; calculating connectivity parameters of sand bodies among wells according to the detection information, wherein the connectivity parameters represent the sand body connectivity performance between the injection well and the monitoring well; determining the sand body communication condition between the injection well and the monitoring well according to the connectivity parameters; and identifying and comparing single channel sand bodies from the composite channel according to the sand body communication condition, wherein the composite channel is formed by splicing and superposing a plurality of single channel sand bodies. This patent determines sand body communication between an injection well and a monitoring well by monitoring tracer velocity. The disadvantage of this method is that there are fewer wells that can be tested for tracers in actual production, especially in offshore fields.
In the application No.: 201710376308.X, relates to a sand body connectivity evaluation method and device. The method comprises the following steps: respectively establishing a transverse connectivity sand body sample, a longitudinal connectivity sand body sample and an internal connectivity sand body sample of a work area, and dividing each type of sand body sample into a training sample and a test sample; training the training samples of each type of sand body samples by using a preset machine learning algorithm, and establishing a corresponding sand body connectivity prediction model; optimizing the corresponding sand body connectivity prediction model according to the test sample of each type of sand body sample so as to enable the prediction result of the corresponding sand body connectivity prediction model to meet the preset condition; and according to the optimized sand body connectivity prediction model, performing sand body connectivity evaluation on sand body data corresponding to sand bodies to be identified in the work area to obtain an evaluation result. The method comprises the steps of training samples of each type of sand body samples through a preset machine learning algorithm, establishing a corresponding sand body connectivity prediction model, and evaluating sand bodies to be recognized by applying the prediction model. The method has the defects of multiple machine learning algorithms, complex process and unfavorable operation.
There is a need for a fast and efficient method for calculating a quantitative connectivity evaluation value between wells. For river facies sand bodies with large work area and more drilled wells, a method which is rapid, efficient and simple to operate is lacked. Therefore, a novel quantitative evaluation method for inter-well plane connectivity of river facies sand bodies is invented to solve the technical problems.
Disclosure of Invention
The invention aims to provide a river facies sand body inter-well plane connectivity quantitative evaluation method which realizes the quantification of inter-well plane connectivity by optimizing evaluation parameters and determining an evaluation method and provides directly-usable data for well pattern adjustment and production allocation and injection allocation of injection and production wells.
The object of the invention can be achieved by the following technical measures: the method for quantitatively evaluating the inter-well plane connectivity of the river facies sand bodies comprises the following steps:
step 1, identifying and dividing sand bodies of single river channels, and establishing a well-to-sample detailed data table;
step 2, setting 5 key parameters closely related to connectivity, and calculating 5 key parameters of a well pair;
step 3, carrying out data processing and weight calculation on the 5 key parameters of the well pair;
and 4, calculating a quantitative evaluation parameter M of the inter-well connectivity, and carrying out quantitative evaluation on the inter-well plane connectivity of the river facies sand body.
The object of the invention can also be achieved by the following technical measures:
before step 1, carrying out conventional geological research work in a target research area according to well completion information, completing stratum contrast division, well logging interpretation and sedimentary facies research work, and determining the sedimentary type on the basis; and (4) if the deposition type is the meandering stream deposition flow, entering the step 1, and if not, ending the flow.
In step 1, on the basis of the deposition of the meandering stream, single-channel division is performed according to geological research results, the specific single-channel sand body in which the sand body of each well is positioned is determined, and the source direction of the single channel is determined.
In step 1, a sample detailed data table of each well pair is established, wherein the sample detailed data table comprises the sand body width of each well, the relative depth of the sand body top surface from the top surface of the stratum, the well distance of the well pair, the included angle between the well pair connecting line and the source direction and the permeability of each well.
In step 2, 5 key parameters X1, X2, X3, X4 and X5 which are closely related to connectivity are set, wherein the width and the well spacing ratio, the depth difference and the sand thickness ratio, the direction deviation degree, the permeability difference and the well spacing ratio and the sand thickness difference and the well spacing ratio are respectively set, and the key parameter values are calculated for each well pair.
In step 2, let
Figure BDA0002814095130000031
Wherein: l is1The width of the sand body of the well A;
L2the sand body width of the well B;
d is the well spacing of the A and B wells;
x1 is calculating the sand width of two wells, summing, and dividing by the well spacing; the geological meaning represented by X1 is to judge whether two sand bodies are isolated; if X1<1, the sand is isolated and not connected, the connection evaluation parameter M of the well pair is 0, and the process ends, otherwise, the process continues.
In step 2, let
Figure BDA0002814095130000041
Wherein: h1 is the depth from the top surface of the sand body of the well A to the top surface of the small layer;
H2the depth from the top surface of the sand body of the well B to the top surface of the small layer;
h1,2the smaller sand thickness of the two sands;
x2 is the relative depth difference of the sand tops of two wells, and the absolute value is taken, and the absolute value is divided by the thickness of the sand body with smaller relative depth of the sand top; the geological meaning represented by X2 is to judge that the smaller the difference in depth between two sand bodies in the longitudinal direction, the better the connectivity.
In step 2, let X3 be cos α
Wherein: alpha is the included angle between the well pair connecting line and the source direction;
x3 is the deviation degree of the well-to-well connecting line and the object source; the geological meaning represented by X3 is that the closer the direction of the injection-production well pair is to the direction of the source, the better the connectivity in the single-channel sand body.
In step 2, let
Figure BDA0002814095130000042
Wherein: k1Is a well AThe sand permeability of (a);
K2the sand permeability of the well B;
d is the well spacing of the A and B wells;
x4 is the difference of the permeability of the two wells and the absolute value is taken, and the absolute value is divided by the well spacing; the geological meaning represented by X4 is that in a single river channel, the smaller the difference between two well permeabilities and the smaller the well spacing, the better the connectivity.
In step 2, let
Figure BDA0002814095130000043
Wherein: h is1The sand body thickness of the well A;
h2the sand body thickness of the well B;
d is the well spacing of the A and B wells;
x5 is the difference of the sand thickness of two wells and the absolute value is taken, and the absolute value is divided by the well spacing; the geological meaning represented by X5 is that in a single river channel, the smaller the difference of the thicknesses of sand bodies encountered by two wells and the smaller the well spacing, the better the connectivity.
In step 3, performing data processing on the 5 parameter results calculated in step 2; firstly, respectively carrying out normalization processing on calculated parameters, normalizing data to be within 0-1, wherein in the normalization process, the smaller the 3 parameters of X2, X4 and X5 are, the better the connectivity is, and the negative correlation is;
Figure BDA0002814095130000051
wherein: e1iNormalized for parameter X1 ith well pair;
X1ithe ith well pair value of the 1 st key parameter X1 in the 5 key parameters set for the previous time;
X1minis the minimum value of the parameter X1;
X1maxis the maximum value of the parameter X1;
i is the well pair count;
Figure BDA0002814095130000052
wherein: e2iNormalized for parameter X2 ith well pair;
X2ithe ith well pair value of the 2 nd key parameter X2 in the 5 key parameters set for the previous time;
X2minis the minimum value of the parameter X2;
X2maxis the maximum value of the parameter X2;
i is the well pair count;
E3i＝X3i
wherein: e3iNormalized for parameter X3 ith well pair;
X3ithe ith well pair value of the 3 rd key parameter X3 in the 5 key parameters set for the previous time;
i is the well pair count;
Figure BDA0002814095130000053
wherein: e4iNormalized for parameter X4 ith well pair;
X4ithe ith well pair value of the 4 th key parameter X4 in the 5 key parameters set for the previous time;
X4minis the minimum value of the parameter X4;
X4maxis the maximum value of the parameter X4;
i is the well pair count;
Figure BDA0002814095130000061
wherein: e5iNormalized for parameter X5 ith well pair;
X5ithe ith well pair value of the 5 th key parameter X5 in the 5 key parameters set in the previous step;
X5minis the minimum value of the parameter X5;
X5maxis the maximum value of the parameter X5;
i is the well pair count.
In step 3, the weight value of each parameter is further calculated, and the weights of 5 parameters are integrated into 1; the formula is as follows: the index weight is the sum of a single index score model coefficient/five index comprehensive score models; and calculating the weight value of each parameter.
In step 4, the quantitative inter-well connectivity evaluation parameter M is defined as:
Figure BDA0002814095130000062
wherein: ei-the normalized parameter value;
zi-weight;
N—1～5；
and calculating the communication evaluation value between each well pair through the formula, and finishing the quantitative communication evaluation of all the well pairs.
The invention discloses a quantitative evaluation method for inter-well plane connectivity of river facies sand bodies, and relates to a quantitative evaluation method for inter-well plane connectivity in river facies sand bodies. The method has clear technical thought and simple application, reduces the complexity of establishing a geological model or a prediction model for quantitative evaluation, has better operability compared with the prior method, has innovation and practicability, and is beneficial to popularization. A simple, quick and feasible method is provided for the quantitative evaluation of the plane connectivity between wells. The method successfully obtains the quantitative connectivity of the well pair and provides basic data for rapid production allocation and injection allocation.
Drawings
FIG. 1 is a flow chart of a method for quantitatively evaluating the inter-well plane connectivity of a river facies sand body according to an embodiment of the present invention;
FIG. 2 is a plan view of two single channel sands identified in an embodiment of the present invention;
FIG. 3 is a schematic representation of an interwell connectivity evaluation parameter in an embodiment of the present invention.
Detailed Description
It is to be understood that the following detailed description is exemplary and is intended to provide further explanation of the invention as claimed. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of the stated features, steps, operations, and/or combinations thereof, unless the context clearly indicates otherwise.
The method for quantitatively evaluating the inter-well plane connectivity of the river-phase sand body comprises the following steps of:
step 1, identifying and dividing sand bodies of single river channels, and establishing a well-to-sample detailed data table;
on the basis of the deposition of the meandering stream, single river channel division is carried out according to geological research results, the specific single river channel sand body in which the sand body of each well is positioned is determined, and the source direction of the single river channel is determined.
And establishing a sample detailed data table of each well pair, wherein the sample detailed data table comprises the sand body width of each well, the relative depth of the top surface of the sand body from the top surface of the stratum, the well distance of the well pair, the included angle between the well pair connecting line and the source direction and the permeability of each well.
In one embodiment, prior to step 1, in the target research area, based on the completed well data, performing conventional geological research, performing stratigraphic correlation and partitioning, well logging interpretation, and sedimentary facies research, and determining the sedimentary type. And (4) if the deposition type is the meandering stream deposition process, entering the step 1, and if not, ending the process.
Step 2, calculating 5 key parameters of the well pair;
5 key parameters (X1, X2, X3, X4 and X5) which are closely related to connectivity are established, namely the width and well spacing ratio, the depth difference and sand thickness ratio, the direction deviation degree, the permeability difference and well spacing ratio and the sand thickness difference and well spacing ratio. The 5 key parameters not only reflect the thickness and the development direction of the sand body, but also reflect the physical property change.
The key parameter values are calculated for each well pair.
Is provided with
Figure BDA0002814095130000081
Wherein: l is1The width of the sand body of the well A;
L2the sand body width of the well B;
d is the well spacing of the A and B wells;
x1 is the calculation of sand width for two wells and summing, then dividing by well spacing. The geological meaning of X1 is to determine whether two sand bodies are isolated. If X1<1, the sand is isolated and not connected, the connection evaluation parameter M of the well pair is 0, and the process ends, otherwise, the step continues.
Is provided with
Figure BDA0002814095130000082
Wherein: h1 is the depth from the top surface of the sand body of the well A to the top surface of the small layer;
H2the depth from the top surface of the sand body of the well B to the top surface of the small layer;
h1,2the smaller sand thickness of the two sands;
and X2 is the relative depth difference of the sand tops of the two wells, and the absolute value is taken, and the absolute value is divided by the thickness of the sand body with the smaller relative depth of the sand top. The geological meaning represented by X2 is to judge that the smaller the difference in depth between two sand bodies in the longitudinal direction, the better the connectivity.
Let X3 be cos alpha
Wherein: alpha is the included angle between the well pair connecting line and the source direction;
x3 is the deviation degree of the well-to-well connecting line and the object source; the geological meaning represented by X3 is that the closer the direction of the injection-production well pair is to the direction of the source, the better the connectivity in the single-channel sand body.
Is provided with
Figure BDA0002814095130000083
Wherein: k1The sand permeability of the well A is shown;
K2the sand permeability of the well B;
d is the well spacing of the A and B wells;
x4 is the difference between the two well permeabilities and the absolute value is taken and divided by the well spacing. The geological meaning represented by X4 is that in a single river channel, the smaller the difference between two well permeabilities and the smaller the well spacing, the better the connectivity.
Is provided with
Figure BDA0002814095130000091
Wherein: h is1The sand body thickness of the well A;
h2the sand body thickness of the well B;
d is the well spacing of the A and B wells;
x5 is the difference between the sand thickness of two wells and the absolute value is taken, and divided by the well spacing. The geological meaning represented by X5 is that in a single river channel, the smaller the difference of the thicknesses of sand bodies encountered by two wells and the smaller the well spacing, the better the connectivity.
For clarity of explanation of the representative meanings of the parameters, see FIG. 3.
Step 3, data processing and weight calculation;
for convenient evaluation, the 5 parameter results calculated in step 2 are subjected to data processing. Firstly, respectively carrying out normalization processing on the calculated parameters, normalizing the data to be within 0-1, and processing by using a conventional normalization method. During the normalization process, the 3 parameters of X2, X4 and X5 are noted to be that the smaller the connectivity is, the better the connectivity is, and the correlation is negative.
Figure BDA0002814095130000092
Wherein: e1iNormalized for parameter X1 ith well pair;
X1ithe ith well pair value of the 1 st key parameter X1 in the 5 key parameters set for the previous time;
X1minis the minimum value of the parameter X1;
X1maxis the maximum value of the parameter X1;
i is the well pair count;
Figure BDA0002814095130000093
wherein: e2iNormalized for parameter X2 ith well pair;
X2ithe ith well pair value of the 2 nd key parameter X2 in the 5 key parameters set for the previous time;
X2minis the minimum value of the parameter X2;
X2maxis the maximum value of the parameter X2;
i is the well pair count;
E3i＝X3i
wherein: e3iNormalized values for parameter X3 for the ith well pair;
X3ithe ith well pair value of the 3 rd key parameter X3 in the 5 key parameters set in the previous step;
i is the well pair count;
Figure BDA0002814095130000101
wherein: e4iNormalized for parameter X4 ith well pair;
X4ithe ith well pair value of the 4 th key parameter X4 in the 5 key parameters set for the previous time;
X4minis the minimum value of the parameter X4;
X4maxis the maximum value of the parameter X4;
i is the well pair count;
Figure BDA0002814095130000102
wherein: e5iNormalized for parameter X5 ith well pair;
X5ithe ith well pair value of the 5 th key parameter X5 in the 5 key parameters set as the previous value;
X5minis the minimum value of the parameter X5;
X5maxis the maximum value of the parameter X5;
i is the well pair count.
Further, the weight value of each parameter is calculated, and the weights of 5 parameters are integrated into 1. The formula is as follows: and the index weight is the sum of the single index score model coefficient and the five index comprehensive score models. And calculating the weight value of each parameter.
And 4, calculating a comprehensive evaluation parameter M.
The quantitative evaluation parameter M for the well connectivity is defined as:
Figure BDA0002814095130000111
wherein: ei-normalized parameter value
Zi weight
N—1～5
The communication evaluation value between each well pair can be calculated through the formula. Thus, the quantitative connectivity evaluation of all the well pairs is completed.
In the specific example 1 applying the invention, the reservoir of the certain oil field at the victory sea is the sediment of the meandering river at the upper section of the building, the sand body communication is complex, the average porosity is 33 percent, and the permeability is 1504 multiplied by 10-3μm2And the development is put into 1993, the number of completed wells is 742, and the water content is 86.7 percent at present. The well spacing is 300-400 m. Chengjiang 11 wells were selected for illustration.
As shown in fig. 1, fig. 1 is a flow chart of the quantitative evaluation method for inter-well plane connectivity of river facies sand body according to the present invention.
And 101, identifying and dividing sand bodies of the single river channel, and establishing a well-to-sample detailed data table.
Ng5 according to the sedimentary phase characteristics of the well region6Two single river course sand bodies can be identified on the layerThe two rivers are basically north-east-south-west, and the source of the substance is from north-west, as shown in fig. 2. Each well forms a different well pair with its surrounding neighbors, and 5 parameters are available for each well pair. As shown in Table 1, the first well pair, well 1 is 26A-3, well 2 is 26A-1, well 1(26A-3) corresponds to parameters of sand width L1 of 214m, relative depth H1 of sand top surface to small layer top surface of 7.2m, sand thickness H1 of 3.6m, and permeability K1 of 1427 × 10-3μm2The parameters corresponding to the well 2(26A-1) are that the sand body width L2 is 388m, the relative depth H2 of the top surface of the sand body from the top surface of the small layer is 8.3m, the sand body thickness H2 is 5.6m, and the permeability K2 is 1861 multiplied by 10-3μm2. The well spacing between the 26A-2 well and the 26A-1 well is 553m, and the included angle between the connecting line of the 26A-2 well and the 26A-1 well and the source direction is 76 degrees. By analogy, a sample detail data table can be built for all well pairs in the zone.
TABLE 1 well-to-sample detailed data Table
Figure BDA0002814095130000112
Figure BDA0002814095130000121
Figure BDA0002814095130000131
Step 102, 5 key parameters of a well pair are calculated.
On the basis of step 101, the parameter X1, well pair 1, is first calculated
Figure BDA0002814095130000132
Since the calculation result X1 is 0.54, which is smaller than 1, the process ends and it can be directly concluded that the overall evaluation value of the well pair is 0. And by analogy, continuously calculating X1 parameters of other well pairs, and finding that X1 parameters of 2 nd, 3 th and 5 th well pairs are all smaller than 1, so that other parameters of the 4 well pairs are not required to be calculated, and a comprehensive evaluation value is directly obtained to be 0. For theWell pairs for which X1 is greater than 1 continue to calculate parameters X2, X3, X4, and X5. The calculation results are written in table 2.
Table 2 comprehensive evaluation value calculation table
Figure BDA0002814095130000133
Figure BDA0002814095130000141
And 103, processing data and calculating parameter weight.
For evaluation, the 5 parameters calculated in step 102 were normalized by a conventional method, and the normalized data are E1, E2, E3, E4 and E5 in table 2. Where E3 ═ X3, since X3 is a number between 0 and 1, it can be used without normalization.
The weights can be obtained by well-established algorithms, such as analytic hierarchy process, fuzzy analytic hierarchy process, coefficient of variation process, entropy weight method, and expert scoring. In this example, principal component analysis was applied to obtain a parameter X1 weight of 0.30, a parameter X2 weight of 0.36, a parameter X3 weight of 0.15, a parameter X4 weight of 0.10, and a parameter X5 weight of 0.09, as shown in table 3.
TABLE 3 parameter weight table
Parameter(s) Each index weight
X1 0.30
X2 0.36
X3 0.15
X4 0.10
X5 0.09
Step 104, calculating a comprehensive evaluation parameter M
From the normalized weights of E1-E5 in Table 2 and the parameters in Table 3, it is easy to calculate the overall evaluation value for each well pair, i.e., obtain quantitative values for the connectivity of each well pair, which are in the last column of Table 2, the well pairs with the highest overall evaluation value among the well pairs listed are the CB22E-4 and CB22A-5 wells, and the M value is 0.71, indicating that the well pair connectivity is the best among the well pairs listed in the Table. There are 4 well pairs with an M value of 0, indicating that the 4 well pairs are not connected. The pair of wells with the worst connectivity in the table are 11K-5 and 11K-2.
The method is applied to a plurality of blocks of the field of the Chengqi field, the coincidence rate with the field of the Chengqi field is more than 90%, and the coincidence rate is high. The method and the parameters are suitable for evaluating the inter-well plane connectivity of the river facies sand body and have reference significance for evaluating the connectivity of other sediment types of sand bodies.
In specific embodiment 2 to which the present invention is applied, the following steps are included:
step 101, identifying and dividing sand bodies of single river channel, and establishing a well-to-sample detailed data table
Ng5 according to the sedimentary phase characteristics of the well region5Two single river channel sand bodies can be identified in the layer, the two rivers are basically in the north-east-south-west direction, the source of the matter is from the north-west direction, as shown in table 4, the first well pair is 11M-1 for the well 1, 11H-4 for the well 2, 469M for the sand body width L1 for the well 1(11M-1) corresponding parameters, 5.2M for the relative depth H1 of the top surface of the sand body from the top surface of the small layer, and 6 for the sand body thickness H14m, a permeability K1 of 1427X 10-3μm2The parameters corresponding to the well 2(11H-4) are that the sand body width L2 is 384m, the relative depth H2 of the top surface of the sand body from the top surface of the small layer is 2.0m, the sand body thickness H2 is 5.6m, and the permeability K2 is 1861 multiplied by 10-3μm2. The well spacing between the 11M-1 well and the 11H-4 well is 230M, and the included angle between the connecting line of the 11M-1 well and the 11H-4 well and the source direction is 24 degrees. By analogy, a sample detail table 4 can be built for all well pairs within a zone.
Table 4 well to sample detail data table
Figure BDA0002814095130000151
Figure BDA0002814095130000161
Step 102, 5 key parameters of a well pair are calculated.
On the basis of step 101, the parameters X1 are first calculated, X1 for well pair 1 is (469+384)/(2 × 230), the result X1 is 1.86, which is greater than 1, so the flow continues, and so on, the X1 parameters for other well pairs continue to be calculated, and the parameters X2, X3, X4 and X5 continue to be calculated for well pairs for which X1 is greater than 1. The calculation results are written in table 5.
Table 5 comprehensive evaluation value calculation table
Figure BDA0002814095130000171
And 103, processing data and calculating parameter weight.
For evaluation, the 5 parameters calculated in step 102 were normalized by a conventional method, and the normalized data are E1, E2, E3, E4 and E5 in table 4. Where E3 ═ X3, X3 can be used as is without normalization, since it is a number between 0 and 1. The weights were taken as in example 1.
Step 104, calculating a comprehensive evaluation parameter M
From the normalized weights of the parameters in table 4, E1 to E5, and table 5, the overall evaluation value for each well pair can be easily calculated.
Finally, it should be noted that: although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that changes may be made in the embodiments and/or equivalents thereof without departing from the spirit and scope of the invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
In addition to the technical features described in the specification, the technology is known to those skilled in the art.

Claims (13)

1. The method for quantitatively evaluating the inter-well plane connectivity of the river facies sand bodies is characterized by comprising the following steps of:
step 1, identifying and dividing sand bodies of single river channels, and establishing a well-to-sample detailed data table;
step 2, setting 5 key parameters closely related to connectivity, and calculating 5 key parameters of a well pair;
step 3, carrying out data processing and weight calculation on the 5 key parameters of the well pair;
and 4, calculating a quantitative evaluation parameter M of the inter-well connectivity, and carrying out quantitative evaluation on the inter-well plane connectivity of the river facies sand body.
2. The method for quantitatively evaluating the inter-well plane connectivity of river facies sand bodies according to claim 1, further comprising, before the step 1, performing conventional geological research work in the target research area according to completed well data, completing stratigraphic comparison and division, well logging interpretation and sedimentary facies research work, and determining the sedimentary type based on the conventional geological research work; and (4) if the deposition type is the meandering stream deposition process, entering the step 1, and if not, ending the process.
3. The method for quantitatively evaluating the inter-well plane connectivity of river facies sand bodies according to claim 1, wherein in step 1, on the basis of the deposition of the meandering river, single-channel division is performed according to geological research results, and the specific single-channel sand body in which the sand body of each well is located is determined, and the source direction of the single channel is determined.
4. The method for quantitatively evaluating the inter-well plane connectivity of river facies sand bodies according to claim 1, wherein in step 1, a detailed data table of samples of each well pair is established, which comprises the width of the sand body where each well is located, the relative depth of the top surface of the sand body from the top surface of the stratum, the well distance of the well pair, the included angle between the connecting line of the well pair and the source direction, and the permeability of each well.
5. The method for quantitatively evaluating the plane connectivity between wells of river-phase sand bodies according to claim 1, wherein in step 2, 5 key parameters X1, X2, X3, X4 and X5 closely related to the connectivity are established, wherein the key parameter values are respectively calculated for each well pair according to the width and the well spacing ratio, the depth difference and the sand thickness ratio, the direction deviation degree, the permeability difference and the well spacing ratio and the sand thickness difference and the well spacing ratio.
6. The method for quantitatively evaluating the inter-well plane connectivity of river facies sand bodies as claimed in claim 5, wherein in step 2, the method is characterized in that
Figure FDA0002814095120000011
Wherein: l is a radical of an alcohol1The width of the sand body of the well A is obtained;
L2the sand body width of the well B;
d is the well spacing of the A and B wells;
x1 is calculating the sand width of two wells, summing, and dividing by the well spacing; the geological meaning represented by X1 is to judge whether two sand bodies are isolated; if X1<1, the sand is isolated and not connected, the connection evaluation parameter M of the well pair is 0, and the process ends, otherwise, the process continues.
7. The method for quantitatively evaluating the inter-well plane connectivity of river facies sand bodies as claimed in claim 5, wherein in step 2, the method is characterized in that
Figure FDA0002814095120000021
Wherein: h1 is the depth from the top surface of the sand body of the well A to the top surface of the small layer;
H2the depth from the top surface of the sand body of the well B to the top surface of the small layer;
h1,2the smaller sand thickness of the two sands;
x2 is the relative depth difference of the sand tops of two wells, and the absolute value is taken, and the absolute value is divided by the thickness of the sand body with smaller relative depth of the sand top; the geological meaning represented by X2 is to judge that the smaller the difference in depth between two sand bodies in the longitudinal direction, the better the connectivity.
8. The method for quantitatively evaluating the inter-well plane connectivity of river-phase sand bodies according to claim 5, wherein in step 2, X3 ═ cos α is defined
Wherein: alpha is the included angle between the well pair connecting line and the source direction;
x3 is the deviation degree of the well-to-well connecting line and the object source; the geological meaning represented by X3 is that the closer the direction of the injection-production well pair is to the direction of the source, the better the connectivity in the single-channel sand body.
9. The method for quantitatively evaluating the inter-well plane connectivity of fluvial facies sand bodies according to claim 5, wherein in step 2, the method is characterized in that
Figure FDA0002814095120000022
Wherein: k1The sand permeability of the well A is shown;
K2the sand permeability of the well B;
d is the well spacing of the A and B wells;
x4 is the difference of the permeability of the two wells and the absolute value is taken, and the absolute value is divided by the well spacing; the geological meaning represented by X4 is that in a single river channel, the smaller the difference between two well permeabilities and the smaller the well spacing, the better the connectivity.
10. The method for quantitatively evaluating the inter-well plane connectivity of river facies sand bodies as claimed in claim 5, wherein in step 2, the method is characterized in that
Figure FDA0002814095120000031
Wherein: h is1The sand body thickness of the well A;
h2the sand thickness of the well B;
d is the well spacing of the A and B wells;
x5 is the difference of the sand thickness of two wells and the absolute value is taken, and the absolute value is divided by the well spacing; the geological meaning represented by X5 is that in a single river channel, the smaller the difference of the thicknesses of sand bodies encountered by two well drilling and the smaller the well spacing, the better the connectivity.
11. The method for quantitatively evaluating the inter-well plane connectivity of river facies sand bodies according to claim 5, wherein in step 3, the 5 parameter results calculated in step 2 are subjected to data processing; firstly, respectively carrying out normalization processing on calculated parameters, normalizing data to be within 0-1, wherein in the normalization process, the smaller the 3 parameters of X2, X4 and X5 are, the better the connectivity is, and the negative correlation is;
Figure FDA0002814095120000032
wherein: e1iNormalized for parameter X1 ith well pair;
X1ithe ith well pair value of the 1 st key parameter X1 in the 5 key parameters set for the previous time;
X1minis the minimum value of the parameter X1;
X1maxis the maximum value of the parameter X1;
i is the well pair count;
Figure FDA0002814095120000033
wherein: e2iNormalized for parameter X2 ith well pair;
X2ithe ith well pair value of the 2 nd key parameter X2 in the 5 key parameters set for the previous time;
X2minis the minimum value of the parameter X2;
X2maxis the maximum value of the parameter X2;
i is the well pair count;
E3i＝X3i
wherein: e3iNormalized for parameter X3 ith well pair;
X3ithe ith well pair value of the 3 rd key parameter X3 in the 5 key parameters set for the previous time;
i is the well pair count;
Figure FDA0002814095120000041
wherein: e4iNormalized for parameter X4 ith well pair;
X4ithe ith well pair value of the 4 th key parameter X4 in the 5 key parameters set for the previous time;
X4minis the minimum value of the parameter X4;
X4maxis the maximum value of the parameter X4;
i is the well pair count;
Figure FDA0002814095120000042
wherein: e5iNormalized for parameter X5 ith well pair;
X5ifor the 5 th of the 5 key parameters setThe ith well pair value of the key parameter X5;
X5minis the minimum value of the parameter X5;
X5maxis the maximum value of the parameter X5;
i is the well pair count.
12. The method for quantitatively evaluating the inter-well plane connectivity of river facies sand bodies according to claim 11, wherein in step 3, the weight value of each parameter is further calculated, and the weight values of 5 parameters are integrated into 1; the formula is as follows: the index weight is the sum of a single index score model coefficient/five index comprehensive score models; and calculating the weight value of each parameter.
13. The method for quantitatively evaluating the inter-well plane connectivity of river facies sand bodies according to claim 1, wherein in step 4, the inter-well connectivity quantitative evaluation parameter M is defined as:
Figure FDA0002814095120000043
wherein: ei-the normalized parameter value;
zi-weight;
N—1～5；
and calculating the communication evaluation value between each well pair through the formula to finish the quantitative connectivity evaluation of all the well pairs.
CN202011397661.4A 2020-12-03 2020-12-03 Quantitative evaluation method for plane connectivity between river-phase sand bodies Active CN114594526B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011397661.4A CN114594526B (en) 2020-12-03 2020-12-03 Quantitative evaluation method for plane connectivity between river-phase sand bodies

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011397661.4A CN114594526B (en) 2020-12-03 2020-12-03 Quantitative evaluation method for plane connectivity between river-phase sand bodies

Publications (2)

Publication Number Publication Date
CN114594526A true CN114594526A (en) 2022-06-07
CN114594526B CN114594526B (en) 2023-07-11

Family

ID=81812474

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011397661.4A Active CN114594526B (en) 2020-12-03 2020-12-03 Quantitative evaluation method for plane connectivity between river-phase sand bodies

Country Status (1)

Country Link
CN (1) CN114594526B (en)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160145994A1 (en) * 2014-11-20 2016-05-26 Petrochina Company Limited Evaluation Method and Evaluation Device for Water Breakthrough Risk of Production Wells in Aquifer Drive Gas Reservoirs
CN107330578A (en) * 2017-05-25 2017-11-07 中国石油天然气股份有限公司 Sand body Connectivity Evaluation method and device
CN107861171A (en) * 2017-11-08 2018-03-30 中国石油大港油田勘探开发研究院 A kind of terrestrial lake basin sedimentary micro flat distribution map preparation method
CN108090656A (en) * 2017-12-04 2018-05-29 中国石油天然气股份有限公司 A kind of method and device of definite sand body connectedness
US20190212460A1 (en) * 2018-01-08 2019-07-11 Dagang Oil Field Of Cnpc Method for secondary exploration of old oil area in fault subsidence basin
CN110376661A (en) * 2019-07-19 2019-10-25 中国地质大学（北京） The control methods of single river channel sand and device between a kind of missile silo
CN110412649A (en) * 2019-06-17 2019-11-05 中国石油化工股份有限公司 A kind of recognition methods of list phase distributary channel
CN111426813A (en) * 2019-01-10 2020-07-17 中国石油天然气股份有限公司 Method and device for determining sand body communication quality and storage medium
AU2020101482A4 (en) * 2020-06-05 2020-09-03 China University Of Geosciences, Beijing Calculation method for quantitatively establishing sea level change and transgression-regression curves
CN111766635A (en) * 2019-04-02 2020-10-13 中国石油天然气股份有限公司 Sand body communication degree analysis method and system

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160145994A1 (en) * 2014-11-20 2016-05-26 Petrochina Company Limited Evaluation Method and Evaluation Device for Water Breakthrough Risk of Production Wells in Aquifer Drive Gas Reservoirs
CN107330578A (en) * 2017-05-25 2017-11-07 中国石油天然气股份有限公司 Sand body Connectivity Evaluation method and device
CN107861171A (en) * 2017-11-08 2018-03-30 中国石油大港油田勘探开发研究院 A kind of terrestrial lake basin sedimentary micro flat distribution map preparation method
CN108090656A (en) * 2017-12-04 2018-05-29 中国石油天然气股份有限公司 A kind of method and device of definite sand body connectedness
US20190212460A1 (en) * 2018-01-08 2019-07-11 Dagang Oil Field Of Cnpc Method for secondary exploration of old oil area in fault subsidence basin
CN111426813A (en) * 2019-01-10 2020-07-17 中国石油天然气股份有限公司 Method and device for determining sand body communication quality and storage medium
CN111766635A (en) * 2019-04-02 2020-10-13 中国石油天然气股份有限公司 Sand body communication degree analysis method and system
CN110412649A (en) * 2019-06-17 2019-11-05 中国石油化工股份有限公司 A kind of recognition methods of list phase distributary channel
CN110376661A (en) * 2019-07-19 2019-10-25 中国地质大学（北京） The control methods of single river channel sand and device between a kind of missile silo
AU2020101482A4 (en) * 2020-06-05 2020-09-03 China University Of Geosciences, Beijing Calculation method for quantitatively establishing sea level change and transgression-regression curves

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
栗亮等: "基于模糊数学方法的砂体连通性评价", 中州煤炭, no. 11 *

Also Published As

Publication number Publication date
CN114594526B (en) 2023-07-11

Similar Documents

Publication Publication Date Title
CN108301825B (en) High-temperature high-pressure reservoir pore pressure analysis method
CN109214026A (en) A kind of shale gas horizontal well initial productivity prediction technique
CN103616731B (en) Altered Volcanic Rocks at some Place Effective Reservoirs defining method and device in a kind of oil-gas exploration
CN112213797A (en) Method for predicting reservoir productivity by using logging-phase combined post-stack seismic attributes
CN107728227B (en) A kind of well pattern maturation zone quick discrimination buried channel sand body method
CN110162851A (en) A kind of data calibration method of cable formation testing pumping numerical simulation and its process
CA3032547A1 (en) System and method for using geological analysis for the designing of stimulation operations
CN110646847B (en) Method for identifying low-order fault breakpoint position based on well logging data of close well pattern area
CN114066084B (en) Method and system for predicting phase permeation curve based on machine learning
CN116910638A (en) Dominant seepage channel intelligent identification method based on big data
CN105626056A (en) Method and device for determining three-dimensional stratum pore pressure of target area
CN109239783B (en) Method for accepting or rejecting earthquake information in well earthquake combination
CN114066666A (en) Method for analyzing connectivity among wells through injection-production profile monitoring data
CN112302594A (en) Method, device and equipment for determining connection structure of heterogeneous water-drive reservoir
CN114594526A (en) Quantitative evaluation method for inter-well plane connectivity of river facies sand body
CN111830558B (en) Fracture zone engraving method
CN117131971A (en) Reservoir dominant seepage channel prediction method based on Xgboost algorithm
CN113325468B (en) Reservoir fracture-cavity distribution range prediction method and device
CN111155980A (en) Water flow dominant channel identification method and device
CN114635682A (en) Multi-layer commingled production capacity prediction method based on shaft multi-permeability coupling
CN110566196B (en) Reservoir connectivity analysis method
RU2119583C1 (en) Method for monitoring development of oil deposits
CN104462792B (en) Log data lithological stratum numerical value reduction method
CN112031740B (en) Shale gas geological dessert evaluation method and device based on stratum elements
CN108875109B (en) Method and system for predicting abnormal formation pressure

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant