CN111596352B

CN111596352B - Method, system and device for analyzing spatial development law of bead bodies and storage medium

Info

Publication number: CN111596352B
Application number: CN202010349125.0A
Authority: CN
Inventors: 朱光有; 田飞; 杨敏; 王萌; 张志遥
Original assignee: Petrochina Co Ltd
Current assignee: Petrochina Co Ltd
Priority date: 2020-04-28
Filing date: 2020-04-28
Publication date: 2023-04-07
Anticipated expiration: 2040-04-28
Also published as: CN111596352A

Abstract

The invention provides a method, a system, a device and a storage medium for analyzing a spatial development rule of a bead body. The method comprises the following steps: acquiring a high-precision three-dimensional seismic amplitude data volume of a work area; determining the development position of each bead body; determining a bead body space development rule comprises determining a bead body longitudinal development rule and/or a bead body transverse development rule; determining a longitudinal development rule comprises respectively determining optimal development windows of the bead bodies corresponding to all the layers of the work area, respectively determining the sum of the first amplitude change rates of the optimal development windows of all the layers, and determining the dominant development layer of the bead bodies according to the numerical value of the sum of the first amplitude change rates of all the layers; determining the transverse development rule comprises respectively determining the optimal development width of the bead body corresponding to each sliding fracture in the work area, respectively determining the sum of the second amplitude change rates in the optimal development width of each sliding fracture, and determining the main control sliding fracture of the bead body according to the numerical value of the sum of the second amplitude change rates of each sliding fracture.

Description

Method, system and device for analyzing spatial development law of bead bodies and storage medium

Technical Field

The invention relates to the technical field of oil exploration, in particular to a method, a system, a device and a storage medium for analyzing a spatial development rule of a bead body by using a data driving method.

Background

The search of oil and gas in the ultra-deep carbonate rock stratum of the Chinese marine basin is one of the trends of future oil and gas exploration in China. In recent years, significant progress has been made successively in deep and ultra-deep oil and gas exploration in the Tarim basin, sichuan basin and Ordos basin. However, the reservoir of the carbonate oil-gas reservoir in the deep layer of China is ancient (mainly ancient world), and the heterogeneity is strong. Due to reservoir depositional environment and the complex geometry inside, quantitative characterization of "beadlets" of carbonate reservoirs presents a significant challenge to geophysicists.

Data mining and exploratory statistics have many practical applications that can take maximum advantage of various types of data, look for previously missed trends and solve problems by analogy. Processing large amounts of data can reduce the degree of uncertainty, which is a common feature of geological interpretation.

Disclosure of Invention

Aiming at the defects of the prior art, the invention aims to provide a method for quantitatively analyzing the spatial development law of a bead body. According to the method, geological interpretation of the bead body is performed by means of data mining and exploratory statistics, so that the method is more beneficial to defining the development advantage zone of the bead body and makes up for the limitation of qualitative analysis.

In order to achieve the above object, the present invention provides a method for analyzing spatial development laws of bead bodies, wherein the method comprises:

acquiring a high-precision three-dimensional seismic amplitude data volume of a work area; the high-precision three-dimensional seismic amplitude data body is a three-dimensional seismic amplitude data body with a main frequency not lower than a preset main frequency, and the preset main frequency is not lower than 50Hz;

determining the development position of each bead body based on the high-precision three-dimensional seismic amplitude data body of the work area;

determining a spatial development rule of the bead bodies based on the high-precision three-dimensional seismic amplitude data body of the work area and the determined development positions of the bead bodies; wherein the determining of the spatial development law of the bead body comprises determining the longitudinal development law of the bead body and/or determining the transverse development law of the bead body;

the determining of the longitudinal development law of the bead body comprises: respectively determining optimal development windows of the bead bodies corresponding to all the layers of the work area, respectively determining the sum of first amplitude change rates of the optimal development windows of all the layers, and determining the dominant development layer of the bead bodies according to the numerical value of the sum of the first amplitude change rates of all the layers;

the determining the transverse development rule of the bead body comprises the following steps: respectively determining the optimal development width of the bead body corresponding to each sliding fracture of the work area, respectively determining the sum of second amplitude change rates of each sliding fracture in the optimal development width, and determining the main control sliding fracture of the bead body according to the numerical value of the sum of the second amplitude change rates of each sliding fracture.

In the above method for analyzing spatial development law of bead bodies, preferably, the preset dominant frequency is 50Hz to 60Hz.

In the above method for analyzing spatial development laws of a beading body, preferably, the acquiring a high-precision three-dimensional seismic amplitude data volume of a work area includes:

acquiring a three-dimensional seismic amplitude data volume of a work area;

judging whether the dominant frequency of the three-dimensional seismic amplitude data body of the work area is lower than the preset dominant frequency or not;

if the dominant frequency of the three-dimensional seismic amplitude data body of the work area is not lower than the preset dominant frequency, the three-dimensional seismic amplitude data body of the work area is the high-precision three-dimensional seismic amplitude data body of the work area;

and if the main frequency of the three-dimensional seismic amplitude data body of the work area is lower than the preset main frequency, performing frequency broadening processing on the three-dimensional seismic amplitude data body of the work area to generate a high-precision three-dimensional seismic amplitude data body of the work area, wherein the main frequency is not lower than the preset main frequency.

In one embodiment, the Frequency broadening processing of the three-dimensional seismic data is performed on the three-dimensional seismic amplitude data volume of the work area by using an ultra High resolution Imaging (HFI) method.

In the above method for analyzing spatial development regularity of a bead body, preferably, the determining the development position of the bead body includes: and calculating the amplitude change rate of different parts of each layer in the work area, wherein the part with the amplitude change rate not lower than a first preset threshold value is the development position of the bead body. More preferably, the calculating the amplitude change rate of different portions of each layer in the work area includes: and calculating the amplitude change rate of the upper part and the lower part, the east part and the west part, the south part and the north part of the earthquake reflection point by taking the earthquake reflection point as a center. In a specific embodiment, the determining the developmental location of a bead comprises: calculating the amplitude change rate of the upper part and the lower part, the east part and the west part, the south part and the north part of the earthquake reflection point by taking the earthquake reflection point as a center; wherein, the upper and lower computing windows are 3-5, the east and west computing windows are 10-20, and the south and north computing windows are 10-20; the part with the amplitude change rate not lower than a first preset threshold value is the bead body development position, wherein the first preset threshold value is 9200-9800.

In the above method for analyzing the spatial development law of a bead body, preferably, the determining optimal development windows of the bead body corresponding to each layer position of the work area respectively includes: aiming at each layer in the work area, respectively calculating the distance between the top of each bead body and the layer, respectively taking each layer as a reference, calculating the number Nwi of the bead bodies within the range from one layer to different depths W ms,

