Seismic wave dispersion and attenuation prediction method and related equipment

Abstract

The invention provides a seismic wave dispersion and attenuation prediction method and related equipment, and relates to the field of geophysical fluctuation induced fluid flow. The method comprises the steps that based on the bubble vibration theory and a pore fluid seepage continuity equation, a target elastic modulus is obtained, and the target elastic modulus is an elastic modulus obtained after bubble vibration correction is added; fusing the target elastic modulus into a strain energy function of a classical pore acoustic elasticity theory to obtain the strain energy function; deducing to obtain a wave equation of the fluid saturated porous medium containing a small amount of bubbles under the action of pressure by combining the strain energy function, the kinetic energy function and the dissipated energy function; and solving the simplified wave equation through the plane wave to obtain the seismic wave speed and attenuation. According to the embodiment of the invention, the frequency dispersion and attenuation of seismic waves can be accurately predicted within a certain effective pressure range.

Description

TECHNICAL FIELD

[0001] The present disclosure relates to the field of geophysical wave-induced fluid flow, and in particular, to a method for predicting seismic wave dispersion and attenuation and related equipment. BACKGROUND

[0002] Fluid-saturated porous media containing a small amount of bubbles in the subsurface formation are very common, but the research on seismic wave propagation in such media under the action of formation pressure is less. The acoustic-elastic theory proposed in the related art quantitatively describes the influence of effective pressure on seismic wave propagation in solid materials. According to the field observation, rock physics experiment and theoretical modeling, the acoustic-elastic effect in fluid-saturated porous media is obviously stronger than that in pure elastic solids, which shows that the fluid-saturated porous media is highly sensitive to effective pressure. In order to further explain the influence of pressure on seismic wave propagation in fluid-saturated porous media, the related art combines the Biot poroelastic theory and the acoustic-elastic theory to establish the poroacoustic-elastic theory. Subsequently, some technologies introduce the influence of nonlinear deformation of cracks on elastic modulus, which provides a good explanation for the nonlinear pressure dependence of seismic wave velocity. However, there are a small amount of bubbles in the actual formation, and according to the bubble vibration theory, even a small amount of bubbles can significantly affect the seismic wave velocity in the fluid.

[0003] At present, the related art does not consider the action of formation pressure and bubble vibration at the same time, and cannot accurately predict the dispersion and attenuation of seismic waves.

[0004] It should be noted that the information disclosed in the above background section is only used to strengthen the understanding of the background of the present disclosure, and therefore can include information that does not constitute prior art known to those of ordinary skill in the art. SUMMARY

[0005] The present disclosure provides a method for predicting seismic wave dispersion and attenuation and related equipment, which at least partially solves the problem that the related art does not consider the action of formation pressure and bubble vibration at the same time, and cannot accurately predict the dispersion and attenuation of seismic waves.

[0006] Other characteristics and advantages of the present disclosure will become apparent from the following detailed description, or will be learned by practice of the present disclosure.

[0007] According to one aspect of the present disclosure, a method for predicting seismic wave dispersion and attenuation is provided, comprising: obtaining a target elastic modulus based on bubble vibration theory and a percolation continuity equation of pore fluid, the target elastic modulus being an elastic modulus with bubble vibration correction; fusing the target elastic modulus into a strain energy function of classical poroacoustic-elastic theory to obtain a corrected strain energy function; deriving a wave equation of fluid-saturated porous media containing a small amount of bubbles under the action of pressure by combining the strain energy function, a kinetic energy function and a dissipation energy function; and obtaining seismic wave velocity and attenuation by simplifying the wave equation through plane wave solution.

[0008] According to another aspect of the present disclosure, there is provided a seismic wave dispersion and attenuation prediction device, comprising an elastic modulus correction module, a strain energy correction module, a wave equation derivation module and a prediction module.

[0009] The elastic modulus correction module is configured to obtain a target elastic modulus based on bubble vibration theory and a pore fluid percolation continuity equation, the target elastic modulus being an elastic modulus corrected by adding bubble vibration; The strain energy correction module is configured to fuse the target elastic modulus into a strain energy function of a classical pore acoustic-elasticity theory to obtain a corrected strain energy function; The wave equation derivation module is configured to derive a wave equation of a fluid-saturated porous medium containing a small amount of bubbles under pressure by combining the strain energy function, a kinetic energy function and a dissipation energy function; The prediction module is configured to obtain seismic wave velocity and attenuation by solving a simplified wave equation by means of a plane wave.

