Novel three-dimensional earthquake travel time simulation method based on PIUFNO

Abstract

The invention is suitable for the technical field of seismic exploration, and provides a PIUFNO-based novel three-dimensional seismic travel time simulation method, which comprises the following specific steps: step 1, establishing a PIUFNO network; step 2, determining the input and the output of the PIUFNO; 3, determining a loss function of PIUFNO, and carrying out network training and result output; and 4, inputting a new model after the network training is finished, and carrying out earthquake travel time simulation. According to the method, the problems of difficulty in three-dimensional model travel time simulation and low efficiency caused by gridding limitation of a traditional numerical method are effectively solved. According to the method, the three-dimensional earthquake travel time of multiple seismic sources can be quickly obtained by inputting different speed models and corresponding background travel time after one-time training, so that the calculation cost is remarkably reduced, and more efficient and practical three-dimensional earthquake travel time simulation is realized. The improvement of the efficiency is of great significance in enhancing the overall effect of three-dimensional seismic exploration and practical application in the fields of oil and gas resource development and the like.

Description

Technical Field

[0001] The invention belongs to the technical field of seismic exploration, and in particular relates to a novel three-dimensional seismic travel time simulation method based on PIUFNO. Background Art

[0002] Oil and gas resources are strategic resources that are essential for national development, so oil and gas exploration is an important task that affects the national economy and people's livelihood. By using the properties of seismic wave propagation, geophysical methods such as seismic imaging can be used to detect the internal structure of the earth. Seismic travel time is an important information transmitted by seismic waves, and seismic travel time simulation has always been an important topic of cross-scale research in geophysics.

[0003] The basis of earthquake travel time simulation is to solve the Eikonal equation, which is a nonlinear first-order partial differential equation. There are two most commonly used numerical solution methods: the fast march method (FMM) and the fast scanning method (FSM). Both methods are based on finite difference approximation and require the underground model to be gridded. In addition, numerical algorithms have theoretical limitations and can usually only perform travel time simulation on a single model and a fixed source. With the increase of dimensions, especially when performing three-dimensional travel time simulation, the computational complexity and cost increase significantly, making it more difficult to achieve efficient and accurate three-dimensional simulation. Therefore, a new three-dimensional earthquake travel time simulation method based on physical information Unet Fourier neural operator (PIUFNO) is proposed. Summary of the invention

[0004] The purpose of the embodiment of the present invention is to provide a novel 3D seismic travel time simulation method based on PIUFNO, aiming to solve the problems raised in the above background technology.

[0005] The embodiment of the present invention is implemented as follows: a novel 3D seismic travel time simulation method based on PIUFNO can obtain seismic travel time information by solving the Eikonal equation in an isotropic medium. The Eikonal equation can be expressed as: (1) in, T It represents the travel time of the first arrival wave of the earthquake. v represents the velocity model, represents the gradient operator, and are the coordinates of the model and the source, respectively. represents the time zero at the earthquake source. In order to avoid the singularity of the earthquake source, T Decompose into additive form T = T 0+ τ , τ is the factorization time factor, T0 is the background travel time corresponding to the velocity value at the source position, and the expression is: (2) in, The seismic waves propagate from Location to The distance is the velocity at the earthquake source. T = T 0+ τ Substituting into formula 1, the factorized form of the eikonal equation can be expressed as: (3) in, It means that the factorization factor at the earthquake source is zero.

[0006] The method comprises the following specific steps: Step 1. Establish PIUFNO network; Step 2: Determine the input and output of PIUFNO; Step 3: Determine the loss function of PIUFNO and perform network training; Step 4: After the network training is completed, input the new model and perform earthquake travel time simulation.

[0007] In a further technical solution, in step 1, PIUFNO converts the velocity model v And background time T 0 as input, the corresponding disturbance travel time τ = T - T 0 as output, and the mapping between them is established through the Unet Fourier Neural Operator (UFNO); In local transformation, and are the input and output of the current module, input After forward Fourier transform, linear high-frequency filtering, inverse Fourier transform, U-net feature transform and local linear transform, the output is finally obtained. .

[0008] In a further technical solution, in step 2, the input includes a velocity model v And background time T 0, speed model v Provides speed information, background travel time T 0 provides the source location information, and directly outputs the corresponding model and the travel time disturbance of the corresponding source after completing PIUFNO training τ .

[0009] In a further technical solution, in step 3, the physical information loss function is as follows: (4) in, N v is the number of velocity models used, iv is the velocity model index value; A further technical solution is that in step 4, after the network training is completed, different velocity models v and corresponding background travel times are input. T 0, output the travel time disturbance of the corresponding model and the corresponding source τ , and finally after addition T = T 0+ τ Get the simulated travel time results.

[0010] The embodiment of the present invention provides a novel 3D seismic travel time simulation method based on PIUFNO, which effectively solves the problem of low efficiency of traditional numerical methods that they can only simulate one model at a time. The present invention can realize 3D seismic travel time simulation of multiple earthquake sources through one training, and is applicable to different velocity models. This method significantly reduces the computational cost and improves the efficiency and practicality of 3D seismic travel time simulation. When testing typical geological velocity models, the method shows excellent performance and can accurately simulate 3D seismic travel time at any earthquake source position. BRIEF DESCRIPTION OF THE DRAWINGS

[0011] Figure 1 It is a 3D travel time simulation framework based on PIUFNO; Figure 2 It is the Unet Fourier neural operator module; Figure 3 The speed model for network input (a and b are two different speed models); Figure 4 is the background travel time input to the network (a and b are the background travel times corresponding to the two velocity models respectively); Figure 5 The velocity model and background travel time input after training and the travel time predicted by PIUFNO simulation (a is the velocity model, b is the corresponding background travel time T 0, c is the travel time simulated by PIUFNO). DETAILED DESCRIPTION

[0012] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0013] The specific implementation of the present invention is described in detail below in conjunction with specific embodiments.

[0014] An embodiment of the present invention provides a novel 3D seismic travel time simulation method based on PIUFNO. In an isotropic medium, seismic travel time information can be obtained by solving the Eikonal equation. The Eikonal equation can be expressed as: (1) in, T It represents the travel time of the first arrival wave of the earthquake. v represents the velocity model, represents the gradient operator, and are the coordinates of the model and the source, respectively. Represents time zero at the earthquake source.

[0015] In order to avoid the singularity of the earthquake source, T Decomposition into multiplication form T = T 0+ τ , τ is the factorization time factor, T 0 is the background travel time corresponding to the velocity value at the earthquake source position, and the expression is: (2) in, The seismic waves propagate from Location to The distance is the velocity at the earthquake source. T = T 0+ τ Substituting into formula 1, the factorized form of the eikonal equation can be expressed as: (3) in, It means that the factorization factor at the earthquake source is zero.

[0016] The method comprises the following specific steps: Step 1. Establish PIUFNO network; The established PIUFNO network is shown in the figure below: Figure 1 As shown, PIUFNO will speed model v And background time T 0 as input, the corresponding disturbance travel time τ = T - T 0 as output, and the mapping between them is established through the Unet Fourier Neural Operator (UFNO).

[0017] The core of PIUFNO is the UFNO block, such as Figure 2 As shown, represents forward fast Fourier transform; represents inverse Fourier transform; Indicates that linear transformation filters out high-frequency components; U Represents U-net feature transformation; represents a local linear transformation; Represents the Relu activation function; Represents input After the three transformations the results are added together. and are the input and output of the current module, input After forward Fourier transform, linear high-frequency filtering, inverse Fourier transform, U-net feature transform and local linear transform, the output is finally obtained. .

