CN111323818A

CN111323818A - Method and device for screening static correction mode of land seismic data

Info

Publication number: CN111323818A
Application number: CN201811522778.3A
Authority: CN
Inventors: 李劲松; 张明; 洪德明; 高银波; 黄荣善
Original assignee: Petrochina Co Ltd
Current assignee: Petrochina Co Ltd
Priority date: 2018-12-13
Filing date: 2018-12-13
Publication date: 2020-06-23
Anticipated expiration: 2038-12-13
Also published as: CN111323818B

Abstract

The application provides a method and a device for screening static correction modes of land seismic data, wherein the method comprises the following steps: acquiring synthetic seismic record groups corresponding to the corrected seismic data of the target area, wherein each synthetic seismic record group comprises calibrated synthetic seismic records of each drilled well; respectively determining the average speed of each drilled single well and the average speed of adjacent wells between every two adjacent drilled wells; and determining the static correction mode corresponding to the synthetic seismic record group meeting the preset conditions as a target static correction mode of the target area according to the values of the average speed of the single well and the average speed of the adjacent wells. The method and the device can more accurately and effectively screen the static correction mode of the land seismic data, the screening process is simple, efficient and high in sensitivity, the most suitable static correction mode of the land seismic data can be selected for different target areas, and therefore the pertinence and the accuracy of the static correction of the land seismic data can be effectively improved.

Description

Method and device for screening static correction mode of land seismic data

Technical Field

The application relates to the technical field of oil development, in particular to a method and a device for screening static correction modes of land seismic data.

Background

The static correction is an essential link in the conventional processing flow of the land seismic data, and the main method of the static correction processing is to correct the seismic data to a uniform reference surface which is generally a horizontal plane, so that the channels participating in stacking are close to the same phase, and the imaging quality is ensured. Therefore, a great deal of research work is put into the research and application of the static correction method, and a plurality of static correction modes are developed. Therefore, how to select a static correction mode suitable for a certain target area also becomes an important research topic in the field of seismic exploration.

In the prior art, the screening method of the static correction mode of the land seismic data generally comprises the steps of firstly determining a plurality of static correction modes, then one of the static correction modes is selected as a target static correction mode aiming at the target area according to the empirical value, or with reference to the difference in travel time of the same stratigraphic interface at the drilled well point and the seismic profile, however, because the basic assumptions of different static correction methods and the requirements on data bases are different, the static correction processing results obtained by adopting different static correction methods are often different, and the target static correction mode is selected only according to the empirical value or the difference of the same stratum interface when traveling on the drilled well point and the seismic profile, so that the accuracy of the static correction processing result aiming at the target area cannot be ensured due to the lack of a more definite screening standard, and the accuracy and the reliability of the whole seismic exploration process of the target area can be influenced.

Therefore, how to design a method capable of effectively screening the static correction mode is an urgent problem to be solved.

Disclosure of Invention

Aiming at the problems in the prior art, the application provides the method and the device for screening the static correction mode of the land seismic data, the static correction mode of the land seismic data can be more accurately and effectively screened, the screening process is simple, efficient and high in sensitivity, the most suitable static correction mode of the land seismic data can be selected according to different target areas, and the pertinence and the accuracy of the static correction of the land seismic data can be effectively improved.

In order to solve the technical problem, the application provides the following technical scheme:

in a first aspect, the present application provides a method for screening static correction modes of land seismic data, comprising:

acquiring synthetic seismic record groups corresponding to all corrected seismic data of a target area, wherein the corrected seismic data are different in static correction mode, and each synthetic seismic record group comprises calibrated synthetic seismic records of all drilled wells in the target area;

respectively determining the average speed of each drilled single well and the average speed of adjacent wells between each two adjacent drilled wells, which are respectively corresponding to each synthetic seismic record group;

and determining the static correction mode corresponding to the synthetic seismic record group meeting the preset condition as a target static correction mode of the target area according to the values of the average speed of each single well and the average speed of adjacent wells corresponding to each synthetic seismic record group.

Further, the separately determining the average velocity of each of the drilled single wells and the average velocity of the adjacent well between each adjacent two of the drilled single wells for each of the synthetic seismic record sets comprises:

determining a single-well average velocity of each of the drilled wells for a target formation from each of the drilled synthetic seismic records;

and determining the average speed of the adjacent wells between every two adjacent drilled wells.

Further, the determining an average velocity of adjacent wells between each adjacent two of the drilled wells comprises:

and determining the average speed of the adjacent well between two adjacent drilled wells according to the difference of the layered altitude depths and the seismic reflection travel time difference of the same section or the sub-section of the two adjacent drilled wells.

Further, the average velocity V of the adjacent well between the two adjacent drilled wells is determined by applying formula one:

V＝(D_w1-X-Bottom—D_w2-X-Bottom)/((T_w1-X-Bottom—T_w2-X-Bottom) /2) formula 1

Wherein, D is_w1-X-BottomAnd T_w1-X-BottomD is the layered altitude depth and seismic reflection travel of x sections of one of two adjacent drilled wells w1_w2-X-BottomAnd T_w2-X-BottomThe other of the two adjacent drilled wells w2 is at the layered elevation depth and seismic reflection travel of section x.

Further, the acquiring a synthetic seismic record group corresponding to each corrected seismic data of the target area includes:

respectively making each of the drilled synthetic seismic records in the synthetic seismic record group corresponding to each of the corrected seismic data;

and calibrating each of the drilled synthetic seismic records according to the wave group similarity of each of the drilled synthetic seismic records and the well-side seismic channel to obtain a synthetic seismic record group corresponding to each of the corrected seismic data.

Further, the determining, according to the values of the average velocity of each single well and the average velocity of adjacent wells corresponding to each synthetic seismic record group, the static correction mode corresponding to the synthetic seismic record group meeting a preset condition as a target static correction mode of the target area includes:

comparing the average velocity of each individual well in each synthetic seismic record group with the average velocity of the corresponding adjacent well;

and if the comparison result is smaller than a preset difference value, judging that the current synthetic seismic record group meets a preset condition, and determining the static correction mode corresponding to the synthetic seismic record meeting the preset condition as a target static correction mode of the target area.

Further, the determining, according to the values of the average velocity of each single well and the average velocity of adjacent wells corresponding to each synthetic seismic record group, the static correction mode corresponding to the synthetic seismic record group that meets the preset condition as the target static correction mode of the target area further includes:

sequentially judging whether the average velocity value of the adjacent wells in each synthetic seismic record group is within a preset stratum lithology velocity range;

if yes, judging that the current synthetic seismic record group meets a preset condition, and determining the static correction mode corresponding to the synthetic seismic record meeting the preset condition as a target static correction mode of the target area.

Further, before the obtaining the synthetic seismic record sets corresponding to the corrected seismic data of the target area, the method further includes:

correcting the seismic data of the target area by applying different static correction modes to obtain each corrected seismic data corresponding to each static correction mode, and

unifying the initial survey points of the plurality of well-drilled well log data within the target area with the initial survey points of the seismic data.

Further, the unifying the initial measurement point of the plurality of well log data of the drilled well within the target area and the initial measurement point of the seismic data includes:

determining the initial measuring point of each drilled well in the depth domain in the target area as the altitude depth according to the ground altitude and the square-complement height during drilling in each drilled well logging data in the target area, and adjusting each drilled well hierarchical data as the altitude depth to unify the initial measuring points of the plurality of drilled well logging data in the target area;

and determining the seismic travel time difference between the initial measuring point of the well depth of each drilled well and a processing datum plane based on the seismic data of each drilled well in the target area and the pre-acquired replacing speed, and determining the initial measuring point of the time domain as a seismic processing datum plane so as to unify the initial measuring points of the plurality of drilled seismic data in the target area.

Further, before the applying the different static correction modes to correct the seismic data of the target area, the method further includes:

stratigraphically dividing each of the drilled wells in the target area using the same stratigraphic division basis;

and performing an environmental correction on each of the well logs in the target zone.

In a second aspect, the present application provides a system for screening stationary corrections for land seismic data, comprising:

a synthetic seismic record group acquisition module, configured to acquire a synthetic seismic record group corresponding to each corrected seismic data of a target area, where static correction manners applied to each corrected seismic data are different, and each synthetic seismic record group includes a calibrated synthetic seismic record of each drilled well in the target area;

the single well and adjacent well average velocity determining module is used for respectively determining the single well average velocity of each drilled well corresponding to each synthetic seismic record group and the adjacent well average velocity between each two adjacent drilled wells;

and the static correction mode screening module is used for determining the static correction mode corresponding to the synthetic seismic record group meeting the preset condition as the target static correction mode of the target area according to the values of the average speed of each single well and the average speed of adjacent wells corresponding to each synthetic seismic record group.

Further, the single well and adjacent well average velocity determination module comprises:

the single-well average velocity obtaining unit is used for determining the single-well average velocity of each drilled well for a target stratum according to each drilled synthetic seismic record;

and the adjacent well average speed acquisition unit is used for determining the adjacent well average speed between every two adjacent drilled wells.

Further, the adjacent well average velocity obtaining unit is specifically configured to:

Further, the synthetic seismic record group acquisition module comprises:

a synthetic seismic record making unit for making each of the drilled synthetic seismic records in the synthetic seismic record group corresponding to each of the corrected seismic data, respectively;

and the synthetic seismic record calibration unit is used for calibrating each drilled synthetic seismic record according to the wave group similarity of each drilled synthetic seismic record and the seismic channel beside the well to obtain the synthetic seismic record group corresponding to each corrected seismic data.

Further, the static correction mode screening module comprises:

a comparison unit for comparing the average velocity of each single well in each synthetic seismic record group with the average velocity of the corresponding adjacent well;

and the first screening unit is used for judging that the current synthetic seismic record group meets a preset condition if the comparison result is smaller than a preset difference value, and determining the static correction mode corresponding to the synthetic seismic record meeting the preset condition as a target static correction mode of the target area.

Further, the static correction mode screening module further comprises:

the judging unit is used for sequentially judging whether the average velocity value of the adjacent wells in each synthetic seismic record group is within a preset stratum lithology velocity range;

and the second screening unit is used for judging that the current synthetic seismic record group meets a preset condition if the average velocity values of the adjacent wells in one synthetic seismic record group are all within a preset stratum lithology velocity range, and determining the static correction mode corresponding to the synthetic seismic record meeting the preset condition as a target static correction mode of the target area.

Further, the system for screening the land seismic data in the static correction mode further comprises:

a static correction module for correcting the seismic data of the target area by applying different static correction modes to obtain each corrected seismic data corresponding to each static correction mode, and

and the preprocessing module is used for unifying the initial measuring point of the plurality of well-drilled well logging data in the target area and the initial measuring point of the seismic data.

Further, the preprocessing module comprises:

the initial measurement point unifying unit of the well logging information is used for determining the initial measurement point of each drilled well in the target area in the depth domain as the altitude depth according to the ground altitude and the square complement height during drilling in each drilled well logging information in the target area, and adjusting each drilled well hierarchical data to the altitude depth so as to unify the initial measurement points of the plurality of drilled well logging information in the target area;

and the initial measuring point unifying unit of the seismic data is used for determining the seismic travel time difference between each drilled well depth initial measuring point and the processing datum plane based on each drilled well seismic data in the target area and the pre-acquired replacing speed, determining the initial measuring point of the time domain as the seismic processing datum plane and unifying the initial measuring points of the plurality of drilled well seismic data in the target area.

the stratum dividing module is used for carrying out stratum division on each drilled well in the target area by applying the same stratum dividing basis;

and the environment correction module is used for performing environment correction on each well-drilled well logging curve in the target area.

In a third aspect, the present application provides an electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the method for screening for statics of land seismic data when executing the program.

In a fourth aspect, the present application provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the method for static correction of land seismic data.

According to the technical scheme, the method and the device for screening the static correction mode of the land seismic data are characterized in that a synthetic seismic record group corresponding to each corrected seismic data of a target area is obtained, wherein the static correction mode applied to each corrected seismic data is different, and each synthetic seismic record group comprises the calibrated synthetic seismic records of each drilled well in the target area, so that accurate and effective data base can be improved for calculation of subsequent data, and the calculation and screening efficiency of the subsequent data can be effectively improved; by respectively determining the average speed of each drilled single well and the average speed of each adjacent well between each two adjacent drilled single wells corresponding to each synthetic seismic record group and introducing the concept of the average speed of the adjacent well stratum, a simple, accurate and effective data basis can be provided for screening of a subsequent static correction mode, so that the screening result of the subsequent static correction mode is accurate and efficient; the static correction mode corresponding to the synthetic seismic record group meeting the preset conditions is determined as the target static correction mode of the target area according to the values of the single-well average velocity and the adjacent-well average velocity corresponding to the synthetic seismic record groups, the target static correction mode most suitable for the target area can be accurately and effectively screened based on the concept of the adjacent-well stratum average velocity, the screening process is simple, efficient and high in sensitivity, the most suitable static correction mode of land seismic data can be selected according to different target areas, the pertinence and the accuracy of static correction of the land seismic data can be effectively improved, and the accuracy and the reliability of the whole seismic exploration process aiming at the target area can be improved.

Drawings

In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.

FIG. 1 is a schematic diagram of the communication between a server A1 and a client device B1 for performing static correction mode screening of land seismic data in accordance with the present invention.

FIG. 2 is a flow chart illustrating a method for screening static correction modes of land seismic data according to an embodiment of the invention.

FIG. 3 is a flowchart illustrating a step 200 of a method for screening static correction modes of land seismic data according to an embodiment of the invention.

FIG. 4 is a flowchart illustrating a step 100 of a method for screening static correction modes of land seismic data according to an embodiment of the present invention.

FIG. 5 is a flowchart illustrating an embodiment of step 300 of the method for screening static correction of land-based seismic data according to the present invention.

FIG. 6 is a flow chart illustrating another embodiment of step 300 of the method for screening static correction of land seismic data according to an embodiment of the present invention.

FIG. 7 is a schematic flow chart of steps A01 and A02 in the method for screening the static correction mode of land seismic data according to the embodiment of the invention.

FIG. 8 is a schematic flow chart of

steps

001 and 002 of the method for screening the stationary calibration mode of the land seismic data according to the embodiment of the present invention.

FIG. 9 is a schematic flow chart of a method for screening static correction modes of land seismic data according to the whole process example of the present invention.

FIG. 10 is a basic diagram of a unified well depth measurement and seismic travel survey initiation survey point in an example of an application of the present invention.

FIG. 11 is a basic diagram of the calculation of the average velocity of the adjacent well by using the difference of the layered altitude and the depth of the same or a sub-segment bottom boundary of the adjacent well and the travel time difference of the seismic reflection in the application example of the invention.

FIG. 12 is a schematic diagram of the division of four drilled strata and interwell formation in a work area in an example of the application of the present invention.

FIG. 13 is a schematic diagram showing the environmental calibration of a log for one of four drilled wells in a work area in accordance with an exemplary embodiment of the present invention.

FIG. 14 is a schematic diagram of a cross-well seismic profile of four drilled wells in a work area when seismic data are processed by a sand dune curve statics correction method in an application example of the invention.

FIG. 15 is a schematic cross-sectional view of a cross-well seismic profile of four drilled wells in a work area when seismic data are processed by a tomographic static correction method without using micro-logging constraints in an application example of the invention.

FIG. 16 is a schematic diagram of a cross-well seismic profile of four drilled wells in a work area when seismic data are processed by a micro-logging constrained tomographic statics method according to an embodiment of the present invention.

FIG. 17 is a schematic diagram of a synthetic seismic record based on seismic data obtained from a sand dune curve statics correction process in an example of an application of the present invention.

Fig. 18 is a schematic cross-sectional view of a well-connected seismic section with four drilled wells in a work area when seismic data processed by a sand dune curve statics correction method in an application example of the invention.

FIG. 19 is a schematic cross-sectional view of a cross-well seismic profile of four drilled wells in a work area when seismic data are processed by a tomographic static correction method without using micro-logging constraints in an application example of the invention.

FIG. 20 is a schematic diagram of a cross-well seismic profile of four drilled wells in a work area when seismic data are processed by a micro-logging constrained tomographic statics method according to an embodiment of the present invention.

FIG. 21 is a schematic structural diagram of a land seismic data static correction type screening apparatus according to an embodiment of the present invention.

Fig. 22 is a schematic structural diagram of an electronic device in an embodiment of the present invention.

Detailed Description

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

The theory of seismic exploration assumes that the excitation point and the receiving point are on a horizontal plane when acquiring seismic data, and the average velocity of the stratum is uniform, i.e. the so-called uniform horizontal laminar medium hypothesis. However, in actual land seismic exploration, the ground is often uneven, the depths of excitation points are different, and the propagation velocity of waves in a low-velocity zone is greatly different from that of waves in a diagenetic stratum, so that the time-distance curve shape of actual seismic acquisition data is influenced, and further the imaging of the seismic data is influenced. To eliminate these effects, it is necessary to perform terrain correction, excitation depth correction, low-velocity band correction, and the like on the original seismic data. The correction is invariable for different seismic reflection interfaces under the same acquisition observation point, so the correction is called static correction in general.

For the same set of seismic acquisition data, the static correction processing results obtained by adopting different static correction methods are different due to different basic assumptions of different static correction methods and different requirements on data bases. In the screening process of seismic data processing, processing results adopting different static correction methods need to be screened. Generally, the long-wavelength static correction problem is reflected remarkably on a seismic section, and mainly causes the seismic reflection layer to show consistent fluctuation from a shallow layer to a middle deep layer on the seismic section, so that the screening of the long-wavelength static correction is very clear and the operation is simple and convenient; however, for screening short-wavelength static correction and even residual static correction processing, the difference of reflection time of the same stratum interface on a drilled well and a seismic section caused by incomplete static correction problem solution is small, so that an effective static correction screening means is always lacked. Particularly, with the rapid development of fine static correction processing technology in recent years, the situation often faced in static correction processing screening is that the processing results obtained by adopting different static correction processing methods have the difference between the reflection time of the same stratum interface on a drilled well and the reflection time of the same stratum interface on a seismic section within an error allowable range, but the processing results of different static correction methods have larger difference in an undrilled well area between well points, and at the moment, due to the lack of a clear judgment standard, the screening of the static correction processing results is more difficult.

In consideration of the problem that the accuracy of a static correction processing result aiming at a target area cannot be ensured in a static correction mode of selecting land seismic data based on empirical values in the prior art, the application provides a screening method for static correction processing of land seismic data. By introducing the concept of the average velocity of the adjacent well stratums, the average velocity of the single well stratums and the average velocity of the adjacent well stratums are respectively calculated for the same section or the sub-section stratums on the basis of fine synthetic seismic record calibration. And (3) screening the static correction processing result of the onshore seismic data by comparing the consistency of the average speed of the stratums of the adjacent wells and the average speed of the stratums of the single well and judging whether the average speed of the stratums of the adjacent wells is within the normal range of the lithological speed of the stratums, and preferably obtaining the optimal static correction method.

Based on the above, the present application further provides a screening apparatus for land seismic data static correction processing, which may be a server a1, see fig. 1, where the server a1 may be in communication connection with a client device B1, a user may input collected relevant data of a target area into the client device B1, the client device B1 may send the relevant data to a server a1 online, the server a1 may receive the relevant data sent by the client device B1 online, and obtain a synthetic seismic record group corresponding to each corrected seismic data of the target area online or offline according to the relevant data, where static correction manners applied to each corrected seismic data are different, and each synthetic seismic record group includes a calibrated synthetic seismic record of each well drilled in the target area, respectively determining the average speed of each drilled well corresponding to each synthetic seismic record group and the average speed of each adjacent well between two adjacent drilled wells, determining the static correction mode corresponding to the synthetic seismic record group meeting preset conditions as a target static correction mode of the target area according to the values of the average speed of each single well and the average speed of the adjacent well corresponding to each synthetic seismic record group, and then sending the target static correction mode of the target area obtained by screening to the client device B1 on line by the server A1, so that a user can know the target static correction mode of the target area obtained by final screening through the client device B1.

It is understood that the client device B1 may include a smart phone, a tablet electronic device, a network set-top box, a portable computer, a desktop computer, a Personal Digital Assistant (PDA), a vehicle-mounted device, a smart wearable device, etc. Wherein, intelligence wearing equipment can include intelligent glasses, intelligent wrist-watch, intelligent bracelet etc..

In practical applications, the screening part for performing the static correction of the land seismic data may be performed on the side of the server a1 as described above, i.e., the architecture shown in fig. 1, or all operations may be performed in the client device B1. Specifically, the selection may be performed according to the processing capability of the client device B1, the limitation of the user usage scenario, and the like. This is not a limitation of the present application. If all of the operations are performed at the client device B1, the client device B1 may further include a processor for performing specific processing of the static correction mode of the land seismic data.

The client device may have a communication module (i.e., a communication unit), and may be communicatively connected to a remote server to implement data transmission with the server. For example, the communication unit may send the relevant data input by the user to the server, so that the server filters the static correction mode of the land seismic data of the target area according to the relevant data. The communication unit can also receive the filtering result returned by the server. The server may include a server on the task scheduling center side, and in other implementation scenarios, the server may also include a server on an intermediate platform, for example, a server on a third-party server platform that is communicatively linked to the task scheduling center server. The server may include a single computer device, or may include a server cluster formed by a plurality of servers, or a server structure of a distributed apparatus.

The server and the client device may communicate using any suitable network protocol, including network protocols not yet developed at the filing date of this application. The network protocol may include, for example, a TCP/IP protocol, a UDP/IP protocol, an HTTP protocol, an HTTPS protocol, or the like. Of course, the network Protocol may also include, for example, an RPC Protocol (Remote Procedure Call Protocol), a REST Protocol (Representational State Transfer Protocol), and the like used above the above Protocol.

In order to more accurately and effectively screen static correction modes of onshore seismic data, an embodiment of the present application provides a specific implementation of a screening method of static correction modes of onshore seismic data, and referring to fig. 2, the screening method of static correction modes of onshore seismic data specifically includes the following contents:

step 100: and acquiring a synthetic seismic record group corresponding to each corrected seismic data of the target area, wherein the corrected seismic data are different in static correction mode, and each synthetic seismic record group comprises calibrated synthetic seismic records of each drilled well in the target area.

In step 100, the screening device for the static correction method of the land seismic data first obtains each corrected seismic data after being corrected by different static correction methods, wherein the static correction method adopted by each corrected seismic data is not repeated.

Then, referring to table 1, assuming that there are 10 drilled wells in the target area, four static correction modes are provided, and the four static correction modes correspond to the corrected seismic data Z1 to Z4, respectively, the land seismic data static correction mode screening apparatus first obtains a set of calibrated synthetic seismic records H11 to H110 corresponding to the drilled wells W1 to W10 in the target area corresponding to the corrected seismic data Z1, respectively, to form a synthetic seismic record group S1. And then, replacing another piece of corrected seismic data Z2, and acquiring calibrated synthetic seismic records H21 to H210 corresponding to the drilled wells W1 to W10 in the target area corresponding to Z2 to form a synthetic seismic record group S2. And repeating the steps until all the synthetic seismic record groups corresponding to the corrected seismic data are obtained.

TABLE 1

It is understood that the static correction method at least includes the following methods:

(1) static correction of field elevation: it is only applicable to areas where the low-dropout belt is not present or where the low-dropout belt does not vary laterally.

(2) Model static correction: the basis of the model method is to establish a model capable of accurately describing the geological geophysical attributes of the near-surface medium. In the traditional model static correction, the description of the geophysical attributes of a near-surface medium is obtained through conventional near-surface investigation methods such as small refraction and micro-logging, and then a near-surface model is obtained through spatial interpolation, so that the correction of a datum plane is completed.

(3) Refraction static correction: the refraction static correction is also a model static correction method strictly speaking, but the model establishment is different from the traditional method, the refraction surface speed and the delay time are mainly obtained by solving equation inversion, then a speed model is established by means of the surface speed, and the calculation of the static correction value is completed on the basis.

(4) First arrival static correction: the first-arrival static correction is the same as the refraction static correction, strictly belongs to model static correction, and neither accurately describes the geophysical attributes of near-surface media nor establishes a laminar refraction model, but obtains the velocity field distribution of the near-surface layer through tomography inversion by utilizing first-arrival time to obtain an optimal model.

(5) Residual static correction: the theoretical assumption and implementation method, which is generally referred to as reflection static correction and reflection residual static correction, determines its inevitable limitations, so its application must be based on good reference plane static correction.

(6) Relative refraction static correction: the method is a method between the static correction of a reference surface and the residual static correction, two conditions which are necessarily met by the refraction static correction are avoided, only partial shot-inspected refracted waves with good quality are concerned, although an accurate model cannot be established, accurate high-frequency components and partial medium-frequency components can be obtained, the method is used as a supplement of a model method, necessary correction can be carried out on the static correction value obtained by the model method, and meanwhile, most high-frequency static correction components caused by a low-speed-reducing zone can be eliminated.

Step 200: and respectively determining the average speed of each drilled single well and the average speed of adjacent wells between each two adjacent drilled wells, which correspond to each synthetic seismic record group.

In step 200, the land seismic data static correction mode screening device determines the average velocity of each drilled single well and the average velocity of each adjacent two drilled neighboring wells corresponding to each synthetic seismic record set. It can be understood that, because the calculated average speed of the adjacent wells has a more definite physical meaning, the calculated average speed can be compared with the average speed of the actual stratum or the lithology speed range of the stratum, so that the quantification degree is effectively improved.

For example, based on table 1, the land seismic data static correction mode screening apparatus needs to obtain 4 × 10 average velocities of the drilled single wells and average velocities of adjacent wells between each two adjacent drilled wells according to calibrated synthetic seismic records H11 to H110, H21 to H210, H31 to H310, and H41 to H410 corresponding to all the drilled wells W1 to W10, respectively.

Step 300: and determining the static correction mode corresponding to the synthetic seismic record group meeting the preset condition as a target static correction mode of the target area according to the values of the average speed of each single well and the average speed of adjacent wells corresponding to each synthetic seismic record group.

In step 300, the screening apparatus for static correction modes of land seismic data respectively determines whether the synthetic seismic record group corresponding to each static correction mode satisfies a preset condition according to the values of the average velocity of each single well and the average velocity of the adjacent well corresponding to each synthetic seismic record group, where the preset condition may be that the difference between the average velocity of each drilled single well and the average velocity of the corresponding adjacent well is smaller than a preset value, or that the average velocity of the adjacent well is within a preset range, and if it is determined that one synthetic seismic record group satisfies the preset condition, the static correction mode uniquely corresponding to the synthetic seismic record group is used as the target static correction mode of the target area.

It can be understood that, if it is determined that more than one synthetic seismic record group satisfies the preset condition, the static correction modes corresponding to the synthetic seismic record groups satisfying the preset condition may all be used as the target static correction modes of the target area.

In addition, if it is determined that any one synthetic seismic record group does not meet the preset rule for the current target area, selecting a new static correction mode which is not used, and executing the steps 100 to 300 again until at least one synthetic seismic record meeting the preset condition is obtained.

As can be seen from the above description, in the method for screening static correction modes of land seismic data provided in the embodiment of the present invention, by obtaining synthetic seismic record groups corresponding to respective corrected seismic data of a target area, where the static correction modes applied to the respective corrected seismic data are different, and each synthetic seismic record group includes a calibrated synthetic seismic record of each drilled well in the target area, an accurate and effective data base can be improved for calculation of subsequent data, and the efficiency of calculation and screening of the subsequent data can be effectively improved; by respectively determining the average speed of each drilled single well and the average speed of each adjacent well between each two adjacent drilled single wells corresponding to each synthetic seismic record group and introducing the concept of the average speed of the adjacent well stratum, a simple, accurate and effective data basis can be provided for screening of a subsequent static correction mode, so that the screening result of the subsequent static correction mode is accurate and efficient; the static correction mode corresponding to the synthetic seismic record group meeting the preset conditions is determined as the target static correction mode of the target area according to the values of the single-well average velocity and the adjacent-well average velocity corresponding to the synthetic seismic record groups, the target static correction mode most suitable for the target area can be accurately and effectively screened based on the concept of the adjacent-well stratum average velocity, the screening process is simple, efficient and high in sensitivity, the most suitable static correction mode of land seismic data can be selected according to different target areas, the pertinence and the accuracy of static correction of the land seismic data can be effectively improved, and the accuracy and the reliability of the whole seismic exploration process aiming at the target area can be improved.

In order to further improve the accuracy of screening the static correction method, in a specific embodiment, referring to fig. 3, the step 200 in the method for screening the static correction method for the land seismic data specifically includes the following steps:

step 201: determining a single-well average velocity for each of the drilled synthetic seismic records for a target formation from each of the drilled synthetic seismic records.

Step 202: determining an average velocity of adjacent wells between each adjacent two of the drilled wells.

In a specific embodiment, the step 202 specifically includes the following steps:

and the screening device of the static correction mode of the land seismic data determines the average speed of the adjacent wells between the two drilled wells according to the layered altitude depth difference and the seismic reflection travel time difference of the same section or the sub-section of the two adjacent drilled wells.

Wherein the average velocity V of the adjacent well between the two adjacent drilled wells is determined by applying the formula I:

Wherein, D is_w1-X-BottomAnd T_w1-X-BottomD is the layered altitude depth and seismic reflection travel of x sections of one of two adjacent drilled wells w1_w2-X-BottomAnd T_w2-X-BottomRespectively in two adjacent drilled wellsAnother well w2 is drilled at the layered elevation depth and seismic reflection travel of section x.

It can be understood that the layered elevation depth difference and the seismic reflection travel time difference of the same section or the sub-section bottom boundary between adjacent wells are obtained, the layered elevation depth difference is taken as a numerator, a half of the seismic reflection travel time difference is taken as a denominator, and the quotient is defined as the adjacent well stratum average velocity, which is hereinafter referred to as the adjacent well average velocity. Meanwhile, the concept of the average velocity of the single-well stratum is used, namely, the quotient of the thickness of one section or sub-section of stratum on the single well and half of the corresponding seismic reflection travel time is called the average velocity of the single well for short on the basis of fine synthetic seismic record calibration.

The traditional static correction screening mainly refers to the difference of the same stratum interface between a drilled well point and a seismic profile during traveling, and the difference is used as a denominator, so that the traveling error caused by improper static correction is amplified, and the sensitivity of the screening method is higher for the condition of smaller static correction error; meanwhile, the calculated average speed of the adjacent well has a more definite physical meaning and can be compared with the actual average speed of the stratum or the lithologic speed range of the stratum, so that the quantification degree is better.

In order to further improve the screening accuracy of the static correction mode by providing a more accurate data base, referring to fig. 4, in a specific embodiment, the step 100 of the method for screening the static correction mode of the land seismic data specifically includes the following steps:

step 101: and respectively making each drilled synthetic seismic record in the synthetic seismic record group corresponding to each corrected seismic data.

It will be appreciated that the apparatus for statics-corrected screening of land-based seismic data utilizes logs, primarily sonic and density logs, to label the elevation data of each section or sub-section of the formation bed boundary onto the logs. For each drilled well, synthetic seismic records are made, the seismic wavelets from which the synthetic seismic records are made being determined by the researcher or extracted from the seismic data.

Step 102: and calibrating each drilled synthetic seismic record according to the wave group similarity of each drilled synthetic seismic record and the seismic channel beside the well to obtain a synthetic seismic record group corresponding to each corrected seismic data.

It can be understood that the land seismic data static correction mode screening device carries out calibration of the synthetic seismic records according to the wave group similarity of the synthetic seismic records and the well-side seismic channels. After calibration of the synthetic seismic record, the conversion relation between the layered elevation depth of each section or sub-section stratum bottom boundary of each well and the seismic reflection travel time is determined. For each well, the layered elevation depth of each segment or sub-segment of the formation bottom boundary corresponds to a unique seismic reflection travel time.

In order to further improve the efficiency and reliability of the static correction mode screening, referring to fig. 5, in an embodiment, the step 300 of the method for screening static correction modes of land seismic data specifically includes the following steps:

step 301A: comparing the average velocity of each individual well in each synthetic seismic record set to the average velocity of the corresponding adjacent well.

Step 302A: and if the comparison result is smaller than a preset difference value, judging that the current synthetic seismic record group meets a preset condition, and determining the static correction mode corresponding to the synthetic seismic record meeting the preset condition as a target static correction mode of the target area.

