CN116931083A - Determination method of azimuth angle gather division scheme - Google Patents
Determination method of azimuth angle gather division scheme Download PDFInfo
- Publication number
- CN116931083A CN116931083A CN202310930152.0A CN202310930152A CN116931083A CN 116931083 A CN116931083 A CN 116931083A CN 202310930152 A CN202310930152 A CN 202310930152A CN 116931083 A CN116931083 A CN 116931083A
- Authority
- CN
- China
- Prior art keywords
- azimuth
- gather
- angle
- dividing
- determining
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000011156 evaluation Methods 0.000 claims abstract description 48
- 238000000638 solvent extraction Methods 0.000 claims description 36
- 238000011158 quantitative evaluation Methods 0.000 abstract description 8
- 238000004590 computer program Methods 0.000 description 15
- 238000004891 communication Methods 0.000 description 8
- 238000001514 detection method Methods 0.000 description 8
- 238000005192 partition Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 239000003795 chemical substances by application Substances 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 239000004973 liquid crystal related substance Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Abstract
The invention discloses a method for determining an azimuth angle gather dividing scheme, which comprises the following steps: determining a plurality of initial dividing angles of the seismic data according to a preset angle dividing step length; dividing the seismic data into a plurality of azimuth angle gather data according to the initial dividing angles and a preset dividing rule for each initial dividing angle, thereby obtaining a plurality of azimuth angle gather dividing schemes to be evaluated; for each azimuth gather dividing scheme to be evaluated, acquiring seismic characteristic data corresponding to each azimuth gather data, and determining azimuth anisotropic crack prediction precision evaluation factors based on the seismic characteristic data; and determining a target azimuth angle gather dividing scheme based on azimuth anisotropic crack prediction precision evaluation factors corresponding to each azimuth angle gather dividing scheme to be evaluated. The invention provides a quantitative evaluation method for an azimuth angle gather dividing scheme of seismic data, which realizes the rapid and efficient determination of the optimal azimuth angle gather dividing scheme of the seismic data, thereby improving the precision of crack prediction.
Description
Technical Field
The embodiment of the invention relates to the technical field of oil and gas geophysics, in particular to a method for determining an azimuth angle gather dividing scheme.
Background
In azimuth anisotropic fracture prediction based on wide azimuth seismic data, the wide azimuth seismic data is first divided into a plurality of azimuth gathers to perform fracture prediction based on the divided azimuth gathers. Because the longitudinal and transverse data acquisition times of the wide-azimuth seismic data are not uniform, the data acquisition times corresponding to azimuth angle gathers determined by different data dividing schemes are also different. Based on this, the division manner of the wide-azimuth seismic data is different, and the accuracy of crack prediction based on the divided seismic data is also different, so it is important to determine a seismic data division scheme with high prediction accuracy.
At present, a plurality of azimuth angle gathers of wide azimuth seismic data are usually obtained according to working experience of professionals, crack prediction is carried out based on the divided azimuth angle gathers, after the completion of the crack prediction, a prediction result is compared with logging data, and finally the rationality of the seismic data division is evaluated according to the precision of the crack prediction result. However, prior to fracture prediction, an objective and effective evaluation method is lacking to determine reasonable seismic data angle gather division schemes, which results in low accuracy of fracture prediction results.
Disclosure of Invention
The embodiment of the invention provides a method for determining an azimuth angle gather dividing scheme, and provides a quantitative evaluation method for the seismic data azimuth angle gather dividing scheme, so that the method for rapidly and efficiently determining the optimal seismic data azimuth angle gather dividing scheme is realized, and the accuracy of crack prediction is improved.
In a first aspect, the present invention provides a method for determining an azimuth gather partitioning scheme, the method comprising:
determining a plurality of initial dividing angles corresponding to the seismic data according to a preset angle dividing step length;
dividing the seismic data into a plurality of azimuth angle gather data according to the initial dividing angles and a preset dividing rule for each initial dividing angle;
determining an angle gather dividing scheme to be evaluated for the multiple azimuth angle gather data corresponding to each initial dividing angle;
obtaining seismic characteristic data corresponding to the azimuth angle gather data for each scheme of dividing the angle gather to be evaluated, and determining azimuth anisotropic fracture prediction precision evaluation factors based on the seismic characteristic data;
and determining a target angle gather dividing scheme based on azimuth anisotropic crack prediction precision evaluation factors corresponding to each angle gather dividing scheme to be evaluated.
In a second aspect, the present invention provides an apparatus for determining an azimuth gather partitioning scheme, the apparatus comprising:
the initial angle determining module is used for determining a plurality of initial dividing angles corresponding to the seismic data according to a preset angle dividing step length;
the angle gather dividing module is used for dividing the seismic data into a plurality of azimuth angle gather data according to the initial dividing angles and a preset dividing rule for each initial dividing angle;
the division scheme determining module is used for determining an angle gather division scheme to be evaluated for the plurality of azimuth angle gather data corresponding to each initial division angle;
the evaluation factor determining module is used for obtaining seismic characteristic data corresponding to the azimuth angle gather data for each to-be-evaluated angle gather division scheme and determining azimuth anisotropic fracture prediction precision evaluation factors based on the seismic characteristic data;
and the target scheme determining unit is used for determining a target angle gather dividing scheme based on azimuth anisotropic crack prediction precision evaluation factors corresponding to each angle gather dividing scheme to be evaluated.
In a third aspect, the present invention provides an electronic device, comprising:
at least one processor; and
a memory communicatively coupled to the at least one processor; wherein, the liquid crystal display device comprises a liquid crystal display device,
the memory stores a computer program executable by the at least one processor to enable the at least one processor to perform the method of determining the azimuth gather partitioning scheme of any one of the embodiments of the present invention.
In a fourth aspect, the present invention provides a computer readable storage medium storing computer instructions for causing a processor to perform a method of determining an azimuth gather partitioning scheme according to any one of the embodiments of the present invention.
According to the technical scheme provided by the embodiment of the invention, a plurality of initial dividing angles corresponding to the seismic data are determined according to the preset angle dividing step length, and then the seismic data are divided into a plurality of azimuth angle gather data according to the initial dividing angles and the preset dividing rules for each initial dividing angle, and further, the azimuth angle gather dividing scheme to be evaluated is determined for the plurality of azimuth angle gather data corresponding to each initial dividing angle. And for each azimuth angle gather division scheme to be evaluated, acquiring seismic characteristic data corresponding to each azimuth angle gather data, determining azimuth anisotropy crack prediction precision evaluation factors based on the seismic characteristic data, and finally determining a target azimuth angle gather division scheme based on the azimuth anisotropy crack prediction precision evaluation factors corresponding to each azimuth angle gather division scheme to be evaluated. The invention provides a quantitative evaluation method for an azimuth angle gather dividing scheme of seismic data, which realizes the rapid and efficient determination of the optimal azimuth angle gather dividing scheme of the seismic data, thereby improving the accuracy of crack prediction.
It should be understood that the description in this section is not intended to identify key or critical features of the embodiments of the invention or to delineate the scope of the invention. Other features of the present invention will become apparent from the description that follows.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of a method for determining an azimuth gather partitioning scheme according to an embodiment of the present invention;
FIG. 2 is a flowchart of a method for determining an azimuth gather partitioning scheme according to a second embodiment of the present invention;
FIG. 3 is a schematic diagram of a determining apparatus for an azimuth gather partitioning scheme according to a third embodiment of the present invention;
fig. 4 is a schematic structural diagram of an electronic device according to a fourth embodiment of the present invention.
Detailed Description
In order that those skilled in the art will better understand the present invention, a technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in which it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present invention without making any inventive effort, shall fall within the scope of the present invention.
It should be noted that, in the description and claims of the present invention and the above figures, the terms "first preset condition", "second preset condition", and the like are used to distinguish similar objects, and are not necessarily used to describe a specific order or precedence. It is to be understood that the data so used may be interchanged where appropriate such that the embodiments of the invention described herein may be implemented in sequences other than those illustrated or otherwise described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
Example 1
Fig. 1 is a flowchart of a method for determining an azimuth gather partitioning scheme according to an embodiment of the present invention, where the embodiment is applicable to a situation of performing quantitative evaluation on an azimuth gather partitioning scheme of seismic data when performing azimuth anisotropic fracture prediction based on the seismic data. The method may be performed by a means of determining the azimuth gather partitioning scheme, which may be implemented in hardware and/or software, which may be configured on a computer device, which may be a notebook, a desktop computer, a smart tablet, etc. As shown in fig. 1, the method includes:
s110, determining a plurality of initial dividing angles corresponding to the seismic data according to the preset angle dividing step length.
The preset angle dividing step length is a preset angle value and is used for representing how many angle values are used for determining an initial dividing angle. For example, the preset angle division step is 5 °. The initial division angle may be understood as a division start point when dividing the azimuth gather. The specific number of the initial dividing angles is not particularly limited, and the initial dividing angles can be adaptively adjusted according to actual requirements. The number of azimuth angle gather dividing schemes corresponds to the number of initial dividing angles, and the number of azimuth angle gather dividing schemes is the number of initial dividing angles.
In this embodiment, an OVT domain gather corresponding to the original seismic data set may be generated by using an existing offset vector slice processing technique (Offset Vector Tile, OVT), and in a subsequent application process, seismic data may be acquired from the OVT domain gather.
In an application scenario of azimuth anisotropy crack prediction based on seismic data, one of a large number of signal detection points is taken as a target detection point. And carrying out azimuth anisotropy crack prediction according to the seismic data corresponding to the target detection point, and representing azimuth anisotropy crack conditions of the work area to be predicted. Therefore, when the quantitative evaluation of the azimuth angle gather dividing scheme is carried out, only the seismic data corresponding to the target detection point is acquired from the OVT domain gather corresponding to the original seismic data set to evaluate the azimuth angle gather dividing scheme. The seismic data includes incident angle, offset, azimuth angle and other data of all seismic traces corresponding to the target detection point.
Optionally, determining the plurality of initial division angles specifically includes the steps of: determining a reference zero angle from the seismic data; and taking the reference zero angle as a dividing starting point, dividing step length at intervals of preset angles, and determining an initial dividing angle to obtain a plurality of initial dividing angles corresponding to the seismic data.
Specifically, since the analog source emits a signal from 360 ° of the target detection point in all directions, in order to clearly distinguish between different angle values, a reference zero angle needs to be determined. For example, the reference zero angle may be the azimuth angle at which the simulated seismic source emits the simulated seismic signal directly north of the target detection point. The initial dividing angle is determined every preset angle dividing step length by taking the reference zero angle as a reference, and the initial dividing angle can be determinedThe initial dividing angle.
For example, if the north direction of the target detection point is taken as a reference zero angle, and the preset angle dividing step is 5 °, the initial dividing angles may be 0 °, 5 °, 10 °, …, 355 °, and 72 initial dividing angles are determined in total.
S120, dividing the seismic data into a plurality of azimuth angle gather data according to the initial dividing angles and preset dividing rules for each initial dividing angle.
The preset dividing rule is a preset rule for dividing the seismic data. For example, the preset division rule may divide the seismic data into a preset number of azimuth gather data, or may divide the seismic data into one azimuth gather data every preset angle with an initial division angle as a division start point.
Optionally, if the preset division rule is to divide the seismic data into a preset number of azimuth gather data. The specific implementation mode for dividing the seismic data into a plurality of azimuth gather data is as follows: determining a unit azimuth gather angle based on the flat angle value and a preset number; and for the seismic data, determining azimuth angle gather data every unit azimuth angle gather angle by taking the initial dividing angle as a dividing starting point to obtain a plurality of azimuth angle gather data.
Illustratively, the flat angle value is 180 °, the preset number is 6, and the unit azimuth gather angle is 180 ° -6=30°. If the initial dividing angle is 0 degrees, determining the seismic data of 0 degrees to 30 degrees azimuth angles, the seismic data of 30 degrees to 60 degrees azimuth angles, the seismic data of 60 degrees to 90 degrees azimuth angles, the seismic data of 90 degrees to 120 degrees azimuth angles, the seismic data of 120 degrees to 150 degrees azimuth angles and the seismic data of 150 degrees to 180 degrees azimuth angles as azimuth angle gather data respectively. If the initial dividing angle is 20 degrees, determining 20-50-degree azimuth seismic data, 50-80-degree azimuth seismic data, 80-110-degree azimuth seismic data, 110-140-degree azimuth seismic data, 140-170-degree azimuth seismic data and 170-200-degree azimuth seismic data as azimuth gather data respectively.
S130, determining an azimuth angle gather dividing scheme to be evaluated for a plurality of azimuth angle gather data corresponding to each initial dividing angle.
In this embodiment, for each initial dividing angle, a dividing manner of a plurality of azimuth gather data corresponding to one initial dividing angle is determined as a scheme for dividing azimuth gathers to be evaluated. Based on this, there are how many initial dividing angles, corresponding to how many azimuth gather dividing schemes to be evaluated.
For example, an initial dividing angle is 0 °, and an azimuth angle gather data is divided every 30 ° to be an azimuth angle gather dividing scheme to be evaluated; and dividing azimuth angle gather data at intervals of 30 degrees by using an initial dividing angle of 20 degrees as another azimuth angle gather dividing scheme to be evaluated.
S140, for each azimuth gather dividing scheme to be evaluated, seismic characteristic data corresponding to each azimuth gather data are obtained, and azimuth anisotropy crack prediction accuracy assessment factors are determined based on the seismic characteristic data.
Wherein the seismic profile data comprises: total number of seismic traces, number of non-anisotropic seismic traces, maximum incidence angle, and maximum incidence angle without anisotropy. Illustratively, taking azimuth gather data of azimuth angles from 0 ° to 30 °, the total number of seismic traces can be understood as how many source shots are co-recorded in the seismic data of azimuth angles from 0 ° to 30 °. The number of non-anisotropic seismic traces is a reference that can be directly calculated based on the prior art. The seismic data of azimuth angles of 0-30 degrees comprise the incidence angle of each seismic channel, and the maximum incidence angle of each seismic channel is determined to be the maximum incidence angle. The maximum incidence angle without anisotropy is also a reference that can be directly calculated based on the prior art.
The azimuth anisotropy crack prediction precision evaluation factor can be understood as a division scheme evaluation index, and is characterized in a numerical form. The larger the numerical value of the azimuth anisotropy fracture prediction precision evaluation factor is, the more ideal the rationality of the azimuth gather division scheme to be evaluated applied to azimuth anisotropy fracture prediction is represented.
Specifically, an azimuth angle gather division scheme to be evaluated is taken as an example for illustration, azimuth angle gather data of a plurality of azimuth angle ranges corresponding to the azimuth angle gather division scheme to be evaluated is obtained, and the total number of seismic traces, the number of non-anisotropic seismic traces, the maximum incident angle and the non-anisotropic maximum incident angle corresponding to each azimuth angle gather data are obtained. For each azimuth angle gather data, determining an effective channel duty ratio according to the total number of the seismic channels and the number of the non-anisotropic seismic channels; determining an effective incident angle ratio from the maximum incident angle and the anisotropic-free maximum incident angle; thereby determining an azimuthal anisotropy fracture prediction accuracy assessment factor based on the effective track occupancy, the effective incident angle occupancy, and the total number of azimuthal gather data.
S150, determining a target azimuth angle gather dividing scheme based on azimuth anisotropy crack prediction precision evaluation factors corresponding to each azimuth angle gather dividing scheme to be evaluated.
In particular, the number of target azimuth gather partitioning schemes may be one or more. If a target azimuth angle gather dividing scheme is determined from the azimuth angle gather dividing schemes to be evaluated, the maximum value of the azimuth anisotropy crack prediction precision evaluation factor corresponding to each azimuth angle gather dividing scheme to be evaluated can be determined, and the azimuth angle gather dividing scheme to be evaluated with the maximum value is determined as the target azimuth angle gather dividing scheme. If multiple target azimuth angle gather partitioning schemes are determined from the azimuth angle gather partitioning schemes to be evaluated, the target azimuth angle gather partitioning schemes can be screened from the multiple azimuth angle gather partitioning schemes to be evaluated according to a preset evaluation factor threshold. For example, an azimuth gather partitioning scheme to be evaluated, in which the azimuth anisotropy crack prediction accuracy evaluation factor is greater than a preset evaluation factor threshold, is determined as a target azimuth gather partitioning scheme.
According to the technical scheme provided by the embodiment of the invention, firstly, the seismic data is acquired, a plurality of initial dividing angles corresponding to the seismic data are determined according to the preset angle dividing step length, then, the seismic data are divided into a plurality of azimuth angle gather data according to the initial dividing angles and the preset dividing rules for each initial dividing angle, and further, the azimuth angle gather dividing scheme to be evaluated is determined for the plurality of azimuth angle gather data corresponding to each initial dividing angle. And for each azimuth angle gather division scheme to be evaluated, acquiring seismic characteristic data corresponding to each azimuth angle gather data, determining azimuth anisotropy crack prediction precision evaluation factors based on the seismic characteristic data, and finally determining a target azimuth angle gather division scheme based on the azimuth anisotropy crack prediction precision evaluation factors corresponding to each azimuth angle gather division scheme to be evaluated. The invention provides a quantitative evaluation method for an azimuth angle gather dividing scheme of seismic data, which realizes the rapid and efficient determination of the optimal azimuth angle gather dividing scheme of the seismic data, thereby improving the accuracy of crack prediction.
Example two
Fig. 2 is a flowchart of a method for determining an azimuth gather partitioning scheme according to a second embodiment of the present invention, where the step of determining an azimuth anisotropic fracture prediction accuracy assessment factor based on seismic feature data is further refined on the basis of the foregoing embodiments, and the embodiments of the present invention may be combined with each of the alternatives in one or more embodiments. As shown in fig. 2, the method for determining the azimuth gather partitioning scheme includes the following steps:
s210, determining a plurality of initial dividing angles corresponding to the seismic data according to the preset angle dividing step length.
S220, dividing the seismic data into a plurality of azimuth angle gather data according to the initial dividing angles and preset dividing rules for each initial dividing angle.
S230, determining an azimuth angle gather dividing scheme to be evaluated for the plurality of azimuth angle gather data corresponding to each initial dividing angle.
S240, acquiring seismic characteristic data corresponding to the azimuth angle gather data for each azimuth angle gather division scheme to be evaluated.
S250, determining the effective channel duty ratio based on the total number of the seismic channels and the number of the non-anisotropic seismic channels.
In this embodiment, since this embodiment is applied to a scenario of azimuth anisotropic fracture prediction, the greater the number of anisotropic seismic traces, the higher the accuracy of fracture prediction. The total number of seismic traces is the sum of the number of non-anisotropic seismic traces and the number of anisotropic seismic traces. The non-anisotropic seismic traces may be determined based on seismic modeling data, while the anisotropic seismic traces are difficult to determine. Thus, the number of anisotropic traces can be determined as long as the number of non-anisotropic traces and the total number of traces are determined. Further, the ratio of the number of anisotropic traces to the total number of traces is determined as the effective trace duty cycle.
S260, determining the effective incidence angle ratio based on the maximum incidence angle and the maximum incidence angle without anisotropy.
In this embodiment, first, an incident angle difference between the maximum incident angle and the maximum incident angle without anisotropy is determined, and a ratio of the incident angle difference to the maximum incident angle is determined as an effective incident angle ratio.
S270, determining an azimuth angle gather coefficient based on the total number of the azimuth angle gather data.
In this embodiment, the total number of azimuth gather data is directly available. If the total number of azimuth gather data is N AG The azimuth gather coefficients can be expressed as
S280, determining an azimuth anisotropic fracture prediction precision evaluation factor based on the effective track duty ratio, the effective incidence angle duty ratio and the azimuth angle gather coefficient.
Specifically, the target expression is adopted to quantify the estimation factor of the accuracy of predicting the anisotropic fracture of the azimuth, and the target expression is:
wherein a is f For the estimation factor of the azimuth anisotropy crack prediction precision, N AG For the total number of azimuth gather data,for azimuth gather coefficients, θ max (i) For the maximum angle of incidence, θ, of the ith azimuth gather data ISO For maximum incidence angle without anisotropy, F (i) is the total number of seismic traces for the ith azimuth gather data, F ISO (i) Number of anisotropic free seismic traces for the ith azimuth gather data, +.>For the effective incident angle ratio, +.>Is the effective track duty cycle.
S290, determining a maximum evaluation factor value from azimuth anisotropic fracture prediction precision evaluation factors corresponding to all azimuth gather division schemes to be evaluated, and determining the azimuth gather division scheme to be evaluated corresponding to the maximum evaluation factor value as a target azimuth gather division scheme.
In this embodiment, each azimuth gather partitioning scheme to be evaluated corresponds to an azimuth anisotropic fracture prediction accuracy evaluation factor, and since the azimuth anisotropic fracture prediction accuracy evaluation factor is a specific quantization index, the target azimuth gather partitioning scheme can be determined from a plurality of azimuth gather partitioning schemes to be evaluated based on the magnitude of the azimuth anisotropic fracture prediction accuracy evaluation factor value. For example, the azimuth anisotropy crack prediction accuracy evaluation factor corresponding to the azimuth gather to be evaluated partition scheme a is 0.43, the azimuth anisotropy crack prediction accuracy evaluation factor corresponding to the azimuth gather to be evaluated partition scheme B is 0.45, and the azimuth anisotropy crack prediction accuracy evaluation factor corresponding to the azimuth gather to be evaluated partition scheme C is 0.46, and then the azimuth gather to be evaluated partition scheme C is determined as the target azimuth gather partition scheme.
According to the technical scheme provided by the embodiment of the invention, when the azimuth anisotropy crack prediction precision evaluation factor is determined based on the seismic characteristic data, the effective channel duty ratio is determined based on the total number of the seismic channels and the number of the non-anisotropy seismic channels, wherein the effective channel duty ratio is the ratio of the number of the anisotropic seismic channels to the total number of the seismic channels; determining an effective incident angle ratio based on the maximum incident angle and the anisotropic-free maximum incident angle, wherein the effective incident angle ratio is the ratio of the anisotropic maximum incident angle to the maximum incident angle; determining an azimuth gather coefficient based on the total number of azimuth gather data; and finally, determining an azimuth anisotropic fracture prediction precision evaluation factor based on the effective track duty ratio, the effective incidence angle duty ratio and the azimuth angle gather coefficient. According to the quantitative evaluation method for the seismic data azimuth angle gather dividing scheme, which is provided by the embodiment, the data characteristics of the total number of azimuth angle gather data, the number of anisotropic seismic channels and the maximum anisotropic incidence angle are comprehensively considered, and the optimal seismic data azimuth angle gather dividing scheme is further determined rapidly and efficiently, so that the accuracy of crack prediction is improved.
Example III
Fig. 3 is a schematic structural diagram of a determining device for an azimuth gather partitioning scheme according to a third embodiment of the present invention, where the determining device can execute the determining method for the azimuth gather partitioning scheme according to the embodiment of the present invention. The device comprises: an initial angle determination module 310, an angle gather partitioning module 320, a partitioning scheme determination module 330, an assessment factor determination module 340, a target scheme determination module 350.
The initial angle determining module 310 is configured to determine a plurality of initial division angles corresponding to the seismic data according to a preset angle division step;
the angle gather dividing module 320 is configured to divide the seismic data into a plurality of azimuth angle gather data for each initial dividing angle according to the initial dividing angle and a preset dividing rule;
the division scheme determining module 330 is configured to determine an azimuth angle gather division scheme to be evaluated for the plurality of azimuth angle gather data corresponding to each initial division angle;
the evaluation factor determining module 340 is configured to obtain, for each azimuth gather to be evaluated partition scheme, seismic feature data corresponding to each azimuth gather data, and determine an azimuth anisotropic fracture prediction accuracy evaluation factor based on the seismic feature data;
the target scheme determining module 350 is configured to determine a target azimuth gather partitioning scheme based on azimuth anisotropic fracture prediction accuracy evaluation factors corresponding to the azimuth gather partitioning schemes to be evaluated.
According to the technical scheme provided by the embodiment of the invention, firstly, the seismic data is acquired, a plurality of initial dividing angles corresponding to the seismic data are determined according to the preset angle dividing step length, then, the seismic data are divided into a plurality of azimuth angle gather data according to the initial dividing angles and the preset dividing rules for each initial dividing angle, and further, the azimuth angle gather dividing scheme to be evaluated is determined for the plurality of azimuth angle gather data corresponding to each initial dividing angle. And for each azimuth angle gather division scheme to be evaluated, acquiring seismic characteristic data corresponding to each azimuth angle gather data, determining azimuth anisotropy crack prediction precision evaluation factors based on the seismic characteristic data, and finally determining a target azimuth angle gather division scheme based on the azimuth anisotropy crack prediction precision evaluation factors corresponding to each azimuth angle gather division scheme to be evaluated. The invention provides a quantitative evaluation method for an azimuth angle gather dividing scheme of seismic data, which realizes the rapid and efficient determination of the optimal azimuth angle gather dividing scheme of the seismic data, thereby improving the accuracy of crack prediction.
Based on the above aspects, the initial angle determining module 310 includes:
a reference angle determining unit for determining a reference zero angle from the seismic data;
the initial angle determining unit is used for determining an initial dividing angle every preset angle dividing step length by taking the reference zero angle as a dividing starting point, and obtaining a plurality of initial dividing angles corresponding to the seismic data.
On the basis of the technical schemes, the angle gather dividing module comprises:
the unit angle determining unit is used for determining a unit azimuth angle gather angle based on the flat angle value and the preset quantity;
the angle gather dividing unit is used for determining azimuth angle gather data of the seismic data at intervals of unit azimuth angle gather angles by taking an initial dividing angle as a dividing starting point to obtain a plurality of azimuth angle gather data.
Based on the above aspects, the evaluation factor determining module 340 includes:
an effective channel duty ratio determining unit for determining an effective channel duty ratio based on the total number of seismic channels and the number of non-anisotropic seismic channels;
an incidence angle duty ratio determining unit for determining an effective incidence angle duty ratio based on the maximum incidence angle and the maximum incidence angle without anisotropy;
a coefficient determination unit for determining an azimuth gather coefficient based on a total number of azimuth gather data;
and the evaluation factor determining unit is used for determining an azimuth anisotropic fracture prediction precision evaluation factor based on the effective track duty ratio, the effective incidence angle duty ratio and the azimuth angle gather coefficient.
In this embodiment, the target expression is used to quantify the estimation factor of the accuracy of predicting the anisotropic fracture of the azimuth, where the target expression is:
wherein a is f For the estimation factor of the azimuth anisotropy crack prediction precision, N AG For the total number of azimuth gather data,for azimuth gather coefficients, θ max (i) For the maximum angle of incidence, θ, of the ith azimuth gather data ISO For maximum incidence angle without anisotropy, F (i) is the total number of seismic traces for the ith azimuth gather data, F ISO (i) Number of anisotropic free seismic traces for the ith azimuth gather data, +.>For the effective incident angle ratio, +.>Is the effective track duty cycle.
Based on the above technical solutions, the target solution determining module 350 includes:
the maximum value determining unit is used for determining the maximum value of the evaluation factors from the azimuth anisotropic fracture prediction precision evaluation factors corresponding to each azimuth gather dividing scheme to be evaluated;
and the target scheme determining unit is used for determining the azimuth angle gather dividing scheme to be evaluated corresponding to the maximum value of the evaluation factor as the target azimuth angle gather dividing scheme.
The determining device for the azimuth angle gather dividing scheme provided by the embodiment of the disclosure can execute the determining method for the azimuth angle gather dividing scheme provided by any embodiment of the disclosure, and has the corresponding functional modules and beneficial effects of the executing method.
It should be noted that each unit and module included in the above apparatus are only divided according to the functional logic, but not limited to the above division, so long as the corresponding functions can be implemented; in addition, the specific names of the functional units are also only for convenience of distinguishing from each other, and are not used to limit the protection scope of the embodiments of the present disclosure.
Example IV
Fig. 4 is a schematic structural diagram of an electronic device according to a fourth embodiment of the present invention. The electronic device 10 is intended to represent various forms of digital computers, such as laptops, desktops, workstations, personal digital assistants, servers, blade servers, mainframes, and other appropriate computers. The electronic device may also represent various forms of mobile devices, such as personal digital processing, cellular telephones, smartphones, wearable electronic devices (e.g., helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions, are meant to be exemplary only, and are not meant to limit implementations of the inventions described and/or claimed herein.
As shown in fig. 4, the electronic device 10 includes at least one processor 11, and a memory, such as a Read Only Memory (ROM) 12, a Random Access Memory (RAM) 13, etc., communicatively connected to the at least one processor 11, in which the memory stores a computer program executable by the at least one processor, and the processor 11 may perform various appropriate actions and processes according to the computer program stored in the Read Only Memory (ROM) 12 or the computer program loaded from the storage unit 18 into the Random Access Memory (RAM) 13. In the RAM13, various programs and data required for the operation of the electronic device 10 may also be stored. The processor 11, the ROM12 and the RAM13 are connected to each other via a bus 13. An input/output (I/O) interface 15 is also connected to bus 13.
Various components in the electronic device 10 are connected to the I/O interface 15, including: an input unit 16 such as a keyboard, a mouse, etc.; an output unit 17 such as various types of displays, speakers, and the like; a storage unit 18 such as a magnetic disk, an optical disk, or the like; and a communication unit 19 such as a network card, modem, wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information/data with other electronic devices via a computer network, such as the internet, and/or various telecommunication networks.
The processor 11 may be a variety of general and/or special purpose processing components having processing and computing capabilities. Some examples of processor 11 include, but are not limited to, a Central Processing Unit (CPU), a Graphics Processing Unit (GPU), various specialized Artificial Intelligence (AI) computing chips, various processors running machine learning model algorithms, digital Signal Processors (DSPs), and any suitable processor, controller, microcontroller, etc. The processor 11 performs the various methods and processes described above, such as the determination of the azimuth gather partitioning scheme.
In some embodiments, the method of determining the azimuth gather partitioning scheme may be implemented as a computer program tangibly embodied on a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program may be loaded and/or installed onto the electronic device 10 via the ROM12 and/or the communication unit 19. When the computer program is loaded into RAM13 and executed by processor 11, one or more steps of the above-described method of determining the azimuth gather partitioning scheme may be performed. Alternatively, in other embodiments, the processor 11 may be configured to perform the method of determining the azimuth gather partitioning scheme in any other suitable way (e.g., by means of firmware).
Various implementations of the systems and techniques described here above may be implemented in digital electronic circuitry, integrated circuit systems, field Programmable Gate Arrays (FPGAs), application Specific Integrated Circuits (ASICs), application Specific Standard Products (ASSPs), systems On Chip (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and/or combinations thereof. These various embodiments may include: implemented in one or more computer programs, the one or more computer programs may be executed and/or interpreted on a programmable system including at least one programmable processor, which may be a special purpose or general-purpose programmable processor, that may receive data and instructions from, and transmit data and instructions to, a storage system, at least one input device, and at least one output device.
A computer program for carrying out methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general purpose computer, special purpose computer, or other programmable apparatus for determining an azimuth gather partitioning scheme, such that the computer programs, when executed by the processor, cause the functions/operations specified in the flowchart and/or block diagram to be implemented. The computer program may execute entirely on the machine, partly on the machine, as a stand-alone software package, partly on the machine and partly on a remote machine or entirely on the remote machine or server.
In the context of the present invention, a computer-readable storage medium may be a tangible medium that can contain, or store a computer program for use by or in connection with an instruction execution system, apparatus, or device. The computer readable storage medium may include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or electronic device, or any suitable combination of the foregoing. Alternatively, the computer readable storage medium may be a machine readable signal medium. More specific examples of a machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage electronic device, a magnetic storage electronic device, or any suitable combination of the foregoing.
To provide for interaction with a user, the systems and techniques described here can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to a user; and a keyboard and a pointing device (e.g., a mouse or a trackball) through which a user can provide input to the electronic device. Other kinds of devices may also be used to provide for interaction with a user; for example, feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form, including acoustic input, speech input, or tactile input.
The systems and techniques described here can be implemented in a computing system that includes a background component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front-end component (e.g., a user computer having a graphical user interface or a web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such background, middleware, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: local Area Networks (LANs), wide Area Networks (WANs), blockchain networks, and the internet.
The computing system may include clients and servers. The client and server are typically remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other. The server can be a cloud server, also called a cloud computing server or a cloud host, and is a host product in a cloud computing service system, so that the defects of high management difficulty and weak service expansibility in the traditional physical hosts and VPS service are overcome. It should be appreciated that various forms of the flows shown above may be used to reorder, add, or delete steps. For example, the steps described in the present invention may be performed in parallel, sequentially, or in a different order, so long as the desired results of the technical solution of the present invention are achieved, and the present invention is not limited herein. The above embodiments do not limit the scope of the present invention. It will be apparent to those skilled in the art that various modifications, combinations, sub-combinations and alternatives are possible, depending on design requirements and other factors. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of the present invention.
Claims (7)
1. A method for determining an azimuth gather partitioning scheme, comprising:
determining a plurality of initial dividing angles corresponding to the seismic data according to a preset angle dividing step length;
dividing the seismic data into a plurality of azimuth angle gather data according to the initial dividing angles and a preset dividing rule for each initial dividing angle;
determining an azimuth angle gather dividing scheme to be evaluated for the azimuth angle gather data corresponding to each initial dividing angle;
for each azimuth gather dividing scheme to be evaluated, acquiring seismic characteristic data corresponding to each azimuth gather data, and determining azimuth anisotropic fracture prediction precision evaluation factors based on the seismic characteristic data;
and determining a target azimuth angle gather dividing scheme based on azimuth anisotropic crack prediction precision evaluation factors corresponding to the azimuth angle gather dividing schemes to be evaluated.
2. The method of claim 1, wherein determining a plurality of initial division angles for the seismic data according to the predetermined angle division step comprises:
determining a reference zero angle from the seismic data;
and taking the reference zero angle as a dividing starting point, dividing step length at intervals of preset angles, and determining an initial dividing angle to obtain a plurality of initial dividing angles corresponding to the seismic data.
3. The method of claim 1, wherein the predetermined division rule is to divide the seismic data into a predetermined number of azimuth gather data, and wherein the dividing the seismic data into a plurality of azimuth gather data according to the initial division angle and the predetermined division rule comprises:
determining a unit azimuth gather angle based on the flat angle value and a preset number;
and determining azimuth angle gather data of the seismic data at intervals of unit azimuth angle gather angles by taking the initial dividing angle as a dividing starting point to obtain a plurality of azimuth angle gather data.
4. The method of claim 1, wherein the seismic profile data comprises: total number of seismic traces, number of non-anisotropic seismic traces, maximum incidence angle, and maximum incidence angle without anisotropy.
5. The method of claim 4, wherein the determining an azimuthal anisotropic fracture prediction accuracy assessment factor based on the seismic feature data comprises:
determining an effective trace duty cycle based on the total number of seismic traces and the number of non-anisotropic seismic traces;
determining an effective incident angle duty cycle based on the maximum incident angle and the anisotropic-free maximum incident angle;
determining an azimuth gather coefficient based on the total number of azimuth gather data;
an azimuthal anisotropy fracture prediction accuracy assessment factor is determined based on the effective trace duty cycle, the effective incident angle duty cycle, and the azimuthal trace set coefficient.
6. The method of claim 5, wherein the azimuthal anisotropy fracture prediction accuracy assessment factor is quantified using a target expression, the target expression being:
wherein a is f For the estimation factor of the azimuth anisotropy crack prediction precision, N AG For the total number of azimuth gather data,for azimuth gather coefficients, θ max (i) For the maximum angle of incidence, θ, of the ith azimuth gather data ISO For maximum incidence angle without anisotropy, F (i) is the total number of seismic traces for the ith azimuth gather data, F ISO (i) Number of anisotropic free seismic traces for the ith azimuth gather data, +.>For the effective incident angle ratio, +.>Is the effective track duty cycle.
7. The method of claim 1, wherein determining the target azimuth gather partitioning scheme based on azimuth anisotropic fracture prediction accuracy assessment factors corresponding to each of the azimuth gather partitioning schemes to be evaluated comprises:
determining the maximum value of the evaluation factors from the azimuth anisotropic fracture prediction precision evaluation factors corresponding to each azimuth gather dividing scheme to be evaluated;
and determining an azimuth angle gather dividing scheme to be evaluated corresponding to the maximum value of the evaluation factor as a target azimuth angle gather dividing scheme.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310930152.0A CN116931083A (en) | 2023-07-27 | 2023-07-27 | Determination method of azimuth angle gather division scheme |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310930152.0A CN116931083A (en) | 2023-07-27 | 2023-07-27 | Determination method of azimuth angle gather division scheme |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116931083A true CN116931083A (en) | 2023-10-24 |
Family
ID=88389437
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310930152.0A Pending CN116931083A (en) | 2023-07-27 | 2023-07-27 | Determination method of azimuth angle gather division scheme |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116931083A (en) |
-
2023
- 2023-07-27 CN CN202310930152.0A patent/CN116931083A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113705628B (en) | Determination method and device of pre-training model, electronic equipment and storage medium | |
| CN113642710B (en) | Quantification method, device, equipment and storage medium of network model | |
| CN114861039A (en) | Parameter configuration method, device, equipment and storage medium of search engine | |
| CN117593115A (en) | Feature value determining method, device, equipment and medium of credit risk assessment model | |
| CN116931083A (en) | Determination method of azimuth angle gather division scheme | |
| CN115438007A (en) | File merging method and device, electronic equipment and medium | |
| CN114999665A (en) | Data processing method and device, electronic equipment and storage medium | |
| CN114692987A (en) | Time sequence data analysis method, device, equipment and storage medium | |
| CN112560987A (en) | Image sample processing method, device, equipment, storage medium and program product | |
| CN114037057B (en) | Pre-training model generation method and device, electronic equipment and storage medium | |
| CN114866437B (en) | Node detection method, device, equipment and medium | |
| CN114037058B (en) | Pre-training model generation method and device, electronic equipment and storage medium | |
| CN117874559A (en) | Determination method, device, equipment and storage medium for lithofacies classification quantity | |
| CN116610453A (en) | Task allocation method and device, electronic equipment and storage medium | |
| CN117235500A (en) | Material category identification method and device, electronic equipment and storage medium | |
| CN117972469A (en) | Igneous rock classification and identification method, device, equipment and storage medium | |
| CN115658826A (en) | Method, device and equipment for determining track stop point and storage medium | |
| CN116361460A (en) | Data integration method and device, storage medium, electronic equipment and product | |
| CN115455019A (en) | Search intention identification method, device and equipment based on user behavior analysis | |
| CN115421696A (en) | Problem code determination method and device | |
| CN116205321A (en) | Method, device, equipment and storage medium for determining carbon consumption | |
| CN116933968A (en) | Method and device for measuring saturation of oil-containing gas, electronic equipment and medium | |
| CN116167519A (en) | Monitoring amount prediction method, device, equipment and medium | |
| CN116631318A (en) | Display panel and debugging method and device thereof | |
| CN115017145A (en) | Data expansion method, device and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |