A random modeling method and device for deep strike-slip fault flower-shaped structure

Abstract

The present application provides a method and device for stochastic modeling of deep strike-slip fault flower-like structures, which relates to the technical fields of structural geology and mathematical geological modeling. The method includes: generating n three-dimensional fractal data bodies based on a fractal noise method; extracting two-dimensional slices along the Z direction from each three-dimensional fractal data body and adding reference depths to obtain multiple sedimentary stratigraphic sequences; adding corresponding dip angles to each sedimentary stratigraphic sequence; superimposing multiple sedimentary stratigraphic sequences to obtain a stratigraphic model; adding a main strike-slip fault, and generating the center point of the dip-slip fault and a mask; adding the same type of simulated faults at the center point of the dip-slip fault, and using a mask to shield the simulated faults that exceed the main strike-slip fault, to obtain a negative flower-like structural model or a positive flower-like structural model. The present invention can effectively generate a realistic strike-slip fault flower-like structural model, and provide training data for finely characterizing strike-slip faults and their flower-like structures using deep learning.

Description

Technical Field

[0001] The present application relates to the technical fields of structural geology and mathematical geological modeling, and in particular to a method and device for stochastic modeling of deep strike-slip fault flower-shaped structures. Background Art

[0002] Strike-slip faults are one of the important manifestations of tectonic movement. Their development process is complex and they are often accompanied by structural deformation features such as flower-shaped structures and fault breccias. Flower-shaped structures are a common structural phenomenon in strike-slip fault zones. They are complex and diverse in form and have significant scale differences. They have an important impact on oil and gas exploration and development, engineering construction and other fields. Therefore, in-depth research on the formation mechanism, spatial distribution characteristics and influence of flower-shaped structures on fluid migration is crucial to understanding the structural evolution history of strike-slip fault zones, predicting the distribution of oil and gas reservoirs and ensuring engineering safety.

[0003] With the rapid development of deep learning technology, the application of deep learning methods for automatic identification and modeling of geological structures has become a hot research field. Deep learning models, especially convolutional neural networks (CNNs), have shown strong recognition capabilities in image recognition and pattern recognition. However, the application of deep learning technology to the identification and modeling of strike-slip fault flower-like structures faces challenges such as data scarcity, inaccurate labels, and insufficient authenticity and diversity, which need to be overcome by artificially synthesizing fault models. In 2020, Wu et al. developed a workflow for automatically constructing a diversified structural model with fold and fault characteristics. The method made some assumptions about typical fold and fault patterns, and then used a set of parameters to simulate the structural features in the three-dimensional model. By randomly selecting parameters from some predefined ranges, a large number of fault models with realistic and diverse structural features can be automatically generated. Subsequent studies mostly use this process to synthesize fault models. Patent CN116338790A discloses a method, device, equipment and storage medium for synthesizing seismic fault training data, in which the synthetic fault model part adopts the workflow developed by Wu et al.

[0004] The above fault modeling methods can only generate folds and dip-slip faults including normal faults and reverse faults, but cannot simulate pinch-out phenomena and strike-slip faults and their flower-like structures. Therefore, a new fault modeling technology is urgently needed, which can effectively generate realistic strike-slip fault flower-like structures on the one hand, and provide high-quality simulation data for deep learning through advanced stochastic modeling technology on the other hand. Summary of the invention

[0005] The purpose of the present invention is to provide a random modeling method and device for deep strike-slip fault flower-like structures in order to solve the problem that the existing fault modeling method can only generate folds and dip-slip faults including normal faults and reverse faults, but cannot simulate pinch-out phenomena and strike-slip faults and their flower-like structures.

[0006] The above-mentioned purpose of the present application is achieved through the following technical solutions:

[0007] S1: Obtain a three-dimensional white noise data volume in the spatial domain; obtain a three-dimensional fractal data volume in the spatial domain by combining Fourier transform and inverse Fourier transform according to the three-dimensional white noise data volume in the spatial domain;

[0008] S2: extracting a preset number of two-dimensional slices along the Z direction from each three-dimensional fractal data volume in the spatial domain, each slice is regarded as a stratigraphic plane, and a reference depth is added to each slice as the depth of the stratigraphic plane. The two-dimensional slices extracted from the same three-dimensional fractal data volume in the spatial domain are regarded as a sedimentary stratigraphic sequence. Applying the above operation to each three-dimensional fractal data volume in the spatial domain obtains n sedimentary stratigraphic sequences;

[0009] S3: adding a dip angle to each sedimentary stratigraphic sequence to obtain n sedimentary stratigraphic sequences containing dip angles; obtaining a stratigraphic model by superimposing n sedimentary stratigraphic sequences containing dip angles;

[0010] S4: Add a main strike-slip fault to the stratigraphic model, and generate a preset number of dip-slip fault center points on the left and right sides of the strike-slip fault and masks of the left and right sides of the main strike-slip fault;

[0011] S5: Add the same type of simulated faults at the center point of the dip-slip fault, use a mask to shield the simulated faults that exceed the main strike-slip fault, and obtain a negative flower structure model or a positive flower structure model.

[0012] Optionally, step S1 includes:

[0013] S11: Based on the fractal noise method, n three-dimensional random number matrices are generated as three-dimensional white noise data volumes in n spatial domains;