and when the maximum value is obtained, the value of W is the optimal development window of the bead body corresponding to the layer. In a specific embodiment, the calculating the number of beads Nwi in the range from a horizon to a different depth, W ms, comprises: setting the initial value of the search window length as W = W ₀ From the horizon to W = W ₀ Number of ms-deep beads Nwi = Nw ₀ And calculates->

Generating a new W value in a certain step size (e.g., 1-5 ms), calculating the number of beads in the depth range from the horizon to W ms Nwi, and determining ^ or ^ R>

Until the depth reaches the top surface of the next horizon.

In the above method for analyzing spatial development laws of beads, preferably, the determining a dominant development horizon of beads according to a value of a sum of the first amplitude change rates of each horizon includes: and determining the first N layers with the maximum value of the sum of the first amplitude change rates as the dominant development layers of the bead bodies. In a specific embodiment, the first 3 horizons with the largest sum of the first amplitude change rates are determined as dominant development horizons of the bead body.

In the above method for analyzing spatial development law of beads, preferably, the determining the longitudinal development law of beads further includes:

and further determining at least one of the top development depth, the bottom development depth, the east-west development width, the north-south development width, the maximum amplitude change rate value, the oil production of the drilling encountering beading body, the gas production of the drilling encountering beading body, the water production of the drilling encountering beading body and the oil-gas-water geochemical parameters of the drilling encountering beading body on the basis of the dominant development horizon of the beading body.

In the above method for analyzing the spatial development law of a bead body, preferably, the determining the optimal development width of the bead body corresponding to each slipping fracture of the work area includes: aiming at each sliding fracture, respectively calculating the distance between each bead body and the sliding fracture, taking the sliding fracture as a reference, calculating the number Nli of the bead bodies extending from one sliding fracture to two sides within the range of L ms,

and the value of L is the optimal development width of the bead body corresponding to the sliding fracture when the value is the maximum value. In one embodiment, the calculating the number Nli of beads extending bilaterally from a walking slip fracture in the range of L ms comprises: setting an initial value of a search window length to L = L ₀ Calculating the extension L = L from one sliding break to both sides ₀ Number of beads in ms range Nli = Nl ₀ And calculates->

Generating a new L value in a certain step length (for example 1-5 ms), calculating the number Nli of the beads extending from a slide break to both sides in the range of L ms, and determining ^ based on the number Nli>

Until it extends sideways from the trip break to another trip break.

In the above method for analyzing a spatial development law of a bead body, preferably, the determining the main control sliding fracture of the bead body according to the magnitude of the sum of the second amplitude change rates of the sliding fractures comprises: and determining the first N sliding fractures with the maximum numerical value of the sum of the second amplitude change rates as the main control sliding fractures of the bead body. In a specific embodiment, the first 5-15 slip fractures with the largest sum of the second amplitude change rates are determined as the main control slip fractures of the bead body.

The invention also provides a system for analyzing the spatial development law of the bead bodies, wherein the system comprises:

the high-precision three-dimensional seismic amplitude data volume acquisition module comprises: the method comprises the steps of obtaining a high-precision three-dimensional seismic amplitude data volume of a work area; the high-precision three-dimensional seismic amplitude data body is a three-dimensional seismic amplitude data body with a main frequency not lower than a preset main frequency, and the preset main frequency is not lower than 50Hz;

position location module is grown to the string of beads body: the system is used for determining the development position of each bead body based on the high-precision three-dimensional seismic amplitude data body of the work area acquired by the acquisition module;

the bead body space development rule analysis module: the system is used for determining the spatial development rule of the bead bodies on the basis of the high-precision three-dimensional seismic amplitude data volume of the work area acquired by the acquisition module and the development positions of the bead bodies determined by the positioning module; wherein the determining of the spatial development law of the bead body comprises determining the longitudinal development law of the bead body and/or determining the transverse development law of the bead body;

the determining of the lateral development law of the bead body comprises: respectively determining the optimal development width of the bead body corresponding to each sliding fracture of the work area, respectively determining the sum of second amplitude change rates of each sliding fracture within the optimal development width, and determining the main control sliding fracture of the bead body according to the numerical value of the sum of the second amplitude change rates of each sliding fracture.

In the above system for analyzing spatial development regularity of a bead body, preferably, the preset dominant frequency is 50Hz to 60Hz.

In the above system for analyzing a spatial development law of a bead body, preferably, the high-precision three-dimensional seismic amplitude data volume acquisition module includes:

a first acquisition unit: the method comprises the steps of obtaining a three-dimensional seismic amplitude data volume of a work area;

a first discrimination unit: judging whether the dominant frequency of the three-dimensional seismic amplitude data of the work area acquired by the first acquisition unit is lower than the preset dominant frequency; if the first judging unit judges that the main frequency of the three-dimensional seismic amplitude data body of the work area is not lower than the preset main frequency, a first processing unit is adopted for processing; if the dominant frequency of the three-dimensional seismic amplitude data body of the work area is lower than the preset dominant frequency, a second processing unit is adopted for processing;

a first processing unit: for performing operations that use the three-dimensional seismic amplitude data volume of the work area as a high-precision three-dimensional seismic amplitude data volume of the work area.

A second processing unit: and the frequency extension processing module is used for carrying out frequency extension processing on the three-dimensional seismic amplitude data body of the work area to generate a high-precision three-dimensional seismic amplitude data body of the work area with the main frequency not lower than the preset main frequency.

In a specific embodiment, the second processing unit performs Frequency broadening processing on the three-dimensional seismic amplitude data volume of the work area by using an ultra High resolution processing method (HFI).

In the above system for analyzing a spatial development law of a bead body, preferably, the module for locating a development position of a bead body includes:

a third processing unit: the system is used for calculating the amplitude change rate of different parts of each layer in the work area based on the high-precision three-dimensional seismic amplitude data volume of the work area acquired by the acquisition module;

a second determination unit: the amplitude change rate of different parts of each layer in the work area calculated by the third processing unit is judged whether to be lower than the first preset threshold value or not; the part with the amplitude change rate not lower than the first preset threshold value is the development position of the bead body.

More preferably, the calculating the amplitude change rate of different portions of each layer in the work area includes: and calculating the amplitude change rate of the upper part and the lower part, the east part and the west part, the south part and the north part of the earthquake reflection point by taking the earthquake reflection point as a center.

In a specific embodiment, the bead body development position locating module comprises:

a third processing unit: the system is used for calculating the amplitude change rate of different parts of each layer in the work area based on the high-precision three-dimensional seismic amplitude data volume of the work area acquired by the acquisition module; calculating the amplitude change rate of different parts of each layer in the work area by taking the seismic reflection point as a center, and calculating the amplitude change rate of three different dimensions of the upper part, the lower part, the east part, the west part, the south part and the north part of the seismic reflection point; wherein, the upper and lower computing windows are 3-5, the east and west computing windows are 10-20, and the south and north computing windows are 10-20;