[0010] According to yet another aspect of the present disclosure, there is provided an electronic device, comprising a memory configured to store instructions, and a processor configured to invoke the instructions stored in the memory to implement the seismic wave dispersion and attenuation prediction method described above.

[0011] According to yet another aspect of the present disclosure, there is provided a computer-readable storage medium having computer instructions stored thereon, the computer instructions being executed by a processor to implement the seismic wave dispersion and attenuation prediction method described above.

[0012] According to yet another aspect of the present disclosure, there is provided a computer program product having instructions stored thereon, the instructions being executed by a computer to cause the computer to implement the seismic wave dispersion and attenuation prediction method described above.

[0013] According to yet another aspect of the present disclosure, there is provided a chip, comprising at least one processor and an interface, the interface being configured to provide program instructions or data for the at least one processor, and the at least one processor being configured to execute the program instructions to implement the seismic wave dispersion and attenuation prediction method described above.

[0014] The method and device for predicting seismic wave dispersion and attenuation provided by the embodiments of the present disclosure can more accurately simulate the seismic wave propagation characteristics in actual geological media by introducing bubble vibration theory and a corrected elastic modulus; the model can better adapt to complex and variable actual geological conditions by combining the classic poroelastic theory with new physical understanding (e.g., bubble vibration), for example, providing more accurate seismic wave propagation parameters for reservoir evaluation; the scheme provides a more comprehensive and accurate framework for predicting seismic wave dispersion and attenuation by comprehensively considering the influence of bubble vibration on seismic wave propagation, and is suitable for complex underground fluid-saturated porous media environment, and has important guiding significance for the field of oil and gas exploration.

[0015] It should be understood that the foregoing general description and the following detailed description are only exemplary and explanatory, and are not limiting to the present disclosure. BRIEF DESCRIPTION OF DRAWINGS

[0016] The accompanying drawings, which are incorporated into and form part of the specification, illustrate embodiments consistent with the present disclosure and, together with the description, serve to explain the principles of the disclosure.

[0017] Obviously, the drawings in the following description are only some embodiments of the present disclosure, and other drawings can be obtained by those skilled in the art without creative labor on the basis of these drawings.

[0018] Figure 1 A flowchart of a method for predicting seismic wave dispersion and attenuation is shown in the embodiments of the present disclosure; Figure 2 A schematic diagram of a fluid-saturated porous medium containing a small amount of bubbles is shown in the embodiments of the present disclosure; Figures 3a-3f A schematic diagram of the dispersion and attenuation of shear waves, fast longitudinal waves and slow longitudinal waves is shown in the embodiments of the present disclosure; Figure 4 A graph showing the change of velocity measured in the laboratory and the predicted seismic wave velocity with effective pressure in the embodiments of the present disclosure is shown; Figure 5 A schematic diagram of a device for predicting seismic wave dispersion and attenuation is shown in the embodiments of the present disclosure; Figure 6 A structural block diagram of an electronic device is shown in the embodiments of the present disclosure. DETAILED DESCRIPTION

[0019] In order to facilitate the understanding of the technical solutions of the present disclosure, the present disclosure will be further described below in conjunction with the drawings.

[0020] The terms "first" and "second" and the like in the description, claims and drawings of the present disclosure are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. Furthermore, the terms "comprises", "comprising", "includes", "including" and the like means "including but not limited to" and are not intended to exclude, for example, additional steps or units, equivalents or additions of the described processes, methods, systems, products or devices. The terms "comprises", "comprising", "includes", "including" and the like mean the following: comprising but not limited to, including but not limited to, and the like.

[0021] Reference herein to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the disclosure. The appearances of the phrase that an embodiment in various places in the specification are not necessarily all referring to the same embodiment, nor are they necessarily mutually exclusive of one another. As will be apparent to those of ordinary skill in the art, embodiments described herein can be combined with one another.

[0022] In the present disclosure, "at least one" means one or more, "multiple" means two or more, "at least two" means two or three and three or more, and "and / or" is used to describe the relationship between associated objects, indicating that there can be three relationships, for example, "A and / or B" can mean: only A, only B, and A and B exist at the same time, where A and B can be singular or plural. "Or" means that there can be two relationships, such as only A, only B; when A and B are not mutually exclusive, it can also mean that there are three relationships, such as only A, only B, and A and B exist at the same time. The character " / " generally represents an "or" relationship between the associated objects before and after it. "At least one of the following" or similar expressions means any combination of these items. For example, at least one of a, b or c can mean: a, b, c, "a and b", "a and c", "b and c", or "a and b and c".

[0023] Fluid-saturated porous media with a small amount of gas bubbles in the subsurface formation are common, but the research on seismic wave propagation in such media under the action of formation pressure is less. The acoustic-elastic theory proposed in the related art quantitatively describes the influence of effective pressure on seismic wave propagation in solid materials. According to the field observation, rock physics experiment and theoretical modeling, the acoustic-elastic effect in the fluid-saturated porous media is obviously stronger than that in the pure elastic solid, which shows that the fluid-saturated porous media is highly sensitive to the effective pressure. Therefore, in order to further explain the influence of pressure on seismic wave propagation in the fluid-saturated porous media, the related art combines the Biot poroelastic theory and the acoustic-elastic theory to establish the poroacoustic-elastic theory. Subsequently, some technologies introduce the influence of nonlinear deformation of cracks on the elastic modulus to provide a better explanation for the nonlinear pressure dependence of seismic wave velocity. However, there are a small amount of gas bubbles in the actual formation, and according to the bubble vibration theory, even a small amount of gas bubbles can significantly affect the seismic wave velocity in the fluid. The related art derives the wave equation of the fluid-saturated porous media with a small amount of gas bubbles by substituting the relationship between the bubble volume fraction and the fluid pressure into the percolation continuity equation, thereby predicting the dispersion and attenuation of seismic waves.

[0024] Although a series of prediction methods for seismic wave dispersion and attenuation are proposed in the related art, the formation pressure action and bubble vibration are not considered at the same time, and the seismic wave dispersion and attenuation cannot be accurately predicted. Therefore, on the basis of the bubble vibration theory and the classical poroacoustic-elastic theory, the wave equation of the fluid-saturated porous media with a small amount of gas bubbles under the action of pressure is derived, and a new theoretical model is established to predict the dispersion and attenuation of seismic waves.

[0025] In view of the defects in the above solutions and the proposed solutions, the inventors have come to the above-mentioned problems after careful research and practice, and therefore the discovery process of the above-mentioned problems and the solutions proposed by the present disclosure to solve the above-mentioned problems should be the contribution of the inventors to the present disclosure in the process of the present disclosure.

[0026] The present example embodiment will be described in detail below with reference to the accompanying drawings and examples.

[0027] First, a seismic wave dispersion and attenuation prediction method is provided in the present embodiment, which can be executed by any electronic device with computing processing capability.

[0028] Figure 1 A flowchart of a seismic wave dispersion and attenuation prediction method in the present embodiment is shown, as shown in Figure 1 The seismic wave dispersion and attenuation prediction method provided in the present embodiment includes S101-S104.

[0029] In S101, based on the bubble vibration theory and the pore fluid percolation continuity equation, a target elastic modulus is obtained, the target elastic modulus being an elastic modulus added with a bubble vibration correction; In S102, the target elastic modulus is fused into a strain energy function of a classical pore acoustic elasticity theory to obtain a corrected strain energy function; In S103, a wave equation of a fluid-saturated porous medium containing a small amount of bubbles under pressure is derived by combining the strain energy function, a kinetic energy function and a dissipation energy function; In S104, the seismic wave velocity and attenuation are obtained by solving the simplified wave equation with a plane wave.

[0030] The embodiments of the present disclosure, based on the bubble vibration theory and the pore fluid percolation continuity equation, calculate an elastic modulus added with a bubble vibration correction, and consider the elastic modulus in a strain energy function of a classical pore acoustic elasticity theory, and then derive a wave equation of a fluid-saturated porous medium containing a small amount of bubbles under pressure by combining a kinetic energy function and a dissipation energy function, so as to more accurately predict the dispersion and attenuation of seismic waves.

[0031] The implementation manners of the above steps are described in detail below.

[0032] In some embodiments, based on the bubble vibration theory and the pore fluid percolation continuity equation, the target elastic modulus is obtained, including: based on the dynamic response of the bubbles under the action of the acoustic wave, a nonlinear vibration equation describing the change of the bubble radius with time is established; the nonlinear vibration equation is simplified under the condition of low frequency to obtain a linear relationship between the bubble radius and the pore pressure; the linear relationship between the bubble radius and the pore pressure is substituted into the corresponding pore fluid percolation continuity equation to obtain the target elastic modulus.

[0033] The bubbles in the pore fluid meet the following five basic assumptions: (1) the bubbles are of the same size and are randomly and uniformly distributed in the pores; (2) the bubbles are not soluble with the liquid in the pores; (3) the gas in the bubbles is assumed to be an ideal gas, and the influence of the steam in the bubbles is ignored; (4) the bubbles always maintain a spherical shape and are not affected by the change of the effective pressure; (5) the amount of the bubbles is very small, and the influence of the bubbles on the modulus of the fluid can be ignored. When the volume fraction of the bubbles is small, the density of the fluid The following formula (1) is used.

[0034] (1) and are the densities of the gas and the liquid.

[0035] The volume fraction of the bubbles can be expressed as: (2) wherein, is the number of bubbles per unit volume, is the instantaneous radius of the bubble, is the volume of the bubble per unit pore, is the pore radius, is the volume per unit pore.

[0036] Assuming that the formation, rupture and coalescence of the bubble are ignored, the bubble will undergo radial forced vibration, and its vibration equation can be expressed as: (3) where . is the pressure at the location where the bubble is located when it does not exist, is the pressure on the liquid side at the gas-liquid boundary, is the pore pressure. is the speed of sound in the fluid, is the initial radius of the bubble, is the density of the liquid.

[0037] When , formula (2) can be linearized as: (4) Omitting terms of the second order and above, the relationship between the time derivative of the bubble volume fraction and the time derivative of the fluid pressure in the pore is: (5) Substituting formula (4) into the corresponding pore fluid seepage continuity equation and simplifying it into a seepage continuity equation similar to the Biot form, the elastic modulus corrected by considering the bubble vibration can be obtained by comparison: (6) where, , , , , . is the internal pressure of the bubble under steady-state conditions, is the static pressure of the liquid, is the surface tension coefficient of the gas-liquid interface, and are the thermal diffusivity and specific heat ratio of the gas. is the viscosity coefficient of the pore fluid, is the angular frequency, is the total porosity, denotes the target elastic modulus, and M denotes the elastic modulus before correction.

[0038] Figure 2A schematic diagram of a fluid-saturated porous medium containing a small amount of bubbles is shown, and a wave equation of the fluid-saturated porous medium containing a small amount of bubbles under pressure is derived by combining a strain energy function, a kinetic energy function, and a dissipation energy function in the above steps. The strain energy function is obtained by fusing a target elastic modulus into a strain energy function of a classical poroelasticity theory, and can specifically be a strain energy function of the classical poroelasticity theory, but the elastic modulus in the strain energy function of the classical poroelasticity theory is replaced by the target elastic modulus after vibration correction of the bubbles.

[0039] In some embodiments, the strain energy function of the fluid-saturated porous medium containing a small amount of bubbles under pressure can be expressed as: (7) , . , and are bulk moduli of the dry rock skeleton, the solid particles, and the pore fluid, respectively. is a shear modulus of the dry rock. is a finite strain of the solid. is a finite strain of the relative fluid flow under Lagrangian description, where , . , , , , , and are third-order elastic moduli.

[0040] In some embodiments, the wave equation of the fluid-saturated porous medium containing a small amount of bubbles under pressure is derived by combining the strain energy function, the kinetic energy function, and the dissipation energy function, including: introducing the strain energy function, the kinetic energy function, and the dissipation energy function into the generalized Lagrangian equation, and deriving the wave equation of the fluid-saturated porous medium containing a small amount of bubbles under pressure by deriving the solid displacement and the fluid relative displacement variables, respectively.

[0041] wherein the kinetic energy function and the dissipation energy function can be written as: (8) (9) wherein , . is tortuosity of the pore space, and are densities of the solid and the fluid, respectively. and ​​where u is the displacement of the solid skeleton and p is the fluid pressure, is the flow of fluid relative to the solid. is the permeability of the rock.

[0042] The general Lagrangian equation considering fluid dissipation effect is shown as follows: (10) (11) Substituting equations (7), (8) and (9) into equations (10) and (11) can obtain: (12) (13) where .

[0043] Assuming that the deformation caused by pressure is a uniform static deformation, which is much larger than the deformation caused by wave disturbance, the corresponding wave equation is (14) (15) wherein the related parameters are as follows: (16) (17) (18) (19) (20) (21) (22) (23) (24) In some embodiments, the seismic wave velocity and attenuation are obtained by solving the simplified wave equation by plane wave, including: setting the seismic wave disturbance as a simple harmonic plane wave form; substituting the plane wave solution into the separated wave equation, and decoupling the equation to obtain the complex wave number; based on the complex wave number, the seismic wave velocity and attenuation are obtained.

[0044] In some embodiments, according to the workflow of the generalized poroelastic theory, the above wave equation can be solved by plane wave to predict the dispersion and attenuation of seismic wave.

[0045] Therefore, the expression of the seismic wave velocity and attenuation is: (25) (26) (27) According to the dual-porosity model, it is known that the influence of the nonlinear deformation of soft pores is considered: (28) (29) (30) (31) (32) (33) wherein, , and are the porosities of the hard pores and the soft pores, respectively, under no pressure. is the volumetric strain of the soft pores, is the effective pressure. and are the volumetric compressibility of the dry rock and the solid particles, respectively. , , , , and are the pressure-sensitive coefficients.

[0046] Substituting equation (6) and equations (30)-(33) into the expression of the coefficients of the wave equation, it is obtained that: (34) (35) (36) (37) wherein, , , , , , . Then, substituting equations (34)-(37) into equations (25)-(27), the dispersion and attenuation of the seismic wave are predicted.

[0047] The embodiment of the present disclosure establishes a new theoretical model by simultaneously considering the classical pore acoustic-elastic theory and the bubble vibration theory, and more accurately predicts the dispersion and attenuation of seismic waves in the actual complex medium underground. The modeling results show that the presence of bubbles will cause significant fluctuations in the dispersion and attenuation of the longitudinal wave, while the dispersion and attenuation of the transverse wave are basically not affected by the bubble content. Through comparison with the experimental data of the rock sample, the rationality and applicability of the method are verified. The embodiment of the present disclosure has important guiding significance for practical applications such as oil and gas exploration and reservoir fluid identification.

[0048] Figures 3a-3f The dispersion and attenuation of the transverse wave, the fast longitudinal wave and the slow longitudinal wave are shown. Figure 4 The comparison of the predicted velocity with the measured velocity in the laboratory is shown. As can be seen from the figure, the predicted seismic wave velocity of the embodiment of the present disclosure is in good agreement with the laboratory measurement data, which shows the rationality and practicability of the embodiment of the present disclosure, and it can accurately predict the dispersion and attenuation of the seismic wave within a certain effective pressure range.

[0049] Based on the same inventive concept, the embodiment of the present disclosure also provides a seismic wave dispersion and attenuation prediction device, as shown in Figure 4 The seismic wave dispersion and attenuation prediction device comprises an elastic modulus correction module 501, a strain energy module 502, a wave equation derivation module 503 and a prediction module 504.

[0050] The elastic modulus correction module 501 is used to obtain a target elastic modulus based on the bubble vibration theory and the pore fluid percolation continuity equation, and the target elastic modulus is the elastic modulus after adding the bubble vibration correction; The strain energy correction module 502 is used to fuse the target elastic modulus into the strain energy function of the classical pore acoustic-elastic theory to obtain a corrected strain energy function; The wave equation derivation module 503 is used to derive the wave equation of the fluid-saturated porous medium containing a small amount of bubbles under the action of pressure by combining the strain energy function, the kinetic energy function and the dissipation energy function; The prediction module 504 is used to obtain the seismic wave velocity and attenuation by solving the simplified wave equation by means of the plane wave.