[0014] S12: Using Fourier transform, the three-dimensional white noise data volume in the spatial domain is converted into a three-dimensional white noise data volume in the frequency domain. The specific formula is as follows:

[0015]

[0016] in, is the length of the three-dimensional white noise data volume in the spatial domain, is the width of the three-dimensional white noise data volume in the spatial domain, is the height of the three-dimensional white noise data volume in the spatial domain; Represents a three-dimensional white noise data volume in the frequency domain; is the spatial frequency in the X direction in the spatial domain, The range is ; is the spatial frequency in the Y direction of the spatial domain, It is the spatial domain The spatial frequency in the direction, The range is , The range is ; represents a three-dimensional white noise data volume in the spatial domain, Represents a three-dimensional white noise data volume in the spatial domain Axis coordinates; Represents a three-dimensional white noise data volume in the spatial domain Axis coordinates, Represents a three-dimensional white noise data volume in the spatial domain Axis coordinates; The range is , The range is , The range is ;

[0017] S13: Construct a frequency domain filter, the formula is:

[0018]

[0019] in, Represents a three-dimensional white noise data volume in the spatial domain The frequency in the direction, express Directional weight; Represents a three-dimensional white noise data volume in the spatial domain The frequency in the direction, express Directional weight; Represents a three-dimensional white noise data volume in the spatial domain The frequency in the direction, express Directional weight; fractal dimension Related constants ;

[0020] The frequency domain filter is multiplied by the three-dimensional white noise data volume in the frequency domain for filtering. The expression is as follows:

[0021]

[0022] in, Represents the original three-dimensional white noise data volume in the frequency domain, represents a three-dimensional white noise data volume in the frequency domain after filtering;

[0023] S14: Apply inverse Fourier transform to the filtered three-dimensional white noise data volume in the frequency domain to obtain n three-dimensional fractal data volumes in the space domain. The calculation formula is:

[0024]

[0025] in, represents the three-dimensional white noise data volume in the frequency domain after filtering, is the spatial frequency in the X direction in the spatial domain, is the spatial frequency in the Y direction of the spatial domain, It is the spatial domain The spatial frequency in the direction, The range is , The range is ; The range is ; represents the filtered spatial domain fractal data volume, Representing 3D fractal data in spatial domain Axis coordinates, Representing 3D fractal data in spatial domain Axis coordinates, Representing 3D fractal data in spatial domain Axis coordinates, The range is , The range is , The range is ; is the length of the three-dimensional fractal data volume in the spatial domain, is the width of the three-dimensional fractal data volume in the spatial domain, is the height of the three-dimensional fractal data volume in the spatial domain.

[0026] Optionally, step S2 includes:

[0027] A preset number of two-dimensional slices are extracted along the Z direction from n three-dimensional fractal data volumes in the spatial domain. Each slice is regarded as a layer plane. Each slice is added with a reference depth as the depth of the layer plane. The calculation formula is:

[0028]

[0029] in, Indicates the depth of the i-th layer extracted from the n-th three-dimensional fractal data volume, represents the i-th slice extracted from the n-th three-dimensional fractal data volume, Indicates the base depth of the i-th layer plane extracted from the n-th three-dimensional fractal data volume;

[0030] The two-dimensional slices extracted from the same three-dimensional fractal data volume in the spatial domain are taken as a sedimentary stratigraphic sequence, and n sedimentary stratigraphic sequences of n three-dimensional fractal data volumes in the spatial domain are obtained.

[0031] Optionally, step S3 includes:

[0032] S31: Add the dip angle to each sedimentary stratigraphic sequence respectively, and the calculation formula is:

[0033]

[0034]

[0035] in, , , The three-dimensional coordinates of the nth sedimentary stratigraphic sequence after adding the dip angle, , , is the original three-dimensional coordinate of the nth sedimentary stratigraphic sequence; is the displacement field, , , As parameters, and Used to control the formation and The slope in the direction, Used to control the offset distance of the center point;

[0036] S32: superimposing n sedimentary stratigraphic sequences with added dip angles in a preset order to obtain a stratigraphic model containing unconformity phenomena.

[0037] Optionally, step S4 includes:

[0038] S41: Generate a strike-slip fault center point on the horizon model , Steering angle And the inclination ;

[0039] S42: Taking a strike-slip fault center point is the origin of coordinates, the strike is the x-axis, the inclination is the y-axis, and the z-axis is the normal direction and points to the upper plate;

[0040] Establish a local coordinate system and transform all points of the layer model from original coordinates to local coordinates:

[0041]

[0042] in, is the x-axis coordinate in the local coordinate system, is the y-axis coordinate in the local coordinate system, is the z-axis coordinate in the local coordinate system, is the rotation matrix:

[0043]

[0044] The plane is the plane where the fault plane is located in the local coordinate system. In the current step, the fault plane is a strike-slip fault;

[0045] Specifically, what is added in step S4 is a strike-slip fault (main strike-slip fault), so this fault plane refers to the strike-slip fault. In step S5, when executing step S42, the fault plane refers to the dip-slip fault. According to the classification, faults can be divided into strike-slip faults and dip-slip faults.

[0046] S43: Calculate the displacement field of the fault plane based on the ellipse formula. The calculation formula is as follows:

[0047]