a second determination unit: the amplitude change rate of different parts of each layer in the work area calculated by the third processing unit is judged whether to be lower than the first preset threshold value or not; the part of which the amplitude change rate is not lower than a first preset threshold value is a bead body development position; wherein the first predetermined threshold is 9200-9800.

In the above system for analyzing spatial development law of bead body, preferably, the module for analyzing spatial development law of bead body comprises:

a fourth processing unit: the system is used for respectively determining optimal development windows of the bead bodies corresponding to all the layers of the work area; aiming at each layer in the work area, respectively calculating the distance between the top of each bead body and the layer, respectively taking each layer as a reference, calculating the number Nwi of the bead bodies within the range from one layer to different depths W ms,

when the maximum value is W, the value of W is the optimal development window of the bead body corresponding to the layer;

a fifth processing unit: the processing unit is used for calculating the sum of the first amplitude change rates of the optimal development windows of the layers based on the bead body optimal development windows corresponding to the layers of the work area determined by the fourth processing unit;

a sixth processing unit: the method is used for determining the dominant development horizon of the bead body based on the numerical value of the sum of the first amplitude change rates of the horizons calculated by the fifth processing unit.

In a specific embodiment, the calculating the number of beads Nwi in a range from a horizon to a different depth of W ms comprises: setting an initial value of search window length as W = W ₀ From the horizon to W = W ₀ Number of ms-deep beads Nwi = Nw ₀ And calculate

Until the depth reaches the top surface of the next horizon.

In the above system for analyzing a spatial development law of a bead body, preferably, the sixth processing unit is configured to determine, based on a sum of the first amplitude change rates of the respective layers calculated by the fifth processing unit, the first N layers with the largest value of the sum of the first amplitude change rates as dominant development layers of the bead body.

In a specific embodiment, the sixth processing unit is configured to determine, based on the sum of the first amplitude change rates of the respective layers calculated by the fifth processing unit, the first 3 layers with the largest value of the sum of the first amplitude change rates as dominant development layers of the bead body.

In the above system for analyzing spatial development laws of beading bodies, preferably, the module for analyzing spatial development laws of beading bodies further comprises:

a seventh processing unit: the method is used for further determining at least one of the top development depth, the bottom development depth, the east-west development width, the north-south development width, the maximum amplitude change rate value, the oil production of the drilling encountering beading body, the gas production of the drilling encountering beading body, the water production of the drilling encountering beading body and the oil-gas-water geochemical parameters of the drilling encountering beading body based on the dominant development horizon of the beading body determined by the sixth processing unit.

an eighth processing unit: the bead body optimal development width corresponding to each sliding fracture of the work area is respectively determined; aiming at each sliding fracture, respectively calculating the distance between each bead body and the sliding fracture, taking the sliding fracture as a reference, calculating the number Nli of the bead bodies extending from one sliding fracture to two sides within the range of L ms,

the value of L is the optimal development width of the bead body corresponding to the sliding fracture when the value is the maximum value;

a ninth processing unit: the second amplitude change rate summation is used for calculating the second amplitude change rate summation in the optimal development width of each slipping fracture based on the optimal development width of the bead body corresponding to each slipping fracture of the work area determined by the eighth processing unit;

a tenth processing unit: and the dominant development horizon of the bead body is determined based on the numerical value of the sum of the second amplitude change rates of the sliding fractures calculated by the ninth processing unit.

In one embodiment, the calculating the number Nli of beads extending bilaterally from a walking slip fracture in the range of L ms comprises: setting an initial value of a search window length to L = L ₀ Calculating the extension L = L from one sliding break to both sides ₀ Number of beads in ms range Nli = Nl ₀ And calculate

Generating new L value according to a certain step length (for example 1-5 ms), calculating the number Nli of the beads extending from one sliding fracture to two sides within the range of L ms, and determining/>

Until it extends sideways from the trip break to another trip break.

In the above system for analyzing a spatial development law of a bead body, preferably, the tenth processing unit is configured to determine, based on a sum of the second amplitude change rates of the respective sliding fractures calculated by the ninth processing unit, the first N sliding fractures with the largest value of the sum of the second amplitude change rates as the main control sliding fractures of the bead body.

In a specific embodiment, the tenth processing unit is configured to determine, based on the sum of the second amplitude change rates of the respective sliding fractures calculated by the ninth processing unit, the first 5 to 15 sliding fractures with the largest value of the sum of the second amplitude change rates as the main control sliding fractures of the bead body.

The invention also provides a device for analyzing the spatial development law of the bead bodies, which comprises a processor and a memory; wherein the content of the first and second substances,

a memory for storing a computer program;

and the processor is used for realizing the steps of the method for analyzing the spatial development law of the bead body when executing the program stored in the memory.

The present invention also provides a computer-readable storage medium storing one or more programs, which are executable by one or more processors to implement the steps of the above-described method for analyzing spatial development regularity of a bead body.

According to the technical scheme provided by the invention, the geological interpretation of the bead body is carried out by utilizing a data mining and exploratory statistics mode, and the plane distribution rule and/or the longitudinal distribution rule of the bead body can be quantitatively determined by analyzing parameters in the technical scheme provided by the invention, so that the main development horizon and/or the main control sliding fracture are determined. The method can make up the limitation of qualitative analysis and dig out the hidden relation of different 'bead bodies'.

Drawings

In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts:

fig. 1 is a schematic flow chart of a method for analyzing spatial development laws of beads according to an embodiment of the present invention;

fig. 2 is a schematic flow chart of a method for analyzing spatial development laws of beads according to an embodiment of the present invention;

fig. 3 is a schematic flow chart of a method for analyzing spatial development laws of beads according to an embodiment of the present invention;

fig. 4 is an optimization schematic diagram of a step of acquiring high-precision three-dimensional seismic amplitude data of a work area in the bead body spatial development law analysis method according to an embodiment of the present invention;

fig. 5 is a schematic structural diagram of a system for analyzing spatial development laws of beads according to an embodiment of the present invention;

fig. 6 is a schematic structural diagram of a high-precision three-dimensional seismic amplitude data volume acquisition module in the system for analyzing a spatial development law of a bead body according to an embodiment of the present invention;

fig. 7 is a schematic structural diagram of a bead development position positioning module in the system for analyzing spatial development laws of beads according to an embodiment of the present invention;

fig. 8 is a schematic structural diagram of a bead spatial development law analysis module in the bead spatial development law analysis system according to an embodiment of the present invention;

fig. 9 is a schematic structural diagram of a bead spatial development law analysis module in the bead spatial development law analysis system according to an embodiment of the present invention;

