CN-115979800-B - Rock elasticity parameter calculation method and device based on full stress strain curve

CN115979800BCN 115979800 BCN115979800 BCN 115979800BCN-115979800-B

Abstract

The application provides a rock elasticity parameter calculation method, device, equipment and storage medium based on a full stress strain curve, which comprises the steps of carrying out filtering treatment on original full stress strain curve data, determining a data calculation interval and obtaining preprocessed data; determining constraint limits of rock physical meaning parameters to obtain Poisson's ratio processing data, adopting a locally optimal self-adaptive seed search method to conduct linear search segment optimization on the preprocessed data to obtain search method processing data, jointly determining a data range of a linear segment of a full stress strain curve according to the obtained Poisson's ratio processing data and the search method processing data, and calculating rock elasticity parameters. Noise interference is reduced by performing operations such as preprocessing and filtering on the curve. Meanwhile, the linear section is determined by the constructed Poisson ratio constraint method and the locally optimal self-adaptive seed search method, so that human selection errors are greatly reduced, and the calculation accuracy of rock elastic parameters is improved.

Inventors

ZHANG FANGNAN
SHEN HUI
LIU WEIHUA
BAI JUN
ZHOU FENG
TAN WENHUI

Assignees

中国石油化工股份有限公司
中国石油化工股份有限公司石油物探技术研究院

Dates

Publication Date: 20260512
Application Date: 20211014

Claims (8)

1. The rock elasticity parameter calculation method based on the full stress strain curve is characterized by comprising the following steps of: The method comprises the steps of S1, preprocessing data, wherein the preprocessing data comprises the steps of carrying out filtering processing on original full-stress strain curve data, determining a data calculation interval according to a filtering processing result to obtain preprocessed data, wherein the filtering processing comprises the steps of carrying out frequency spectrum analysis on waveforms of the original full-stress strain curve data to obtain analysis parameters, filtering high-frequency information of the analysis parameters through a low-pass filter, and retaining effective low-frequency information to obtain low-frequency data; S2, determining constraint limits of the petrophysical meaning parameters according to the preprocessed data to obtain Poisson ratio processed data, wherein the constraint limits of the petrophysical meaning parameters are determined by solving the circumferential axial change rate of the preprocessed data, namely, the instantaneous Poisson ratio, wherein the position where the instantaneous Poisson ratio is smaller than 0 is used as the lower boundary of the constraint limits of the petrophysical meaning parameters, and the position where the instantaneous Poisson ratio is larger than 0.5 is used as the upper boundary of the constraint limits of the petrophysical meaning parameters; S3, according to the preprocessed data, adopting a local optimal self-adaptive seed search method to conduct linear search segment optimization on the preprocessed data to obtain search method processed data; S4, determining the data range of the linear section of the full stress strain curve together according to the obtained Poisson ratio processing data and the searching method processing data; and S5, calculating the rock elasticity parameter according to the data range of the determined full stress strain curve linear section.
2. The method of claim 1, wherein determining the data calculation interval comprises detecting a maximum position of an axial and lateral extensometer of the low frequency data and deleting all data after the maximum position to obtain the data calculation interval.
3. The method according to claim 2, wherein the specific method of using the locally optimal adaptive seed search method is: Setting a search seed, namely, considering that the initial value, namely, the position of a half of the stress peak value is in a linear section of the full stress strain curve of the preprocessed data, searching in two directions from the position of the half of the stress peak value, and considering that the slope difference at two ends is smaller than the constraint value and still is in the linear section until the slope change rate exceeds the constraint value after one edge is taken one step, considering that the edge has taken the end of the linear section, starting to take the other edge to the end, and considering that the two ends at the moment are the two ends of the linear section, so as to determine the range of the linear section of the full stress strain curve of the preprocessed data, and obtaining the search method processed data.
4. The method of claim 3, wherein the determining the data range of the linear section of the full stress strain curve according to the obtained poisson's ratio processing data and the searching method processing data comprises comparing the upper and lower boundaries of the poisson's ratio processing data and the searching method processing data, and taking the value, which is close to the initial searching seed direction, between the upper and lower boundaries as two end points of the data range of the linear section of the full stress strain curve.
5. The method of claim 4, wherein the slope of the linear segment of the full stress-strain curve is the elastic modulus of the rock, and the specific method for calculating the elastic modulus of the rock comprises determining the stress value at the starting point of the linear segment And longitudinal strain Endpoint stress value And longitudinal strain Calculating the rock elastic modulus according to the formula (1), and calculating the rock elastic poisson ratio according to the formula (2): (1) (2) In the formula: -rock elastic modulus, MPa; -rock elastic poisson ratio; -stress value of the starting point of the straight line segment on the stress-axial strain relationship curve; -stress value of the end point of the straight line segment on the stress-axial strain relationship curve; Stress is A longitudinal strain value at that time; Stress is A longitudinal strain value at that time; Stress is The circumferential strain value at that time; Stress is Circumferential strain values at that time.
6. A rock elasticity parameter calculation device based on a full stress strain curve, comprising: the data preprocessing module is used for carrying out filtering processing on original full-stress strain curve data, determining a data calculation interval according to a filtering processing result to obtain preprocessed data, wherein the filtering processing comprises the steps of carrying out frequency spectrum analysis on waveforms of the original full-stress strain curve data to obtain analysis parameters, filtering high-frequency information of the analysis parameters through a low-pass filter, and retaining effective low-frequency information to obtain low-frequency data; The Poisson ratio processing data module is used for determining a constraint limit of the petrophysical meaning parameter according to the preprocessed data to obtain Poisson ratio processing data, wherein the constraint limit of the petrophysical meaning parameter is determined by solving the circumferential axial change rate of the preprocessed data, namely, the instantaneous Poisson ratio, wherein the position of the instantaneous Poisson ratio is smaller than 0 is used as the lower boundary of the constraint limit of the petrophysical meaning parameter, and the position of the instantaneous Poisson ratio is larger than 0.5 is used as the upper boundary of the constraint limit of the petrophysical meaning parameter; the search method processing data module is used for carrying out linear search segment optimization on the preprocessed data by adopting a local optimal self-adaptive seed search method according to the preprocessed data to obtain search method processing data; the data processing module of the full stress strain curve linear section jointly determines the data range of the full stress strain curve linear section according to the obtained Poisson ratio processing data and the searching method processing data; and the rock elasticity parameter calculation module is used for calculating rock elasticity parameters according to the data range of the determined full stress strain curve linear section.
7. A rock elasticity parameter calculation device based on a full stress strain curve, comprising a memory and a processor, the memory having stored thereon a computer program which, when executed by the processor, performs the rock elasticity parameter calculation method based on a full stress strain curve as claimed in any one of claims 1 to 5.
8. A storage medium storing a computer program executable by one or more processors for implementing the rock elasticity parameter calculation method based on a full stress strain curve as claimed in any one of claims 1 to 5.

