CN113156494B - Method and device for calculating equivalent coverage times of surface elements - Google Patents
Method and device for calculating equivalent coverage times of surface elements Download PDFInfo
- Publication number
- CN113156494B CN113156494B CN202010013587.5A CN202010013587A CN113156494B CN 113156494 B CN113156494 B CN 113156494B CN 202010013587 A CN202010013587 A CN 202010013587A CN 113156494 B CN113156494 B CN 113156494B
- Authority
- CN
- China
- Prior art keywords
- energy
- target layer
- excitation
- calculating
- seismic
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000005284 excitation Effects 0.000 claims abstract description 172
- 238000012360 testing method Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000013461 design Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 7
- 239000002360 explosive Substances 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000002699 waste material Substances 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
-
- 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 provides a method and a device for calculating the equivalent coverage times of a surface element, wherein the method comprises the following steps: acquiring a first excitation parameter, and acquiring a first seismic record of the standard gun according to the first excitation parameter; the excitation parameters are determined according to the excitation mode; calculating first target layer energy of the standard cannon according to the first seismic record; acquiring second seismic records of a plurality of second excitation parameters of the standard cannon, and calculating second target layer energy of the standard cannon according to the second seismic records; the energy of the second target layer and the energy of the first target layer are obtained according to the same time window data and offset data; determining a curve of the energy ratio changing along with the excitation parameters according to the second excitation parameters, the first target layer energy and the second target layer energy; and calculating the equivalent coverage times according to the curve. The invention can optimize the design scheme of the observation system and ensure the quality of the seismic data by obtaining more accurate effective coverage times of the surface elements.
Description
Technical Field
The invention relates to the technical field of seismic data acquisition, in particular to a method and a device for calculating equivalent coverage times of a surface element.
Background
Because of dense barriers, the shot points are difficult to be uniformly distributed in the complicated ground surface area. When the seismic data acquisition is carried out in the field, in order to acquire the whole data, the safety distance limit of the obstacle is avoided, and various excitation modes and excitation parameters are adopted, so that the energy difference of adjacent single shots is larger. At present, three-dimensional acquisition design software is used for counting offset pairs in CMP (chemical mechanical polishing) surface elements to calculate the coverage times, and the influence of initial excitation energy on actual reflected waves is not considered, so that the situation that the CMP surface elements are relatively uniform in coverage times and the transverse energy of actual data is very large (influence of underground structures is not considered), the quality of seismic data cannot be ensured, and further the later seismic attribute analysis and research are influenced.
Disclosure of Invention
The invention provides a method and a device for calculating the equivalent coverage times of a surface element, which can obtain more accurate equivalent coverage times of the surface element, and optimize the design scheme of an observation system and ensure the quality of seismic data by obtaining the effective coverage times of the surface element.
In a first aspect, an embodiment of the present invention provides a method for calculating a number of coverage times of a bin, where the method includes: acquiring a first excitation parameter, and acquiring a first seismic record of a standard gun according to the first excitation parameter; the excitation parameters are determined according to an excitation mode; calculating first target layer energy of the standard cannon according to the first seismic record; acquiring second seismic records of a plurality of second excitation parameters of the standard cannon, and calculating second target layer energy of the standard cannon according to the second seismic records; the second target layer energy and the first target layer energy are obtained according to the same time window data and offset data; determining a curve of the energy ratio changing along with the excitation parameters according to the second excitation parameters, the first target layer energy and the second target layer energy; calculating equivalent coverage times according to the curve; calculating a first destination layer energy of the standard shot from the first seismic record, comprising:
determining target time window data and target offset data at a target layer according to the first seismic record;
calculating first target layer energy of the standard cannon according to the target time window data and the target offset data;
calculating the energy of a first target layer of the standard gun by using the following formula:
wherein A is t For the energy of the target layer, N is the number of time window samples, and A is the amplitude of the seismic record.
In a second aspect, an embodiment of the present invention further provides a device for calculating a number of coverage equivalent to a bin, where the device includes: the acquisition module is used for acquiring a first excitation parameter and acquiring a first seismic record of the standard gun according to the first excitation parameter; the excitation parameters are determined according to an excitation mode; the excitation parameters are determined according to an excitation mode; the first energy module is used for calculating first target layer energy of the standard cannon according to the first seismic record; the second energy module is used for acquiring second seismic records of a plurality of second excitation parameters of the standard cannon and calculating second target layer energy of the standard cannon according to the second seismic records; the second target layer energy and the first target layer energy are obtained according to the same time window data and offset data; the fitting module is used for determining a curve of the energy ratio changing along with the excitation parameters according to the second excitation parameters, the first target layer energy and the second target layer energy; the calculation module is used for calculating the equivalent coverage times according to the curve;
the first energy module is specifically used for:
determining target time window data and target offset data at a target layer according to the first seismic record;
calculating first target layer energy of the standard cannon according to the target time window data and the target offset data;
the first energy module is specifically used for: calculating the energy of a first target layer of the standard gun by using the following formula:
wherein A is t For the purpose of layer energy, N isThe number of window samples, A, is the amplitude of the seismic record.
In a third aspect, an embodiment of the present invention further provides a computer device, including a memory, and a processor, where the memory stores a computer program that can run on the processor, and when the processor executes the computer program, the computing method for the number of equivalent coverage times of the bin is implemented.
In a fourth aspect, embodiments of the present invention also provide a computer readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the above-described method of calculating the number of equivalent coverage of a bin.
The embodiment of the invention has the following beneficial effects: the embodiment of the invention provides a calculation scheme of the equivalent coverage times of a surface element, which comprises the steps of obtaining a first excitation parameter which is determined in advance according to an excitation mode, obtaining a first seismic record of a standard gun corresponding to the first excitation parameter, and calculating the first target layer energy of the standard gun in a target layer according to the information of the first seismic record; then, second seismic records of the standard cannon at a plurality of second excitation parameters are obtained, and under the condition that time window data and offset data are the same, second target layer energy corresponding to each second excitation parameter is calculated; and finally, determining a curve of the energy ratio changing along with the excitation parameters according to the second excitation parameters, the first target layer energy and the second target layer energy, and calculating the equivalent coverage times under different excitation parameters according to the curve. According to the embodiment of the invention, the earthquake record is obtained according to the excitation parameters, the energy of the target layer under different excitation parameters is calculated according to the earthquake record, and a relation curve between the energy ratio and the excitation parameters is fitted so as to calculate the equivalent coverage times based on the actual energy contribution condition of a single gun. The embodiment of the invention can optimize the design scheme of the observation system by obtaining more accurate effective coverage times of the surface elements and ensure the quality of the seismic data.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
In order to make the above objects, features and advantages of the present invention more comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the description of the embodiments or the prior art will be briefly described, and it is obvious that the drawings in the description below are some embodiments of the present invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of a method for calculating the number of equivalent coverage times of a bin according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a calculation flow of the number of equivalent coverage of a bin according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a fitted curve of the energy ratio of a single gun to a standard gun with different excitation parameters according to the variation of the excitation parameters provided by the embodiment of the invention;
FIG. 4 is a block diagram of a computing device for the number of equivalent coverage of a bin according to an embodiment of the present invention;
FIG. 5 is a block diagram of another embodiment of a device for calculating the number of equivalent coverage of a bin;
fig. 6 is a block diagram of a computer device according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
At present, when seismic data are acquired, the problem that the number of times of CMP surface element coverage obtained by the existing method for calculating the number of times of coverage is relatively uniform and has great difference with the energy of a data body is caused because single shot energy imbalance and large difference of the number of times of coverage are not considered in the design stage, and serious offset noise is generated during offset. Even if the amplitude-based equalization processing technique is adopted, the offset arcing problem cannot be completely eliminated.
Based on the above, the method and the device for calculating the equivalent coverage times of the surface element provided by the embodiment of the invention can calculate the effective energy ranges of the reflected waves of different excitation modes and excitation parameters according to the actual conditions of the work area, calculate the coverage times based on the actual contribution degree, obtain the effective coverage times of the surface element, optimize the design scheme of an observation system, ensure the quality of data and reduce the production cost.
For the convenience of understanding the present embodiment, a method for calculating the number of equivalent coverage of a bin disclosed in the present embodiment is first described in detail.
The embodiment of the invention provides a method for calculating the equivalent coverage times of a surface element, which is shown in a flow chart of the method for calculating the equivalent coverage times of the surface element in fig. 1, and comprises the following steps:
step S102, obtaining a first excitation parameter, and obtaining a first seismic record of the standard shot according to the first excitation parameter.
In the embodiment of the invention, a work area may need to adopt different excitation sources due to the limitation of the surface condition, and the first excitation parameter is determined according to the excitation mode. For example, when an explosive is used as the excitation source, the first excitation parameter may be the amount of excitation; when a vibroseis is used as the excitation source, the first excitation parameter may be the number of shocks, or a combination of both. The value of the first excitation parameter may be predetermined. And taking the single cannon generating the first excitation parameter as a standard cannon, and acquiring a first seismic record generated by excitation of the standard cannon according to the value of the first excitation parameter.
Step S104, calculating the first target layer energy of the standard cannon according to the first seismic record.
In the embodiment of the invention, the first seismic record comprises time window data, offset data and other information, the time window data, offset data and other information of a target layer in the first seismic record are obtained, and the energy generated by a standard gun is calculated to obtain the energy of the first target layer.
It should be noted that, the target layer of exploration can be determined according to the requirement of the project geological task, and the corresponding time window and offset range can be determined on the single shot record.
Step S106, obtaining second seismic records of a plurality of second excitation parameters of the standard cannon, and calculating second target layer energy of the standard cannon according to the second seismic records.
In the embodiment of the invention, the second excitation parameter and the first excitation parameter are the same excitation parameter of a standard gun, but have different values. For example, when the excitation source is an explosive, if the first excitation parameter is 6kg, the second excitation parameter may be 3kg, 5kg, 8kg, 9.6kg, or the like. And a plurality of different second excitation parameters are required to be acquired, and second address records are acquired respectively according to the second excitation parameter values. And calculating the energy of the standard gun in the target layer according to the information in the second seismic record to obtain the energy of the second target layer. The values of the second objective layer energies obtained are different due to the different values of the second excitation parameters.
It should be noted that, since the first excitation parameter corresponds to the same excitation source as the second excitation parameter, the type of the second excitation parameter is identical to the type of the first excitation parameter, for example, for an explosive excitation source, if the first excitation parameter is a dose, then the second excitation parameter is also a dose, but each second excitation parameter is different from the value of the first excitation parameter.
In addition, when calculating the energy of the second target layer, it is necessary to select the same time window data and offset data as when calculating the energy of the first target layer.
And S108, determining a curve of the energy ratio changing along with the excitation parameters according to the second excitation parameters, the first target layer energy and the second target layer energy.
In the embodiment of the invention, a plurality of energy ratios can be obtained according to the energy of the first target layer and the energy of the second target layer, and the energy ratios can be obtained by fitting according to the corresponding relation between the energy ratios and the second excitation parameter values. A plot of energy ratio as a function of excitation parameters.
In the work area using the explosive as the main excitation source, the controllable source construction is needed in larger towns, and the excitation energy of the controllable source construction has a larger difference relative to the explosive. The same excitation source also needs to adopt different excitation parameters under different surface conditions, and the excitation energy of the excitation source also has great difference. Therefore, referring to the schematic diagram of the fitted curve of the energy ratio of the single cannon to the standard cannon with the variation of the excitation parameters, shown in fig. 3, a certain project in a certain area adopts explosive as a main excitation source, and a controllable source is adopted for excitation in a town area. It is necessary to fit curves of different energy ratios as a function of excitation parameters for different types of excitation parameters.
Step S110, calculating the equivalent coverage times according to the curve.
In the embodiment of the invention, the energy ratio of the standard cannon corresponding to different excitation parameters can be obtained according to the curve and used as the equivalent coverage times of the excitation parameters, so that the equivalent coverage times of the whole work area can be calculated. In the process of determining the equivalent coverage times, the influence of different excitation sources on energy when different excitation parameter values are taken into consideration, the determination process of the coverage times is optimized, the design scheme of an observation system is optimized, and the quality of seismic data is ensured.
In the embodiment of the invention, before the excitation is performed in the work area, the excitation parameter values corresponding to different equivalent coverage times can be obtained according to the curve, so that the numerical value determination process of the excitation parameter is optimized, the method can be used for guiding the excitation process, and the waste of production cost caused by multiple tests or unreasonable values of the excitation parameter is reduced.
The embodiment of the invention provides a calculation scheme of the equivalent coverage times of a surface element, which comprises the steps of obtaining a first excitation parameter which is determined in advance according to an excitation mode, obtaining a first seismic record of a standard gun corresponding to the first excitation parameter, and calculating the first target layer energy of the standard gun in a target layer according to the information of the first seismic record; then, second seismic records of the standard cannon at a plurality of second excitation parameters are obtained, and under the condition that time window data and offset data are the same, second target layer energy corresponding to each second excitation parameter is calculated; and finally, determining a curve of the energy ratio changing along with the excitation parameters according to the second excitation parameters, the first target layer energy and the second target layer energy, and calculating the equivalent coverage times under different excitation parameters according to the curve. According to the embodiment of the invention, the earthquake record is obtained according to the excitation parameters, the energy of the target layer under different excitation parameters is calculated according to the earthquake record, and a relation curve between the energy ratio and the excitation parameters is fitted so as to calculate the equivalent coverage times based on the actual energy contribution condition of a single gun. The embodiment of the invention can optimize the design scheme of the observation system by obtaining more accurate effective coverage times of the surface elements and ensure the quality of the seismic data.
In view of the fact that the time window data and offset data are to be determined in order to obtain a more accurate number of equivalent coverage, the following steps may be performed before the first excitation parameters are obtained:
acquiring a plurality of seismic data corresponding to a plurality of test excitation parameters; values of the first excitation parameters are determined from the plurality of seismic data among a plurality of test excitation parameters.
In the embodiment of the invention, the parameter values of the plurality of test excitation parameters are different. The test excitation parameters can be artificially set excitation parameters, or excitation parameters selected according to historical data, seismic data excited by the test excitation parameters are obtained, and the value of the first excitation parameter is determined in a plurality of test excitation parameters according to the quality of the seismic data. By determining the value of the first excitation parameter of better quality, a clearer target layer can be obtained, and the time window data and offset data meeting the requirements can be conveniently selected at the target layer.
The seismic record comprises time window data, offset data and the like, and the energy of the reflected wave is mainly influenced by the excitation parameters relative to a specific stratum, so that the energy of a target stratum in the time window and offset range determined by a standard gun can be calculated. The method can be specifically carried out according to the following steps:
determining target time window data and target offset data at a target layer according to the first seismic record; and calculating the first target layer energy of the standard cannon according to the target time window data and the target offset data.
In the embodiment of the invention, the target time window data and the target offset data are determined at the target layer of the first seismic record according to the requirements, and the range of the energy to be calculated in the seismic record is determined through the target time window data and the target offset data.
According to the determined time window of the target layer, the root mean square amplitude energy is calculated to obtain the energy of the first target layer, and the energy can be calculated by using the following formula:wherein A is t For the energy of the target layer, N is the number of time window samples, and A is the amplitude of the seismic record.
In the embodiment of the invention, the number N of time window samples is determined according to the time window data, and the corresponding range of the first target layer energy in the seismic record is determined according to N. Amplitude values of the seismic records may be obtained from the first seismic record. The target layer, the target time window data, the target offset data and the formula used in calculating the second target layer energy are identical to those used in calculating the first target layer energy.
The present solution is described below in connection with a specific embodiment. Referring to the schematic diagram of the calculation flow of the number of equivalent coverage of the bin shown in fig. 2, a basis is provided for implementing design of the observation system, so as to ensure quality of the collected data, and the method can be specifically implemented according to the following steps:
step S201, determining a typical earthquake single shot record of an exploration work area and excitation parameters thereof, and taking the record as a standard shot;
in seismic acquisition, the observation system has a great relationship with excitation sources in consideration of signal-to-noise ratio, energy, cost, efficiency and the like of a single shot, and the observation systems have great differences among different excitation sources in the same work area. The standard gun excitation parameters are the excitation parameters of the main excitation source of the current acquisition project. And obtaining the optimal excitation parameters of the project according to the early-stage test, namely the standard gun excitation parameters. Taking a certain three-dimensional in the Jizhong exploration area as an example, taking well cannon excitation as a main part, and performing early-stage test, wherein the optimal excitation dosage is 6kg, and then a single cannon record excited by 6kg of dosage is the standard cannon of the work area.
Step S202, determining a time window of a target layer on a standard gun seismic record and a gun offset range;
step S203, the energy of a target layer of the standard gun in the determined time window and the determined offset range is obtained;
step S204, the energy ratio of the energy of the single gun recording target layer with different excitation parameters in the determined time window and offset range to the energy of the standard gun is obtained;
step S205, fitting the energy ratio with the change of the excitation parameter to form a curve of the energy ratio with the change of the excitation parameter;
step S206, the energy ratio of different excitation parameters is obtained according to the curve to obtain the equivalent coverage times.
The embodiment of the invention provides a method and a device for calculating the equivalent coverage times of a surface element, which can calculate the effective energy ranges of reflection waves of different excitation modes and excitation parameters according to the actual conditions of a work area, calculate the coverage times based on the actual contribution degree, obtain the effective coverage times of the surface element, further optimize the design scheme of an observation system, ensure the quality of data and reduce the production cost. In the embodiment of the invention, the energy factor is considered in the method for calculating the coverage times, so that the field coverage times condition can be reflected more truly, the field offset design is guided, the data quality is improved, and the production cost is reduced.
The embodiment of the invention also provides a device for calculating the equivalent coverage times of the surface element, referring to a structural block diagram of the device for calculating the equivalent coverage times of the surface element shown in fig. 4, the device comprises:
an acquisition module 71, configured to acquire a first excitation parameter, and acquire a first seismic record of the standard shot according to the first excitation parameter; the excitation parameters are determined according to the excitation mode; the excitation parameters are determined according to the excitation mode;
a first energy module 72 for calculating a first target zone energy of the standard shot from the first seismic record; a second energy module 73, configured to obtain a second seismic record of a plurality of second excitation parameters of the standard shot, and calculate a second target layer energy of the standard shot according to the second seismic record; the energy of the second target layer and the energy of the first target layer are obtained according to the same time window data and offset data; a fitting module 74 for determining a curve of the energy ratio as a function of the excitation parameters from the second excitation parameters, the first energy of the target layer and the second energy of the target layer; a calculating module 75, configured to calculate the equivalent coverage number according to the curve.
In one embodiment, referring to the structural block diagram of the computing device for another number of equivalent coverage of a bin shown in fig. 5, the device further includes a determining module 76 for: acquiring a plurality of seismic data corresponding to a plurality of test excitation parameters; values of the first excitation parameters are determined from the plurality of seismic data among a plurality of test excitation parameters.
In one embodiment, the first energy module is specifically configured to: determining target time window data and target offset data at a target layer according to the first seismic record; and calculating the first target layer energy of the standard cannon according to the target time window data and the target offset data.
In one embodiment, the first energy module is specifically configured to: the energy of the first target layer of the standard cannon is calculated by using the following formula:wherein A is t For the energy of the target layer, N is the number of time window samples, and A is the amplitude of the seismic record.
The embodiment of the present invention further provides a computer device, referring to the schematic block diagram of the structure of the computer device shown in fig. 6, where the computer device includes a memory 81 and a processor 82, and the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, the processor implements the steps of any of the methods described above.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described computer device may refer to corresponding procedures in the foregoing method embodiments, which are not repeated here
Embodiments of the present invention also provide a computer readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the steps of any of the methods described above.
It will be appreciated by those skilled in the art that 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 and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations 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.
Finally, it should be noted that: the above examples are only specific embodiments of the present invention, and are not intended to limit the scope of the present invention, but it should be understood by those skilled in the art that the present invention is not limited thereto, and that the present invention is described in detail with reference to the foregoing examples: any person skilled in the art may modify or easily conceive of the technical solution described in the foregoing embodiments, or perform equivalent substitution of some of the technical features, while remaining within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention, and are intended to be included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (6)
1. The method for calculating the number of equivalent coverage times of the surface element is characterized by comprising the following steps of:
acquiring a first excitation parameter, and acquiring a first seismic record of a standard gun according to the first excitation parameter; the excitation parameters are determined according to an excitation mode;
calculating first target layer energy of the standard cannon according to the first seismic record;
acquiring second seismic records of a plurality of second excitation parameters of the standard cannon, and calculating second target layer energy of the standard cannon according to the second seismic records; the second target layer energy and the first target layer energy are obtained according to the same time window data and offset data;
determining a curve of the energy ratio changing along with the excitation parameters according to the second excitation parameters, the first target layer energy and the second target layer energy;
calculating equivalent coverage times according to the curve;
calculating a first destination layer energy of the standard shot from the first seismic record, comprising:
determining target time window data and target offset data at a target layer according to the first seismic record;
calculating first target layer energy of the standard cannon according to the target time window data and the target offset data;
calculating the energy of a first target layer of the standard gun by using the following formula:
wherein A is t For the energy of the target layer, N is the number of time window samples, and A is the amplitude of the seismic record.
2. The method of claim 1, further comprising, prior to acquiring the first excitation parameter:
acquiring a plurality of seismic data corresponding to a plurality of test excitation parameters;
determining a value of a first excitation parameter from the plurality of test excitation parameters based on the plurality of seismic data.
3. A computing device for a number of equivalent coverage of a bin, comprising:
the acquisition module is used for acquiring a first excitation parameter and acquiring a first seismic record of the standard gun according to the first excitation parameter; the excitation parameters are determined according to an excitation mode;
the first energy module is used for calculating first target layer energy of the standard cannon according to the first seismic record;
the second energy module is used for acquiring second seismic records of a plurality of second excitation parameters of the standard cannon and calculating second target layer energy of the standard cannon according to the second seismic records; the second target layer energy and the first target layer energy are obtained according to the same time window data and offset data;
the fitting module is used for determining a curve of the energy ratio changing along with the excitation parameters according to the second excitation parameters, the first target layer energy and the second target layer energy;
the calculation module is used for calculating the equivalent coverage times according to the curve;
the first energy module is specifically used for:
determining target time window data and target offset data at a target layer according to the first seismic record;
calculating first target layer energy of the standard cannon according to the target time window data and the target offset data;
the first energy module is specifically used for: calculating the energy of a first target layer of the standard gun by using the following formula:
wherein A is t For the energy of the target layer, N is the number of time window samples, and A is the amplitude of the seismic record.
4. The apparatus of claim 3, further comprising a determination module configured to:
acquiring a plurality of seismic data corresponding to a plurality of test excitation parameters;
determining a value of a first excitation parameter from the plurality of test excitation parameters based on the plurality of seismic data.
5. A computer device comprising a memory, a processor, the memory having stored therein a computer program executable on the processor, characterized in that the processor, when executing the computer program, implements the steps of the method of any of the preceding claims 1 to 2.
6. A computer readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the method of any one of the preceding claims 1 to 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010013587.5A CN113156494B (en) | 2020-01-07 | 2020-01-07 | Method and device for calculating equivalent coverage times of surface elements |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010013587.5A CN113156494B (en) | 2020-01-07 | 2020-01-07 | Method and device for calculating equivalent coverage times of surface elements |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113156494A CN113156494A (en) | 2021-07-23 |
| CN113156494B true CN113156494B (en) | 2023-10-31 |
Family
ID=76881343
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010013587.5A Active CN113156494B (en) | 2020-01-07 | 2020-01-07 | Method and device for calculating equivalent coverage times of surface elements |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113156494B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103543465A (en) * | 2012-07-12 | 2014-01-29 | 中国石油天然气集团公司 | Method for determining effective excitation points on basis of illumination energy of target strata |
| CN103576213A (en) * | 2012-07-27 | 2014-02-12 | 中国石油天然气集团公司 | Method for determining uniformity of properties of shot-geophone distances |
| CN105487117A (en) * | 2015-12-21 | 2016-04-13 | 中国石油天然气集团公司 | Three-dimensional earthquake observation system optimization method and apparatus |
| CN106066491A (en) * | 2016-08-02 | 2016-11-02 | 中国石油天然气集团公司 | A kind of high density wide-azimuth 3-d seismic exploration observation system optimization method |
| CN109884698A (en) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | Seismic prospecting observation system quantitative evaluation method based on target zone |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6545944B2 (en) * | 2001-05-30 | 2003-04-08 | Westerngeco L.L.C. | Method for acquiring and processing of data from two or more simultaneously fired sources |
-
2020
- 2020-01-07 CN CN202010013587.5A patent/CN113156494B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103543465A (en) * | 2012-07-12 | 2014-01-29 | 中国石油天然气集团公司 | Method for determining effective excitation points on basis of illumination energy of target strata |
| CN103576213A (en) * | 2012-07-27 | 2014-02-12 | 中国石油天然气集团公司 | Method for determining uniformity of properties of shot-geophone distances |
| CN105487117A (en) * | 2015-12-21 | 2016-04-13 | 中国石油天然气集团公司 | Three-dimensional earthquake observation system optimization method and apparatus |
| CN106066491A (en) * | 2016-08-02 | 2016-11-02 | 中国石油天然气集团公司 | A kind of high density wide-azimuth 3-d seismic exploration observation system optimization method |
| CN109884698A (en) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | Seismic prospecting observation system quantitative evaluation method based on target zone |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113156494A (en) | 2021-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (en) | Dolomite oil and gas reservoir distribution prediction method and device | |
| CN105425292A (en) | Oil and gas prediction method and oil and gas prediction device | |
| CN111399050B (en) | High-angle crack prediction method and device | |
| CN110442834B (en) | Method, device, computer equipment and readable storage medium for constructing logging curve | |
| US11209563B2 (en) | Data optimization method and integral prestack depth migration method | |
| CN104570066A (en) | Method for building seismic inversion low-frequency models | |
| CN110095811B (en) | Method and device for constructing and processing velocity model of paste rock stratum | |
| CN104330823A (en) | Method for determining observational parameters of vertical seismic profile | |
| CN109884698B (en) | Quantitative evaluation method of seismic exploration observation system based on target layer | |
| CN113156494B (en) | Method and device for calculating equivalent coverage times of surface elements | |
| CN103513279B (en) | A kind of illumination analysis computing method based on seismic wave equation and calculation element | |
| CN104614762A (en) | Loose sandstone gas reservoir boundary determination method and device | |
| CN111435174B (en) | Method and device for compensating amplitude of seismic data in strong reflection area | |
| CN104570115A (en) | Surface wave attenuating method and device | |
| CN109991663B (en) | Work area seismic velocity body correction method and device | |
| CN113126155B (en) | Sandstone reservoir prediction method aiming at strong reflection influence distributed among coal rocks | |
| CN104133239A (en) | Method for defining earthquake single-shot record target stratum time window | |
| CN113126146A (en) | Detection method, control device and storage medium for complex geology of canyon region | |
| CN108897046B (en) | Far detection sound wave reverse time migration imaging method and device | |
| CN113267810B (en) | Seismic exploration full-depth velocity modeling method and device | |
| CN111812716A (en) | Pre-stack quantitative prediction method, device and equipment for total organic carbon content of shale gas reservoir | |
| CN112630832A (en) | Method and device for predicting gas content based on variation of attenuation factor with incident angle | |
| CN108318936A (en) | A kind of STRATIGRAPHIC DIVISION treating method and apparatus | |
| CN111781661A (en) | Method and device for predicting sedimentary microphase plane spread in small-well area | |
| CN112285782B (en) | Near-surface seismic wave absorption attenuation investigation method and device |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |