A fitting modeling method and system for full-band rock physical measurement data

Abstract

This invention relates to the field of petroleum and natural gas geophysics, and discloses a method and system for fitting and modeling full-band rock physical measurement data. The method includes: acquiring full-band rock velocity data and rock parameter data; determining the number of attenuation mechanisms and the frequency band range of each attenuation mechanism based on the full-band rock velocity data; obtaining the relaxation time data of each attenuation mechanism based on the frequency band range and rock parameter data; obtaining a fitting expression based on the relaxation time data of all attenuation mechanisms, and constructing a full-band rock physical measurement data fitting model. Based on full-band rock physical measurement data and combined with linear elastic body theory, this invention achieves rapid fitting and modeling of full-band rock velocity attenuation and dispersion characteristics, quickly obtaining a concise and efficient full-band rock physical measurement data fitting model that can uniformly simulate the velocity attenuation and dispersion characteristics of rocks under multiple attenuation mechanisms across the entire frequency band.

Description

Technical Field

[0001] This invention relates to the field of petroleum and natural gas geophysics, and in particular to a fitting modeling method and system for full-band rock physical measurement data. Background Technology

[0002] Rock physics, a crucial branch of geophysics, serves as a vital bridge for transforming geophysical data into geological parameters for reservoirs, fluids, and even oil and gas deposits. Existing rock physics measurement equipment falls into two categories: one involves stress-strain measurements conducted within the seismic frequency band, typically the low-frequency range of 5-100 Hz, whose rock physics characteristics best reflect the realities of seismic acquisition; the other utilizes ultrasonic equipment for direct rock physics velocity measurements, generally within the 10 Hz range. 6 Hz. Currently, the industry generally conducts two measurements simultaneously when performing measurements at different frequency bands. Rocks exhibit different elastic characteristics in different measurement frequency bands, requiring a unified model that can characterize the full-band attenuation and dispersion characteristics of rocks.

[0003] Existing theoretical models describing velocity decay and dispersion in porous rocks can be broadly categorized into two types: First, the BISQ model, designed for high frequencies, which can simulate the macroscopic and microscopic jet flow characteristics of rocks in the high-frequency band; second, mesoscale models, such as the White and Johnson models, designed for seismic frequency bands. These models can effectively describe the decay and dispersion phenomena caused by compressive flow in the seismic frequency pore fluid band. However, neither of these models can simultaneously characterize the velocity decay and dispersion characteristics of rocks across a wide frequency range. Furthermore, these modeling methods are computationally inefficient in numerical simulations due to their complex construction, making them unsuitable for efficient industrial applications.

[0004] Therefore, there is an urgent need for a rapid fitting and modeling method for full-band rock physical measurement data. Summary of the Invention

[0005] This invention provides a fitting modeling method and system for full-band rock physics measurement data to address the shortcomings of existing technologies that cannot simultaneously characterize the velocity attenuation and dispersion characteristics of rocks across large frequency bands.

[0006] This invention provides a fitting and modeling method for full-band rock physical measurement data, comprising:

[0007] Acquire rock velocity data and rock parameter data across the entire frequency band, where the entire frequency band is 10. -1 -10 6 The frequency band range of Hz;

[0008] Based on the rock velocity data across the entire frequency band, determine the number of attenuation mechanisms and the frequency band range data for each attenuation mechanism;

[0009] Based on the frequency band range and rock parameter data of each attenuation mechanism, the relaxation time data of each attenuation mechanism is obtained;

[0010] Based on the relaxation time data of all attenuation mechanisms, a fitting expression is obtained, and a fitting model for full-band rock physics measurement data is constructed.

[0011] In one embodiment, the rock density parameter data includes: rock porosity data, rock density data, fluid composition data, fluid density data, fluid saturation data, and dry rock shear modulus.

[0012] In one embodiment, the frequency band range data includes frequency band range data and frequency band range velocity data, wherein the frequency band range data includes the maximum frequency value and the minimum frequency value of the frequency band range, and the frequency band range velocity data includes the maximum velocity value and the minimum velocity value of the frequency band range.

[0013] In one implementation, obtaining the relaxation time data for each attenuation mechanism based on its frequency band range and rock parameter data includes:

[0014] Based on the frequency data of the frequency band range of each attenuation mechanism, the characteristic angular frequency of each attenuation mechanism is obtained;

[0015] Based on the frequency band velocity data and rock density data of each attenuation mechanism, the extreme value of the rock bulk modulus of each attenuation mechanism is obtained.

[0016] The relaxation time data for each decay mechanism are obtained based on the characteristic angular frequency and the extreme value of the rock bulk modulus for each decay mechanism.

[0017] In one implementation, obtaining the characteristic angular frequency of each attenuation mechanism based on the frequency data of the bandwidth range of each attenuation mechanism includes:

[0018] Based on the frequency data of the bandwidth of each attenuation mechanism, the characteristic angular frequency of each attenuation mechanism is obtained through the first expression;

[0019] The first expression is:

[0020]

[0021] In the first expression, ω represents the characteristic angular frequency. max ω represents the maximum frequency value within the bandwidth. min This indicates the minimum frequency value within the bandwidth.

[0022] In one embodiment, the extreme values ​​of rock bulk modulus include low-frequency extreme values ​​of rock bulk modulus and high-frequency extreme values ​​of rock bulk modulus.

[0023] In one implementation, obtaining the rock bulk modulus extreme value for each attenuation mechanism based on the frequency band velocity data and rock density data for each attenuation mechanism includes:

[0024] Based on the frequency band velocity data and rock density data of each attenuation mechanism, the low-frequency extreme value of the rock bulk modulus of each attenuation mechanism is obtained through the second expression.

[0025] Based on the frequency band velocity data and rock density data of each attenuation mechanism, the high-frequency extreme value of the rock bulk modulus of each attenuation mechanism is obtained through the third expression.

[0026] The second expression is:

[0027]

[0028] In the second expression, K min This represents the low-frequency extreme value of the bulk modulus of rock, v min ρ represents the minimum velocity within the frequency band, and ρ represents the rock density data.

[0029] The third expression is:

[0030]

[0031] In the third expression, K max This represents the low-frequency extreme value of the bulk modulus of rock, v max ρ represents the maximum velocity within the frequency band, and ρ represents the rock density data.

[0032] In one implementation, obtaining the relaxation time data for each attenuation mechanism based on its characteristic angular frequency and the extreme value of the rock bulk modulus includes:

[0033] Based on the extreme values ​​of the bulk modulus of rock for each attenuation mechanism and the shear modulus of dry rock, the plane wave modulus of each attenuation mechanism is obtained.

[0034] Based on the plane wave modulus and characteristic angular frequency of each attenuation mechanism, the expression for the complex plane wave modulus of rock for each attenuation mechanism is obtained;

[0035] Based on the expression for the complex plane wave modulus of rock for each attenuation mechanism, the expression for the high-frequency limit value of the bulk modulus of rock for each attenuation mechanism is obtained.

[0036] Based on the characteristic angular frequency of each attenuation mechanism, the relaxation time data of each attenuation mechanism is obtained and the expression is solved.

[0037] The relaxation time data for each attenuation mechanism is obtained by solving the expressions for the high-frequency limit value of the rock bulk modulus and the relaxation time data of the characteristic angular frequency.

[0038] In one embodiment, the expression for the plane wave modulus includes an expression for the maximum plane wave modulus and an expression for the minimum plane wave modulus, wherein the expression for the maximum plane wave modulus is:

[0039]

[0040] In the expression for the maximum modulus of a plane wave, E max K represents the maximum plane wave modulus. max μ represents the high-frequency extreme value of the bulk modulus of rock. m Indicates the shear modulus of dry rock;

[0041] The expression for the minimum plane wave modulus is:

[0042]

[0043] In the expression for the minimum plane wave modulus, E min K represents the minimum plane wave modulus. min μ represents the low-frequency extreme value of the bulk modulus of rock. m This represents the shear modulus of dry rock.

[0044] In one implementation, the expression for the complex plane wave modulus of rock is:

[0045]

[0046] In the expression for the complex plane wave modulus of rock, E ω E represents the complex plane wave modulus of rock. min τ represents the minimum plane wave modulus, ω represents the characteristic angular frequency, and τ represents the minimum plane wave modulus. ε τ σ This represents a set of relaxation times to be determined.

[0047] In one implementation scheme, the expression for the high-frequency limit value of the rock bulk modulus is:

[0048]

[0049] In the expression for the high-frequency limit value of the bulk modulus of rock, E min τ represents the minimum plane wave modulus. ε τ σ This represents a set of relaxation times to be determined.

[0050] In one implementation, the expression for solving the relaxation time data of the characteristic angular frequency is:

[0051]

[0052] In the expression for solving the relaxation time data of the characteristic angular frequency, ω represents the characteristic angular frequency, and τ ε τ σ This represents a set of relaxation times to be determined.

[0053] In one implementation, the step of obtaining a fitting expression based on the relaxation time data of all attenuation mechanisms and constructing a full-band rock physics measurement data fitting model includes:

[0054] Based on the characteristic angular frequency and relaxation time data of each attenuation mechanism, the quality factor of each attenuation mechanism is obtained.

[0055] Based on the characteristic angular frequency, relaxation time, and preset elastic modulus of each attenuation mechanism, the longitudinal wave dynamic modulus of each attenuation mechanism is obtained.

[0056] Based on the rock density data and preset elastic modulus of each attenuation mechanism, the longitudinal wave complex velocity of each attenuation mechanism is obtained.

[0057] Based on the P-wave complex velocity of each attenuation mechanism, the P-wave velocity fitting curve of each attenuation mechanism is obtained, and the P-wave velocity fitting curves of all attenuation mechanisms constitute a fitting expression.

[0058] Based on the quality factors of all attenuation mechanisms and the P-wave velocity fitting curves, a fitting model for full-band rock physics measurement data is constructed.

[0059] In one implementation, obtaining the quality factor of each attenuation mechanism based on its characteristic angular frequency and relaxation time data includes:

[0060] Based on the characteristic angular frequency and relaxation time data of each decay mechanism, the quality factor of each decay mechanism is obtained through the fourth expression.

[0061] The fourth expression is:

[0062]

[0063] In the fourth expression, Q represents the quality factor, ω represents the characteristic angular frequency, and τ′ ε τ′ σ This represents a set of relaxation times after the solution is obtained.

[0064] In one implementation, obtaining the longitudinal wave dynamic modulus of each attenuation mechanism based on its characteristic angular frequency, relaxation time, and preset elastic modulus includes:

[0065] Based on the characteristic angular frequency, relaxation time, and preset elastic modulus of each attenuation mechanism, the longitudinal wave dynamic modulus of each attenuation mechanism is obtained through the fifth expression.

[0066] The fifth expression is:

[0067]

[0068] In the fifth expression, E (ω) The dynamic modulus of the longitudinal wave is represented by i, the imaginary unit is represented by L, the number of attenuation mechanisms is represented by ω, and the characteristic angular frequency is represented by τ′. ε τ′ σ This represents a set of relaxation times after the solution is obtained.

[0069] In one implementation, obtaining the P-wave complex velocity of each attenuation mechanism based on the rock density data and a preset elastic modulus of each attenuation mechanism includes:

[0070] Based on the rock density data and preset elastic modulus of each attenuation mechanism, the longitudinal wave complex velocity of each attenuation mechanism is obtained through the sixth expression.

[0071] The sixth expression is:

[0072]

[0073] In the sixth expression, V (ω) E represents the complex velocity of the longitudinal wave. R This represents the preset elastic modulus, and ρ represents the rock density data.

[0074] In one implementation, obtaining the P-wave velocity fitting curve for each attenuation mechanism based on the P-wave complex velocity of each attenuation mechanism includes:

[0075] Based on the complex P-wave velocity of each attenuation mechanism, the P-wave velocity fitting curve for each attenuation mechanism is obtained through the seventh expression.

[0076] The seventh expression is:

[0077]

[0078] In the seventh expression, V p V represents the longitudinal wave velocity fitting curve. (ω) This represents the complex velocity of the longitudinal wave.

[0079] This invention also provides a fitting and modeling system for full-band rock physical measurement data, comprising:

[0080] The data acquisition module is used to acquire rock velocity data and rock parameter data across the entire frequency band.

[0081] The first data processing module is used to: determine the number of attenuation mechanisms and the frequency band range data of each attenuation mechanism based on the rock velocity data of the entire frequency band;

[0082] The second data processing module is used to: obtain the relaxation time data of each attenuation mechanism based on the frequency band range and rock parameter data of each attenuation mechanism;

[0083] The fitting modeling module is used to: obtain fitting expressions based on the relaxation time data of all decay mechanisms, and construct a fitting model for full-band rock physics measurement data.

[0084] The present invention also provides an electronic device, including a processor and a memory storing a computer program, wherein the processor executes the computer program to implement the fitting modeling method for full-band rock physical measurement data described above.

[0085] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the fitting and modeling method for full-band rock physical measurement data as described above.

[0086] The present invention also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer is able to execute any of the above-described fitting modeling methods for full-band rock physical measurement data.

[0087] This invention provides a fitting modeling method and system for full-band rock physics measurement data. Based on full-band rock physics measurement data and combined with linear elastic body theory, it enables rapid fitting modeling of the velocity attenuation and dispersion characteristics of rocks across the entire frequency band. Based on this method, a concise and efficient fitting model for full-band rock physics measurement data can be quickly obtained. This model can uniformly simulate the velocity attenuation and dispersion characteristics of rocks under various attenuation mechanisms across the entire frequency band. It has high numerical simulation efficiency and high fitting accuracy, which can meet the needs of efficient industrial applications and also contribute to the research of oil and gas geophysics. Attached Figure Description

[0088] To more clearly illustrate the technical solutions in this invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0089] Figure 1This is a flowchart illustrating a fitting and modeling method for full-band rock physical measurement data provided by the present invention.

[0090] Figure 2 The graph shows the measured rock velocity data as a function of frequency.

[0091] Figure 3 The fitting results of modeling under four decay mechanisms are shown in the figure.

[0092] Figure 4 The figure shows the fitting effect of modeling under one decay mechanism.

[0093] Figure 5 This is a schematic diagram of the structure of a fitting and modeling system for full-band rock physical measurement data provided by the present invention.

[0094] Figure 6 This is a schematic diagram of the structure of the electronic device provided by the present invention. Detailed Implementation

[0095] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, embodiments of this invention, and should not be construed as limiting the invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention. In the description of this invention, it should be understood that the terminology used is for descriptive purposes only and should not be construed as indicating or implying relative importance.

[0096] The following is combined with Figures 1-6 This invention describes the fitting and modeling method and system for full-band rock physics measurement data.

[0097] Figure 1 This is a flowchart illustrating the fitting and modeling method for full-band rock physical measurement data provided by the present invention. (Refer to...) Figure 1 The present invention provides a fitting and modeling method for full-band rock physical measurement data, which may include:

[0098] Step S110: Obtain rock velocity data and rock parameter data across the entire frequency band, where the entire frequency band is 10. -1 -10 6 The frequency band range of Hz;

[0099] Step S120: Based on the rock velocity data across the entire frequency band, determine the number of attenuation mechanisms and the frequency band range data for each attenuation mechanism;

[0100] Step S130: Based on the frequency band range and rock parameter data of each attenuation mechanism, obtain the relaxation time data of each attenuation mechanism;

[0101] Step S140: Based on the relaxation time data of all attenuation mechanisms, obtain the fitting expression and construct a fitting model for full-band rock physics measurement data.

[0102] It should be noted that the execution subject of the fitting and modeling method for full-band rock physical measurement data provided by the present invention can be any network-side device / terminal-side device that meets the technical requirements, such as a fitting and modeling device for full-band rock physical measurement data.

[0103] like Figure 2 As shown, the continuous sample points on the left represent low- to medium-frequency rock velocity data, while the individual sample points on the right represent high-frequency rock velocity data. The measured velocity of the rock sample will fluctuate within a certain range according to a certain pattern, which is called dispersion. The main reason for this is that the pore fluid is pressure-driven and flows at different scales. Step S110 can obtain the full-frequency velocity data of the rock sample (also known as full-frequency velocity sample point data) through high-frequency and low-frequency rock physics experiments, and measure basic rock parameters such as rock porosity, rock density, fluid composition, fluid density, fluid saturation, and dry rock shear modulus.

