METHOD FOR THREE-DIMENSIONAL VELOCITY GEOLOGICAL MODELING WITH STRUCTURES AND VELOCITIES RANDOMLY ARRANGED
20230384470 · 2023-11-30
Assignee
Inventors
- Peng JIANG (Jinan, CN)
- Yuxiao REN (Jinan, CN)
- Qifeng WANG (Jinan, CN)
- Zhiwu ZUO (Jinan, CN)
- Xinji Xu (JiNan, CN)
- Kai WANG (Jinan, CN)
- Lei CHEN (Jinan, CN)
- Chuanyi MA (Jinan, CN)
- Shuai CAO (Jinan, CN)
- Senlin YANG (Jinan, CN)
- Qingyang WANG (Jinan, CN)
- Xianglong MENG (Jinan, CN)
Cpc classification
International classification
Abstract
A method for three-dimensional velocity geological modeling with structures and velocities randomly arranged, including determining base points in three-dimensional space, building equation according to the base points to determine planar layered model, complicating a tilt layer of planar layered model, and building a fold layer model of a surface in three-dimensional space; building three-dimensional fault folded model based on the three-dimensional surface fold layer model combined with a fault plane of a random reference point and displacement of each point in a global coordinate system; building a velocity model containing a salt body based on the three-dimensional fault folded model, and simulating salt body intrusion in a geological body of a certain depth; and performing a random velocity amplitude to realize three-dimensional velocity modeling according to the layered type which has been set and according to the set velocity range and the velocity difference range between each layer of geology.
Claims
1. A method for three-dimensional velocity geological modeling with structures and velocities randomly arranged, comprising: determining a base point in a three-dimensional space, building an equation according to the base point to determine a planar layered model, and complicating a tilt layer of the planar layered model to build a fold layer model of a curved surface in the three-dimensional space; building a three-dimensional fault folded model based on the three-dimensional surface fold layer model combined with a fault plane of a random reference point and displacement of each point in a global coordinate system; building a velocity model containing salt body based on the three-dimensional fault folded model, and simulating salt body intrusion in a geological body of certain depth; and performing a random velocity amplitude to realize three-dimensional velocity modeling according to the layered type which has been set and according to the set velocity range and the velocity difference range between each layer of geology; the complicating a tilt layer of the planar layered model specifically comprising: the planar layered model is determined by building an equation according to the base points, and the different layer models are divided into different categories; a fluctuation function is built for each point based on the plane model; by adjusting a period and amplitude of a trigonometric function in the fluctuation function, a tilt term is built for the surface, the tilt layer is further complicated, and the fold layer model of the surface in the three-dimensional space is built.
2. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1, wherein the specific process of building an equation to determine a planar layered model according to the base points comprises: an equation is built according to the base points (X.sub.ref,Y.sub.ref,Z.sub.ref) to determine the planar layered model, and the calculation formula is as follows:
H(X,Y)=(X−X.sub.ref)+(Y−Y.sub.ref)tan φ wherein φ represents an angle of tilt.
3. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1, wherein the specific process comprises: a folded model is built based on the planar layered model, and the calculation formula is as follows:
D(X,Y)=b.sub.1(X−X.sub.ref)+b.sub.2(Y−Y.sub.ref), wherein X.sub.refY.sub.ref are base point coordinates, and the value of b.sub.1b.sub.2 is randomly taken.
4. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1, wherein the specific process of building a three-dimensional fault folded model comprises: the formula for adding faults to the folded model is as follows:
c.sub.1(X−X.sub.ref)+c.sub.2(Y−Y.sub.ref)+c.sub.3(Z−Z.sub.ref)=0, wherein c.sub.1c.sub.2c.sub.3 are calculated by a rotation matrix:
5. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1, wherein the specific process of simulating an upward salt body intrusion in a geological body of a certain depth comprises: fitting the intrusion by a two-dimensional Gaussian function, defining the height of vertical intrusion by an amplitude, determining the size by variances, and the direction determining by a clockwise rotation angle; setting an affected area with a certain thickness, the maximum intrusion height being at the bottom stratum, the closer to the surface in the affected area, the smaller the influence, while the stratum above the affected area remains unchanged, completing the addition of a salt body.
6. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 5, wherein the formula for building the salt body is as follows:
G(X,Y)=A exp(−(d.sub.1(X−X.sub.ref).sup.2+d.sub.3(Y−Y.sub.ref).sup.2+2d.sub.2(X−X.sub.ref)(Y−Y.sub.ref))), wherein
7. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1, wherein the specific process of performing a random velocity amplitude to achieve three-dimensional velocity modeling comprises: generating a vector V′ of n+1 elements randomly according to the number of layers n; generating a vector V.sub.1 by accumulating the vector V′; taking a random velocity reference value M, M∈[x.sub.1,x.sub.2] and x.sub.1,x.sub.2 as the upper and lower bounds of velocity; taking the last element v.sub.endV.sub.1; evaluating the velocity as V=(V.sub.1|v.sub.end).Math.M; and taking velocity V.sub.salt body, V.sub.salt body∈[M,M+Δv] randomly for the salt body, wherein Δv is a randomly added velocity value.
8. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1, wherein in the process of three-dimensional velocity modeling, the earth surface folded model is rotated by 90° counterclockwise along a central axis to determine fault strike, stratum thickness and velocity distribution range planned according to geological conditions from a geological survey report before tunneling, and the weight of the modeling parameters in different ranges is set.
9. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1, wherein the method for three-dimensional velocity geological modeling with structures and velocities randomly arranged further comprises acquiring a tunnel seismic record, performing feature extraction processing on the tunnel seismic record using a convolutional neural network, and adding tunnel receiver position information on an additional channel to complete data encoding.
10. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 9, wherein a convolutional neural network is used to decode the encoded data, multi-objective learning is performed, and a three-dimensional velocity model and three-dimensional velocity modeling parameters are obtained by respectively processing the decoding results of the decoder using a convolutional neural network and a fully connected neural network.
11. The method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1, wherein in the process of building a velocity model, the loss function used comprises a velocity model loss function and a modeling parameter loss function, wherein the velocity model loss function is used for fitting a real three-dimensional velocity model corresponding to the observation data and a network modeling three-dimensional velocity model; the modeling parameter loss function is used to fit the real three-dimensional velocity model modeling parameters corresponding to the observation data with the network modeling parameters.
12. A computer-readable storage medium, comprising: the computer-readable storage medium storing therein a plurality of instructions adapted to be loaded by a processor of a terminal device and to perform the method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1.
13. A terminal device, comprising: a processor and a computer-readable storage medium, the processor using for implementing a plurality of instructions; the computer-readable storage medium using for storing a plurality of instructions adapted to be loaded by a processor and to perform the method for three-dimensional velocity geological modeling with structures and velocities randomly arranged according to claim 1.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0062] The accompanying drawings constituting a part of the present disclosure are used to provide further understanding of the present disclosure. Exemplary embodiments of the present disclosure and descriptions thereof are used to explain the present disclosure, and do not constitute an improper limitation to the present disclosure.
[0063]
[0064]
[0065]
[0066]
[0067]
[0068]
[0069]
[0070]
[0071]
DETAILED DESCRIPTION
[0072] The present disclosure is further described below with reference to the accompanying drawings and embodiments.
[0073] It should be noted that the following detailed descriptions are all exemplary and are intended to provide a further description of the present disclosure. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by a person of ordinary skill in the technical field to which the present disclosure belongs.
[0074] It should be noted that terms used herein are only for describing specific implementations and are not intended to limit exemplary implementations according to the present disclosure. As used herein, the singular form is intended to include the plural form, unless the context clearly indicates otherwise. In addition, it should further be understood that terms “include” and/or “include” used in this specification indicate that there are features, steps, operations, devices, components, and/or combinations thereof.
[0075] As shown in
[0078] In this example, the planar layered model is determined by building an equation according to the base points (X.sub.ref,Y.sub.ref,Z.sub.ref), and the calculation formula is as follows:
H(X,Y)=(X−X.sub.ref)+(Y−Y.sub.ref)tan φ [0079] in this example, a folded model is built based on the planar layered model, and the calculation formula is as follows:
D(X,Y)=b.sub.1(X−X.sub.ref)+b.sub.2(Y−Y.sub.ref). [0081] Step S2: Based on the built three-dimensional surface fold layer model, building an equation via a determined random reference point, determining a fault plane passing through the reference point, and then determining displacement of each point in a global coordinate system via a rotation matrix, and building a three-dimensional fault folded model, as shown in
c.sub.1(X−X.sub.ref)+c.sub.2(Y−Y.sub.ref)+c.sub.3(Z−Z.sub.ref)=0,
where
[0083] c.sub.1c.sub.2c.sub.3 are calculated by a rotation matrix:
a rotation
[0084] matrix
where Φθ are randomly taken from [0,2π]. [0085] Step S3: Building a velocity model containing a salt body based on a folded model, as shown in
[0086] The formula for building the salt body in this example is as follows:
G(X,Y)=A exp(−(d.sub.1(X−X.sub.ref).sup.2+d.sub.3(Y−Y.sub.ref).sup.2+2d.sub.2(X−X.sub.ref)(Y−Y.sub.ref))),
where
[0087] A represent the height of a vertical intrusion of the salt body, and the size of the salt body is controlled by σ.sub.X.sup.2σ.sub.Y.sup.2; an affected area of the salt body is set as [A.sub.max+5,A.sub.max+15], where A.sub.max represents the maximum intrusion height; in the affected area, the shallower the layer, the smaller the amplitude A of the corresponding Gaussian function, and the stratum above the affected area remains unchanged. As shown in
[0089] The velocity evaluating procedure is as follows: [0090] 1. generating a vector V′ of n+1 elements randomly according to the number of layers n; [0091] 2. generating a vector V.sub.1 by accumulating the vector V′; [0092] 3. taking a random velocity reference value M, M∈[x.sub.1,x.sub.2] (x.sub.1x.sub.2 is the upper and lower bounds of velocity, which is taken as M∈[2000 m/s,4000 m/s] in this example) [0093] 4. taking the last element V.sub.1v.sub.end; [0094] 5 evaluating the velocity as V=(V.sub.1/v.sub.end).Math.M; [0095] 6. taking velocity V.sub.saltbody, V.sub.saltbody∈[M,M=Δv] randomly for the salt body, where Δv is a randomly added velocity value, which is taken as 300 m/s to 500 m/s in this example.
[0096] Based on the above solution, it is possible to automatically build a three-dimensional velocity model. Further, the method can be extended to a method for building a tunnel seismic velocity model based on deep learning, which includes: [0097] a tunnel velocity model database building module configured to generate a three-dimensional velocity model in front of a tunnel based on an on-site geological exploration report in a large number and constitute a tunnel velocity model database; [0098] a finite difference forward modeling module: the acoustic wave fluctuation equation can be represented as follows:
[0099] A finite difference method is used to forward simulate the elastic wave fluctuation equation, and the receivers arranged on the tunnel wall are used to receive the amplitude information of seismic waves for further deep learning modeling.
[0100] A three-dimensional seismic velocity module built based on deep learning is configured to build a tunnel inversion deep neural network, where an input of the network is seismic observation data under a three-dimensional observation mode, and an output is a predicted velocity model.
[0101] Based on the above three-dimensional velocity model, the folded model can be used as a tunnel three-dimensional velocity model by rotating 90° counterclockwise along the central axis.
[0102] The tunnel model of this embodiment is shown in
[0105] In other examples, the geological model may be built by other such parameters. Receivers and sources may be substituted.
[0106] A geological velocity model in this example database is shown in
[0107] The batch database built in this example includes 10000 tunnel seismic velocity models, which are randomly divided into a training set, a validation set and a test set by a ratio of 8:1:1. [0108] Step S3: A tunnel velocity model built neural network is built, where an input of the network is seismic observation data, the network is input according to Batchsize=8, the size of the input data is [8, 1, 12, 2000], and the output is a three-dimensional tunnel velocity model and three-dimensional velocity modeling parameters; the entire neural network includes one encoder and one decoder. Seismic observation information is input into an encoder module which is a global feature encoder composed of a six-layer convolutional neural network, where the encoder converts an input into a feature vector containing velocity model information, and the size of the output feature vector is [8, 63, 12, 6]; a channel is added to the feature vector, and a receiver corresponding to seismic data and source information is added, i.e., receiver position information (x.sub.n,y.sub.n,z.sub.n,x,y,z), where x.sub.n,y.sub.n and z.sub.n represent, the coordinate of the n.sup.th source point, and x, y and z represent the coordinate of the geophone. The size of the processed feature vector is [8, 64, 12, 6], and this vector is input into a decoder for decoding; a decoder adding position encoding information includes a five-layer convolutional layer followed by two-layer convolutional layers and a three-layer full-connection layer respectively, and a velocity model and modeling parameters are respectively output.
[0109] Regarding the velocity model, after obtaining the corresponding velocity model of the data, the velocity model built by the neural network is compared with the velocity model corresponding to the input seismic data, and the least squares loss function is used to calculate the difference between the two velocity models, and the network is transmitted back to optimize the network parameters.
[0110] Regarding the modeling parameters, after obtaining the modeling parameters, the modeling parameters built by the neural network are compared with the modeling parameters corresponding to the input seismic data, and the differences are respectively calculated by using the least squares loss function and the minimum absolute error loss function, and the sum of the two loss functions is added and the gradient return transmission is performed.
[0111] Because there are two outputs to optimize the network, in the process of network learning, the weight of the loss function generated by the two outputs is adjusted. In the initial stage of network training, the weights of the modeling parameters are larger and the weights of the velocity model are smaller, and the weights of the velocity model gradually increase as the training progresses.
[0112] The least squares loss function can be represented as
L.sub.m=||m.sub.est−m.sub.tru||.sup.2
[0113] The minimum absolute error loss function can be represented as:
L.sub.m=||m.sub.est−m.sub.tru||.sup.1
[0114] m.sub.est represents a model wave of the model output by the deep neural network for modeling the velocity of the tunnel, and m.sub.tru represents an actual geological model in the database of the velocity model of the tunnel.
[0115] In the network training process of S3 stage, an Adam optimizer is used, the learning rate is kept constant by 5×10.sup.−5, the BatchSize of network training stage is 8, and the total number of iterations is 150.
[0116] In Step S4, the trained funnel velocity constructed neural network is used to test the velocity constructing effect on a test set, and the test result is shown in
[0117] The primary network parameters and hardware conditions in this embodiment are: calculations are performed using NVIDIARTX3090 GPU*4. Algorithms are written using MATLAB and PYTHON.
[0118] The following product examples are further provided: [0119] a computer-readable storage medium having stored therein a plurality of instructions adapted to be loaded by a processor of a terminal device and to perform the method for three-dimensional velocity geological modeling with structures and velocities randomly arranged.
[0120] A terminal device including a processor and a computer-readable storage medium, where the processor is used for implementing various instructions; the computer-readable storage medium is used to store a plurality of instructions adapted to be loaded by a processor and to perform the method for three-dimensional velocity geological modeling with structures and velocities randomly arranged.
[0121] A person skilled in the art should understand that the embodiments of the present disclosure may be provided as a method, a system, or a computer program product. Therefore, the present disclosure may use a form of hardware-only embodiments, software-only embodiments, or embodiments combining software and hardware. in addition, the present disclosure may use a form of a computer program product implemented on one or more computer-usable storage media (including but not limited to a disk memory, a CD-ROM, an optical memory, and the like) including computer-usable program code.
[0122] The present disclosure is described with reference to flowcharts and/or block diagrams of the method, device (system), and computer program product in the embodiments of the present disclosure. Computer program instructions can implement each procedure and/or block in the flowcharts and/or block diagrams and a combination of procedures and/or blocks in the flowcharts and/or block diagrams. These computer program instructions may be provided to a general-purpose computer, a special-purpose computer, an embedded processor, or a processor of another programmable data processing device to generate a machine, so that an apparatus configured to implement functions specified in one or more procedures in the flowcharts and/or one or more blocks in the block diagrams is generated by using instructions executed by the computer or the processor of another programmable data processing device.
[0123] These computer program instructions may alternatively be stored in a computer-readable memory that can instruct a computer or another programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory generate an artifact that includes an instruction apparatus. The instruction apparatus implements a specific function in one or more procedures in the flowcharts and/or in one or more blocks in the block diagrams.
[0124] These computer program instructions may further be loaded onto a computer or another programmable data processing device, so that a series of operations and steps are performed on the computer or another programmable device, thereby generating computer-implemented processing. Therefore, the instructions executed on the computer or another programmable device provide steps for implementing a specific function in one or more processes in the flowcharts and/or in one or more blocks in the block diagrams.
[0125] The foregoing descriptions are merely exemplary embodiments of the present disclosure, but are not intended to limit the present disclosure. The present disclosure may include various modifications and changes for a person skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present disclosure shall fall within the protection scope of the present disclosure.