Description

Rock elasticity parameter calculation method and device based on full stress strain curve Technical Field The application relates to the field of geophysical exploration, in particular to a rock elasticity parameter calculation method, device, equipment and storage medium based on a full stress strain curve. Background The invention relates to the field of rock physical experiment data processing, in particular to a calculation method for obtaining rock elasticity parameters by preprocessing full-stress strain curves based on rock mechanics uniaxial and triaxial compression deformation tests, filtering and the like, and jointly determining line elasticity segments of the stress strain curves by adopting a Poisson ratio constraint method and a locally optimal self-adaptive seed search method. Namely a rock elasticity parameter calculation method based on a full stress strain curve. The rock elasticity parameter calculation method and the rock elasticity parameter calculation device can calculate rock elasticity parameters in rock mechanical tests more accurately. The deformation, strength and failure characteristics exhibited by rock under force are known as rock mechanical properties. The basic mechanical properties of the rock are known, and mechanical tests are generally required to be carried out on the rock to obtain the mechanical parameters of the rock. The rock mechanical test is mainly to simulate the original stress state of the rock, and a special loading device is used for applying load to the rock and simultaneously testing the deformation of the rock under the action of the load. And drawing a load-deformation relation curve by using the measured load value and the corresponding deformation value. According to the change trend of the stress-strain curve, the mechanical properties of the rock can be analyzed and researched, and deformation and strength parameters reflecting the mechanical properties of the rock, such as elastic modulus, poisson ratio, stress threshold value and the like, can be calculated. In the process of implementing the present invention, the inventor finds that at least the following problems exist in the prior art: At present, the Young modulus and Poisson's ratio are usually determined according to data obtained when the peak strength is 40% -60% based on calculation of the rock elasticity parameters of the full stress strain curve. But this approach does not represent a linear elastic range well for curves of different rocks. In the rock mechanics test, noise of the sensor can interfere with full stress strain curve data, and further, errors can be caused in the obtained result. In addition to the peak strength method, there is also a method of manually picking up the elastic section of the stress strain line. The manual picking-up method is time-consuming and labor-consuming, and meanwhile, error interference caused by sensor noise cannot be well solved. For this reason, it is required to develop a rock elasticity parameter calculation method capable of solving the problem of determining the line elasticity range of the full stress strain curve and reducing the sensor disturbance error. . Disclosure of Invention Aiming at the problems, the application provides a rock elasticity parameter calculation method, a device, equipment and a storage medium based on a full stress strain curve. The application provides a rock elasticity parameter calculation method based on a full stress strain curve, which comprises the following steps: The data preprocessing comprises the steps of carrying out filtering processing on original full-stress strain curve data, and determining a data calculation interval according to a filtering processing result to obtain preprocessed data; s2, determining constraint limits of rock physical meaning parameters according to the preprocessed data to obtain Poisson ratio processing data; S3, according to the preprocessed data, adopting a local optimal self-adaptive seed search method to conduct linear search segment optimization on the preprocessed data to obtain search method processed data; S4, determining the data range of the linear section of the full stress strain curve together according to the obtained Poisson ratio processing data and the searching method processing data; and S5, calculating the rock elasticity parameter according to the data range of the determined full stress strain curve linear section. In some embodiments, the filtering process includes performing spectrum analysis on the waveform of the original full stress strain curve data to obtain analysis parameters, filtering out high frequency information of the analysis parameters by a low pass filter, and retaining effective low frequency information to obtain low frequency data. In some embodiments, the determining the data calculation interval includes detecting a maximum position of the axial and lateral extensometers of the low frequency data and deleting all data after the ma