[0104] The velocity dispersion characteristics of rocks are controlled by multiple mechanisms at multiple scales, often exhibiting multi-stage dispersion phenomena. Step S120 can determine the appropriate number of attenuation mechanisms and the frequency band range data of each attenuation mechanism based on experimentally measured rock velocity data, according to the rock velocity variation characteristics and simulation needs. Specifically, the variation law of rock velocity data can be divided into multiple attenuation mechanisms. Each attenuation mechanism is characterized by a process in which the rock velocity value changes from relatively stable to linearly increasing, and then back to relatively stable within a certain frequency range. The frequency band range data includes frequency data and velocity data. The frequency data includes the maximum and minimum frequency values ​​within the frequency band range, and the velocity data includes the maximum and minimum velocity values ​​within the frequency band range. For example, the frequency band range measured in this embodiment is 10. -1 -10 6 Within the frequency band, there are 5 relatively stable speed steps. The speed change between each two steps can be simulated by an attenuation mechanism. Therefore, in this embodiment, it is preferred to use a model containing 4 sets of attenuation mechanisms to carry out the modeling work, that is, the number of attenuation mechanisms is determined to be 4.

[0105] In one embodiment, step S130 may include:

[0106] Based on the frequency data of the bandwidth of each attenuation mechanism, the characteristic angular frequency of each attenuation mechanism is obtained through the first expression, where the first expression is:

[0107]

[0108] In the first expression, ω represents the characteristic angular frequency. max ω represents the maximum frequency value within the bandwidth. min Indicates the minimum frequency value within the bandwidth;

[0109] Based on the frequency band velocity data and rock density data for each attenuation mechanism, the low-frequency extreme value of the rock bulk modulus for each attenuation mechanism is obtained through the second expression, where the second expression is:

[0110]

[0111] In the second expression, K min This represents the low-frequency extreme value of the bulk modulus of rock, v min ρ represents the minimum velocity within the frequency band, and ρ represents the rock density data.

[0112] Based on the frequency band velocity data and rock density data for each attenuation mechanism, the high-frequency extreme values ​​of the rock bulk modulus for each attenuation mechanism are obtained through the third expression, where the third expression is:

[0113]

[0114] In the third expression, K max This represents the low-frequency extreme value of the bulk modulus of rock, v max ρ represents the maximum velocity within the frequency band, and ρ represents the rock density data.

[0115] The relaxation time data for each decay mechanism are obtained based on the characteristic angular frequency and the extreme value of the rock bulk modulus for each decay mechanism.

[0116] Furthermore, rock density data can be obtained both experimentally and by calculation, and the calculation formula is as follows:

[0117]

[0118] In the formula, ρ s ρ represents the density of the solid skeleton of rock. f Represents the equivalent density of the fluid. Representing rock porosity, when it is determined from fluid composition data that two fluids are present, the equivalent fluid density can be obtained using the following formula:

[0119] ρ f =S1ρ f1+S2ρ f2 ,

[0120] Where, ρ f1 ρ f2 S1 and S2 represent the densities of the two fluids, respectively, and S1 and S2 represent the saturation levels of the two fluids, respectively.

[0121] In one embodiment, obtaining the relaxation time data for each attenuation mechanism based on the characteristic angular frequency and the extreme value of the rock bulk modulus for each attenuation mechanism may include:

[0122] Based on the extreme values ​​of the bulk modulus of rock for each attenuation mechanism and the shear modulus of dry rock, the plane wave modulus of each attenuation mechanism is obtained. The shear modulus of dry rock does not change significantly with frequency and can be obtained through rock physics experiments. The expression for the plane wave modulus includes expressions for the maximum and minimum values. The expression for the maximum plane wave modulus is:

[0123]

[0124] In the expression for the maximum modulus of a plane wave, E max K represents the maximum plane wave modulus. max μ represents the high-frequency extreme value of the bulk modulus of rock. m This represents the shear modulus of dry rock.

[0125] The expression for the minimum plane wave modulus is:

[0126]

[0127] In the expression for the minimum plane wave modulus, E min K represents the minimum plane wave modulus. min μ represents the low-frequency extreme value of the bulk modulus of rock. m Indicates the shear modulus of dry rock;

[0128] Based on the plane wave modulus and characteristic angular frequency of each attenuation mechanism, the expression for the complex plane wave modulus of the rock for each attenuation mechanism is obtained. Specifically, according to the linear volume model theory, the expression for the complex plane wave modulus of the rock for each relaxation mechanism is:

[0129]

[0130] In the expression for the complex plane wave modulus of rock, E ω E represents the complex plane wave modulus of rock. min τ represents the minimum plane wave modulus, ω represents the characteristic angular frequency, and τ represents the minimum plane wave modulus. ε τ σ This represents a set of relaxation times to be determined.

[0131] Based on the expression for the complex plane wave modulus of rock for each attenuation mechanism, the expression for the high-frequency limit value of the bulk modulus of rock for each attenuation mechanism is obtained, where, when ω→ω max At that time, the expression for the high-frequency limit value of the bulk modulus of rock can be approximated as:

[0132]

[0133] In the expression for the high-frequency limit value of the bulk modulus of rock, E min τ represents the minimum plane wave modulus. ε τ σ This represents a set of relaxation times to be determined.

[0134] Based on the characteristic angular frequency of each attenuation mechanism, the relaxation time data solution expression for the characteristic angular frequency of each attenuation mechanism is obtained, where the relaxation time data solution expression for the characteristic angular frequency is:

[0135]

[0136] In the expression for solving the relaxation time data of the characteristic angular frequency, ω represents the characteristic angular frequency, and τ ε τ σ This represents a set of relaxation times to be determined.

[0137] The relaxation time data for each attenuation mechanism is obtained by solving the expressions for the high-frequency limit value of the rock bulk modulus and the relaxation time data of the characteristic angular frequency.

[0138] In one embodiment, step S140 may include:

[0139] Based on the characteristic angular frequency and relaxation time data of each attenuation mechanism, the quality factor of each attenuation mechanism is obtained through the fourth expression, whereby:

[0140]