fig. 10 is a schematic structural diagram of a bead spatial development law analysis module in the bead spatial development law analysis system according to an embodiment of the present invention;

fig. 11 is a schematic structural diagram of a device for analyzing spatial development laws of beads according to an embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be described in detail and completely with reference to the drawings in the embodiments of the present invention. It is to be understood that the embodiments described are only a few embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

The beads described in the present invention refer to beads in the conventional sense in the art. In the area of study, collapsed caverns and their debris appear as strong amplitude "beaded" reflections on seismic sections, mainly due to the large wave impedance difference between the cavern debris and the aodovician surrounding rock. These beaded reflections are interpreted as collapsed paleo-caverns and the well bore proved to be a good production zone. On an earthquake section, bead reflection is associated with micro fracture, a sunken structure is shown on a plane, and an area which is not influenced by cave collapse belongs to a cave undeveloped area. In order to image and identify the beaded reflection related to the collapse cause of the paleo-karst and the associated micro-fracture and crack, the beaded reflection characteristic and the linear structure characteristic of the micro-fracture and crack are extracted by the seismic reflection inclination angle, curvature, coherence and amplitude gradient attributes, so that the aim of detecting and identifying the fracture-cave system is fulfilled.

In the present invention, each horizon of the work area is not limited to each of all horizons in the work area, and may be each of the main horizons, each of the target study horizons, and the like.

The slip fractures of the work area are not limited to all slip fractures in the work area, and may be slip fractures in one or more layers, large slip fractures or medium-large slip fractures, and the like.

The principles and spirit of the present invention are explained in detail below with reference to several representative embodiments of the invention.

Referring to fig. 1, an embodiment of the present invention provides a method for analyzing spatial development laws of beads, including:

step S1: acquiring a high-precision three-dimensional seismic amplitude data volume of a work area; the high-precision three-dimensional seismic amplitude data body is a three-dimensional seismic amplitude data body with a main frequency not lower than a preset main frequency, and the preset main frequency is not lower than 50Hz;

step S2: determining the development position of each bead body based on the high-precision three-dimensional seismic amplitude data volume of the work area;

and step S3: determining a longitudinal development rule of the bead bodies based on the high-precision three-dimensional seismic amplitude data body of the work area and the determined development positions of the bead bodies;

wherein, confirm the longitudinal development law of string of beads body includes: determining optimal development windows of the bead bodies corresponding to all the layers of the work area respectively, determining the sum of first amplitude change rates of the optimal development windows of all the layers respectively, and determining the dominant development layer of the bead bodies according to the numerical value of the sum of the first amplitude change rates of all the layers.

Referring to fig. 2, an embodiment of the present invention provides a method for analyzing spatial development laws of beads, including:

step S2: determining the development position of each bead body based on the high-precision three-dimensional seismic amplitude data body of the work area;

and step S4: determining a transverse development rule of the bead bodies based on the high-precision three-dimensional seismic amplitude data body of the work area and the determined development positions of the bead bodies;

wherein, confirm the horizontal development law of string of beads body includes: respectively determining the optimal development width of the bead body corresponding to each sliding fracture of the work area, respectively determining the sum of second amplitude change rates of each sliding fracture within the optimal development width, and determining the main control sliding fracture of the bead body according to the numerical value of the sum of the second amplitude change rates of each sliding fracture.

Referring to fig. 3, an embodiment of the present invention provides a method for analyzing spatial development laws of beads, including:

and step S3: determining a longitudinal development rule of the bead bodies based on the high-precision three-dimensional seismic amplitude data volume of the work area and the determined development positions of the bead bodies;

wherein, confirm the longitudinal development law of string of beads body includes: respectively determining optimal development windows of the bead bodies corresponding to all the layers of the work area, respectively determining the sum of first amplitude change rates of the optimal development windows of all the layers, and determining the dominant development layer of the bead bodies according to the numerical value of the sum of the first amplitude change rates of all the layers;

wherein, confirm the transverse development law of string of beads body includes: respectively determining the optimal development width of the bead body corresponding to each sliding fracture of the work area, respectively determining the sum of second amplitude change rates within the optimal development width of each sliding fracture, and determining the main control sliding fracture of the bead body according to the numerical value of the sum of the second amplitude change rates of each sliding fracture;

wherein, the sequence of the steps S3 and S4 can be adjusted.

An embodiment of the invention provides a method for analyzing spatial development laws of bead bodies, which comprises the following steps:

and step S4: determining a transverse development rule of the bead bodies based on the high-precision three-dimensional seismic amplitude data volume of the work area and the determined development positions of the bead bodies;

wherein, confirm the horizontal development law of string of beads body includes: respectively determining the optimal development width of the bead body corresponding to each sliding fracture in the dominant development layer of the bead body, respectively determining the sum of the second amplitude change rates of each sliding fracture in the optimal development width, and determining the main control sliding fracture of the bead body according to the numerical value of the sum of the second amplitude change rates of each sliding fracture.

At present, the seismic data body in the field of oil exploration adopts the spacing of 25m multiplied by 25m, collects 0-6s of reflected wave amplitude data, and adopts the vertical sampling interval of 2 ms. In the bead body space development rule analysis method provided by the embodiment of the invention, according to the technical requirements for identifying the three-dimensional transmission and conduction system, the main frequency of the seismic data amplitude data body in the oil and gas reservoir development interval should exceed the preset main frequency, and the lowest preset main frequency in the technology should be more than 50Hz, and the preset main frequency can be usually set to be 50Hz-60Hz.

Referring to FIG. 4, in one embodiment, acquiring a high-precision three-dimensional seismic amplitude data volume for a work area in step S1 may further comprise the steps of:

s11: acquiring a three-dimensional seismic amplitude data volume of a work area;

s12: judging whether the dominant frequency of the three-dimensional seismic amplitude data body of the work area is lower than the preset dominant frequency or not;

s13: if the dominant frequency of the three-dimensional seismic amplitude data body of the work area is not lower than the preset dominant frequency, the three-dimensional seismic amplitude data body of the work area is the high-precision three-dimensional seismic amplitude data body of the work area;

s14: if the dominant frequency of the three-dimensional seismic amplitude data body of the work area is lower than the preset dominant frequency, performing frequency extension processing on the three-dimensional seismic amplitude data body of the work area to generate a high-precision three-dimensional seismic amplitude data body of the work area, of which the dominant frequency is not lower than the preset dominant frequency;

the Frequency broadening processing of the three-dimensional seismic data on the three-dimensional seismic amplitude data volume of the work area may be performed by using a method such as an ultra High resolution Imaging (HFI) method.

In step S2, determining the development position of each bead body based on the high-precision three-dimensional seismic amplitude data volume of the work area may be performed by a method commonly used in the art. In addition, the more preferable scheme is to determine the development position of the bead body by adopting a quantitative analysis mode, and the method specifically comprises the following steps:

calculating the amplitude change rate of different parts of each layer in the work area, wherein the part of which the amplitude change rate is not lower than a first preset threshold value is a bead body development position;

wherein, the calculating the amplitude change rate of different parts of each layer in the work area may further include the steps of: and calculating the amplitude change rate of the upper part and the lower part, the east part and the west part, the south part and the north part of the earthquake reflection point by taking the earthquake reflection point as a center.

For example, the position of bead development is determined by: calculating the amplitude change rate of the upper part and the lower part, the east part and the west part, the south part and the north part of the earthquake reflection point by taking the earthquake reflection point as a center; wherein, the upper and lower computing windows are 3-5, the east and west computing windows are 10-20, and the south and north computing windows are 10-20; the part with the amplitude change rate not lower than a first preset threshold value is the bead body development position, wherein the first preset threshold value is 9200-9800.

In an embodiment, in step S3, the step of determining the optimal development windows of the bead bodies corresponding to each layer of the work area includes the following steps: aiming at each layer in the work area, respectively calculating the distance between the top of each bead body and the layer, respectively taking each layer as a reference, calculating the number Nwi of the bead bodies within the range from one layer to different depths W ms,

and when the maximum value is obtained, the value of W is the optimal development window of the bead body corresponding to the layer.

Wherein the calculating the number of beads Nwi ranging from one horizon to different depths W ms may be performed as follows: setting an initial value of search window length as W = W ₀ From the horizon to W = W ₀ Number of ms-deep beads Nwi = Nw ₀ And calculate

Until the depth reaches the top surface of the next horizon.

In an embodiment, in step S3, the step of determining the dominant development horizon of the bead body according to the magnitude of the sum of the first amplitude change rates of the horizons specifically includes the following steps: and determining the first N layers with the maximum value of the sum of the first amplitude change rates as the dominant development layers of the bead bodies.

In one embodiment, the bead body space development law analysis is performed on the raised Ordovician carbonate fracture-cave reservoir in the Tarim basin tower, and the step S3 specifically comprises the following steps: aiming at the characteristic that the bead body is obviously strongly reflected by earthquake, longitudinal development law analysis is carried out on non-integration surfaces and related secondary non-integration surfaces of positions such as a heavy point position (a good interior lattice group (a top surface of the good interior lattice group, a top surface of a good section 3, a top surface of a good section 5), a room group, a eagle mountain group (a top surface of the eagle mountain group, a top surface of a eagle section 2, a top surface of a eagle section 4), a Paglai dam group and the like in the Ordovician system, and 8 positions are totally included);

respectively determining optimal development windows of the bead bodies corresponding to the various layers by taking the 8 main layers as a reference; taking the main horizon H1 as an example for explanation, calculating the distance between the top of each bead and the main horizon H1, setting the length W ms of a search window by taking the three-dimensional depth of H1 as a reference, calculating the number Nwi of beads in the W ms depth range of the H1 stratum, and calculating different Nwi obtained along with the increase of W according to a certain step length (1-5 ms) until the depth reaches the top surface of the next horizon; determining

Is greater than or equal to>

When the maximum value is W, the value of W is the optimal development window of the bead body corresponding to the layer; correspondingly, 8 main layers are calculated to obtain the optimal development window of the bead body corresponding to each layer; />

And calculating the sum Ei of the amplitude change rates of the optimal development windows of all the layers, and preferably selecting 3 layers with the largest Ei value as the dominant development layers of the bead body to serve as key layers of oil-gas exploration, thereby completing the analysis of the longitudinal development rule of the bead body.

In one embodiment, in step S3, determining the longitudinal development rule of the bead body further includes: and further determining at least one of the top development depth, the bottom development depth, the east-west development width, the north-south development width, the maximum amplitude change rate value, the oil production of the drilling encountering beading body, the gas production of the drilling encountering beading body, the water production of the drilling encountering beading body and the oil-gas-water geochemical parameters of the drilling encountering beading body on the basis of the dominant development horizon of the beading body.

In an embodiment, in step S4, the step of determining the optimal development width of the bead body corresponding to each slipping fracture specifically includes: aiming at each sliding fracture, respectively calculating the distance between each bead body and the sliding fracture, taking the sliding fracture as a reference, calculating the number Nli of the bead bodies extending from one sliding fracture to two sides within the range of L ms,

wherein, the number Nli of the bead bodies extending from one sliding fracture to two sides within the range of L ms can be calculated by the following method: setting an initial value of a search window length to L = L ₀ Calculating the extension L = L from one sliding break to both sides ₀ Number of beads in ms range Nli = Nl ₀ And calculate

Until it extends sideways from the trip break to another trip break.

In an embodiment, in step S4, the determining the main control sliding fracture of the bead body according to the magnitude of the sum of the second amplitude change rates of the sliding fractures specifically includes: and determining the first N sliding fractures with the largest numerical value of the sum of the second amplitude change rates as the main control sliding fractures of the bead body.

In one embodiment, the bead body space development law analysis is performed on the raised Ordovician carbonate fracture-cave reservoir in the Tarim basin tower, and the step S4 specifically comprises the following steps:

respectively calculating the optimal development width of the bead body corresponding to each sliding fracture by taking each sliding fracture in the work area as a reference; taking a certain sliding fracture F1 as an example, calculating each bead body and slidingSetting the length L ms of a search window by taking the glide fracture F1 as a reference, calculating the number Nli of bead strings in the length range of L ms extending from the glide fracture F1 to two sides, calculating different Nli obtained along with the increase of L according to a certain step length (10-50 ms), and determining the nearest distance of the fracture F1

The maximum value of (a) is,

the value of L is the optimal development width of the bead body corresponding to the sliding fracture when the value is the maximum value; correspondingly, all the sliding fractures are calculated to obtain the optimal development width of the bead body corresponding to each sliding fracture;

and calculating the sum Ei of the amplitude change rates of the bead bodies within the optimal development width of the sliding fractures, and preferably selecting 5-15 sliding fractures with the largest Ei value as the main control sliding fractures of the bead bodies so as to finish the analysis of the longitudinal development rule of the bead bodies.

In an embodiment, in step S4, the optimal development width of the bead body corresponding to each large-scale sliding fracture in the dominant development layer of the bead body may be respectively determined, the sum of the second amplitude change rates within the optimal development width of each large-scale sliding fracture is respectively determined, and the main control sliding fracture of the bead body is determined according to the numerical value of the sum of the second amplitude change rates of each large-scale sliding fracture; wherein, the large-scale sliding refers to that the first M sliding parts with the largest length in all sliding parts are broken.