[0048] in Describe the fault displacement distribution on the fault plane, represents the maximum displacement of the center point, is the normalized radial distance of the center point, and the calculation formula is:

[0049]

[0050] in, is the radius along the strike direction, is the radius along the inclination direction, It is used to control the rotation angle of the displacement attenuation ellipse;

[0051] Calculate the strike displacement component separately and the tilt displacement component , the calculation formula is as follows:

[0052]

[0053]

[0054] S44: Fault plane in local coordinate system Generate a preset number of disturbance points within the preset range of :

[0055]

[0056] in, , and are the x-coordinate, y-coordinate and z-coordinate of the random point in the local coordinate system, is a random number; random points are disturbance points;

[0057] Use disturbance points to perturb the fault plane :

[0058]

[0059] S45: Calculate the three-dimensional fault displacement field on the fault plane. The calculation formula is as follows:

[0060]

[0061]

[0062]

[0063]

[0064] in, Indicates left disk; Indicates right plate; Indicates upper plate; Indicates the lower plate; , It represents the displacement ratio between the left and right disks; ( ) represents the displacement ratio of the upper plate and the lower plate; It is used to control the rotation angle of the displacement attenuation ellipse. When , the fault is a dip-slip fault. When , the fault is a strike-slip fault. When , the fault contains displacement in both strike and dip; , and They represent the fault displacement along the strike, dip and perpendicular to the fault plane in the local coordinate system respectively;

[0065] In the current step, , generating a strike-slip fault, i.e. a pure strike-slip fault;

[0066] S46: applying the three-dimensional fault displacement field in the local coordinate system to obtain a simulated fault;

[0067] The simulated fault is a simulated fault of a normal fault type or a simulated fault of a reverse fault type;

[0068] Convert the simulated fault back to the original coordinate system, and the calculation formula is as follows:

[0069]

[0070] in, is the rotation matrix The inverse matrix of , and They are the x-axis, y-axis, and z-axis coordinates of the original coordinate system after applying the three-dimensional fault displacement field in the local coordinate system. The model is recorded as , It is a strike-slip fault model;

[0071] S47: Select n points within the preset range of the strike-slip fault center point in the local coordinate system As the center point of the dip-slip fault, the left disk of the main strike-slip fault in the local coordinate system is generated. and right disk The mask is calculated as follows:

[0072]

[0073]

[0074] Will , and Convert to the original coordinate system and get , and , the calculation formula is as follows:

[0075]

[0076]

[0077] .

[0078] Optionally, step S5 includes:

[0079] S51: The center point of the dip-slip fault in is the center point of the simulated fault, and the dip angle is generated. and the angle of travel , Strike angle with the main strike-slip fault The calculation formula is as follows:

[0080]

[0081]

[0082] in, It indicates that the center point of the simulated fault is located on the left disk of the main strike-slip fault; It indicates that the center point of the simulated fault is located on the right disk of the main strike-slip fault;

[0083] S52: Generate the inclination angle by combining the calculation formulas from step S42 to step S46 and the angle of travel , generate a simulated fault. In this step, the simulated fault generated is a dip-slip fault;

[0084] Determine the type of simulated fault to be generated; the types of simulated faults include: normal fault simulated fault and reverse fault simulated fault;

[0085] Add the same type of simulated fault at the center of the dip-slip fault to obtain the Iteration Model, , The model adds Model of dip-slip faults;

[0086] S53: Use a mask to shield the simulated faults that exceed the main strike-slip fault. The calculation formula is as follows:

[0087]

[0088]

[0089] in, It indicates that the center point of the simulated fault is located on the left disk of the main strike-slip fault; It indicates that the center point of the simulated fault is located on the right disk of the main strike-slip fault;

[0090] S54: Repeat steps S51 to S53 until Each point in the model corresponds to a simulated fault, resulting in a positive structural model or a negative structural model.

[0091] An electronic device comprises a processor, a memory, a user interface and a network interface, wherein the memory is used to store instructions, the user interface and the network interface are used to communicate with other devices, and the processor is used to execute the instructions stored in the memory so that the electronic device executes a random modeling method for deep strike-slip fault flower-like structures.

[0092] A computer-readable storage medium stores instructions, and when the instructions are executed, a random modeling method of deep strike-slip fault flower-shaped structure is executed.

[0093] The beneficial effects of the technical solution provided by this application are:

[0094] Generate n three-dimensional fractal data bodies based on the fractal noise method; extract two-dimensional slices from each three-dimensional fractal data body along the Z direction and add reference depth to obtain multiple sedimentary stratigraphic sequences; add corresponding dip angles to each sedimentary stratigraphic sequence; superimpose multiple sedimentary stratigraphic sequences to obtain a stratigraphic model; add a main strike-slip fault, and randomly generate multiple dip-slip fault center points and masks of the left and right plates of the strike-slip fault; add simulated faults of the same type, use masks to shield simulated faults that exceed the strike-slip fault, and finally obtain a negative flower or positive flower structural model. The present invention can effectively generate a realistic strike-slip fault flower-shaped structural model, provide training data for using deep learning to finely characterize strike-slip faults and their flower-shaped structures, not only can it alleviate the problems of insufficient data and inaccurate labels, improve the recognition accuracy of the model, but also provide more reliable technical support for geological research and resource development. BRIEF DESCRIPTION OF THE DRAWINGS

[0095] The present application will be further described below with reference to the accompanying drawings and embodiments, in which:

[0096] Figure 1 is a step diagram in an embodiment of the present application;

[0097] Figure 2 is an exemplary schematic diagram of three-dimensional white noise in the spatial domain in an embodiment of the present invention;

[0098] Figure 3 is an exemplary schematic diagram of a filter in an embodiment of the present invention;

[0099] Figure 4 is an exemplary schematic diagram of white noise frequency domain filtering results in an embodiment of the present invention;

[0100] Figure 5 is an exemplary schematic diagram of white noise spatial domain filtering results in an embodiment of the present invention;

[0101] Figure 6 is an exemplary schematic diagram of three sedimentary stratigraphic sequences in an embodiment of the present invention;

[0102] Figure 7 is an exemplary schematic diagram of three added dip angle sedimentary stratigraphic sequences in the embodiments of the present application;

[0103] Figure 8 is an exemplary schematic diagram of a stratigraphic model in an embodiment of the present application;

[0104] Fig. 9 This is an exemplary schematic diagram of adding a main strike-slip fault to the stratigraphic model in the embodiment of the present application;

[0105] Fig.10 is an exemplary schematic diagram of the left and right disk masks of the main strike-slip fault in the embodiment of the present application;

[0106] Fig.11 This is an exemplary schematic diagram of the positive flower and negative flower structures finally generated in the embodiment of the present application;

[0107] Fig.12 It is a schematic diagram of the structure of an electronic device in an embodiment of the present application. DETAILED DESCRIPTION

[0108] In order to have a clearer understanding of the technical features, purposes and effects of the present application, the specific implementation methods of the present application are now described in detail with reference to the accompanying drawings.

[0109] An embodiment of the present application provides a random modeling method for deep strike-slip fault flower-like structures.

[0110] Please refer to Figure 1 , Figure 1 This is a step diagram of a random modeling method for deep strike-slip fault flower-shaped structure in an embodiment of the present application, comprising:

[0111] S1: Obtain a three-dimensional white noise data volume in the spatial domain; obtain a three-dimensional fractal data volume in the spatial domain by combining Fourier transform and inverse Fourier transform according to the three-dimensional white noise data volume in the spatial domain;

[0112] Step S1 includes:

[0113] S11: Based on the fractal noise method, n three-dimensional random number matrices are generated as three-dimensional white noise data volumes in n spatial domains;

[0114] Specifically, in this application, n is set to 3, and the size of the three-dimensional random number matrix is ​​set to ,like Figure 2 shown.

[0115] S12: Using Fourier transform, the three-dimensional white noise data volume in the spatial domain is converted into a three-dimensional white noise data volume in the frequency domain. The specific formula is as follows:

[0116]

[0117] in, is the length of the three-dimensional white noise data volume in the spatial domain, is the width of the three-dimensional white noise data volume in the spatial domain, is the height of the three-dimensional white noise data volume in the spatial domain; Represents a three-dimensional white noise data volume in the frequency domain; is the spatial frequency in the X direction in the spatial domain, The range is ; is the spatial frequency in the Y direction of the spatial domain, It is the spatial domain The spatial frequency in the direction, The range is , The range is ; represents a three-dimensional white noise data volume in the spatial domain, Represents a three-dimensional white noise data volume in the spatial domain Axis coordinates; Represents a three-dimensional white noise data volume in the spatial domain Axis coordinates, Represents a three-dimensional white noise data volume in the spatial domain Axis coordinates; The range is , The range is , The range is ;

[0118] S13: Construct a frequency domain filter, the formula is:

[0119]

[0120] in, Represents a three-dimensional white noise data volume in the spatial domain The frequency in the direction, express Directional weight; Represents a three-dimensional white noise data volume in the spatial domain The frequency in the direction, express Directional weight; Represents a three-dimensional white noise data volume in the spatial domain The frequency in the direction, express Directional weight; fractal dimension Related constants ;

[0121] Specifically, in this application is 1, is 1, is 2, D is 1.5.

[0122] The frequency domain filter is multiplied by the three-dimensional white noise data volume in the frequency domain for filtering. The expression is as follows:

[0123]

[0124] in, Represents the original three-dimensional white noise data volume in the frequency domain, represents a three-dimensional white noise data volume in the frequency domain after filtering;

[0125] S14: Apply inverse Fourier transform to the filtered three-dimensional white noise data volume in the frequency domain to obtain n three-dimensional fractal data volumes in the space domain. The calculation formula is:

[0126]

[0127] in, represents the three-dimensional white noise data volume in the frequency domain after filtering, is the spatial frequency in the X direction in the spatial domain, is the spatial frequency in the Y direction of the spatial domain, It is the spatial domain The spatial frequency in the direction, The range is , The range is ; The range is ; represents the filtered spatial domain fractal data volume, Representing 3D fractal data in spatial domain Axis coordinates, Representing 3D fractal data in spatial domain Axis coordinates, Representing 3D fractal data in spatial domain Axis coordinates, The range is , The range is , The range is ; is the length of the three-dimensional fractal data volume in the spatial domain, is the width of the three-dimensional fractal data volume in the spatial domain, is the height of the three-dimensional fractal data volume in the spatial domain.