[0141] In the fourth expression, Q represents the quality factor, ω represents the characteristic angular frequency, and τ′ ε τ′ σ This represents a set of relaxation times after the solution is obtained;

[0142] Based on the characteristic angular frequency, relaxation time, and preset elastic modulus (which is generally a fixed value set according to actual needs) of each attenuation mechanism, the longitudinal wave dynamic modulus of each attenuation mechanism is obtained through the fifth expression, where the fifth expression is:

[0143]

[0144] In the fifth expression, E (ω) The dynamic modulus of the longitudinal wave is represented by ω, where L represents the number of attenuation mechanisms (l = 1, 2, ..., L represents each attenuation mechanism), ω represents the characteristic angular frequency, and τ′ represents the dynamic modulus of the longitudinal wave. ε τ′ σ This represents a set of relaxation times after the solution is obtained, where i represents the imaginary unit;

[0145] Based on the rock density data and preset elastic modulus of each attenuation mechanism, the P-wave complex velocity of each attenuation mechanism is obtained through the sixth expression, where the sixth expression is:

[0146]

[0147] In the sixth expression, V (ω) E represents the complex velocity of the longitudinal wave. R ρ represents the preset elastic modulus, and ρ represents the rock density data.

[0148] Based on the P-wave complex velocity of each attenuation mechanism, the P-wave velocity fitting curve for each attenuation mechanism is obtained through the seventh expression. The P-wave velocity fitting curves of all attenuation mechanisms constitute the fitting expression, where the seventh expression is:

[0149]

[0150] In the seventh expression, V p V represents the longitudinal wave velocity fitting curve. (ω) Represents the complex velocity of the longitudinal wave;

[0151] Based on the quality factors of all attenuation mechanisms and the P-wave velocity fitting curves, a fitting model for full-band rock physics measurement data is constructed.

[0152] Figure 3 The image shows the fitting effect of a model that uses four attenuation mechanisms to fit full-band rock physical measurement data. The model shows a high degree of agreement with experimentally measured velocity samples. Based on this model, subsequent 2D and 3D numerical simulations can be carried out efficiently. It should be noted that the more attenuation mechanisms selected, the better the fitting effect, but the complexity of the model and computational efficiency will decrease. In most cases, using a single attenuation mechanism can basically meet the velocity attenuation dispersion characteristics of rocks across the entire frequency band. Figure 4 As shown. Therefore, the number of attenuation mechanisms can be set according to actual needs.

[0153] The fitting and modeling method for full-band rock physics measurement data provided by this invention, based on full-band rock physics measurement data and combined with linear elastic body theory, enables rapid fitting and modeling of the velocity attenuation and dispersion characteristics of rocks across the entire frequency band. Based on this method, a concise and efficient fitting model for full-band rock physics measurement data can be quickly obtained. This model can uniformly simulate the velocity attenuation and dispersion characteristics of rocks under various attenuation mechanisms across the entire frequency band. It has high numerical simulation efficiency and high fitting accuracy, which can meet the needs of efficient industrial applications and also contribute to the research of oil and gas geophysics.

[0154] The following describes the fitting and modeling system for full-band rock physical measurement data provided by the present invention. The fitting and modeling system for full-band rock physical measurement data described below can be referred to in correspondence with the fitting and modeling method for full-band rock physical measurement data described above.

[0155] Reference Figure 5 The present invention provides a fitting and modeling system for full-band rock physical measurement data, which may include:

[0156] The data acquisition module is used to acquire rock velocity data and rock parameter data across the entire frequency band, where the entire frequency band is 10. -1 -10 6 The frequency band range of Hz;

[0157] The first data processing module is used to: determine the number of attenuation mechanisms and the frequency band range data of each attenuation mechanism based on the rock velocity data of the entire frequency band;

[0158] The second data processing module is used to: obtain the relaxation time data of each attenuation mechanism based on the frequency band range and rock parameter data of each attenuation mechanism;

[0159] The fitting modeling module is used to: obtain fitting expressions based on the relaxation time data of all decay mechanisms, and construct a fitting model for full-band rock physics measurement data.

[0160] In one implementation, the second data processing module may include:

[0161] The characteristic angular frequency acquisition submodule is used to: obtain the characteristic angular frequency of each attenuation mechanism based on the frequency data of the frequency band range of each attenuation mechanism, using a first expression, where the first expression is:

[0162]

[0163] In the first expression, ω represents the characteristic angular frequency. max ω represents the maximum frequency value within the bandwidth. min Indicates the minimum frequency value within the bandwidth;

[0164] The submodule for obtaining the low-frequency extremum of rock bulk modulus is used to: obtain the low-frequency extremum of rock bulk modulus for each attenuation mechanism based on the frequency band velocity data and rock density data of each attenuation mechanism, using a second expression, where the second expression is:

[0165]

[0166] In the second expression, K min This represents the low-frequency extreme value of the bulk modulus of rock, v min ρ represents the minimum velocity within the frequency band, and ρ represents the rock density data.

[0167] The submodule for obtaining the high-frequency extremum of rock bulk modulus is used to: obtain the high-frequency extremum of rock bulk modulus for each attenuation mechanism based on the frequency band velocity data and rock density data of each attenuation mechanism, using a third expression, where the third expression is:

[0168]

[0169] In the third expression, K max This represents the low-frequency extreme value of the bulk modulus of rock, v max ρ represents the maximum velocity within the frequency band, and ρ represents the rock density data.

[0170] The relaxation time data acquisition submodule is used to: obtain the relaxation time data of each decay mechanism based on the characteristic angular frequency and the extreme value of the rock bulk modulus of each decay mechanism.

[0171] In one implementation, the relaxation time data acquisition submodule may include:

[0172] The plane wave modulus acquisition submodule is used to: obtain the plane wave modulus of each attenuation mechanism based on the rock bulk modulus extremum and the shear modulus of dry rock for each attenuation mechanism.

[0173] The rock complex plane wave modulus acquisition submodule is used to: obtain the expression for the rock complex plane wave modulus of each attenuation mechanism based on the plane wave modulus and characteristic angular frequency of each attenuation mechanism;

[0174] The high-frequency limit value of rock bulk modulus is obtained by sub-module, which is used to: obtain the expression of the high-frequency limit value of rock bulk modulus for each attenuation mechanism based on the expression of the complex plane wave modulus of rock for each attenuation mechanism.

[0175] The relaxation time data solution expression is obtained by a submodule, which is used to: obtain the relaxation time data solution expression for the characteristic angular frequency of each attenuation mechanism based on the characteristic angular frequency of each attenuation mechanism;

[0176] The relaxation time solution submodule is used to: solve the expression based on the expression of the high-frequency limit value of the rock bulk modulus and the relaxation time data of the characteristic angular frequency for each attenuation mechanism, and obtain the relaxation time data for each attenuation mechanism.

[0177] In one embodiment, the expression for the plane wave modulus includes an expression for the maximum plane wave modulus and an expression for the minimum plane wave modulus, wherein the expression for the maximum plane wave modulus is:

[0178]

[0179] In the expression for the maximum modulus of a plane wave, E max K represents the maximum plane wave modulus. max μ represents the high-frequency extreme value of the bulk modulus of rock. m Indicates the shear modulus of dry rock;

[0180] The expression for the minimum plane wave modulus is:

[0181]

[0182] In the expression for the minimum plane wave modulus, E min K represents the minimum plane wave modulus. min μ represents the low-frequency extreme value of the bulk modulus of rock. m This represents the shear modulus of dry rock.

[0183] In one implementation, the expression for the complex plane wave modulus of rock is:

[0184]

[0185] In the expression for the complex plane wave modulus of rock, E ω E represents the complex plane wave modulus of rock. min τ represents the minimum plane wave modulus, ω represents the characteristic angular frequency, and τ represents the minimum plane wave modulus. ε τ σ This represents a set of relaxation times to be determined.

[0186] In one implementation scheme, the expression for the high-frequency limit value of the rock bulk modulus is:

[0187]

[0188] In the expression for the high-frequency limit value of the bulk modulus of rock, E min τ represents the minimum plane wave modulus. ε τ σ This represents a set of relaxation times to be determined.

[0189] In one implementation, the expression for solving the relaxation time data of the characteristic angular frequency is:

[0190]

[0191] In the expression for solving the relaxation time data of the characteristic angular frequency, ω represents the characteristic angular frequency, and τ ε τ σ This represents a set of relaxation times to be determined.

[0192] In one implementation, the fitting modeling module may include:

[0193] The quality factor acquisition submodule is used to: obtain the quality factor of each attenuation mechanism based on its characteristic angular frequency and relaxation time data, using a fourth expression, where the fourth expression is:

[0194]

[0195] In the fourth expression, Q represents the quality factor, ω represents the characteristic angular frequency, and τ′ ε τ′ σ This represents a set of relaxation times after the solution is obtained;

[0196] The longitudinal wave dynamic modulus acquisition submodule is used to: obtain the longitudinal wave dynamic modulus of each attenuation mechanism based on its characteristic angular frequency, relaxation time, and preset elastic modulus, using the fifth expression, where the fifth expression is:

[0197]

[0198] In the fifth expression, E (ω) Let L represent the longitudinal wave dynamic modulus, L represent the number of attenuation mechanisms, ω represent the characteristic angular frequency, and τ′ represent the characteristic angular frequency. ε τ′ σ This represents a set of relaxation times after the solution is obtained, where i represents the imaginary unit;

[0199] The P-wave complex velocity acquisition submodule is used to: obtain the P-wave complex velocity for each attenuation mechanism based on the rock density data and preset elastic modulus of each attenuation mechanism, using the sixth expression, where the sixth expression is:

[0200]

[0201] In the sixth expression, V (ω) E represents the complex velocity of the longitudinal wave. R ρ represents the preset elastic modulus, and ρ represents the rock density data.

[0202] The P-wave velocity fitting curve acquisition submodule is used to: obtain the P-wave velocity fitting curve for each attenuation mechanism based on the complex P-wave velocity of each attenuation mechanism, using the seventh expression. The P-wave velocity fitting curves of all attenuation mechanisms constitute the fitting expression, where the seventh expression is:

[0203]

[0204] In the seventh expression, V p V represents the longitudinal wave velocity fitting curve. (ω) Represents the complex velocity of the longitudinal wave;

[0205] The modeling submodule is used to construct a fitting model for full-band rock physics measurement data based on the quality factors and P-wave velocity fitting curves of all attenuation mechanisms.

[0206] Figure 6 An example is a schematic diagram of the physical structure of an electronic device, such as... Figure 6 As shown, the electronic device may include: a processor 810, a communication interface 820, a memory 830, and a communication bus 840, wherein the processor 810, the communication interface 820, and the memory 830 communicate with each other via the communication bus 840. The processor 810 can call logical instructions in the memory 830 to execute a fitting modeling method for full-band rock physical measurement data, the method including:

[0207] Acquire rock velocity data and rock parameter data across the entire frequency band, where the entire frequency band is 10. -1 -10 6 The frequency band range of Hz;

[0208] Based on the rock velocity data across the entire frequency band, determine the number of attenuation mechanisms and the frequency band range data for each attenuation mechanism;

[0209] Based on the frequency band range and rock parameter data of each attenuation mechanism, the relaxation time data of each attenuation mechanism is obtained;

[0210] Based on the relaxation time data of all attenuation mechanisms, a fitting expression is obtained, and a fitting model for full-band rock physics measurement data is constructed.

[0211] Furthermore, the logical instructions in the aforementioned memory 830 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, essentially, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0212] On the other hand, the present invention also provides a computer program product, the computer program product comprising a computer program that can be stored on a non-transitory computer-readable storage medium, wherein when the computer program is executed by a processor, the computer is capable of executing the fitting modeling method for full-band rock physical measurement data provided by the above methods, the method comprising:

[0213] Acquire rock velocity data and rock parameter data across the entire frequency band, where the entire frequency band is 10. -1 -10 6 The frequency band range of Hz;

[0214] Based on the rock velocity data across the entire frequency band, determine the number of attenuation mechanisms and the frequency band range data for each attenuation mechanism;

[0215] Based on the frequency band range and rock parameter data of each attenuation mechanism, the relaxation time data of each attenuation mechanism is obtained;

[0216] Based on the relaxation time data of all attenuation mechanisms, a fitting expression is obtained, and a fitting model for full-band rock physics measurement data is constructed.

[0217] In another aspect, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements a fitting modeling method for full-band rock physical measurement data provided by the methods described above, the method comprising:

[0218] Acquire rock velocity data and rock parameter data across the entire frequency band, where the entire frequency band is 10. -1 -10 6 The frequency band range of Hz;

[0219] Based on the rock velocity data across the entire frequency band, determine the number of attenuation mechanisms and the frequency band range data for each attenuation mechanism;

[0220] Based on the frequency band range and rock parameter data of each attenuation mechanism, the relaxation time data of each attenuation mechanism is obtained;

[0221] Based on the relaxation time data of all attenuation mechanisms, a fitting expression is obtained, and a fitting model for full-band rock physics measurement data is constructed.

[0222] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0223] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0224] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A fitting and modeling method for full-band rock physical measurement data, characterized in that, include: Acquire rock velocity data and rock parameter data across the entire frequency band, where the entire frequency band is 10. -1 -10 6 The frequency band range of Hz; Based on the rock velocity data across the entire frequency band, determine the number of attenuation mechanisms and the frequency band range data for each attenuation mechanism; Based on the frequency band range and rock parameter data of each attenuation mechanism, the relaxation time data of each attenuation mechanism is obtained; Based on the relaxation time data of all attenuation mechanisms, a fitting expression is obtained, and a fitting model for full-band rock physics measurement data is constructed.

2. The fitting and modeling method for full-band rock physical measurement data according to claim 1, characterized in that, Rock parameter data include: rock porosity data, rock density data, fluid composition data, fluid density data, fluid saturation data, and dry rock shear modulus; The frequency band range data includes frequency band range data and velocity band range data. The frequency band range data includes the maximum frequency value and the minimum frequency value of the frequency band range, and the velocity band range data includes the maximum velocity value and the minimum velocity value of the frequency band range.

3. The fitting and modeling method for full-band rock physical measurement data according to claim 2, characterized in that, The relaxation time data for each attenuation mechanism is obtained based on its frequency band range and rock parameter data, including: Based on the frequency data within the bandwidth of each attenuation mechanism, the characteristic angular frequency of each attenuation mechanism is obtained; based on the velocity data within the bandwidth of each attenuation mechanism and the rock density data, the characteristic angular frequency of each attenuation mechanism is obtained. The extreme values ​​of the bulk modulus of the rock; The relaxation time data for each decay mechanism are obtained based on the characteristic angular frequency and the extreme value of the rock bulk modulus for each decay mechanism.

4. The fitting and modeling method for full-band rock physical measurement data according to claim 3, characterized in that, The step of obtaining the characteristic angular frequency of each attenuation mechanism based on the frequency data of the frequency band range of each attenuation mechanism includes: Based on the frequency data of the bandwidth of each attenuation mechanism, the characteristic angular frequency of each attenuation mechanism is obtained through the first expression; The first expression is: ， In the first expression, ω represents the characteristic angular frequency. mαx ω represents the maximum frequency value within the bandwidth. min This indicates the minimum frequency value within the bandwidth.

5. The fitting and modeling method for full-band rock physical measurement data according to claim 4, characterized in that, The extreme values ​​of rock bulk modulus include low-frequency extreme values ​​and high-frequency extreme values. The process of obtaining the extreme values ​​of rock bulk modulus for each attenuation mechanism based on the frequency band velocity data and rock density data for each attenuation mechanism includes: Based on the frequency band velocity data and rock density data of each attenuation mechanism, the low-frequency extreme value of the rock bulk modulus of each attenuation mechanism is obtained through the second expression. Based on the frequency band velocity data and rock density data of each attenuation mechanism, the high-frequency extreme value of the rock bulk modulus of each attenuation mechanism is obtained through the third expression. The second expression is: ， In the second expression, K min This represents the low-frequency extreme value of the bulk modulus of rock, v min ρ represents the minimum velocity within the frequency band, and ρ represents the rock density data. The third expression is: ， In the third expression, K mαx This represents the low-frequency extreme value of the bulk modulus of rock, v mαx ρ represents the maximum velocity within the frequency band, and ρ represents the rock density data.