[0051] In some embodiments, the elastic modulus correction module 501 is used to establish a nonlinear vibration equation describing the change of the bubble radius with time based on the dynamic response of the bubble under the action of the acoustic wave; the nonlinear vibration equation is simplified under the condition of low frequency to obtain a linear relationship between the bubble radius and the pore pressure; and the linear relationship between the bubble radius and the pore pressure is substituted into the corresponding pore fluid percolation continuity equation to obtain the target elastic modulus.

[0052] In some embodiments, the strain energy correction module 502 is configured to retain the form of the strain energy function of the classical poroelasticity theory, but replace the elastic modulus in it with the target elastic modulus corrected by the bubble vibration.

[0053] In some embodiments, the wave equation derivation module 503 is configured to derive the strain energy function, kinetic energy function and dissipation energy function into the generalized Lagrangian equation, and derive the wave equation of the fluid-saturated porous medium containing a small amount of bubbles under pressure with respect to the solid displacement and fluid relative displacement variables, respectively.

[0054] In some embodiments, the prediction module 504 is configured to set the seismic wave disturbance in the form of a simple harmonic plane wave; substitute the plane wave solution into the separated wave equation, and decouple the equation to obtain the complex wave number; and obtain the seismic wave velocity and attenuation based on the complex wave number.

[0055] The terms "first", "second", and the like in the present disclosure are merely used to distinguish different devices, modules or units, and do not imply the order or interdependence of the functions performed by these devices, modules or units.

[0056] As to the seismic wave dispersion and attenuation prediction device in the above embodiments, the specific manner in which each module performs operations has been described in detail in the embodiments of the seismic wave dispersion and attenuation prediction method, and will not be described in detail here.

[0057] It should be noted that, although several modules or units of the device for action execution are mentioned in the above detailed description, such division is not mandatory. In fact, according to the embodiments of the present disclosure, the features and functions of two or more modules or units described above can be embodied in one module or unit. Conversely, the features and functions of one module or unit described above can be further divided into several modules or units.

[0058] Some of the block diagrams shown in the drawings are functional entities, which do not necessarily correspond to physically or logically independent entities. These functional entities can be implemented in the form of software, or in one or more hardware modules or integrated circuits, or in different network and / or processor devices and / or microcontroller devices.

[0059] The electronic device provided by the embodiments of the present disclosure will be described below with reference to Figure 6 The electronic device 600 shown is merely an example, and should not impose any limitation on the functions and use range of the embodiments of the present disclosure. Figure 6 The electronic device 600 shown is merely an example, and should not impose any limitation on the functions and use range of the embodiments of the present disclosure.

[0060] Figure 6An architecture schematic of an electronic device 600 is shown. As shown, the electronic device 600 includes, but is not limited to, at least one processor 610, at least one memory 620. Figure 6

[0061] The memory 620 is configured to store instructions.

[0062] In some embodiments, the memory 620 can include a readable medium in the form of a volatile storage unit such as a random access memory (RAM) 6201 and / or a cache 6202, and can further include a non-volatile storage unit (ROM) 6203.

[0063] In some embodiments, the memory 620 can further include a program / utility 6204 having a set (at least one) of program modules 6205, including but not limited to an operating system, one or more application programs, other program modules, and program data, each of which or a combination can include implementation of a network environment.

[0064] In some embodiments, the memory 620 can store an operating system. The operating system can be a real-time operating system (RTX), LINUX, UNIX, WINDOWS, or OS X.

[0065] In some embodiments, the memory 620 can also store data.

[0066] As an example, the processor 610 can read data stored in the memory 620, which can be stored in the same storage address as the instructions, or in a different storage address from the instructions.

[0067] The processor 610 is configured to invoke the instructions stored in the memory 620 to implement the steps of various exemplary embodiments according to the present disclosure described in the above "Exemplary Methods" section of the present specification. For example, the processor 610 can perform each step of the above-described seismic wave dispersion and attenuation prediction method embodiments.

[0068] It should be noted that the above processor 610 can be a general-purpose processor or a special-purpose processor. The processor 610 can include one or more processing cores, and the processor 610 performs various functional applications and data processing by running instructions.

[0069] In some embodiments, the processor 610 can include a central processing unit (CPU) and / or a baseband processor.

[0070] ​In some embodiments, the processor 610 can determine one instruction according to the priority identification and / or the function category information carried in the respective control instructions.