[0128] In one embodiment provided in the present application, a three-dimensional white noise data volume in the spatial domain is converted into a three-dimensional white noise data volume in the frequency domain by using Fourier transform. The frequency domain filter, such as Figure 3 As shown. Frequency domain three-dimensional white noise data volume filtering, such as Figure 4 As shown. The three-dimensional fractal data volume in the spatial domain is Figure 5 As shown. Use the standard_normal method in numpy to directly generate a three-dimensional random number matrix.

[0129] S2: extracting a preset number of two-dimensional slices along the Z direction from each three-dimensional fractal data volume in the spatial domain, each slice is regarded as a stratigraphic plane, and a reference depth is added to each slice as the depth of the stratigraphic plane. The two-dimensional slices extracted from the same three-dimensional fractal data volume in the spatial domain are regarded as a sedimentary stratigraphic sequence. Applying the above operation to each three-dimensional fractal data volume in the spatial domain obtains n sedimentary stratigraphic sequences;

[0130] Step S2 includes:

[0131] A preset number of two-dimensional slices are extracted along the Z direction from n three-dimensional fractal data volumes in the spatial domain. Each slice is regarded as a layer plane. Each slice is added with a reference depth as the depth of the layer plane. The calculation formula is:

[0132]

[0133] in, Indicates the depth of the i-th layer extracted from the n-th three-dimensional fractal data volume, represents the i-th slice extracted from the n-th three-dimensional fractal data volume, Indicates the base depth of the i-th layer plane extracted from the n-th three-dimensional fractal data volume;

[0134] The two-dimensional slices extracted from the same three-dimensional fractal data volume in the spatial domain are taken as a sedimentary stratigraphic sequence, and n sedimentary stratigraphic sequences of n three-dimensional fractal data volumes in the spatial domain are obtained.

[0135] Specifically, Figure 6 As shown, in the present application, 8 slices are extracted from the first fractal, 2 slices are extracted from the second fractal, and 17 slices are extracted from the third fractal, totaling 27 slices. The reference depths of each stratigraphic surface are 1, 7, 14, 22, 27, 33, 39, 45, 53, 58, 64, 69, 75, 80, 92, 104, 113, 124, 129, 138, 147, 152, 159, 170, 181, 187, and 192, respectively, and finally 3 sedimentary stratigraphic sequences are obtained.

[0136] S3: adding a dip angle to each sedimentary stratigraphic sequence to obtain n sedimentary stratigraphic sequences containing dip angles; obtaining a stratigraphic model by superimposing n sedimentary stratigraphic sequences containing dip angles;

[0137] Step S3 includes:

[0138] S31: Add the dip angle to each sedimentary stratigraphic sequence respectively, and the calculation formula is:

[0139]

[0140]

[0141] in, , , The three-dimensional coordinates of the nth sedimentary stratigraphic sequence after adding the dip angle, , , is the original three-dimensional coordinate of the nth sedimentary stratigraphic sequence; is the displacement field, , , As parameters, and Used to control the formation and The slope in the direction, Used to control the offset distance of the center point;

[0142] In one embodiment provided in the present application, Figure 7 As shown, the first sedimentary stratigraphic sequence , , are 0.145, 0.071 and -21.57 respectively; the second sedimentary stratigraphic sequence , , are 0.045, 0.127 and -17.21 respectively; the third sedimentary stratigraphic sequence , , They are 0.0096, 0.06697 and -7.657 respectively.

[0143] S32: superimposing n sedimentary stratigraphic sequences with added dip angles in a preset order to obtain a stratigraphic model containing unconformity phenomena.

[0144] In one embodiment provided in the present application, a stratigraphic model containing unconformity phenomena, such as Figure 8 In this example, n sedimentary stratigraphic sequences are stacked from top to bottom with the added dip angle.

[0145] S4: adding a main strike-slip fault to the stratigraphic model, and generating a preset number of dip-slip fault center points on the left and right sides of the strike-slip fault and masks of the left and right sides of the main strike-slip fault;

[0146] Step S4 includes:

[0147] S41: Generate a strike-slip fault center point on the horizon model , Steering angle And the inclination ;

[0148] S42: Taking a strike-slip fault center point is the origin of coordinates, the strike is the x-axis, the inclination is the y-axis, and the z-axis is the normal direction and points to the upper plate;

[0149] Establish a local coordinate system and transform all points of the layer model from original coordinates to local coordinates:

[0150]

[0151] in, is the x-axis coordinate in the local coordinate system, is the y-axis coordinate in the local coordinate system, is the z-axis coordinate in the local coordinate system, is the rotation matrix:

[0152]

[0153] The plane is the plane where the fault plane is located in the local coordinate system. In the current step, the fault plane is a strike-slip fault;

[0154] In an embodiment provided in the present application, a strike-slip fault (main strike-slip fault) is added in step S4, so the fault plane refers to the strike-slip fault. In step S5, when executing step S42, the fault plane refers to the dip-slip fault. According to the classification, faults can be divided into strike-slip faults and dip-slip faults.

[0155] S43: Calculate the displacement field of the fault plane based on the ellipse formula. The calculation formula is as follows:

[0156]

[0157] in Describe the fault displacement distribution on the fault plane, represents the maximum displacement of the center point, is the normalized radial distance of the center point, and the calculation formula is:

[0158]

[0159] in, is the radius along the strike direction, is the radius along the inclination direction, It is used to control the rotation angle of the displacement attenuation ellipse;

[0160] Calculate the strike displacement component separately and the tilt displacement component , the calculation formula is as follows:

[0161]

[0162]

[0163] S44: Fault plane in local coordinate system Generate a preset number of disturbance points within the preset range of :

[0164]

[0165] in, , and are the x-coordinate, y-coordinate and z-coordinate of the random point in the local coordinate system, is a random number; random points are disturbance points;

[0166] Use disturbance points to perturb the fault plane :

[0167]

[0168] S45: Calculate the three-dimensional fault displacement field on the fault plane. The calculation formula is as follows:

[0169]

[0170]

[0171]

[0172]

[0173] in, Indicates left disk; Indicates right plate; Indicates upper plate; Indicates the lower plate; , It represents the displacement ratio between the left and right disks; ( ) represents the displacement ratio of the upper plate and the lower plate; It is used to control the rotation angle of the displacement attenuation ellipse. When , the fault is a dip-slip fault. When , the fault is a strike-slip fault. When , the fault contains displacement in both strike and dip; , and They represent the fault displacement along the strike, dip and perpendicular to the fault plane in the local coordinate system respectively;

[0174] In the current step, , generating a strike-slip fault, i.e. a pure strike-slip fault;

[0175] S46: applying the three-dimensional fault displacement field in the local coordinate system to obtain a simulated fault;

[0176] The simulated fault is a simulated fault of a normal fault type or a simulated fault of a reverse fault type;

[0177] Convert the simulated fault back to the original coordinate system, and the calculation formula is as follows:

[0178]

[0179] in, is the rotation matrix The inverse matrix of , and They are the x-axis, y-axis, and z-axis coordinates of the original coordinate system after applying the three-dimensional fault displacement field in the local coordinate system. The model is recorded as , It is a strike-slip fault model;

[0180] In one embodiment provided in this application, Models such as Fig. 9 The simulated fault is a strike-slip fault or a dip-slip fault.

[0181] S47: Select n points within the preset range of the strike-slip fault center point in the local coordinate system As the center point of the dip-slip fault, the left disk of the main strike-slip fault in the local coordinate system is generated. and right disk The mask is calculated as follows:

[0182]

[0183]

[0184] Will , and Convert to the original coordinate system and get , and , the calculation formula is as follows:

[0185]

[0186]

[0187] .

[0188] In one embodiment provided in the present application, the left disk of the strike-slip fault in the local coordinate system and right disk The mask, such as Fig.10shown.

[0189] S5: Add the same type of simulated faults at the center point of the dip-slip fault, use a mask to shield the simulated faults that exceed the main strike-slip fault, and obtain a negative flower structure model or a positive flower structure model.

[0190] Step S5 includes:

[0191] S51: The center point of the dip-slip fault in is the center point of the simulated fault, and the dip angle is generated. and the angle of travel , Strike angle with the main strike-slip fault The calculation formula is as follows:

[0192]

[0193]

[0194] in, It indicates that the center point of the simulated fault is located on the left disk of the main strike-slip fault; It indicates that the center point of the simulated fault is located on the right disk of the main strike-slip fault;

[0195] S52: Generate the inclination angle by combining the calculation formulas from step S42 to step S46 and the angle of travel , generate a simulated fault. In this step, the simulated fault generated is a dip-slip fault;

[0196] Determine the type of simulated fault to be generated; the types of simulated faults include: normal fault simulated fault and reverse fault simulated fault;

[0197] Add the same type of simulated fault at the center of the dip-slip fault to obtain the Iteration Model, , The model adds Model of dip-slip faults;

[0198] S53: Use a mask to shield the simulated faults that exceed the main strike-slip fault. The calculation formula is as follows:

[0199]

[0200]

[0201] in, It indicates that the center point of the simulated fault is located on the left disk of the main strike-slip fault; It indicates that the center point of the simulated fault is located on the right disk of the main strike-slip fault;

[0202] S54: Repeat steps S51 to S53 until Each point in the model corresponds to a simulated fault, resulting in a positive structural model or a negative structural model.

[0203] In one embodiment provided in the present application, a positive flower structure model or a negative flower structure model, such as Fig.11 shown.

[0204] The present application also discloses an electronic device. Fig.12 , Fig.12 The electronic device 500 may include: at least one processor 501 , at least one network interface 504 , a user interface 503 , a memory 505 , and at least one communication bus 502 .

[0205] The communication bus 502 is used to realize the connection and communication between these components.

[0206] The user interface 503 may include a display screen, and the optional user interface 503 may also include a standard wired interface or a wireless interface.

[0207] The network interface 504 may optionally include a standard wired interface or a wireless interface (such as a WI-FI interface).

[0208] The present application also discloses a computer-readable storage medium storing a plurality of instructions suitable for loading by a processor to execute the above-mentioned random modeling method of deep strike-slip fault flower-shaped structure.

[0209] The above are only exemplary embodiments of the present disclosure and cannot be used to limit the scope of the present disclosure. That is, any equivalent changes and modifications made according to the teachings of the present disclosure are still within the scope of the present disclosure. After considering the disclosure of the specification and practice, those skilled in the art will easily think of other embodiments of the present disclosure.