6. The fitting and modeling method for full-band rock physical measurement data according to claim 5, characterized in that, The relaxation time data for each attenuation mechanism is obtained based on its characteristic angular frequency and the extreme value of the rock bulk modulus, including: Based on the extreme values ​​of the bulk modulus of rock for each attenuation mechanism and the shear modulus of dry rock, the plane wave modulus of each attenuation mechanism is obtained. Based on the plane wave modulus and characteristic angular frequency of each attenuation mechanism, the expression for the complex plane wave modulus of rock for each attenuation mechanism is obtained; Based on the expression for the complex plane wave modulus of rock for each attenuation mechanism, the expression for the high-frequency limit value of the bulk modulus of rock for each attenuation mechanism is obtained. Based on the characteristic angular frequency of each attenuation mechanism, the relaxation time data of each attenuation mechanism is obtained and the expression is solved. The relaxation time data for each attenuation mechanism is obtained by solving the expression for the high-frequency limit value of the rock bulk modulus and the relaxation time data of the characteristic angular frequency. The expression for the plane wave modulus includes expressions for the maximum and minimum values ​​of the plane wave modulus. The expression for the maximum plane wave modulus is: ， In the expression for the maximum modulus of a plane wave, E max K represents the maximum plane wave modulus. max μ represents the high-frequency extreme value of the bulk modulus of rock. m Indicates the shear modulus of dry rock; The expression for the minimum plane wave modulus is: , In the expression for the minimum plane wave modulus, E min K represents the minimum plane wave modulus. min μ represents the low-frequency extreme value of the bulk modulus of rock. m Indicates the shear modulus of dry rock; The expression for the complex plane wave modulus of rock is: , In the expression for the complex plane wave modulus of rock, E ω E represents the complex plane wave modulus of rock. min This represents the minimum plane wave modulus, and ω represents the characteristic angular frequency. This represents a set of relaxation times to be determined. The expression for the high-frequency limit value of the bulk modulus of rock is: , In the expression for the high-frequency limit value of the bulk modulus of rock, E min This represents the minimum plane wave modulus. This represents a set of relaxation times to be determined. The expression for solving the relaxation time data of the characteristic angular frequency is: , In the expression for solving the relaxation time data of the characteristic angular frequency, ω represents the characteristic angular frequency. This represents a set of relaxation times to be determined.

7. The fitting and modeling method for full-band rock physical measurement data according to claim 6, characterized in that, The process involves obtaining a fitting expression based on the relaxation time data of all attenuation mechanisms and constructing a fitting model for full-band rock physics measurement data, including: Based on the characteristic angular frequency and relaxation time data of each attenuation mechanism, the quality factor of each attenuation mechanism is obtained. Based on the characteristic angular frequency, relaxation time, and preset elastic modulus of each attenuation mechanism, the longitudinal wave dynamic modulus of each attenuation mechanism is obtained. Based on the rock density data and preset elastic modulus of each attenuation mechanism, the longitudinal wave complex velocity of each attenuation mechanism is obtained. Based on the P-wave complex velocity of each attenuation mechanism, the P-wave velocity fitting curve of each attenuation mechanism is obtained, and the P-wave velocity fitting curves of all attenuation mechanisms constitute a fitting expression. Based on the quality factors of all attenuation mechanisms and the P-wave velocity fitting curves, a fitting model for full-band rock physics measurement data is constructed. The quality factor of each attenuation mechanism is obtained based on its characteristic angular frequency and relaxation time data, including: Based on the characteristic angular frequency and relaxation time data of each decay mechanism, the quality factor of each decay mechanism is obtained through the fourth expression. The fourth expression is: , In the fourth expression, Q represents the quality factor, and ω represents the characteristic angular frequency. This represents a set of relaxation times after the solution is obtained; The step of obtaining the longitudinal wave dynamic modulus of each attenuation mechanism based on its characteristic angular frequency, relaxation time, and preset elastic modulus includes: Based on the characteristic angular frequency, relaxation time, and preset elastic modulus of each attenuation mechanism, the longitudinal wave dynamic modulus of each attenuation mechanism is obtained through the fifth expression. The fifth expression is: , In the fifth expression, E (ω) The dynamic modulus of the longitudinal wave is represented by L, the number of attenuation mechanisms is represented by ω, and the characteristic angular frequency is represented by ω. This represents a set of relaxation times after the solution is obtained, where i represents the imaginary unit; The step of obtaining the P-wave complex velocity for each attenuation mechanism based on the rock density data and preset elastic modulus for each attenuation mechanism includes: Based on the rock density data and preset elastic modulus of each attenuation mechanism, the longitudinal wave complex velocity of each attenuation mechanism is obtained through the sixth expression. The sixth expression is: , In the sixth expression, V (ω) E represents the complex velocity of the longitudinal wave. R ρ represents the preset elastic modulus, and ρ represents the rock density data; the longitudinal wave velocity of each attenuation mechanism is obtained based on the complex velocity of the longitudinal wave for each attenuation mechanism. The fitted curves include: Based on the complex P-wave velocity of each attenuation mechanism, the P-wave velocity fitting curve for each attenuation mechanism is obtained through the seventh expression. The seventh expression is: ， In the seventh expression, V p V represents the longitudinal wave velocity fitting curve. (ω) This represents the complex velocity of the longitudinal wave.

8. A fitting and modeling system for full-band rock physical measurement data, characterized in that, include: The data acquisition module is used to acquire rock velocity data and rock parameter data across the entire frequency band, where the entire frequency band is 10. -1 -10 6 The frequency range in Hz; The first data processing module is used to: determine the number of attenuation mechanisms and the frequency band range data of each attenuation mechanism based on the rock velocity data of the entire frequency band; The second data processing module is used to: obtain the relaxation time data of each attenuation mechanism based on the frequency band range and rock parameter data of each attenuation mechanism; The fitting modeling module is used to: obtain fitting expressions based on the relaxation time data of all decay mechanisms, and construct a fitting model for full-band rock physics measurement data.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the fitting and modeling method for full-band rock physical measurement data as described in any one of claims 1 to 7.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the fitting and modeling method for full-band rock physical measurement data as described in any one of claims 1 to 7.

Images