CN112711918A - Air bubble size distribution inversion method based on Gaussian function fitting - Google Patents
Air bubble size distribution inversion method based on Gaussian function fitting Download PDFInfo
- Publication number
- CN112711918A CN112711918A CN202011464262.5A CN202011464262A CN112711918A CN 112711918 A CN112711918 A CN 112711918A CN 202011464262 A CN202011464262 A CN 202011464262A CN 112711918 A CN112711918 A CN 112711918A
- Authority
- CN
- China
- Prior art keywords
- size distribution
- bubble
- coefficient
- gaussian function
- inversion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Databases & Information Systems (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a Gaussian function fitting-based bubble size distribution inversion method. And then establishing a cost function required by inversion by combining with the equivalent density fluid approximate model corrected by bubble vibration as a forward model, and optimizing by utilizing a mixed optimization algorithm DEPSO to obtain a corrected Gaussian function coefficient to be inverted so as to obtain the size distribution of bubbles in the gas-containing marine sediments. The method considers the existence form of the size distribution of the bubbles in the sediment, can reduce the number of parameters to be inverted to a certain degree, can realize inversion of the size distribution of the bubbles by only utilizing attenuation information, and has good coincidence between the attenuation data fitted by an inversion result and the measured attenuation data. The method is an innovative method applied in the field of submarine shallow gas exploration and identification.
Description
Technical Field
The invention belongs to the technical field of marine acoustics and seabed shallow gas exploration and identification, and particularly relates to a gas bubble size distribution inversion method based on Gaussian function fitting.
Background
In modern engineering construction of submarine tunnels, cross-sea bridges and the like, unnecessary engineering accidents are often caused by negligence of distribution of shallow gas in submarine sediments. The sediments in the domestic long triangle and Hangzhou bay areas have a large amount of shallow formation gas distribution, and the dynamic balance of factors such as internal sediments, fluid, pore water, an upper coating, gas, temperature, pressure and the like is very weak, so that accident disasters are more and more serious.
For engineering applications, in addition to knowing the geographical location of the airborne deposits, more information needs to be known: the content of gas and the size distribution of the bubbles. For example, meteorologists need to estimate how much the ocean source contributes to atmospheric methane by the content of gas in the sediment; the sediment pore pressure and sediment strength are sensitive to the bubble content and size distribution, and knowledge of these parameters allows the oil exploration industry to more reliably locate offshore structures and avoid blowouts during drilling operations. In addition, the acoustic properties of gas bearing sediments are sensitive to bubble content and bubble size distribution, making this information important for shallow sea sound fields. The acoustic inversion technique for bubbles in marine sediments is to infer the content of bubbles and the distribution of bubble sizes from the measured compressional wave acoustic properties by corresponding acoustic theory. Since the bubble sizes in the natural marine sediment medium often show a complex distribution, a gaussian distribution, a rayleigh distribution, and the like. Generally, the attenuation coefficient of the gas-containing sediment is measured in a frequency band corresponding to the size distribution of the bubbles, the result of the attenuation coefficient often presents a plurality of obvious peak values, and the number of fitted parameters can be reduced to a certain extent by considering that the corrected Gaussian function is used for fitting the size distribution of the bubbles.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provides a bubble size distribution inversion method based on Gaussian function fitting.
The purpose of the invention is achieved by the following technical scheme: a bubble size distribution inversion method based on Gaussian function fitting comprises the following steps: firstly, aiming at the attenuation coefficient peak value of gas-containing sediments observed in experiments, representing the size distribution function of bubbles by using a modified Gaussian function, then setting a cost function as the square sum of the simulated attenuation coefficient and the measured attenuation coefficient error, and inverting the modified Gaussian function coefficient by using an optimization algorithm; and finally fitting the size distribution of the bubbles by using the Gaussian function coefficient obtained by inversion. The method considers the existence form of the size distribution of the bubbles in the sediment, can reduce the number of parameters to be inverted to a certain degree, can realize inversion of the size distribution of the bubbles by only utilizing attenuation information, and has good coincidence between the attenuation data fitted by an inversion result and the measured attenuation data. The method is an innovative method applied in the field of submarine shallow gas exploration and identification.
Further, the method specifically comprises the following steps:
(1) constructing the bubble size distribution of the model input, wherein the first term represents the attenuation coefficient of the resonance peak of the fitting bubble realized by using a plurality of modified Gaussian distributions, and the second term adjusts the attenuation coefficient of the non-resonance peak by using a plurality of linear functions, in the initially constructed bubble size distribution, the parameter ai3Other parameters are randomly generated corresponding to the bubble radius of the attenuation peak value;
(a) distribution of bubble radius size, a being the variation of bubble radius, Ni、ai1、ai2、ai3For the coefficient of the Gaussian function to be solved, kj、bjIs the linear function coefficient to be solved;
(2) setting initial model input environment parameters:
(3) substituting the equivalent density fluid approximate model corrected by bubble vibration to calculate the attenuation coefficient of the gas-containing marine sediment;
ρefffor deposit equivalent fluid density, KeffFor the equivalent bulk modulus of the deposit,ω is the angular frequency, b (a) is the damping coefficient for the bubble size vibration, i is the imaginary unit, ω is the equivalent fluid density for bubble vibration correction0Is the bubble resonance frequency, beta is the sediment porosity;
(4) calculating an objective function of an inversion setWhereinAs a measure of the attenuation coefficient of the sound wave at different frequencies,predicting a value for the forward model at the corresponding frequency;
(5) repeating steps (1) - (4) using a two-stage hybrid optimization algorithm (DEPSO: differential evolution algorithm and particle swarm optimization).
Further, the model input environment parameters include: physical parameters of the sediment: particle size, particle density, particle bulk modulus, pore water density, pore water bulk modulus, pore water viscosity coefficient, porosity, permeability, tortuosity; gas physical parameters: gas density, static pressure, thermal diffusivity of the gas, surface tension, specific heat ratio.
The invention has the beneficial effects that: according to the invention, the corrected Gaussian function fitted bubble size distribution is used, compared with the traditional bubble inversion method, the possible distribution in nature is considered, the number of the parameters to be inverted is reduced to a certain extent, and the attenuation coefficient fitted by utilizing the inverted bubble size distribution is well matched with the measured attenuation coefficient.
Drawings
FIG. 1 is a schematic flow chart of the present invention.
FIG. 2 is a graph showing the comparison of observed attenuation coefficients to model fitted attenuation coefficients.
FIG. 3 is a schematic diagram of the optimal bubble size distribution obtained by model fitting.
Detailed Description
The invention will be described in detail below with reference to the following drawings:
the invention discloses a Gaussian function fitting-based bubble size distribution inversion method, which comprises the following steps of: and fitting the unknown bubble size distribution function in the marine sediments by using the corrected Gaussian function, and solving the coefficient of the corrected Gaussian function by using an optimization algorithm to invert the bubble size distribution in the marine sediments. The method specifically comprises the following steps: firstly, aiming at the attenuation coefficient peak value of the gas-containing sediment observed in the experiment, the size distribution function of the bubbles is expressed by a modified Gaussian function. And then setting the cost function as the square sum of the simulated attenuation coefficient and the measured attenuation coefficient error, and inverting the corrected Gaussian function coefficient by using an optimization algorithm. And finally fitting the size distribution of the bubbles by using the Gaussian function coefficient obtained by inversion. The method can accurately invert the size distribution of the bubbles under the condition that the acoustic velocity and the attenuation information of the gas-containing sediments are not completely known, and the fitted attenuation coefficient is well matched with the experimentally measured attenuation coefficient by taking the inversion only using the attenuation coefficient as an example.
(1) Constructing the bubble size distribution of the model input, wherein the first term represents the attenuation coefficient of fitting the bubble formant with a plurality of modified Gaussian distributions, and the second term adjusts the non-formant with a plurality of linear functionsThe attenuation coefficient of (2). In the initially constructed bubble size distribution, the parameter ai3Other parameters were randomly generated corresponding to the bubble radius of the decay peak.
(a) distribution of bubble radius size, a being the variation of bubble radius, Ni、ai1、ai2、ai3For the coefficient of the Gaussian function to be solved, kj、bjIs the linear function coefficient to be solved;
(2) setting initial model input environment parameters: physical parameters of the sediment: particle size, particle density, particle bulk modulus, pore water density, pore water bulk modulus, pore water viscosity coefficient, porosity, permeability, tortuosity, and the like; gas physical parameters: gas density, static pressure, thermal diffusivity of the gas, surface tension, specific heat ratio, and the like.
(3) Substituting the equivalent density fluid approximate model corrected by bubble vibration to calculate the attenuation coefficient of the aerated marine sediment.
ρeffFor deposit equivalent fluid density, KeffFor the equivalent bulk modulus of the deposit,for the equivalent fluid density of bubble vibration correction, ω is the angular frequency, b (a) is the corresponding bubble sizeDamping coefficient of vibration, i being imaginary unit, ω0Is the bubble resonance frequency, beta is the sediment porosity;
(4) calculating an objective function of an inversion setWhereinAs a measure of the attenuation coefficient of the sound wave at different frequencies,and predicting values for the forward model at the corresponding frequency.
(5) Repeating steps (1) - (4) using a two-stage hybrid optimization algorithm (DEPSO: differential evolution algorithm and particle swarm optimization).
Fig. 1 is a diagram illustrating inversion of bubble size distribution by an inversion method for fitting bubble size distribution based on a modified gaussian function, and comparison of an attenuation coefficient obtained by fitting with an attenuation coefficient predicted by a model. In the figure, the solid line represents the measured attenuation coefficient, and the dotted line represents the attenuation coefficient predicted by the model. 4 attenuation peaks caused by bubble resonance can be obviously observed from the attenuation data of 500-3000Hz, the inversion method can accurately fit the four peaks A-D in the attenuation curve, and the attenuation of 500-3000Hz frequency is consistent with the measured attenuation data.
Fig. 2 is the optimal bubble size distribution obtained by model fitting, and it can be clearly seen that the bubble size distribution has 3 peaks of 1.236mm, 1.614mm and 2.025mm, and the peak of 3.62mm is not obvious, which corresponds to the smaller amplitude of the attenuation coefficient peak a.
Fig. 3 is a bubble size distribution inversion flow chart, and the whole inversion processing process is described in detail in the "implementation method".
The frequencies of the four formants A-D and their corresponding attenuation coefficients are given in Table 1. The attenuation of the formant A reaches 109dB/m at 900Hz, the attenuation of the formant B reaches 296dB/m at 1600Hz, the attenuation of the formant C reaches 365dB/m at 2000Hz, and the attenuation of the formant D reaches 476dB/m at 2600 Hz. The bubble radii predicted by the corresponding decay peak model were 3.62mm, 2.025mm, 1.614mm, and 1.236mm, respectively.
TABLE 1 model prediction of observed formants
It should be understood that equivalent substitutions and changes to the technical solution and the inventive concept of the present invention should be made by those skilled in the art to the protection scope of the appended claims.
Claims (3)
1. A bubble size distribution inversion method based on Gaussian function fitting is characterized by comprising the following steps: the method comprises the following steps: firstly, aiming at the attenuation coefficient peak value of gas-containing sediments observed in experiments, representing the size distribution function of bubbles by using a modified Gaussian function, then setting a cost function as the square sum of the simulated attenuation coefficient and the measured attenuation coefficient error, and inverting the modified Gaussian function coefficient by using an optimization algorithm; and finally fitting the size distribution of the bubbles by using the Gaussian function coefficient obtained by inversion.
2. The bubble size distribution inversion method based on gaussian function fitting according to claim 1, wherein: the method specifically comprises the following steps:
(1) constructing the bubble size distribution of the model input, wherein the first term represents the attenuation coefficient of the resonance peak of the fitting bubble realized by using a plurality of modified Gaussian distributions, and the second term adjusts the attenuation coefficient of the non-resonance peak by using a plurality of linear functions, in the initially constructed bubble size distribution, the parameter ai3Other parameters are randomly generated corresponding to the bubble radius of the attenuation peak value;
(a) distribution of bubble radius size, a being the variation of bubble radius, Ni、ai1、ai2、ai3For the coefficient of the Gaussian function to be solved, kj、bjIs the linear function coefficient to be solved;
(2) setting initial model input environment parameters:
(3) substituting the equivalent density fluid approximate model corrected by bubble vibration to calculate the attenuation coefficient of the gas-containing marine sediment;
ρefffor deposit equivalent fluid density, KeffFor the equivalent bulk modulus of the deposit,ω is the angular frequency, b (a) is the damping coefficient for the bubble size vibration, i is the imaginary unit, ω is the equivalent fluid density for bubble vibration correction0Is the bubble resonance frequency, beta is the sediment porosity;
(4) calculating an objective function of an inversion setWhereinAs a measure of the attenuation coefficient of the sound wave at different frequencies,predicting a value for the forward model at the corresponding frequency;
(5) repeating steps (1) - (4) using a two-stage hybrid optimization algorithm.
3. The bubble size distribution inversion method based on gaussian function fitting according to claim 2, wherein: the model input environment parameters comprise: physical parameters of the sediment: particle size, particle density, particle bulk modulus, pore water density, pore water bulk modulus, pore water viscosity coefficient, porosity, permeability, tortuosity; gas physical parameters: gas density, static pressure, thermal diffusivity of the gas, surface tension, specific heat ratio.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011464262.5A CN112711918B (en) | 2020-12-14 | 2020-12-14 | Bubble size distribution inversion method based on Gaussian function fitting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011464262.5A CN112711918B (en) | 2020-12-14 | 2020-12-14 | Bubble size distribution inversion method based on Gaussian function fitting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112711918A true CN112711918A (en) | 2021-04-27 |
CN112711918B CN112711918B (en) | 2023-09-01 |
Family
ID=75541932
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011464262.5A Active CN112711918B (en) | 2020-12-14 | 2020-12-14 | Bubble size distribution inversion method based on Gaussian function fitting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112711918B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5784160A (en) * | 1995-10-10 | 1998-07-21 | Tsi Corporation | Non-contact interferometric sizing of stochastic particles |
WO2015094307A1 (en) * | 2013-12-19 | 2015-06-25 | Halliburton Energy Services, Inc. | Pore size classification in subterranean formations based on nuclear magnetic resonance (nmr) relaxation distributions |
CN110530765A (en) * | 2019-09-26 | 2019-12-03 | 哈尔滨工程大学 | Underwater bubble group's size distribution parameter inversion method based on measuring non-linear parameters |
CN110648342A (en) * | 2019-09-30 | 2020-01-03 | 福州大学 | Foam infrared image segmentation method based on NSST significance detection and image segmentation |
-
2020
- 2020-12-14 CN CN202011464262.5A patent/CN112711918B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5784160A (en) * | 1995-10-10 | 1998-07-21 | Tsi Corporation | Non-contact interferometric sizing of stochastic particles |
WO2015094307A1 (en) * | 2013-12-19 | 2015-06-25 | Halliburton Energy Services, Inc. | Pore size classification in subterranean formations based on nuclear magnetic resonance (nmr) relaxation distributions |
CN110530765A (en) * | 2019-09-26 | 2019-12-03 | 哈尔滨工程大学 | Underwater bubble group's size distribution parameter inversion method based on measuring non-linear parameters |
CN110648342A (en) * | 2019-09-30 | 2020-01-03 | 福州大学 | Foam infrared image segmentation method based on NSST significance detection and image segmentation |
Non-Patent Citations (1)
Title |
---|
黎章龙;屈科;张薇;何树斌;: "水中气泡分布的共振估计法的误差分析及修正", 海洋技术学报, no. 06 * |
Also Published As
Publication number | Publication date |
---|---|
CN112711918B (en) | 2023-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Raaen et al. | FORMEL: A step forward in strength logging | |
CN108957542B (en) | Method for establishing seismic wave attenuation rock physical drawing board | |
CN111859632B (en) | Rock physical model construction method and processing terminal for hydrate reservoir | |
Hovem | Viscous attenuation of sound in suspensions and high‐porosity marine sediments | |
CN113640119A (en) | Method for determining stress-related rock dynamic Biot coefficient | |
CN111381292B (en) | Logging interpretation method and device for predicting sandstone hydrocarbon-bearing reservoir | |
CN112711918B (en) | Bubble size distribution inversion method based on Gaussian function fitting | |
CN106597546B (en) | A method of detection river bed surface layer silt physical property | |
CN117055114A (en) | Quantitative analysis method for free gas saturation of reservoir sediment | |
CN109709610B (en) | Rock crack detection method and system | |
CN104483702B (en) | A kind of seismic forward simulation method being applicable to nonuniform motion water body | |
CN115906715A (en) | Method and system for calculating movement speed of compressible silt seabed soil under wave action | |
Hovem | The nonlinearity parameter of saturated marine sediments | |
CN114859420A (en) | Shallow sea target sorting and underwater target motion situation and depth judgment method | |
Noye | Tide-well systems II: The frequency response of a linear tide-well system | |
Kulik et al. | Using two-layer compliant coatings to control turbulent boundary layer | |
Zong et al. | An exact expression for the effective bulk modulus for acoustic wave propagation in cylindrical patchy-saturation rocks | |
Wang et al. | Comparison of sound speed and attenuation measurements to the corrected effective density fluid model for gassy sediments | |
Lalita et al. | Underwater acoustic propagation using Monterey-Miami parabolic equation in shallow water Kayeli Bay Buru Distric | |
CN107764697A (en) | Gas potential detection method based on the progressive equation non-linear inversion of pore media | |
Hovem | Attenuation of sound in marine sediments | |
Shaikh et al. | Acoustic propagation in ocean sediments using biot model | |
CN112799127B (en) | Seismic wave frequency dispersion and attenuation numerical calculation method considering non-uniform difference of seepage of fractured pore rock | |
CN113358466B (en) | Method and system for determining transfer coefficient of dynamic stress of layered foundation soil layer interface | |
Singh et al. | Shear-wave propagation in a loosely bonded inhomogeneous fiber-reinforced stratum under the influence of horizontal initial stress growing with depth |
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 |