[0018] Step 2: Determine the input and output of PIUFNO; The purpose is to use PIUFNO to simulate the travel time of multiple earthquake sources and arbitrary velocity models at the same time, and to establish input [ v ; T 0] and output τ The mapping relationship between v Provides speed information, T 0 provides the source location information. After completing PIUFNO training, the corresponding model and the travel time disturbance of the corresponding source can be directly output. τ , and finally after a simple addition operation T = T 0+ τ The final travel time result can be obtained.

[0019] Step 3: Determine the loss function of PIUFNO and perform network training; The deep learning framework based on physical information does not need to generate a large amount of training data and manual labels, which reduces the cost of data generation and enhances the interpretability of the deep learning framework. The physical information loss function is as follows: (4) in, N v is the number of velocity models used, and iv is the velocity model index value.

[0020] Step 4: After the network training is completed, input the new model and perform earthquake travel time simulation.

[0021] After the training is completed, by inputting different speed models v and the corresponding background travel time T0, output the travel time disturbance of the corresponding model and the corresponding source τ , and finally after the addition operation T = T 0+τ to get the simulated travel time result.

[0022] To verify the effectiveness of this method, Figure 3 a and Figure 3 b shows the two velocity models and the background travel times corresponding to the two velocity models T 0 (respectively Figure 4 a and 4b) are input into the network for training. This model contains 60×60×60 grid points, the spatial sampling interval is 25 m, the number of training times is 500, and the Adam optimizer is used for network optimization training. The learning rate is 0.0001, and after 250 training times, the learning rate becomes half of the original.

[0023] By training PIUFNO, we can predict the travel time of earthquakes with different source locations and different velocity models. Figure 5 The velocity model shown in a and Figure 5 The corresponding background travel time is shown in b T 0 is input into the network, and the source position is set at the center (0.75, 0.75, 0.75). The travel time simulated by PIUFNO is as follows Figure 5 As shown in c.

[0024] According to experiments, the time required for this method to simulate the travel time using the displayed model is 0.054 seconds, while the traditional numerical method takes 0.095 seconds to simulate the travel time of the same model, which greatly improves the overall calculation efficiency. This efficiency improvement is of great significance to enhancing the overall effect of 3D seismic exploration and its practical application in fields such as oil and gas resource development.

[0025] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A new 3D seismic travel time simulation method based on PIUFNO, characterized in that: The specific steps include: Step 1. Establish PIUFNO network; Step 2: Determine the input and output of PIUFNO; Step 3: Determine the loss function of PIUFNO and perform network training; Step 4: After the network training is completed, input the new model and perform earthquake travel time simulation.

2. The novel 3D seismic travel time simulation method based on PIUFNO according to claim 1, characterized in that: In step 1, PIUFNO converts the velocity model v And background time T 0 as input, the corresponding disturbance travel time τ = T - T 0 as output, and the mapping between them is established through UFNO; In local transformation, and are the input and output of the current module, input After forward Fourier transform, linear high-frequency filtering, inverse Fourier transform, U-net feature transform and local linear transform, the output is finally obtained. .

3. The novel 3D seismic travel time simulation method based on PIUFNO according to claim 2 is characterized in that: In step 2, the input includes the velocity model v And background time T 0, speed model v Provides speed information, background travel time T 0 provides the source location information, and directly outputs the corresponding model and the travel time disturbance of the corresponding source after completing PIUFNO training τ .

4. The novel 3D seismic travel time simulation method based on PIUFNO according to claim 3 is characterized in that: In step 3, the physical information loss function is as follows: ； in, N v is the number of velocity models used, and iv is the velocity model index value.

5. The novel 3D seismic travel time simulation method based on PIUFNO according to claim 4 is characterized in that: In step 4, after the training is completed, the new velocity model v and the background travel time are input. T 0, output the travel time disturbance of the corresponding model and the corresponding source τ , and finally after addition operation T = T 0+ τ Get the final simulated travel time result.

Images