[0071] The processor 610 and the memory 620 in the present disclosure can be separately arranged or integrated together. As an example, the processor 610 and the memory 620 can be integrated on a single board or a system on chip (SOC).

[0072] As shown in Figure 6 The electronic device 600 is in the form of a general computing device. The electronic device 600 can also include a bus 630.

[0073] The bus 630 can be in the form of one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, a graphics acceleration port, a processor or a local bus using any of a variety of bus structures.

[0074] The electronic device 600 can also communicate with one or more external devices 640 (such as a keyboard or pointing device, a Bluetooth device, etc.) that can be used to interact with the electronic device 600, and / or one or more devices that enable a user to interact with the electronic device 600 and / or any devices (such as a router, a modem, etc.) that enable the electronic device 600 to communicate with one or more other computing devices. Such communication can occur via an input / output (I / O) interface 650.

[0075] In addition, the electronic device 600 can communicate with one or more networks (such as a local area network (LAN), a wide area network (WAN), and / or a public network, such as the Internet) via a network adapter 660.

[0076] As shown in Figure 6 The network adapter 660 communicates with other modules of the electronic device 600 via the bus 630.

[0077] It should be understood that although not shown in the figure, other hardware and / or software modules can be used in connection with the electronic device 600, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data archival storage systems, etc.

[0078] It can be understood that the structure shown in the embodiments of the present disclosure does not constitute a specific limitation on the electronic device 600. In other embodiments of the present disclosure, the electronic device 600 can include more or fewer components than those shown in Figure 6 or combine certain components, or split certain components, or different arrangement of components. Figure 6 The components shown can be implemented in hardware, software, or a combination of software and hardware.

[0079] The present disclosure also provides a computer readable storage medium, having stored thereon computer instructions, which when executed by a processor implement the method for predicting seismic wave dispersion and attenuation described in the above method embodiments.

[0080] The computer readable storage medium in the embodiments of the present disclosure is a computer readable medium that can send, propagate or transfer a computer program for use by or in connection with the instruction execution system, apparatus or device. As an example, the computer readable storage medium is a non-volatile storage medium.

[0081] In some embodiments, more specific examples of the computer readable storage medium in the present disclosure can include, but are not limited to, an electrical connection with one or more wires, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, a U disk, a mobile hard disk or any appropriate combination of the above.

[0082] In the embodiments of the present disclosure, the computer readable storage medium can include a data signal propagating in the baseband or as a part of a carrier wave, in which computer instructions (readable program code) are carried.

[0083] In some examples, the computer instructions contained on the computer readable storage medium can be transmitted by any appropriate medium, including but not limited to wireless, wired, optical fiber, RF, etc., or any appropriate combination of the above.

[0084] The embodiments of the present disclosure also provide a computer program product, which stores instructions, which when executed by a computer, cause the computer to implement the method for predicting seismic wave dispersion and attenuation described in the above method embodiments. The above instructions can be program codes. In specific implementation, the program codes can be written in any combination of one or more programming languages. The program codes can be executed entirely on a user computing device, partially on a user device, as an independent software package, partially on a user computing device and partially on a remote computing device, or entirely on a remote computing device or server.

[0085] The embodiments of the present disclosure also provide a chip, which includes at least one processor and an interface; the interface is used to provide program instructions or data for the at least one processor; the at least one processor is used to execute the program instructions to implement the method for predicting seismic wave dispersion and attenuation described in the above method embodiments.

[0086] In some embodiments, the chip can further include a memory for holding program instructions and data, the memory being located within the processor or external to the processor.

[0087] Those of skill in the art would understand that information and signals can be represented using any of a variety of different technologies and techniques. For example, data, instructions, commands, information, signals, bits, symbols, and chips that can be referenced throughout the above description can be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, optical fields or particles, or any combination thereof.

[0088] Those of skill would further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the disclosure herein can be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components, blocks, modules, circuits, and steps have been described above generally in terms of their functionality, whereas implementation of the combination of functional operations into hardware and software

[0089] It is intended to encompass any variations, uses, or adaptations of the disclosure following, in general, the principles of the disclosure and including such departures from the present disclosure that come within known The description and examples are to be regarded as illustrative only, and not restrictive, of the scope and spirit of the disclosure, which is indicated by the appended claims.

Claims