In order to further improve the efficiency and reliability of the static correction mode screening, referring to fig. 6, in another embodiment, the step 300 of the method for screening static correction modes of land seismic data specifically includes the following steps:

step 301B: and sequentially judging whether the average velocity value of the adjacent wells in each synthetic seismic record group is within a preset stratum lithology velocity range.

Step 302B: if yes, judging that the current synthetic seismic record group meets a preset condition, and determining the static correction mode corresponding to the synthetic seismic record meeting the preset condition as a target static correction mode of the target area.

It is understood that the two specific embodiments of step 300 may be performed simultaneously, and the specific screening principle is to compare the average velocity of the adjacent well with the average velocity of the single well, and if there is better consistency between the two, or the numerical range of the average velocity of the adjacent well is within the normal lithologic velocity range of the stratum, it is determined that the fluctuation of the seismic reflection interface between the two adjacent wells is the seismic travel time difference caused by the average velocity of the section or the sub-section, which indicates that the static correction process does not change the geological geophysical characteristics of the section or the sub-section of the stratum, and the static correction method is reasonably applied; if the difference between the average speed of the adjacent well and the average speed of the single well is large, or the numerical range of the average speed of the adjacent well is not in the normal lithological speed range of the stratum, the static correction treatment destroys the geological and geophysical characteristics of the stratum of the section or the sub-section, and the application of the static correction method is unreasonable. Namely, the screening device of the static correction mode of the land seismic data compares the average velocity of each single well in each synthetic seismic record group with the average velocity of the corresponding adjacent well, and then sequentially judges whether the average velocity of the adjacent well in each synthetic seismic record group is within a preset range of stratigraphic lithology velocities, if the synthetic seismic record group simultaneously meets the two judgment conditions, the static correction mode corresponding to the synthetic seismic record is determined as the target static correction mode of the target area.

In order to further improve the screening accuracy of the static correction mode by providing a more accurate data base, referring to fig. 7, in a specific embodiment, step 100 of the method for screening the static correction mode of the land seismic data of the present application further includes the following preprocessing contents:

step A01: correcting the seismic data of the target area by applying different static correction modes to obtain each corrected seismic data corresponding to each static correction mode, and

step A02: unifying the initial survey points of the plurality of well-drilled well log data within the target area with the initial survey points of the seismic data.

In a specific embodiment, the step a02 specifically includes the following steps:

(1) according to the ground elevation and the square-of-complement height during drilling in each well-drilled well-logging data in the target area, determining the initial measuring point of each well-drilled well in the depth domain in the target area as the elevation depth, and adjusting each well-drilled layered data to the elevation depth to unify the initial measuring points of the well-drilled well-logging data in the target area.

(2) And determining the seismic travel time difference between the initial measuring point of the well depth of each drilled well and a processing datum plane based on the seismic data of each drilled well in the target area and the pre-acquired replacing speed, and determining the initial measuring point of the time domain as a seismic processing datum plane so as to unify the initial measuring points of the plurality of drilled seismic data in the target area.

It is understood that the execution sequence of step a01 and step a02 is not sequential.

Based on the above, referring to fig. 8, the following is also specifically included before step a 01:

step 001: stratigraphically partitioning each of the drilled wells in the target zone using the same stratigraphic partitioning basis.

Step 002: performing an environmental correction on each of the well logs in the target zone.

From the above description, it can be known that the method for screening the static correction mode of the land seismic data provided by the embodiment of the application can more accurately and effectively screen the static correction mode of the land seismic data, has a simple, efficient and high sensitivity screening process, can select the most suitable static correction mode of the land seismic data for different target areas, and further can effectively improve the pertinence and the accuracy of the static correction of the land seismic data.

To further explain the scheme, the present application further provides a full-flow example of the method for screening static correction modes of land seismic data, and referring to fig. 9, the full-flow example of the method for screening static correction modes of land seismic data specifically includes the following contents:

s1, researching the drilled stratum division and the inter-well system layer in the work area, wherein the specific contents are as follows:

and S101, according to the difference of research precision, the precision of stratum division is also different. Generally, stratigraphic divisions will be made at different levels of groups, segments, layers, etc. One well has completed the depth correction, stratigraphic division to group level, etc. of the underlying geological work after geological completion. In the oil and gas exploration stage, researchers mainly divide each well-drilled stratum into stratum units such as groups and sections according to well logging curves (mainly including natural gamma, natural potential and resistivity) and information such as rock debris logging and according to the lithology and electrical characteristics of the stratum, and then subdivide the sections into sub-sections.

S102, because each drilled well is divided into strata independently, and because the completion time is different, the phenomenon that the strata division among wells is not uniform often exists. During the overall research, the new stratigraphic division comparison needs to be carried out on each drilled group, section and sub-section so as to ensure the consistency of the stratigraphic division basis among wells.

S2, correcting the environment of the logging curve, specifically comprising the following steps:

synthetic seismic records are an important means of establishing the conversion relationship between the depth of the earth strata and the travel time of seismic reflection, and well logs are an important bridge in the synthetic seismic records. Before making a synthetic seismic record, the environmental correction needs to be performed on the well-drilled well logging curve to eliminate the influence on the well logging curve due to objective factors such as different well logging times, different well logging instruments, well bore collapse and the like.

S3, unifying the initial measuring points of the well depth measurement and the seismic travel measurement, wherein the specific contents are as follows:

and unifying the well depth measurement and the initial measurement point of the seismic travel time measurement. To ensure correct depth and seismic travel time conversion, the initial survey point of the logging data and the initial survey point of the seismic data need to be unified in the depth domain and the seismic reflection travel time domain.

S301, determining the initial measuring point of the depth domain as the sea level. The depth values of the logging curve and the drilled well stratification are based on the drilling bushing position as an initial measurement point, and when the initial measurement point is adjusted to the sea level, the ground elevation of the well point and the square bushing height during drilling need to be utilized. If the ground elevation of the well point is recorded as Ds and the bushing height of the drilling machine is recorded as KB, the elevation depth D of the initial measurement point of the well logging is measured_w＝D_s+ KB, accordingly, all of the drilled layered data is also adjusted to elevation depth. The seismic data processing reference plane is an equal-depth plane, and the elevation depth is recorded as Dd. The unification of the depth field starting measurement points is now completed.

S302, determining a starting measuring point of a time domain as a seismic processing reference plane, namely a plane with 0 millisecond seismic reflection travel time. Knowing the seismic travel time difference between the well depth starting measurement point and the processing datum, the well depth time measurement starting point can be adjusted to the seismic processing datum. In the actual seismic data processing, the travel time difference is calculated as (D)_w—D_d) Replacement speed. The replacement speed is a basic parameter for seismic data processing, can be obtained by calculation after field small refraction acquisition or micro-logging measurement, and can also be obtained by calculation by using the relation between the offset distance and the first arrival time after a processing person picks up the first arrival. And after the travel time difference between the well depth initial measuring point and the seismic datum plane is considered in the synthetic record, the initial measuring points are unified in a time domain.

S4, calibrating the drilled synthetic seismic record, which comprises the following specific contents:

s401, making a synthetic seismic record. And marking the layered elevation data of each section or sub-section of the stratum bottom boundary of the stratum system layer subjected to stratum division and stratum system layer in the S1 on the logging curve by using the logging curve of the elevation depth, mainly a sonic logging curve and a density logging curve, which is subjected to environment correction and depth conversion in the S2. Selecting a synthetic seismic record making method considering Ds, KB, Dd and replacement speed, making a synthetic seismic record for each drilled well, and making seismic wavelets of the synthetic seismic record determined by researchers or extracted from seismic data.

S402, calibrating the synthetic seismic record according to the wave group similarity of the synthetic seismic record and the well-side seismic channel. After calibration of the synthetic seismic record, the conversion relation between the layered elevation depth of each section or sub-section stratum bottom boundary of each well and the seismic reflection travel time is determined. For each well, the layered elevation depth of each segment or sub-segment of the formation bottom boundary corresponds to a unique seismic reflection travel time.

S5, calculating the average speed of the single well and the average speed of the adjacent wells according to the time-depth relation determined by the synthetic seismic records, wherein the specific contents are as follows:

s501, for all drilled wells in the work area, calculating the average speed of the single well according to the single well synthetic seismic record and the stratum near the target layer. Taking the well in FIG. 10 as an example, the elevation depth D of the top and bottom strata of the Y-section stratum is read₁And D₂And corresponding seismic reflection travel time T₁And T₂The average speed of a single well of the Y-section stratum can be obtained:

V_Yave＝(D₂-D₁)/((T₂-T₁)/2)

and S502, defining the average speed of the adjacent wells. Taking two adjacent wells w1 and w2 in FIG. 11 as an example, for the X-section stratum, the depth of bottom boundary elevation of well w1 is D_w1X-BottomCorresponding seismic reflection travel time T_w1-X-Bottom(ii) a The bottom boundary altitude depth of the well w2 is D_w2-X-BottomCorresponding seismic reflection travel time T_w2-X-Bottom. Then for the two adjacent wells, the average velocity V of the adjacent wells of the X section stratum_X12Is defined as:

V_X12＝(D_w1-X-Bottom—D_w2-X-Bottom)/((T_w1-X-Bottom—T_w2-X-Bottom)/2)

s503, calculating the average speed of the adjacent wells for all the adjacent wells in the work area according to the stratum near the target layer.

S6, performing static correction processing screening, specifically including the following contents:

s601, repeating S4-S5 on the seismic data obtained by adopting different static correction processing methods, and completing the calibration of the drilled synthetic seismic record on different static correction data.

S602, based on the same set of seismic data, aiming at the stratum near the target layer, calculating the average speed of a single well and the average speed of adjacent wells;

and S603, comparing the average speed of the single well with the average speed of the adjacent wells for the same section or the sub-section stratum based on the same set of seismic data. If the two are relatively good in consistency or the numerical range of the average velocity of the adjacent wells is in the normal lithological velocity range of the stratum, judging that the fluctuation of the seismic reflection interface between the two adjacent wells is the seismic travel time difference caused by the average velocity of the section or the sub-section, and indicating that the geological geophysical characteristics of the section or the sub-section stratum are not changed by static correction processing, wherein the static correction method is reasonably applied; if the difference between the average speed of the adjacent well and the average speed of the single well is large, or the numerical range of the average speed of the adjacent well is not in the normal lithological speed range of the stratum, the static correction treatment destroys the geological and geophysical characteristics of the stratum of the section or the sub-section, and the application of the static correction method is unreasonable.

S604, repeat S501 to S603 for different static correction processing data. And preferably selecting the optimal static correction method according to the judgment result of the S603. At this time, the screening of the static correction process is completed.

From the above description, it can be seen that the method for screening the static correction mode of the land seismic data provided by the full-flow example of the present application can improve an accurate and effective data base for the calculation of the subsequent data and can effectively improve the efficiency of the calculation and screening of the subsequent data by obtaining the synthetic seismic record groups corresponding to the respective corrected seismic data of the target area, wherein the static correction modes applied to the respective corrected seismic data are different, and each synthetic seismic record group includes the calibrated synthetic seismic record of each drilled well in the target area; by respectively determining the average speed of each drilled single well and the average speed of each adjacent well between each two adjacent drilled single wells corresponding to each synthetic seismic record group and introducing the concept of the average speed of the adjacent well stratum, a simple, accurate and effective data basis can be provided for screening of a subsequent static correction mode, so that the screening result of the subsequent static correction mode is accurate and efficient; the static correction mode corresponding to the synthetic seismic record group meeting the preset conditions is determined as the target static correction mode of the target area according to the values of the single-well average velocity and the adjacent-well average velocity corresponding to the synthetic seismic record groups, the target static correction mode most suitable for the target area can be accurately and effectively screened based on the concept of the adjacent-well stratum average velocity, the screening process is simple, efficient and high in sensitivity, the most suitable static correction mode of land seismic data can be selected according to different target areas, the pertinence and the accuracy of static correction of the land seismic data can be effectively improved, and the accuracy and the reliability of the whole seismic exploration process aiming at the target area can be improved.

In one specific example of use, the study target zone is a dwarfism sandstone formation. There were 4 drilled wells in the study area at one coverage, see figure 12, where only 2 wells were in full coverage, QianShao1 and M16, respectively.

(1) Dividing the drilled stratum and making the interwell layer. And carrying out stratum division and interwell formation on 4 drilled wells in the primary coverage range of the research work area. According to the logging curve (mainly natural gamma, natural potential and resistivity) and the data of rock debris logging, each well-drilled stratum is divided into stratum units such as groups and sections according to the lithological and electrical characteristics of the stratum, and the sections are subdivided into sub-sections. The study target layer is a Jurassic sandstone-shale stratum, and the work of stratum layering and interwell systematical layers is mainly aimed at the stratum near the target layer and is subdivided into the grades of sections and subsections. Based on the original layering of predecessors, according to a unified standard, K of Qianshao1₁Q₁Adjusting the altitude to-3477 to-3512M, and adjusting the J of M011₁S₂ ¹Adjusting the altitude to-4008 to-4035M, and adjusting the K of M003₁Q₁Adjusted to-3612 m by altitude-3629.

(2) And performing logging curve quality inspection on four drilled wells in the work area, and performing environmental correction on logging curves of 2 wells. FIG. 13 is a comparison of log curves before and after calibration of the environment of 1 well therein. Wherein P-Wave refers to a curve obtained after a horizontal scale is converted into speed by an acoustic logging curve, Den represents a density logging curve, GR represents a natural gamma logging curve, and Cal represents a hole diameter curve. Within the range defined by the black ellipse, the density logging has abnormal values due to the collapse of the well diameter, and the precision of the synthetic seismic record is seriously influenced. And (5) performing borehole collapse correction on the density curve by using a statistical method. The blue density curve is the original measurement curve and the red density curve is the corrected density curve.

(3) 4 drilled wells are in the range of the primary coverage area in the work area, and fig. 14, 15 and 16 are superposed sections obtained after processing by a sand dune curve static correction method, a chromatography static correction method without using micro-logging constraint and a chromatography static correction method using micro-logging constraint respectively. The red single trace seismic trace at the well location in the figure is the synthetic seismic record for that well, and the orange color on the curve is labeled as the stratigraphic hierarchy name and the location of the bottom boundary on the synthetic record. From the comparison of the synthetic seismic records with the actual seismic sections, the single-well synthetic seismic record of the three sections is calibrated within a reasonable error range; the blue lines are the same seismic picking horizon. According to the contrast of the blue line and the seismic strong amplitude reflection interface, the three sections are mainly shown as the stratum fluctuation difference among wells. It can be seen that the processing results obtained by using different static correction processing methods have small difference in processing results at each well point, and have large difference in processing results of different static correction methods in the non-drilling area.

(4) And (4) giving a seismic processing datum plane, replacing the speed, the ground elevation of the well point and the bushing height of the well point, and making a synthetic seismic record. A synthetic record was made for one of the wells based on seismic data from the dune curve statics process, see figure 17.

(5) Because the inter-well average speed calculation sensitivity adopted by the invention is very high, the inter-well average speed is calculated only for two wells with a full coverage area with stable processing results, and is compared with the single-well average speed. FIG. 18 is a cross-section of four well seismic data processed by a sand dune curve statics method in a work area. After the synthetic seismic records of the drilled wells are calibrated, elevation values of data of section or sub-section bottom boundary layering of the single well, corresponding seismic travel time values and the calculated average speed of the single well are marked beside the wells respectively, and the average speed of adjacent wells is marked between the two wells. By comparison, it can be seen that the upper and lower values of the average velocity of the adjacent wells are stable, are closer to the average velocity of the single well, and are within a reasonable range of the average velocity of the sand shale formation, fig. 18 is a graph obtained by calibrating a fine synthetic seismic record, and for two wells in a full coverage area, the layered elevation depth and seismic reflection travel time of the bottom boundary of each section or sub-section are marked on a seismic section, and the calculated average velocity of the adjacent wells and the calculated average velocity of the single well are also marked on the seismic section.