[0210] This application is intended to cover any variation, use or adaptation of the present disclosure, which follows the general principles of the present disclosure and includes common knowledge or customary technical means in the art not described in the present disclosure. The description and examples are to be regarded as exemplary only, and the scope and spirit of the present disclosure are defined by the claims.

Claims

1. A random modeling method for deep strike-slip fault flower-shaped structures, characterized in that: The method comprises the following steps: S1: Obtain three-dimensional white noise data volume in the spatial domain; According to the three-dimensional white noise data volume in the spatial domain, a three-dimensional fractal data volume in the spatial domain is obtained by combining Fourier transform and inverse Fourier transform; Step S1 includes: S11: Based on the fractal noise method, n three-dimensional random number matrices are generated as three-dimensional white noise data volumes in n spatial domains; S12: Using Fourier transform, the three-dimensional white noise data volume in the spatial domain is converted into a three-dimensional white noise data volume in the frequency domain. The specific formula is as follows: in, is the length of the three-dimensional white noise data volume in the spatial domain, is the width of the three-dimensional white noise data volume in the spatial domain, is the height of the three-dimensional white noise data volume in the spatial domain; Represents a three-dimensional white noise data volume in the frequency domain; is the spatial frequency in the X direction in the spatial domain, The range is ; is the spatial frequency in the Y direction of the spatial domain, It is the spatial domain The spatial frequency in the direction, The range is , The range is ; represents a three-dimensional white noise data volume in the spatial domain, Represents a three-dimensional white noise data volume in the spatial domain Axis coordinates; Represents a three-dimensional white noise data volume in the spatial domain Axis coordinates, Represents a three-dimensional white noise data volume in the spatial domain Axis coordinates; The range is , The range is , The range is ; S13: Construct a frequency domain filter, the formula is: in, Represents a three-dimensional white noise data volume in the spatial domain The frequency in the direction, express Directional weight; Represents a three-dimensional white noise data volume in the spatial domain The frequency in the direction, express Directional weight; Represents a three-dimensional white noise data volume in the spatial domain The frequency in the direction, express Directional weight; fractal dimension Related constants ; The frequency domain filter is multiplied by the three-dimensional white noise data volume in the frequency domain for filtering. The expression is as follows: in, Represents the original three-dimensional white noise data volume in the frequency domain, represents a three-dimensional white noise data volume in the frequency domain after filtering; S14: Apply inverse Fourier transform to the filtered three-dimensional white noise data volume in the frequency domain to obtain n three-dimensional fractal data volumes in the space domain. The calculation formula is: in, represents the three-dimensional white noise data volume in the frequency domain after filtering, is the spatial frequency in the X direction in the spatial domain, is the spatial frequency in the Y direction of the spatial domain, It is the spatial domain The spatial frequency in the direction, The range is , The range is ; The range is ; represents the filtered spatial domain fractal data volume, Representing 3D fractal data in spatial domain Axis coordinates, Representing 3D fractal data in spatial domain Axis coordinates, Representing 3D fractal data in spatial domain Axis coordinates, The range is , The range is , The range is ; is the length of the three-dimensional fractal data volume in the spatial domain, is the width of the three-dimensional fractal data volume in the spatial domain, is the height of the three-dimensional fractal data volume in the spatial domain; S2: extracting a preset number of two-dimensional slices along the Z direction from each three-dimensional fractal data volume in the spatial domain, each slice is regarded as a stratigraphic plane, and a reference depth is added to each slice as the depth of the stratigraphic plane. The two-dimensional slices extracted from the same three-dimensional fractal data volume in the spatial domain are regarded as a sedimentary stratigraphic sequence, and n sedimentary stratigraphic sequences are obtained; S3: adding a dip angle to each sedimentary stratigraphic sequence to obtain n sedimentary stratigraphic sequences containing dip angles; obtaining a stratigraphic model by superimposing n sedimentary stratigraphic sequences containing dip angles; S4: adding a main strike-slip fault to the stratigraphic model, and generating a preset number of dip-slip fault center points on the left and right sides of the strike-slip fault and masks of the left and right sides of the main strike-slip fault; S5: Add the same type of simulated faults at the center point of the dip-slip fault, use a mask to shield the simulated faults that exceed the main strike-slip fault, and obtain a negative flower structure model or a positive flower structure model.

2. A method for stochastic modeling of deep strike-slip fault flower-shaped structures according to claim 1, characterized in that: Step S2 includes: A preset number of two-dimensional slices are extracted along the Z direction from n three-dimensional fractal data volumes in the spatial domain. Each slice is regarded as a layer plane. A reference depth is added to each slice as the depth of the layer plane. The calculation formula is: in, Indicates the depth of the i-th layer extracted from the n-th three-dimensional fractal data volume, represents the i-th slice extracted from the n-th three-dimensional fractal data volume, Indicates the base depth of the i-th layer plane extracted from the n-th three-dimensional fractal data volume; The two-dimensional slices extracted from the same three-dimensional fractal data volume in the spatial domain are taken as a sedimentary stratigraphic sequence, and n sedimentary stratigraphic sequences of n three-dimensional fractal data volumes in the spatial domain are obtained.

3. A stochastic modeling method for deep strike-slip fault flower-shaped structures according to claim 1, characterized in that: Step S3 includes: S31: Add the dip angle to each sedimentary stratigraphic sequence respectively, and the calculation formula is: in, , , The three-dimensional coordinates of the nth sedimentary stratigraphic sequence after adding the dip angle, , , is the original three-dimensional coordinate of the nth sedimentary stratigraphic sequence; is the displacement field, , , As parameters, and Used to control the formation and The slope in the direction, Used to control the offset distance of the center point; S32: superimposing n sedimentary stratigraphic sequences with added dip angles in a preset order to obtain a stratigraphic model containing unconformity phenomena.

4. A method for stochastic modeling of deep strike-slip fault flower-shaped structures according to claim 1, characterized in that: Step S4 includes: S41: Generate a strike-slip fault center point on the horizon model , Steering angle And the inclination ; S42: Taking a strike-slip fault center point is the origin of coordinates, the strike is the x-axis, the inclination is the y-axis, and the z-axis is the normal direction and points to the upper plate; Establish a local coordinate system and transform all points of the layer model from original coordinates to local coordinates: in, is the x-axis coordinate in the local coordinate system, is the y-axis coordinate in the local coordinate system, is the z-axis coordinate in the local coordinate system, is the rotation matrix: The plane is the plane where the fault plane is located in the local coordinate system. In the current step, the fault plane is a strike-slip fault; S43: Calculate the displacement field of the fault plane based on the ellipse formula. The calculation formula is as follows: in Describe the fault displacement distribution on the fault plane, represents the maximum displacement of the center point, is the normalized radial distance of the center point, and the calculation formula is: in, is the radius along the strike direction, is the radius along the inclination direction, It is used to control the rotation angle of the displacement attenuation ellipse; Calculate the strike displacement component separately and the tilt displacement component , the calculation formula is as follows: S44: Fault plane in local coordinate system Generate a preset number of disturbance points within the preset range of : in, , and are the x-coordinate, y-coordinate and z-coordinate of the random point in the local coordinate system, is a random number; random points are disturbance points; Use disturbance points to perturb the fault plane : S45: Calculate the three-dimensional fault displacement field on the fault plane. The calculation formula is as follows: in, Indicates left disk; Indicates right plate; Indicates upper plate; Indicates the lower plate; , It represents the displacement ratio between the left and right disks; ( ) represents the displacement ratio of the upper plate and the lower plate; It is used to control the rotation angle of the displacement attenuation ellipse. When , the fault is a dip-slip fault. When , the fault is a strike-slip fault. When , the fault contains displacement in both strike and dip; , and They represent the fault displacement along the strike, dip and perpendicular to the fault plane in the local coordinate system respectively; In the current step, , generating a strike-slip fault, i.e. a pure strike-slip fault; S46: applying the three-dimensional fault displacement field in the local coordinate system to obtain a simulated fault; The simulated fault is a simulated fault of a normal fault type or a simulated fault of a reverse fault type; Convert the simulated fault back to the original coordinate system, and the calculation formula is as follows: in, is the rotation matrix The inverse matrix of , and They are the x-axis, y-axis, and z-axis coordinates of the original coordinate system after applying the three-dimensional fault displacement field in the local coordinate system. The model is recorded as , It is a strike-slip fault model; S47: Select n points within the preset range of the strike-slip fault center point in the local coordinate system As the center point of the dip-slip fault, the left disk of the main strike-slip fault in the local coordinate system is generated. and right disk The mask is calculated as follows: Will , and Convert to the original coordinate system and get , and , the calculation formula is as follows: 。 5. A method for stochastic modeling of deep strike-slip fault flower-shaped structures according to claim 4, characterized in that: Step S5 includes: S51: The center point of the dip-slip fault in is the center point of the simulated fault, and the dip angle is generated. and the angle of travel , Strike angle with the main strike-slip fault The calculation formula is as follows: in, It indicates that the center point of the simulated fault is located on the left disk of the main strike-slip fault; It indicates that the center point of the simulated fault is located on the right disk of the main strike-slip fault; S52: Generate the inclination angle by combining the calculation formulas from step S42 to step S46 and the angle of travel , generate a simulated fault. In this step, the simulated fault generated is a dip-slip fault; Determine the type of simulated fault to be generated; the types of simulated faults include: normal fault simulated fault and reverse fault simulated fault; Add the same type of simulated fault at the center of the dip-slip fault to obtain the Iteration Model, , The model adds Model of dip-slip faults; S53: Use a mask to shield the simulated faults that exceed the main strike-slip fault. The calculation formula is as follows: in, It indicates that the center point of the simulated fault is located on the left disk of the main strike-slip fault; It indicates that the center point of the simulated fault is located on the right disk of the main strike-slip fault; S54: Repeat steps S51 to S53 until Each point in the model corresponds to a simulated fault, resulting in a positive structural model or a negative structural model.

6. An electronic device, characterized in that: The electronic device comprises a processor (501), a memory (505), a user interface (503) and a network interface (504), wherein the memory (505) is used to store instructions, the user interface (503) and the network interface (504) are used to communicate with other devices, and the processor (501) is used to execute the instructions stored in the memory (505) so that the electronic device executes the method according to any one of claims 1 to 5.

7. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores instructions, and when the instructions are executed by a computer, the method according to any one of claims 1 to 5 is executed.

Images