CN112765785B - Multi-scale rock mechanical layer well logging division method - Google Patents
Multi-scale rock mechanical layer well logging division method Download PDFInfo
- Publication number
- CN112765785B CN112765785B CN202011634468.8A CN202011634468A CN112765785B CN 112765785 B CN112765785 B CN 112765785B CN 202011634468 A CN202011634468 A CN 202011634468A CN 112765785 B CN112765785 B CN 112765785B
- Authority
- CN
- China
- Prior art keywords
- rock
- rock mechanical
- mechanical
- logging
- layers
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Devices Affording Protection Of Roads Or Walls For Sound Insulation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to the field of reservoir geomechanics, in particular to a multi-scale rock mechanics layer logging division method. Constructing a rock mechanical layer discrimination index by calculating dynamic and static rock mechanical parameters and adopting an equal frequency conversion method; dividing a rock mechanical layer with a single scale and calculating the density of the rock through the clustering analysis of the mechanical parameters of the logging cast points; and changing the threshold value, circularly performing logging cast-point rock mechanical parameter clustering analysis, and completing conventional logging division of rock mechanical layers with different scales. The invention provides a multi-scale rock mechanics layer logging division method which has high practical value and has practical value for reservoir geomechanical modeling and stress field numerical simulation.
Description
Technical Field
The invention relates to the field of reservoir geomechanics, in particular to a multi-scale rock mechanics layer logging division method.
Background
The rock mechanical layer refers to a set of rock layers with consistent rock mechanical properties or similar rock mechanical behaviors, and is related to deformation characteristics of the rock after stress, and the rock mechanical layers can be further divided into discrete intervals. The rock mechanical behavior is dimensional; rock mechanical behavior may vary at different scales, and current research limits rock mechanical layer research to a single scale. How to establish a reasonable discrimination index and realize division of rock mechanical layers is a practical problem faced by reservoir geomechanical research.
Disclosure of Invention
The invention aims to solve the problems and provides a multi-scale rock mechanical layer logging dividing method which can quantitatively divide rock mechanical layers with different scales based on conventional logging information.
The technical scheme of the invention is as follows: a multi-scale rock mechanical layer logging dividing method comprises the following specific steps:
firstly, calculating mechanical parameters of logging rocks;
calculating rock mechanical parameters by using logging information, and converting the dynamic and static rock mechanical parameters to obtain rock static mechanical parameters, wherein the rock mechanical parameters comprise Young modulus, Poisson ratio, internal friction angle, cohesion and compressive strength; according to the data such as the sound wave time difference, the rock density, the mud content, the rock porosity and the like, the rock mechanics parameters are calculated, and the related calculation formula is as follows:
Sc＝Ed[0.008Vsh+0.0045(1-Vsh)] (4)
in formulae (1) to (5), EdDynamic Young's modulus, MPa; mu.sdThe Poisson ratio is dynamic and dimensionless; c is cohesion, MPa; scCompressive strength, MPa; vshThe mud content is percentage, and has no dimension; rhobIs rock density, kg/m3；ΔtpAnd Δ tsRespectively longitudinal wave time difference and transverse wave time difference, mu s/ft;
is the internal friction angle, °; phi is the porosity,%;
secondly, constructing a rock mechanical layer discrimination index;
converting n groups of rock mechanical parameters of different sizes and different units of a single well into units and scales corresponding to the Young modulus of the rock by adopting an equal frequency conversion method, averaging the n groups of rock mechanical parameters corresponding to each converted logging point to construct a new logging mechanical parameter discrimination index G (E1, E2.. Em), wherein m is the number of the logging points;
the method for equal frequency conversion is that the mechanical parameters of different types of rocks of a single well are marked with serial numbers 1, 2 and 3 from small to large in sequence; according to the numerical value corresponding to the marking serial number of the Young modulus of the rock, respectively endowing the Young modulus of the rock corresponding to the mechanical parameters of different types of rocks, and unifying the mechanical parameters of different types of rocks to the unit and the scale corresponding to the Young modulus of the rock;
thirdly, logging, casting rock mechanical parameter clustering analysis;
setting a threshold T in unit GPa; firstly, starting from an E1 data point, preliminarily setting the data point as a rock mechanical layer, dividing E2 and E1 into the same rock mechanical layer when the absolute value of the difference between the discrimination indexes G of E2 and E1 is smaller than T for the next data point E2, otherwise, dividing the two rock mechanical layers into different rock mechanical layers, and finishing the rock mechanical layer division of the E2 logging data point; for the next data point E3, when the absolute value of the discrimination index G between the E3 and the rock mechanical layer on the upper portion is smaller than T, dividing the E3 and the rock mechanical layer on the upper portion into the same rock mechanical layer, otherwise, dividing the two rock mechanical layers into different rock mechanical layers, and completing the rock mechanical layer division of the E3 logging data point; so as to finish the division of the rock mechanical layer of E1 and E2.
Fourthly, dividing the mechanical layers of the rocks with different scales and calculating the density;
setting the difference iteration step length d of the rock mechanical layers with different scales, wherein the unit is as follows: GPa; and changing the threshold T (T + d), completing the division of the rock mechanical layers of the scale according to the clustering analysis of the rock mechanical parameters of the logging throw point in the third step, sequentially changing the threshold T, dividing the rock mechanical layers of different scales, and simultaneously calculating the density of the rock mechanical layers of different scales by using the number of the rock mechanical layers in the unit length in the vertical direction.
The invention has the beneficial effects that: constructing a rock mechanical layer discrimination index by calculating dynamic and static rock mechanical parameters and adopting an equal frequency conversion method; and dividing rock mechanical layers with different scales and calculating the density of the rock through the clustering analysis of the mechanical parameters of the logging cast-point rock. The invention provides a multi-scale rock mechanics layer well logging division method which has high practical value and low prediction cost, and the prediction result has practical value for reservoir geomechanical modeling and stress field numerical simulation.
Drawings
FIG. 1 is a flow chart of a multi-scale rock mechanics layer well logging partitioning method.
FIG. 2X 1 well different scale rock mechanics layer division scheme.
Detailed Description
The following description of the embodiments of the present invention refers to the accompanying drawings:
the invention takes the Chinese Qiongsoutheast basin X1 well as an example to illustrate the specific implementation process of the invention.
Firstly, calculating rock mechanical parameters by using formulas (1) to (5);
secondly, constructing a rock mechanical layer discrimination index G;
setting T to be 3GPa, and dividing the rock mechanical layer with the size by adopting a logging cast point method;
and d set in the fourth step is 3GPa, dividing the rock mechanical layers with different scales, and calculating corresponding density (figure 2).
The present invention has been described above by way of example, but the present invention is not limited to the above-described specific embodiments, and any modification or variation made based on the present invention is within the scope of the present invention as claimed.
Claims (1)
1. A multi-scale rock mechanical layer well logging division method comprises the following implementation steps:
firstly, calculating mechanical parameters of logging rocks;
calculating rock mechanical parameters by using logging information, and converting the dynamic and static rock mechanical parameters to obtain rock static mechanical parameters, wherein the rock mechanical parameters comprise Young modulus, Poisson ratio, internal friction angle, cohesion and compressive strength; according to the data such as the sound wave time difference, the rock density, the mud content, the rock porosity and the like, the rock mechanics parameters are calculated, and the related calculation formula is as follows:
Sc＝Ed[0.008Vsh+0.0045(1-Vsh)] (4)
in formulae (1) to (5), EdDynamic Young's modulus, MPa; mu.sdThe dynamic Poisson ratio is dimensionless; c is cohesion, MPa; scCompressive strength, MPa; vshThe mud content is the percentage of the mud without dimension; rhobIs rock density, kg/m3；ΔtpAnd Δ tsRespectively longitudinal wave time difference and transverse wave time difference, mu s/ft;
is the internal friction angle, °; phi is the logging porosity,%;
secondly, constructing a rock mechanical layer discrimination index;
converting n groups of rock mechanical parameters of different sizes and different units of a single well into units and scales corresponding to the Young modulus of the rock by adopting an equal frequency conversion method, and averaging the n groups of rock mechanical parameters corresponding to each converted logging point to construct a new logging mechanical parameter discrimination index G (E1, E2.. Em), wherein m is the number of the logging points;
the method for equal frequency conversion is that the mechanical parameters of different types of rocks of a single well are marked with serial numbers 1, 2 and 3 from small to large in sequence; according to the numerical value corresponding to the marking serial number of the Young modulus of the rock, respectively endowing the Young modulus of the rock corresponding to the mechanical parameters of different types of rocks, and unifying the mechanical parameters of different types of rocks to the unit and the scale corresponding to the Young modulus of the rock;
thirdly, logging, casting rock mechanical parameter clustering analysis;
setting a threshold T in unit GPa; firstly, starting from an E1 data point, preliminarily setting the data point as a rock mechanical layer, dividing E2 and E1 into the same rock mechanical layer when the absolute value of the difference between the discrimination indexes G of E2 and E1 is smaller than T for the next data point E2, otherwise, dividing the two rock mechanical layers into different rock mechanical layers, and finishing the rock mechanical layer division of the E2 logging data point; for the next data point E3, when the absolute value of the discrimination index G between the E3 and the rock mechanical layer on the upper portion is smaller than T, dividing the E3 and the rock mechanical layer on the upper portion into the same rock mechanical layer, otherwise, dividing the two rock mechanical layers into different rock mechanical layers, and completing the rock mechanical layer division of the E3 logging data point; so as to finish the division of the rock mechanical layer of E1 and E2.
Fourthly, dividing mechanical layers of rocks and stones with different scales and calculating density;
setting the difference iteration step length d of the rock mechanical layers with different scales, wherein the unit is as follows: GPa; and changing the threshold T to be T + d, completing the division of the rock mechanical layers of the scale according to the third step of logging cast rock mechanical parameter cluster analysis, sequentially changing the threshold T, realizing the division of the rock mechanical layers of different scales, and simultaneously calculating the density of the rock mechanical layers of different scales by utilizing the number of the rock mechanical layers in the unit length in the vertical direction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011634468.8A CN112765785B (en) | 2020-12-31 | 2020-12-31 | Multi-scale rock mechanical layer well logging division method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011634468.8A CN112765785B (en) | 2020-12-31 | 2020-12-31 | Multi-scale rock mechanical layer well logging division method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112765785A CN112765785A (en) | 2021-05-07 |
| CN112765785B true CN112765785B (en) | 2022-05-06 |
Family
ID=75699729
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011634468.8A Active CN112765785B (en) | 2020-12-31 | 2020-12-31 | Multi-scale rock mechanical layer well logging division method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112765785B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5416697A (en) * | 1992-07-31 | 1995-05-16 | Chevron Research And Technology Company | Method for determining rock mechanical properties using electrical log data |
| CN110779795A (en) * | 2019-11-04 | 2020-02-11 | 中国石油大学(华东) | Method for determining size of geomechanical modeling grid unit of fractured reservoir |
| CN111425193A (en) * | 2020-01-21 | 2020-07-17 | 东北石油大学 | Reservoir compressibility evaluation method based on clustering analysis logging rock physical facies division |
-
2020
- 2020-12-31 CN CN202011634468.8A patent/CN112765785B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5416697A (en) * | 1992-07-31 | 1995-05-16 | Chevron Research And Technology Company | Method for determining rock mechanical properties using electrical log data |
| CN110779795A (en) * | 2019-11-04 | 2020-02-11 | 中国石油大学(华东) | Method for determining size of geomechanical modeling grid unit of fractured reservoir |
| CN111425193A (en) * | 2020-01-21 | 2020-07-17 | 东北石油大学 | Reservoir compressibility evaluation method based on clustering analysis logging rock physical facies division |
Non-Patent Citations (3)
| Title |
|---|
| Effects of perforation fluid movement on downhole packer with shock loads;Qiao Deng .etc;《Journal of Petroleum Science and Engineering》;20201215;全文 * |
| 中国南方海相与陆相页岩裂缝发育特征及主控因素对比;王濡岳 等;《石油与天然气地质》;20180831;第39卷(第4期);631-640 * |
| 超剥带残留地层测井划分方法研究;陈钢花 等;《测井技术》;20150630;第39卷(第3期);368-372 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112765785A (en) | 2021-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110644980B (en) | Comprehensive classification evaluation method for ultra-low permeability oil reservoir | |
| CN110779795B (en) | Method for determining size of geomechanical modeling grid unit of fractured reservoir | |
| CN112502701B (en) | Classification evaluation method for comprehensive geology-engineering of low-permeability reservoir | |
| CN113901681B (en) | Three-dimensional compressibility evaluation method for dual desserts of shale gas reservoir in whole life cycle | |
| CN112394397B (en) | Shale gas reservoir three-dimensional rock mechanical parameter field modeling method | |
| CN105221141A (en) | A kind of mud shale brittleness index Forecasting Methodology | |
| CN112304754B (en) | Shale brittleness logging evaluation method considering diagenesis and pressure change | |
| CN111832227B (en) | Shale gas saturation determination method, device and equipment based on deep learning | |
| CN111027818B (en) | Shale oil classification evaluation method | |
| CN108104805B (en) | Method for acquiring transverse wave time difference profile by using well logging and core data | |
| CN112765785B (en) | Multi-scale rock mechanical layer well logging division method | |
| CN112282742B (en) | Prediction method for shale oil high-quality reservoir | |
| Sornette et al. | Linking short-timescale deformation to long-timescale tectonics | |
| CN110847887A (en) | Method for identifying and evaluating cracks of fine-grain sedimentary continental facies shale | |
| CN115586569A (en) | Stratum horizontal ground stress calculation method based on data driving under theoretical model constraint | |
| CN113419278B (en) | Well-seismic joint multi-target simultaneous inversion method based on state space model and support vector regression | |
| CN109975189A (en) | Porous sandstone Reservoir Productivity Prediction Method and device | |
| CN111101914B (en) | Horizontal well fracturing segment cluster optimization method and equipment | |
| CN106291690B (en) | Accumulate pattern the Forward Modeling in the stratum of orbital period tuning | |
| CN111159869A (en) | Automatic calculation method for quantitative parameters between sandstone uranium deposit and ancient river channel | |
| CN110989034A (en) | Method for inverting logging transverse wave time difference by regression-fractal interpolation method | |
| CN114109374A (en) | Method and device for determining position of target window of shale gas reservoir | |
| Wang et al. | Impacts of the base-level cycle on pore structure of mouth bar sand bars: A case study of the Paleogene Kongdian formation, Bohai bay basin, China | |
| CN113051261B (en) | Dictionary learning-based high-resolution sound velocity profile sparse coding and storage method | |
| CN117077446B (en) | Unconventional reservoir comprehensive brittleness index evaluation method and application system |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20220422 Address after: 430074 No. 388 Lu Lu, Hongshan District, Hubei, Wuhan Applicant after: CHINA University OF GEOSCIENCES (WUHAN CITY) Address before: 266580 No. 66 Changjiang West Road, Huangdao District, Qingdao, Shandong. Applicant before: CHINA University OF PETROLEUM (EAST CHINA) |
|
| TA01 | Transfer of patent application right |