(6) Fig. 19 and 20 are cross-well seismic profiles of four drilled wells in a work area, and the seismic data are processed by tomographic statics without using micro-logging constraints and tomographic statics with micro-logging constraints, respectively. After the synthetic seismic records of the drilled wells are calibrated, elevation values of data of section or sub-section bottom boundary layering of the single well, corresponding seismic travel time values and the calculated average speed of the single well are marked beside the wells respectively, and the average speed of adjacent wells is marked between the two wells. By comparison, the average speed of the adjacent wells has large up-and-down change, and the difference with the average speed of a single well is large, and the average speed of the adjacent wells is far beyond the reasonable range of the average speed of the sand-shale stratum, so that the average speed level of the compact carbonate rock is reached. Wherein, fig. 19 is the two wells calibrated by the fine synthetic seismic record, and for the full coverage area, the layered elevation depth and seismic reflection travel time of each section or sub-section bottom boundary are marked on the seismic profile, and the calculated average velocity of the adjacent well and the average velocity of the single well are also marked on the seismic profile. FIG. 20 is a diagram of two wells in a full coverage area, wherein the layered elevation depth and seismic reflection travel time of the bottom boundary of each section or sub-section are marked on a seismic section, and the calculated average velocities of adjacent wells and the calculated average velocity of a single well are also marked on the seismic section.

(7) Through comparison of fig. 18, 19 and 20, it is judged that the seismic data of the sand dune curve static correction process best maintains the geological geophysical characteristics of the stratum, while the other two static correction processes destroy the geological geophysical characteristics of the stratum, and the sand dune curve static correction process is preferably the static correction method of the embodiment. The static correction process screening for the examples is completed at this point.

In order to more accurately and effectively screen the static correction method of the onshore seismic data, an embodiment of the present application provides a specific implementation of a screening apparatus for the static correction method of the onshore seismic data, which can implement all the contents in the screening method for the static correction method of the onshore seismic data, and referring to fig. 21, the screening apparatus for the static correction method of the onshore seismic data specifically includes the following contents:

a synthetic seismic record group obtaining module 10, configured to obtain a synthetic seismic record group corresponding to each corrected seismic data of a target area, where static correction manners applied to each corrected seismic data are different, and each synthetic seismic record group includes a calibrated synthetic seismic record of each drilled well in the target area.

And a single well and adjacent well average velocity determination module 20, configured to determine an average velocity of each drilled single well and an average velocity of each adjacent two drilled adjacent wells corresponding to each synthetic seismic record group, respectively.

And the static correction mode screening module 30 is configured to determine, according to the values of the average velocity of each single well and the average velocity of adjacent wells corresponding to each synthetic seismic record group, the static correction mode corresponding to the synthetic seismic record group that meets a preset condition as a target static correction mode of the target area.

The embodiment of the device for screening static correction modes of land seismic data provided by the application can be specifically used for executing the processing flow of the embodiment of the method for screening static correction modes of land seismic data in the above embodiment, and the function of the device is not described herein again, and reference may be made to the detailed description of the method embodiment.

As can be seen from the above description, in the screening apparatus for static correction of land seismic data provided in the embodiment of the present invention, the synthetic seismic record group acquisition module 10 is used to acquire the synthetic seismic record group corresponding to each corrected seismic data of the target area, where the static correction applied to each corrected seismic data is different, and each synthetic seismic record group includes the calibrated synthetic seismic record of each well drilled in the target area, so that an accurate and effective data basis can be improved for calculation of subsequent data, and the efficiency of calculation and screening of subsequent data can be effectively improved; the average speed of each drilled single well and the average speed of each adjacent two drilled adjacent wells which correspond to each synthetic seismic record group are respectively determined by the average speed determination module 20 of the single well and the average speed of the adjacent well, and a simple, accurate and effective data basis can be provided for screening of a subsequent static correction mode by introducing the concept of the average speed of the adjacent well stratum, so that the screening result of the subsequent static correction mode is accurate and efficient; the static correction mode corresponding to the synthetic seismic record group meeting the preset conditions is determined as the target static correction mode of the target area through the static correction mode screening module 30 according to the values of the single-well average velocity and the adjacent-well average velocity corresponding to the synthetic seismic record groups, the target static correction mode most suitable for the target area can be accurately and effectively screened based on the concept of the adjacent-well stratum average velocity, the screening process is simple, efficient and high in sensitivity, the most suitable static correction mode of land seismic data can be selected for different target areas, the pertinence and the accuracy of static correction of the land seismic data can be effectively improved, and the accuracy and the reliability of the whole seismic exploration process for the target area can be improved.

In order to further improve the screening accuracy of the static correction method, in an embodiment, the average velocity determination module 20 of the single well and the adjacent wells in the static correction method screening apparatus for onshore seismic data of the present application specifically includes the following contents:

the single-well average velocity obtaining unit 21 is configured to determine a single-well average velocity of each drilled well for the target formation according to each drilled synthetic seismic record;

and an adjacent well average speed obtaining unit 22, configured to determine an adjacent well average speed between each adjacent two of the drilled wells.

The adjacent well average speed obtaining unit 22 is specifically configured to:

and determining the average speed of the adjacent well between two adjacent drilled wells according to the difference of the layered altitude depths and the seismic reflection travel time difference of the same section or the sub-section of the two adjacent drilled wells. And the average speed V of the adjacent well between the two adjacent drilled wells is determined by applying the formula I:

In order to further improve the screening accuracy of the static correction mode by providing a more accurate data base, in a specific embodiment, the synthetic seismic record set acquisition module 10 in the static correction mode screening device for land seismic data of the application specifically includes the following contents:

a synthetic seismic record making unit 11, configured to make each of the drilled synthetic seismic records in the synthetic seismic record group corresponding to each of the corrected seismic data, respectively, according to each of the corrected seismic data;

and the synthetic seismic record calibration unit 12 is configured to calibrate each of the drilled synthetic seismic records according to the wave group similarity of each of the drilled synthetic seismic records and the well-side seismic channel, so as to obtain a synthetic seismic record group corresponding to each of the corrected seismic data.

In order to further improve the efficiency and reliability of the static correction mode screening, in an embodiment, the static correction mode screening module 30 in the device for screening static correction modes of land seismic data of the present application specifically includes the following contents:

a comparison unit 31A for comparing the average velocity of each single well in each of the synthetic seismic record groups with the average velocity of the corresponding adjacent well;

the first screening unit 32A is configured to determine that the current synthetic seismic record group meets a preset condition if the comparison result is smaller than a preset difference, and determine the static correction mode corresponding to the synthetic seismic record meeting the preset condition as a target static correction mode of the target area.

In order to further improve the efficiency and reliability of the static correction mode screening, in another embodiment, the static correction mode screening module 30 in the device for screening static correction modes of land seismic data of the present application specifically includes the following contents:

the judging unit 31B is configured to sequentially judge whether the value of the average velocity of the adjacent wells in each synthetic seismic record group is within a preset stratigraphic lithology velocity range;

the second screening unit 32B is configured to determine that the current synthetic seismic record group meets a preset condition if the average velocity values of the adjacent wells in a certain synthetic seismic record group are within a preset stratigraphic lithology velocity range, and determine the static correction mode corresponding to the synthetic seismic record meeting the preset condition as a target static correction mode of the target area.

In order to further improve the screening accuracy of the static correction mode by providing a more accurate data base, in a specific embodiment, the screening device of the static correction mode of the land seismic data further specifically comprises the following preprocessing contents:

a static correction module A01, configured to apply different static correction modes to correct the seismic data of the target area, to obtain each corrected seismic data corresponding to each static correction mode, and

a preprocessing module A02 for unifying the initial measurement point of the plurality of well log data and the initial measurement point of the seismic data within the target area.

The preprocessing module a02 specifically includes the following contents:

(1) and the initial measurement point unifying unit of the well logging information is used for determining the initial measurement point of each drilled well in the target area in the depth domain as the altitude depth according to the ground altitude and the square complement height during drilling in each drilled well logging information in the target area, and adjusting each drilled layered data to the altitude depth so as to unify the initial measurement points of the plurality of drilled well logging information in the target area.

(2) And the initial measuring point unifying unit of the seismic data is used for determining the seismic travel time difference between each drilled well depth initial measuring point and the processing datum plane based on each drilled well seismic data in the target area and the pre-acquired replacing speed, determining the initial measuring point of the time domain as the seismic processing datum plane and unifying the initial measuring points of the plurality of drilled well seismic data in the target area.

Based on the above, the system for screening the static correction mode of the land seismic data further comprises the following contents:

the stratum dividing module 01 is used for carrying out stratum division on each drilled well in the target area by applying the same stratum dividing basis;

and the environment correction module 02 is used for performing environment correction on each well-drilled logging curve in the target area.

From the above description, the screening device for the static correction mode of the land seismic data provided by the embodiment of the application can more accurately and effectively screen the static correction mode of the land seismic data, has a simple and efficient screening process and high sensitivity, can select the most suitable static correction mode of the land seismic data for different target areas, and further can effectively improve the pertinence and accuracy of the static correction of the land seismic data.

The embodiment of the present application further provides a specific implementation manner of an electronic device, which is capable of implementing all steps in the method for screening static correction modes of land seismic data in the foregoing embodiment, and referring to fig. 22, the electronic device specifically includes the following contents:

a processor (processor)601, a memory (memory)602, a communication interface (communications interface)603, and a bus 604;

the processor 601, the memory 602 and the communication interface 603 complete mutual communication through the bus 604; the communication interface 603 is used for realizing information transmission among a screening device of a static correction mode of land seismic data, a client terminal and other participating mechanisms;

the processor 601 is used to call the computer program in the memory 602, and the processor executes the computer program to implement all the steps in the method for screening static correction mode of land seismic data in the above embodiment, for example, the processor executes the computer program to implement the following steps:

As can be seen from the above description, in the electronic device provided in the embodiment of the present invention, by obtaining the synthetic seismic record groups corresponding to the corrected seismic data of the target area, where the static correction manners applied to the corrected seismic data are different, and each synthetic seismic record group includes the calibrated synthetic seismic record of each drilled well in the target area, an accurate and effective data base can be improved for the calculation of the subsequent data, and the efficiency of the calculation and screening of the subsequent data can be effectively improved; by respectively determining the average speed of each drilled single well and the average speed of each adjacent well between each two adjacent drilled single wells corresponding to each synthetic seismic record group and introducing the concept of the average speed of the adjacent well stratum, a simple, accurate and effective data basis can be provided for screening of a subsequent static correction mode, so that the screening result of the subsequent static correction mode is accurate and efficient; the static correction mode corresponding to the synthetic seismic record group meeting the preset conditions is determined as the target static correction mode of the target area according to the values of the single-well average velocity and the adjacent-well average velocity corresponding to the synthetic seismic record groups, the target static correction mode most suitable for the target area can be accurately and effectively screened based on the concept of the adjacent-well stratum average velocity, the screening process is simple, efficient and high in sensitivity, the most suitable static correction mode of land seismic data can be selected according to different target areas, the pertinence and the accuracy of static correction of the land seismic data can be effectively improved, and the accuracy and the reliability of the whole seismic exploration process aiming at the target area can be improved.

Embodiments of the present application also provide a computer-readable storage medium capable of implementing all steps in the method for screening static correction mode of land seismic data in the above embodiments, where the computer-readable storage medium stores a computer program, and the computer program is executed by a processor to implement all steps of the method for screening static correction mode of land seismic data in the above embodiments, for example, the processor implements the following steps when executing the computer program:

As can be seen from the above description, the computer-readable storage medium according to the embodiment of the present invention obtains the synthetic seismic record groups corresponding to the corrected seismic data of the target area, where the static correction manners applied to the corrected seismic data are different, and each synthetic seismic record group includes the calibrated synthetic seismic record of each drilled well in the target area, so as to improve an accurate and effective data base for subsequent data calculation and improve the efficiency of subsequent data calculation and screening; by respectively determining the average speed of each drilled single well and the average speed of each adjacent well between each two adjacent drilled single wells corresponding to each synthetic seismic record group and introducing the concept of the average speed of the adjacent well stratum, a simple, accurate and effective data basis can be provided for screening of a subsequent static correction mode, so that the screening result of the subsequent static correction mode is accurate and efficient; the static correction mode corresponding to the synthetic seismic record group meeting the preset conditions is determined as the target static correction mode of the target area according to the values of the single-well average velocity and the adjacent-well average velocity corresponding to the synthetic seismic record groups, the target static correction mode most suitable for the target area can be accurately and effectively screened based on the concept of the adjacent-well stratum average velocity, the screening process is simple, efficient and high in sensitivity, the most suitable static correction mode of land seismic data can be selected according to different target areas, the pertinence and the accuracy of static correction of the land seismic data can be effectively improved, and the accuracy and the reliability of the whole seismic exploration process aiming at the target area can be improved.

The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the hardware + program class embodiment, since it is substantially similar to the method embodiment, the description is simple, and the relevant points can be referred to the partial description of the method embodiment.

The foregoing description has been directed to specific embodiments of this disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.

Although the present application provides method steps as described in an embodiment or flowchart, additional or fewer steps may be included based on conventional or non-inventive efforts. The order of steps recited in the embodiments is merely one manner of performing the steps in a multitude of orders and does not represent the only order of execution. When an actual apparatus or client product executes, it may execute sequentially or in parallel (e.g., in the context of parallel processors or multi-threaded processing) according to the embodiments or methods shown in the figures.

The systems, devices, modules or units illustrated in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions. One typical implementation device is a computer. In particular, the computer may be, for example, a personal computer, a laptop computer, a vehicle-mounted human-computer interaction device, a cellular telephone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

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.

In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.

The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.

Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.

As will be appreciated by one skilled in the art, embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, embodiments of the present description may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects.

The embodiments of this specification may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The described embodiments may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.

The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment. In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of an embodiment of the specification. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.

The above description is only an example of the present specification, and is not intended to limit the present specification. Various modifications and variations to the embodiments described herein will be apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the embodiments of the present specification should be included in the scope of the claims of the embodiments of the present specification.

Claims

1. A method for screening static correction modes of land seismic data is characterized by comprising the following steps:

2. The method of claim 1, wherein determining the average velocity of each of the drilled single wells and the average velocity of each of the neighboring two drilled wells between the respective single well and the neighboring two drilled wells for each of the synthetic seismic record sets comprises:

3. The method of claim 2, wherein determining the average velocity of adjacent wells between each adjacent two of the drilled wells comprises:

4. The method of claim 3, wherein the average velocity V between two adjacent wells is determined by applying formula one:

5. The method of claim 1, wherein the obtaining a synthetic seismic record set corresponding to each corrected seismic data for the target area comprises:

6. The method of claim 1, wherein determining the static correction mode corresponding to the synthetic seismic record group satisfying a predetermined condition as the target static correction mode for the target area according to the average velocity of the single well and the average velocity of the neighboring well corresponding to the synthetic seismic record group comprises:

7. The method of claim 6, wherein the determining the static correction mode corresponding to the synthetic seismic record group satisfying the predetermined condition as the target static correction mode of the target area according to the average velocity of the single well and the average velocity of the neighboring well corresponding to the synthetic seismic record group further comprises:

8. The method of claim 1, further comprising, prior to the obtaining the synthetic seismic record sets corresponding to the respective corrected seismic data for the target area:

9. The method of claim 8, wherein unifying the initial measurement point of the plurality of well log data and the initial measurement point of the seismic data within the target area comprises:

10. The method of claim 8, wherein prior to applying the different static correction modes to correct the seismic data for the target area, further comprising:

11. A system for static correction screening of land seismic data, comprising:

12. The system of claim 11, wherein the single well and neighboring well average velocity determination module comprises:

13. The system of claim 12, wherein the neighboring average velocity acquisition unit is further configured to:

14. The system of claim 13, wherein the average velocity V between two adjacent wells is determined using the equation one:

15. The system for static correction of land seismic data as recited in claim 11, wherein said synthetic seismic record set acquisition module comprises:

16. The system of claim 11, wherein the static correction mode screening module comprises:

17. The system of claim 16, wherein the static correction mode screening module further comprises:

18. The system for static correction screening of land seismic data as recited in claim 11, further comprising:

19. The system for static correction screening of land seismic data as recited in claim 18, wherein the preprocessing module comprises:

20. The system for static correction screening of land seismic data as recited in claim 18, further comprising:

21. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program performs the steps of the method of static correction screening of land seismic data as claimed in any one of claims 1 to 10.

22. A computer-readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the method for static correction style screening of land seismic data as claimed in any one of claims 1 to 10.