1. A method for predicting seismic wave dispersion and attenuation, characterized in that: include: Based on the bubble vibration theory and the pore fluid flow continuity equation, a target elastic modulus is obtained, wherein the target elastic modulus is the elastic modulus after adding bubble vibration correction; Integrating the target elastic modulus into the strain energy function of the classical poroacoustic elasticity theory to obtain a modified strain energy function; Combined with the strain Energy function, kinetic energy function and dissipated energy function are used to derive the wave equation of a fluid-saturated porous medium containing a small amount of bubbles under pressure; The wave equation is simplified by plane wave solution to obtain the seismic wave velocity and attenuation.

2. The method according to claim 1, characterized in that The target elastic modulus is obtained based on the bubble vibration theory and the pore fluid seepage continuity equation, including: Based on the dynamic response of bubbles under the action of sound waves, a nonlinear vibration equation describing the change of bubble radius with time is established; The nonlinear vibration equation is simplified under low-frequency conditions to obtain a linear relationship between bubble radius and pore pressure. The linear relationship between the bubble radius and the pore pressure is substituted into the corresponding pore fluid flow continuity equation to obtain the target elastic modulus.

3. The method according to claim 1, characterized in that The target elastic modulus is integrated into the strain energy function of the classical poroacoustic elasticity theory, including: The strain energy function form of the classical poroacoustic elasticity theory is retained, but the elastic modulus therein is replaced by the target elastic modulus after adding the bubble vibration correction.

4. The method according to claim 1, wherein Combining the strain energy function, kinetic energy function, and dissipated energy function, the wave equation for a fluid-saturated porous medium containing a small amount of bubbles under pressure is derived, including: The strain energy function, the kinetic energy function and the dissipated energy function are introduced into the generalized Lagrangian equation, and the solid displacement and fluid relative displacement variables are differentiated respectively to obtain the wave equation of the fluid-saturated porous medium containing a small amount of bubbles under the pressure.

5. The method according to claim 1, wherein The wave equation is simplified by plane wave solution to obtain seismic wave velocity and attenuation, including: Set the seismic wave disturbance to the form of simple harmonic plane wave; Substitute the plane wave solution into the separated wave equation and decouple the equation to obtain the complex wave number; Based on the complex wave number, seismic wave velocity and attenuation are obtained.

6. The method according to any one of claims 1 to 5, characterized in that: The target elastic modulus is calculated by the following formula: ； ； ； ； ； ； in, represents the target elastic modulus, represents the elastic modulus before correction, represents the density of the liquid, represents the pore radius, represents the number of bubbles per unit volume, represents the initial radius of the bubble, represents the internal pressure of the bubble under steady-state conditions, represents the static pressure of the liquid, represents the surface tension coefficient of the gas-liquid interface, represents the thermal diffusivity of the gas, represents the specific heat ratio of the gas, represents the viscosity coefficient of the pore fluid, represents the angular frequency, Represents the total porosity.

7. A seismic wave dispersion and attenuation prediction device, characterized in that: include: An elastic modulus correction module is used to obtain a target elastic modulus based on bubble vibration theory and the continuity equation of pore fluid flow. The target elastic modulus is the elastic modulus after bubble vibration correction is added. a strain energy correction module, configured to integrate the target elastic modulus into the strain energy function of the classical poroacoustic elasticity theory to obtain a corrected strain energy function; A wave equation derivation module is used to combine the strain energy function, kinetic energy function and dissipated energy function to derive a wave equation for a fluid-saturated porous medium containing a small amount of bubbles under pressure; The prediction module is used to simplify the wave equation by solving the plane wave to obtain the seismic wave velocity and attenuation.

8. An electronic device, characterized in that: include: a memory for storing instructions; A processor is used to call the instructions stored in the memory to implement the seismic wave dispersion and attenuation prediction method according to any one of claims 1 to 6.

9. A computer-readable storage medium having computer instructions stored thereon, characterized in that: When the computer instructions are executed by a processor, the seismic wave dispersion and attenuation prediction method according to any one of claims 1 to 6 is implemented.

10. A computer program product, characterized in that The computer program product stores instructions, and when the instructions are executed by a computer, the computer implements the seismic wave dispersion and attenuation prediction method according to any one of claims 1 to 6.

Images