Another embodiment of the present invention provides a method for analyzing a bead spatial development law of a seam-cave reservoir of an aoto-series carbonate rock raised in a talimu basin tower, the method including:

(1) Acquiring a high-precision three-dimensional seismic amplitude data volume of a work area:

acquiring a three-dimensional seismic amplitude data volume of a work area;

judging whether the dominant frequency of the three-dimensional seismic amplitude data body of the work area is lower than the preset dominant frequency or not; the preset main frequency is 50Hz-60Hz;

if the dominant frequency of the three-dimensional seismic amplitude data body of the work area is lower than the preset dominant frequency, performing frequency extension processing on the three-dimensional seismic amplitude data body of the work area to generate a high-precision three-dimensional seismic amplitude data body of the work area, wherein the dominant frequency is not lower than the preset dominant frequency;

(2) Determining the development position of each bead body based on the high-precision three-dimensional seismic amplitude data body of the work area:

calculating the amplitude change rate of the upper part and the lower part, the east part and the west part, the south part and the north part of the earthquake reflection point by taking the earthquake reflection point as a center; wherein, the upper and lower computing windows are 3-5, the east and west computing windows are 10-20, and the south and north computing windows are 10-20; the part with the amplitude change rate not lower than a first preset threshold value is a bead body development position, wherein the first preset threshold value is 9200-9800;

(3) Based on the high-precision three-dimensional seismic amplitude data volume of the work area and the determined development positions of the bead bodies, determining the longitudinal development rule of the bead bodies:

the development key position of the engineering region of the raised Ordovician carbonate rock fracture-cave reservoir in the Tarim basin tower (8 positions including the unconformity surfaces of the good Rier grid group (the top surface of the good Rier grid group, the top surface of the good 3 sections and the top surface of the good 5 sections) in the Ordovician, a room group, the eagle mountain group (the top surface of the eagle mountain group, the top surface of the eagle 2 sections and the top surface of the eagle 4 sections), the Paglai dam group and other positions and the related secondary unconformity surfaces); analyzing longitudinal development rules of the bead bodies according to the 8 layers;

respectively determining optimal development windows of the bead bodies corresponding to the various layers by taking the 8 main layers as a reference; taking the main horizon H1 as an example for explanation, calculating the distance between the top of each bead and the main horizon H1, setting the length W ms of a search window by taking the three-dimensional depth of H1 as a reference, calculating the number Nwi of beads in the W ms depth range of the H1 stratum, and counting according to a certain step length (1-5 ms is adopted)Calculating different Nwi obtained along with the increase of W until the depth reaches the top surface of the next horizon; determining

Is greater than or equal to>

When the maximum value is W, the value of W is the optimal development window of the bead body corresponding to the layer; correspondingly, 8 main layers are calculated to obtain the optimal development window of the bead body corresponding to each layer;

calculating the sum Ei of the amplitude change rates of the optimal development windows of all the layers, preferably selecting 3 layers with the largest Ei value as the dominant development layers of the bead body as key layers of oil-gas exploration, and thus completing the analysis of the longitudinal development rule of the bead body;

(4) Determining the transverse development rule of the bead bodies in the dominant development layer determined in the step (3) based on the high-precision three-dimensional seismic amplitude data body of the work area and the determined development positions of the bead bodies:

respectively calculating the optimal development width of the bead body corresponding to each large-scale gliding fracture by taking each large-scale gliding fracture in the dominant development horizon determined in the step (3) in the work area as a reference; taking a certain large-scale sliding fracture F1 as an example, calculating the closest distance between each bead body and the large-scale sliding fracture F1, setting the length L ms of a search window based on the large-scale sliding fracture F1, calculating the number Nli of bead bodies within the length range of extending L ms from the large-scale sliding fracture F1 to two sides, calculating different Nli obtained along with the increase of L according to a certain step length (10-50 ms), and determining

Is greater than or equal to>

The value of L is the optimal development width of the bead body corresponding to the large-scale sliding fracture when the value is the maximum value; correspondingly, all large-scale sliding fractures are calculated to obtain the corresponding large-scale sliding fractureOptimal development width of the bead body; wherein the large-scale sliding refers to the first 15 sliding fractures with the largest length in all sliding in the dominant development horizon determined in the step (3) in the work area.

And calculating the sum Ei of the amplitude change rates of the bead bodies within the optimal development width of the large-scale sliding fractures, and preferably selecting 5 large-scale sliding fractures with the largest Ei value as the main control sliding fractures of the bead bodies so as to complete the analysis of the longitudinal development rule of the bead bodies.

The embodiment of the invention also provides a system for analyzing the spatial development law of the bead bodies, and preferably the system is used for realizing the method embodiment.

Fig. 5 is a block diagram of a system for analyzing spatial development laws of beads according to an embodiment of the present invention, and as shown in fig. 5, the apparatus includes: high-precision three-dimensional seismic amplitude data volume acquisition module 51, bead body development position positioning module 52 and bead body space development rule analysis module 53

The high-precision three-dimensional seismic amplitude data volume acquisition module 51: the method comprises the steps of obtaining a high-precision three-dimensional seismic amplitude data volume of a work area; the high-precision three-dimensional seismic amplitude data body is a three-dimensional seismic amplitude data body with a main frequency not lower than a preset main frequency, and the preset main frequency is not lower than 50Hz;

the bead body development position positioning module 52: the system is used for determining the development position of each bead body based on the high-precision three-dimensional seismic amplitude data body of the work area acquired by the acquisition module;

the bead body space development rule analysis module 53: for determining the development position of each bead body based on the high-precision three-dimensional seismic amplitude data body of the work area acquired by the acquisition module and the development position of each bead body determined by the positioning module, determining a spatial development rule of the bead body; wherein the determining of the spatial development law of the bead body comprises determining the longitudinal development law of the bead body and/or determining the transverse development law of the bead body;

the method for determining the longitudinal development rule of the bead body comprises the following steps: respectively determining optimal development windows of the bead bodies corresponding to all the layers of the work area, respectively determining the sum of first amplitude change rates of the optimal development windows of all the layers, and determining the dominant development layer of the bead bodies according to the numerical value of the sum of the first amplitude change rates of all the layers;

In the high-precision three-dimensional seismic amplitude data volume acquisition module 51, the preset main frequency may be 50Hz to 60Hz.

With continued reference to fig. 6, the high-precision three-dimensional seismic amplitude data volume acquisition module 51 may include:

the first acquisition unit 511: the method comprises the steps of obtaining a three-dimensional seismic amplitude data volume of a work area;

first discrimination unit 512: judging whether the dominant frequency of the three-dimensional seismic amplitude data of the work area acquired by the first acquisition unit 511 is lower than the preset dominant frequency; if the first judging unit 512 judges that the main frequency of the three-dimensional seismic amplitude data body of the work area is not lower than the preset main frequency, a first processing unit 513 is adopted for processing; if the dominant frequency of the three-dimensional seismic amplitude data body of the work area is lower than the preset dominant frequency, processing by using a second processing unit 514;

the first processing unit 513: for performing operations that use the three-dimensional seismic amplitude data volume of the work area as a high-precision three-dimensional seismic amplitude data volume of the work area.

The second processing unit 514: and the frequency extension processing module is used for carrying out frequency extension processing on the three-dimensional seismic amplitude data body of the work area to generate a high-precision three-dimensional seismic amplitude data body of the work area with the main frequency not lower than the preset main frequency.

The second processing unit 514 may perform Frequency broadening processing on the three-dimensional seismic amplitude data volume of the work area by using an ultra High resolution processing method (HFI).

With continued reference to fig. 7, the bead development position locating module 52 may include:

the third processing unit 521: the system is used for calculating the amplitude change rate of different parts of each layer in the work area based on the high-precision three-dimensional seismic amplitude data volume of the work area acquired by the acquisition module;

second determination section 522: the amplitude change rate of different parts of each layer in the work area calculated by the third processing unit 521 is lower than the first predetermined threshold; the part with the amplitude change rate not lower than the first preset threshold value is the development position of the bead body.

The third processing unit 521 is configured to calculate amplitude change rates of different positions of each layer in the work area based on the high-precision three-dimensional seismic amplitude data volume of the work area acquired by the acquisition module; and the calculating the amplitude change rate of different parts of each layer in the work area comprises: and calculating the amplitude change rate of the upper part and the lower part, the east part and the west part, the south part and the north part of the earthquake reflection point by taking the earthquake reflection point as a center.

In one embodiment, the bead development position locating module 52 includes:

the third processing unit 521: the system is used for calculating the amplitude change rate of different parts of each layer in the work area based on the high-precision three-dimensional seismic amplitude data volume of the work area acquired by the acquisition module; calculating the amplitude change rate of different parts of each layer in the work area by taking the seismic reflection point as a center, and calculating the amplitude change rate of three different dimensions of the upper part, the lower part, the east part, the west part, the south part and the north part of the seismic reflection point; wherein, the upper and lower computing windows are 3-5, the east and west computing windows are 10-20, and the south and north computing windows are 10-20;

second determination section 522: the amplitude change rate of different parts of each layer in the work area calculated by the third processing unit 521 is lower than the first predetermined threshold; the part of which the amplitude change rate is not lower than a first preset threshold value is a bead body development position; wherein the first predetermined threshold is 9200-9800.

With continued reference to fig. 8, the bead spatial developmental law analysis module 53 may include:

the fourth processing unit 531: for dividingDetermining optimal development windows of the bead bodies corresponding to all the layers of the work area; aiming at each layer in the work area, respectively calculating the distance between the top of each bead body and the layer, respectively taking each layer as a reference, calculating the number Nwi of the bead bodies within the range from one layer to different depths W ms,

the fifth processing unit 532: the fourth processing unit 531 is configured to calculate a sum of first amplitude change rates of the optimal development windows of the layers based on the optimal development windows of the beads corresponding to the layers of the work area determined by the fourth processing unit 531;

the sixth processing unit 533: for determining dominant developmental horizons of beads based on the magnitude of the sum of said first amplitude change rates of the horizons calculated by the fifth processing unit 532.

Further, the calculating a number of beads Nwi in a range of W ms from a horizon to a different depth may comprise: setting an initial value of search window length as W = W ₀ From the horizon to W = W ₀ Number of beads of ms depth Nwi = Nw ₀ And calculate

Until the depth reaches the top surface of the next horizon.

Further, the sixth processing unit 533 may be configured to determine, based on the sum of the first amplitude change rates of the respective layers calculated by the fifth processing unit 532, the first N layers with the largest value of the sum of the first amplitude change rates as dominant development layers of the bead body.

Further, referring to fig. 9, the bead spatial development law analysis module 53 may include: the seventh processing unit 534: for further determining at least one of the top development depth, the bottom development depth, the east-west development width, the north-south development width, the maximum amplitude change rate value, the drilling bead oil production, the drilling bead gas production, the drilling bead water production, and the drilling bead geochemical parameters of the bead based on the dominant development horizon of the bead determined by the sixth processing unit 533.

With continued reference to fig. 10, the bead spatial developmental law analysis module 53 may include:

the eighth processing unit 535: the bead body optimal development width corresponding to each sliding fracture of the work area is respectively determined; aiming at each sliding fracture, respectively calculating the distance between each bead body and the sliding fracture, taking the sliding fracture as a reference, calculating the number Nli of the bead bodies extending from one sliding fracture to two sides within the range of L ms,

the ninth processing unit 536: for calculating a sum of second amplitude change rates within the optimal development width of each slipping fracture based on the optimal development width of the bead body corresponding to each slipping fracture of the work area determined by the eighth processing unit 535;

the tenth processing unit 537: for determining a dominant development horizon for the bead based on the magnitude of the sum of said second amplitude change rates for each gliding break calculated by the ninth processing unit 536.

Further, the calculating the number Nli of beads extending from a walking slip fracture to both sides in a range of L ms includes: setting an initial value of a search window length to L = L ₀ Calculating the extension L = L from one sliding break to both sides ₀ Number of beads in ms range Nli = Nl ₀ And calculate

Generated in certain steps (e.g. 1-5 ms)New L value, calculating the number Nli of beads extending from a glide break to both sides in the range of L ms, and determining ^ H>

Until it extends sideways from the trip break to another trip break.

Further, the tenth processing unit 537 is configured to determine, based on the sum of the second amplitude change rates of the respective sliding fractures calculated by the ninth processing unit 536, the first N sliding fractures with the largest value of the sum of the second amplitude change rates as the main sliding fractures of the bead body.

For specific implementation processes of the units and the modules, reference may be made to the description of the method embodiment, and details are not described herein again.

Fig. 11 is a schematic diagram of a device for analyzing spatial development laws of beads according to an embodiment of the present invention. The device for analyzing spatial development rules of beads shown in fig. 11 is a general-purpose data processing device, which includes a general-purpose computer hardware structure including at least a processor 1000 and a memory 1111; the processor 1000 is configured to execute the molecular structure generating program stored in the memory to implement the method for analyzing the spatial development law of a bead body according to the embodiments of the method (for a specific method, refer to the description of the embodiments of the method, and no further description is given here).

The embodiment of the present invention further provides a computer-readable storage medium, where one or more programs are stored in the storage medium, and the one or more programs may be executed by one or more processors to implement the method for analyzing spatial development laws of bead bodies according to the embodiments of the methods (for a specific method, refer to the description of the above method embodiments, and are not described herein again).

The preferred embodiments of the present invention have been described above with reference to the accompanying drawings. The many features and advantages of the embodiments are apparent from the detailed specification, and thus, it is intended by the appended claims to cover all such features and advantages of the embodiments which fall within the true spirit and scope thereof. Further, since numerous modifications and changes will readily occur to those skilled in the art, it is not desired to limit the embodiments of the invention to the exact construction and operation illustrated and described, and accordingly, all suitable modifications and equivalents may be resorted to, falling within the scope thereof.

As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.

Claims

1. A method for analyzing spatial development law of bead bodies, wherein the method comprises the following steps:

determining the development position of each bead body based on the high-precision three-dimensional seismic amplitude data volume of the work area;

the determining of the lateral development law of the bead body comprises: respectively determining the optimal development width of the bead body corresponding to each sliding fracture of the work area, respectively determining the sum of second amplitude change rates within the optimal development width of each sliding fracture, and determining the main control sliding fracture of the bead body according to the numerical value of the sum of the second amplitude change rates of each sliding fracture;

the optimal development window of the bead bodies respectively determining each layer position of the work area to correspond comprises: aiming at each layer in the work area, respectively calculating the distance between the top of each bead body and the layer, respectively taking each layer as a reference, calculating the number Nwi of the bead bodies within the range from one layer to different depths W ms,

the step of respectively determining the optimal development width of the bead body corresponding to each sliding fracture of the work area comprises the following steps: aiming at each sliding fracture, respectively calculating the distance between each bead body and the sliding fracture, taking the sliding fracture as a reference, calculating the number Nli of the bead bodies extending from one sliding fracture to two sides within the range of L ms,

and the value of L is the optimal development width of the bead body corresponding to the sliding fracture when the value is the maximum value.

2. The method of claim 1, wherein the acquiring a high-precision three-dimensional seismic amplitude data volume for a work zone comprises:

acquiring a three-dimensional seismic amplitude data volume of a work area;

and if the dominant frequency of the three-dimensional seismic amplitude data body of the work area is lower than the preset dominant frequency, performing frequency extension processing on the three-dimensional seismic amplitude data body of the work area to generate a high-precision three-dimensional seismic amplitude data body of the work area, wherein the dominant frequency is not lower than the preset dominant frequency.

3. The method of claim 1, wherein the determining a developmental location of a bead comprises:

and calculating the amplitude change rate of different parts of each layer in the work area, wherein the part with the amplitude change rate not lower than a first preset threshold value is the development position of the bead body.

4. The method of claim 3, wherein said calculating the rate of change of amplitude at different locations of each level in the work area comprises: and calculating the amplitude change rate of the upper part and the lower part, the east part and the west part, the south part and the north part of the earthquake reflection point by taking the earthquake reflection point as a center.

5. The method of claim 1, wherein the determining a dominant developmental horizon for a bead as a function of a numerical magnitude of a sum of the first rates of amplitude change for the horizons comprises: and determining the first N layers with the maximum value of the sum of the first amplitude change rates as the dominant development layers of the bead bodies.

6. The method of claim 1 or 5, wherein the determining longitudinal development regularity of beads further comprises:

7. The method of claim 1, wherein the determining the dominant slip fracture of a bead from the magnitude of the sum of the second rate of amplitude change of each slip fracture comprises: and determining the first N sliding fractures with the maximum numerical value of the sum of the second amplitude change rates as the main control sliding fractures of the bead body.

8. A system for analyzing spatial development laws of a beading body, wherein the system comprises:

the bead body space development rule analysis module: the system comprises an acquisition module, a positioning module, a processing module and a control module, wherein the acquisition module is used for acquiring a high-precision three-dimensional seismic amplitude data volume of the work area and a development position of each bead body determined by the positioning module; wherein the determining of the spatial development law of the bead body comprises determining the longitudinal development law of the bead body and/or determining the transverse development law of the bead body;

the optimal development window of the bead bodies respectively determining each layer position of the work area to correspond comprises: aiming at each layer in the work area, the distance between the top of each bead body and the layer is respectively calculated, and the distance between the top of each bead body and the layer is respectively calculated by taking each layer as a referenceThe number Nwi,

and the value of L is the optimal development width of the bead body corresponding to the walking and sliding fracture when the value is the maximum value.

9. The system of claim 8, wherein the high-precision three-dimensional seismic amplitude data volume acquisition module comprises:

a first processing unit: performing an operation of using the three-dimensional seismic amplitude data volume of the work area as a high-precision three-dimensional seismic amplitude data volume of the work area;

10. The system of claim 8, wherein the bead development position location module comprises:

a second determination unit: the amplitude change rate of different parts of each layer in the work area calculated by the third processing unit is judged whether to be lower than a first preset threshold value or not; the part with the amplitude change rate not lower than the first preset threshold value is the development position of the bead body.

11. The system of claim 10, wherein said calculating the rate of change of amplitude at different locations of each level in the work area comprises: and calculating the amplitude change rate of the upper part and the lower part, the east part and the west part, the south part and the north part of the earthquake reflection point by taking the earthquake reflection point as a center.

12. The system of claim 8, wherein the bead spatial developmental law analysis module comprises:

a sixth processing unit: the dominant development horizon of the bead bodies is determined based on the numerical value of the sum of the first amplitude change rates of the horizons calculated by the fifth processing unit.

13. The system of claim 12, wherein the sixth processing unit is configured to determine, based on the sum of the first amplitude change rates of the respective horizons calculated by the fifth processing unit, the first N horizons with the largest value of the sum of the first amplitude change rates as dominant developmental horizons of a bead.

14. The system of claim 8, wherein the bead spatial developmental law analysis module further comprises:

15. The system of claim 8, wherein the bead spatial developmental law analysis module comprises:

an eighth processing unit: the optimal development width of the bead body corresponding to each sliding fracture of the work area is respectively determined; aiming at each sliding fracture, respectively calculating the distance between each bead body and the sliding fracture, taking the sliding fracture as a reference, calculating the number Nli of the bead bodies extending from one sliding fracture to two sides within the range of L ms,

16. The system according to claim 15, wherein the tenth processing unit is configured to determine, based on the sum of the second amplitude change rates of the respective skid fractures calculated by the ninth processing unit, the first N skid fractures with the largest value of the sum of the second amplitude change rates as the main skid fractures of the bead body.

17. A device for analyzing the spatial development law of a bead body comprises a processor and a memory; wherein the content of the first and second substances,

a memory for storing a computer program;

a processor for implementing the steps of the method for analyzing spatial laws of bead body according to any one of claims 1-7 when executing the program stored in the memory.

18. A computer readable storage medium storing one or more programs, the one or more programs being executable by one or more processors to perform the steps of the method of spatial law analysis of bead bodies according to any one of claims 1 to 7.