CN112162322A - Synchronous property control method and device of seismic exploration and acquisition system - Google Patents
Synchronous property control method and device of seismic exploration and acquisition system Download PDFInfo
- Publication number
- CN112162322A CN112162322A CN202011043503.9A CN202011043503A CN112162322A CN 112162322 A CN112162322 A CN 112162322A CN 202011043503 A CN202011043503 A CN 202011043503A CN 112162322 A CN112162322 A CN 112162322A
- Authority
- CN
- China
- Prior art keywords
- wire harness
- arrival
- arrival time
- seismic
- distance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 27
- 230000005540 biological transmission Effects 0.000 claims abstract description 49
- 230000005284 excitation Effects 0.000 claims abstract description 20
- 239000000523 sample Substances 0.000 claims abstract description 9
- 238000003908 quality control method Methods 0.000 claims description 33
- 238000004590 computer program Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 19
- 230000002159 abnormal effect Effects 0.000 description 6
- 239000013598 vector Substances 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000007689 inspection Methods 0.000 description 4
- 230000005856 abnormality Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005422 blasting Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 230000029305 taxis Effects 0.000 description 1
- 230000007723 transport mechanism Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3808—Seismic data acquisition, e.g. survey design
Abstract
The embodiment of the application discloses a synchronous property control method and device of a seismic exploration and acquisition system. The method comprises the following steps: obtaining the equivalent speed of the first line beam on each shot line on the transmission distance from the excitation point to the demodulator probe; obtaining a full-offset first-arrival wave transmission model corresponding to the first wire harness according to the equivalent speed of the first wire harness on each gun line; if the collected second wire harness is in the same offset distance domain as the first wire harness, calculating a theoretical value of first arrival time of the second wire harness according to a first arrival wave transmission model of the full offset distance; and executing synchronous property control operation of the seismic exploration and acquisition system according to the theoretical value of the first arrival time and the pre-acquired actual value of the first arrival time of the second wire harness.
Description
Technical Field
The embodiment of the application relates to the field of information processing, in particular to a synchronous property control method and device of a seismic exploration and acquisition system.
Background
Marine seismic exploration has single-ship operation and multi-ship cooperative operation, the source excitation and recording systems under a single ship and a single clock are possibly asynchronous, and the operating systems under multiple ships and multiple clocks are also possibly asynchronous. Therefore, the synchronous quality control of the seismic exploration and acquisition system is the key point of field quality control.
The synchronism quality control is the synchronism of the time of the excitation of the quantitative quality control seismic source ship on the seismic data and the time of the receiving of the instrument ship.
Fig. 1(a) is a diagram illustrating the operation result of quality control synchronization using a time-break signal in the related art. As shown in fig. 1(a), the Time Break (TB) signal check in the auxiliary track indicates that the synchronization between the systems is abnormal, but this method is suitable for the synchronization between the systems under the same clock or synchronous clocks, and cannot control the synchronization between multiple sets of clocks. In addition, in practical application, the auxiliary track TB is normal, which does not mean that the synchronization of the seismic data is always normal, as shown in fig. 1(a), in a certain work area, the TB signal is normal, but the synchronization is abnormal through time slice inspection. Therefore, in the field quality control work, the problems are difficult to find in time by the existing method.
FIG. 1(b) is a schematic diagram showing the operation results of quality control synchronization in the related art using the seismic data volume time slicing technique. As shown in fig. 1(b), if the spatial discontinuity on the time slice is not reasonable, it can be determined that the data synchronization is abnormal. However, a 3D stacked data body with normal seismic data and synchronization abnormal data is needed, a large amount of data needs to be collected, so that the quality control timeliness is poor, when a problem of synchronization abnormality is found, a large amount of data DNP exists, and the experience requirements on quality control personnel are high.
Disclosure of Invention
In order to solve any technical problem, the embodiments of the present application provide a method and an apparatus for controlling synchronization properties of a seismic exploration acquisition system.
To achieve the purpose of the embodiments of the present application, the embodiments of the present application provide a method for controlling the synchronization property of a seismic exploration acquisition system, including:
obtaining the equivalent speed of the first line beam on each shot line on the transmission distance from the excitation point to the demodulator probe;
obtaining a full-offset first-arrival wave transmission model corresponding to the first wire harness according to the equivalent speed of the first wire harness on each gun line;
if the collected second wire harness is in the same offset distance domain as the first wire harness, calculating a theoretical value of first arrival time of the second wire harness according to a first arrival wave transmission model of the full offset distance;
and executing synchronous property control operation of the seismic exploration and acquisition system according to the theoretical value of the first arrival time and the pre-acquired actual value of the first arrival time of the second wire harness.
A synchronous property control device for a seismic survey acquisition system, comprising:
the acquisition module is used for acquiring the equivalent speed of the first line beam on each shot line in the transmission distance from the excitation point to the demodulator probe;
the establishing module is used for obtaining a first arrival wave transmission model of the full offset distance corresponding to the first wire harness according to the equivalent speed of the first wire harness on each gun line;
the calculation module is used for calculating a theoretical value of the first arrival time of the second wire harness according to the first arrival wave transmission model of the full offset distance if the acquired second wire harness and the first wire harness are in the same offset distance domain;
and the quality control module is used for executing synchronous quality control operation of the seismic exploration and acquisition system according to the theoretical value of the first arrival time and the pre-acquired actual value of the first arrival time of the second wire harness.
A storage medium having a computer program stored therein, wherein the computer program is arranged to perform the method as described above when executed.
An electronic device comprising a memory having a computer program stored therein and a processor arranged to execute the computer program to perform the method as described above.
One of the above technical solutions has the following advantages or beneficial effects:
the method comprises the steps of obtaining the equivalent speed of a first wire harness on each shot line in the transmission distance from an excitation point to a demodulator probe, obtaining a first arrival wave transmission model of a full offset distance corresponding to the first wire harness according to the equivalent speed of the first wire harness on each shot line, calculating a theoretical value of first arrival time of the second wire harness according to the first arrival wave transmission model of the full offset distance if the collected second wire harness and the first wire harness are in the same offset distance domain, executing synchronous property control operation of a seismic exploration collection system according to the theoretical value of the first arrival time and a pre-obtained actual value of the first arrival time of the second wire harness, achieving the purpose of quantitatively and qualitatively controlling the synchronism among systems, and improving the accuracy and the efficiency of quality control.
Additional features and advantages of the embodiments of the application will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by the practice of the embodiments of the application. The objectives and other advantages of the embodiments of the application may be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The accompanying drawings are included to provide a further understanding of the embodiments of the present application and are incorporated in and constitute a part of this specification, illustrate embodiments of the present application and together with the examples of the embodiments of the present application do not constitute a limitation of the embodiments of the present application.
FIG. 1(a) is a diagram illustrating the operation result of using a time-break signal to control the synchronization of quality control in the related art;
FIG. 1(b) is a schematic diagram showing the operation results of the quality control synchronization of the related art using the seismic data volume time slicing technique;
FIG. 2 is a flow chart of a method for controlling the synchronization properties of a seismic survey acquisition system according to an embodiment of the present application;
FIG. 3 is a schematic illustration of first-arrival wave propagation in seismic exploration for an ocean bottom cable according to embodiments of the present application;
fig. 4(a) is a schematic diagram of an observation system of a work area 8S4L according to an embodiment of the present disclosure;
FIG. 4(b) is a schematic diagram of the scanning results of the equivalent velocity provided by the embodiments of the present application;
FIG. 5(a) is a schematic diagram of the principle of snell's law provided by an embodiment of the present application;
FIG. 5(b) is a vector diagram of direct and taxis waves provided by an embodiment of the present application;
FIG. 5(c) is a diagram of the scanning result of the time constant provided by the embodiment of the present application;
FIG. 6 is a schematic diagram of theoretical fitting positions of a first-break curve provided in an embodiment of the present application;
fig. 7 is a schematic diagram of a time synchronization quality control result according to an embodiment of the present disclosure;
FIG. 8 is a schematic diagram of a synchronous nature control device of a seismic survey acquisition system according to an embodiment of the present disclosure.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application more apparent, the embodiments of the present application will be described in detail below with reference to the accompanying drawings. It should be noted that, in the embodiments of the present application, features in the embodiments and the examples may be arbitrarily combined with each other without conflict.
In the process of implementing the present application, the inventor conducts technical analysis on the related art, and finds that the related art has at least the following problems, including:
the existing quality control synchronization technology is time-break signal inspection and time slicing inspection. Practice proves that the auxiliary track is normally interrupted, and the data synchronism cannot be guaranteed to be normal, because the system synchronism under the same clock or each synchronous clock is assumed. The time slice inspection can control data to be abnormal synchronously, but needs a three-dimensional stacked data body of seismic data, and needs a period of time for acquiring the seismic data of a continuous slice and a period of time for generating the stacked data body, so that the quality control timeliness is poor, and the requirement of rapid quality control in the field of seismic exploration cannot be met.
The method provided by the embodiment of the application aims to quickly and quantitatively control the synchronism under various marine seismic exploration modes during field seismic exploration, and comprises the synchronism among seismic exploration systems under a single-ship single clock and multiple-ship multiple clocks.
Fig. 2 is a flowchart of a method for controlling a synchronization property of a seismic exploration acquisition system according to an embodiment of the present disclosure. As shown in fig. 2, the method shown in fig. 2 includes:
in one exemplary embodiment, the first beam may be a beam of seismic source excitation, such as the first beam.
By adopting the first wire harness, the influence of other wire harnesses on the accuracy of the acquisition result can be reduced, and the accuracy of the acquired data is ensured.
the inventors have found that the equivalent velocity field for each line bundle is the same in the common offset domain. That is, the equivalent velocity fields corresponding to different offsets in the first arrival wave transmission model are applicable to other beams in the common offset domain.
May be the equivalent velocity field of the first beam. In the practical application of the method, the material is,
in an exemplary embodiment, based on the same equivalent velocity field of each line beam in the common offset distance domain, seismic trace data according to a second acquired line beam can be determined, equivalent velocity information of the second line beam at a corresponding offset distance is determined, and calculation of a theoretical value is completed based on the equivalent velocity information.
204, executing synchronous property control operation of the seismic exploration and acquisition system according to the theoretical value of the first arrival time and the pre-acquired actual value of the first arrival time of the second wire harness;
the aim of quantitatively controlling the synchronism among the systems can be achieved by comparing the ratio of the theoretical value with the actual value.
The method comprises the steps of obtaining the equivalent speed of a first wire harness on each gun line in the transmission distance from an excitation point to a demodulator probe, obtaining a first arrival wave transmission model of a full offset distance corresponding to the first wire harness according to the equivalent speed of the first wire harness on each gun line, calculating a theoretical value of first arrival time of the second wire harness according to the first arrival wave transmission model of the full offset distance if the collected second wire harness and the first wire harness are in the same offset distance domain, executing synchronous property control operation of a seismic exploration and collection system according to the theoretical value of the first arrival time and a pre-obtained actual value of the first arrival time of the second wire harness, achieving the purpose of quantitatively and qualitatively controlling the synchronism among the systems, and improving the accuracy and efficiency of quality control.
The method provided by the embodiments of the present application is explained as follows:
in an exemplary embodiment, the full offset first arrival wave transmission model includes:
T＝(X/V)+C；
wherein T is the first-arrival wave transmission time, X is the equivalent propagation distance, C is the time constant, and V is the first-arrival wave propagation equivalent velocity.
FIG. 3 is a schematic diagram of first-arrival wave propagation in seismic exploration for an ocean bottom cable according to an embodiment of the present application. As shown in fig. 3, the first beam has a propagation distance X from the excitation point to the detection point.
According to the pythagorean theorem, the first arrival propagation distance X is obtained by the following method, including:
X＝SQRT{(Hr-Hs)2+offset2}；
wherein Hr represents the water depth of the position of the wave detection point; hs represents the depth of water at the excitation point, and offset represents the distance between the pickup point and the excitation point in the horizontal direction.
The full offset first arrival wave transmission model is obtained by the following method:
if the direct wave is adopted for calculation, setting the equivalent speed as a fixed value, and acquiring the numerical value of the time constant under each offset distance;
if the refracted wave is adopted for calculation, the time constant is set to be zero, and the value of the equivalent speed at each offset distance is obtained.
When a direct wave model is adopted for picking, V is 1500m/S, and C value needs to be scanned; when a refracted wave model is adopted for pickup, the V value needs to be scanned, and the C value is 0.
The scanning of the C value of the direct wave model is explained as follows:
as the offset distance increases, the refracted wave component gradually mixes into the first-motion wave, and the velocity significantly changes. Therefore, when a direct wave model is adopted, a constant speed parameter can introduce a constant error. Under a stable observation system, the error is unchanged.
Fig. 4(a) is a schematic diagram of an observation system of a work area 8S4L according to an embodiment of the present application. The observation system has 16 offset conditions on a single section, constant errors C under different conditions are scanned by a computer, and specific scanning results are shown in fig. 4 (b).
Description will be given of the scanning of the velocity values of the refracted wave model:
fig. 5(a) is a schematic diagram of snell's law provided in an embodiment of the present application. As shown in fig. 5(a), when the incident angle is larger than the critical angle, the refracted wave will propagate along the reflective interface.
In actual acquisition construction, the real speed of seismic waves propagating in a sea water layer and a seabed stratum cannot be mastered, and the size of a critical angle and the real position of a refraction point cannot be confirmed. In view of the above problems, the inventors propose that considering a refracted wave as two components, a direct wave and a gliding wave, the vector sum of the direct wave and the gliding wave is equal to the equivalent propagation distance X, regardless of the position of the refracted point.
FIG. 5(b) is a vector diagram of direct waves and gliding waves provided by an embodiment of the present application. As shown in fig. 5(b), the vector a may be expressed as a direct wave, the vector b may be expressed as a gliding wave, and the direction in which the sum of the two vectors is located is the propagation distance X direction. Therefore, in the direction of the equivalent propagation distance X, an equivalent velocity V' is necessarily present, and the negative jump-off time of the refracted wave on the seismic channel can be fitted.
The equivalent velocity can be obtained by scanning, and is used for eliminating constant errors generated by a direct wave model.
Under the same observation system model, 16 offset conditions are shared on a single section, and the scanning result is shown in fig. 5(c) specifically by scanning the equivalent velocity V' under different conditions through a computer.
In an exemplary embodiment, said calculating a theoretical value of first arrival time of said second beam according to said first arrival wave propagation model at full offset distance comprises:
acquiring a single-channel section result of the seismic data of the second wire harness;
determining a critical location in the profile result below a first negative jump and above a first positive jump;
determining a target equivalent speed corresponding to the critical position from the first-motion wave transmission model of the full offset distance;
and calculating a theoretical value of the first arrival time of the second wire harness based on the target equivalent speed information.
Fitting can be performed by using a near-path profile, and the purpose of quality control synchronism is achieved. The method comprises the steps of drawing the seismic data after observation into a single-channel section, and fitting at a critical position below a first negative jump and above a first positive jump.
Fig. 6 is a schematic diagram of theoretical fitting positions of a first-break curve according to an embodiment of the present disclosure. As shown in fig. 6, the dotted line in the left graph indicates the fitting position of the synchronization curve, and the dotted line in the right graph indicates the fitting position of the curve in which the synchronization abnormality exists.
In an exemplary embodiment, the performing a synchronous property control operation of the seismic survey acquisition system according to the theoretical value of the first arrival time and the pre-acquired actual value of the first arrival time of the second line beam comprises:
comparing the actual value of the first arrival time in each seismic channel of the second wire harness with the theoretical value of the first arrival time;
determining the number of seismic channels of which the actual value conforming to the first arrival time is greater than the theoretical value of the first arrival time;
judging whether the proportion between the number of the earthquakes and the total number of the seismic channels meets a preset synchronous judgment condition or not;
and if the judgment result is that the synchronization judgment condition is met, determining that the synchronization quality control of the seismic exploration and acquisition system is passed.
Fig. 7 is a schematic diagram of a time synchronization quality control result according to an embodiment of the present disclosure. As shown in fig. 7, the area selected by the box in fig. 7 is the theoretical fitting position. Under the same work area, the time synchronization normal first arrival fitting curve should be always kept on the theoretical fitting position. The left side is the fitting condition of a single section of a certain gun line which is synchronous and normal, and the right side is the fitting condition of a single section of a certain abnormal running line.
The purpose of quantitative quality control of statistical properties is realized by adopting a first-break wave automatic pickup technology. Due to the existence of observation errors, each path cannot accord with the theoretical position, but more than 90% of the first-motion wave positions of the seismic channels can be fitted. And if the first arrival position of more than 20% of the seismic traces is not matched with the fitting position, judging that the data possibly has a synchronization problem.
And when the project is started, recording the equivalent speed or the empirical constant C of each blasting line of the first wire harness corresponding to each cable single-channel section. And under the condition that the observation system is stable and unchanged, the equivalent speed or the empirical constant C of the first wire harness is taken as a standard to control the data of each wire harness collected at the back. And after the abnormality is found, further checking in time and taking corresponding measures.
As can be seen from the above, in the method provided in the embodiment of the present application, a full-offset first-arrival wave transmission model is established according to the equivalent speed of each gun line of the first wire harness corresponding to each cable single-path section, and under the condition that the observation system is relatively stable, the equivalent speed in the full-offset first-arrival wave transmission model is used for deduction, so as to control the synchronism of the acquired data of each wire harness. The full-offset first-arrival wave transmission model is constructed by collecting prior information such as offset, cable depth, air gun depth and equivalent speed, theoretical first-arrival curve fitting can be carried out on collected data, the fitting result and the actual first arrival are subjected to difference finding, and the synchronism among the systems is quantitatively controlled through a threshold value of model precision.
Through the test and application of a plurality of marine seismic exploration projects, the quantitative and rapid quality control of the synchronism under the condition of multiple clocks is realized. In addition, the method can quickly and quantitatively control the synchronism of various marine seismic exploration modes, including the synchronism of seismic exploration systems under a single ship and a single clock and multiple ships and multiple clocks, and the quality control result can be given after half an hour of a line.
Compared with the prior art, the method provided by the embodiment of the application can be used for controlling the synchronism of a single ship and a single clock and controlling the synchronism of each seismic exploration system under multiple clocks of the ship, and has the advantages of short time consumption and no need of manual experience.
FIG. 8 is a schematic diagram of a synchronous nature control device of a seismic survey acquisition system according to an embodiment of the present disclosure. As shown in fig. 8, the apparatus shown in fig. 8 includes:
the acquisition module is used for acquiring the equivalent speed of the first line beam on each shot line in the transmission distance from the excitation point to the demodulator probe;
the establishing module is used for obtaining a first arrival wave transmission model of the full offset distance corresponding to the first wire harness according to the equivalent speed of the first wire harness on each gun line;
the calculation module is used for calculating a theoretical value of the first arrival time of the second wire harness according to the first arrival wave transmission model of the full offset distance if the acquired second wire harness and the first wire harness are in the same offset distance domain;
and the quality control module is used for executing synchronous quality control operation of the seismic exploration and acquisition system according to the theoretical value of the first arrival time and the pre-acquired actual value of the first arrival time of the second wire harness.
In an exemplary embodiment, the full offset first arrival wave transmission model includes:
T＝(X/V)+C；
wherein T is the first-arrival wave transmission time, X is the equivalent propagation distance, C is the time constant, and V is the first-arrival wave propagation equivalent velocity.
In an exemplary embodiment, the transmission distance is calculated by:
X＝SQRT{(Hr-Hs)2+offset2}
wherein Hr represents the water depth of the position of the wave detection point; hs represents the depth of water at the excitation point, and offset represents the distance between the pickup point and the excitation point in the horizontal direction.
In an exemplary embodiment, the full offset first arrival wave transmission model is obtained by:
if the direct wave is adopted for calculation, setting the equivalent speed as a fixed value, and acquiring the numerical value of the time constant under each offset distance;
if the refracted wave is adopted for calculation, the time constant is set to be zero, and the value of the equivalent speed at each offset distance is obtained.
In an exemplary embodiment, the calculation module is specifically configured to obtain a single-pass section result of the seismic data of the second beam; determining a critical location in the profile result below a first negative jump and above a first positive jump; determining a target equivalent speed corresponding to the critical position from the first-motion wave transmission model of the full offset distance; and calculating a theoretical value of the first arrival time of the second wire harness based on the target equivalent speed information.
In an exemplary embodiment, the quality control module is specifically configured to compare an actual value of a first arrival time in each seismic channel of the second line bundle with a theoretical value of the first arrival time; determining the number of seismic channels of which the actual value conforming to the first arrival time is greater than the theoretical value of the first arrival time; judging whether the proportion between the number of the earthquakes and the total number of the seismic channels meets a preset synchronous judgment condition or not; and if the judgment result is that the synchronization judgment condition is met, determining that the synchronization quality control of the seismic exploration and acquisition system is passed.
The device provided by the embodiment of the application obtains the equivalent speed of a first wire harness on each gun line on the transmission distance from an excitation point to a demodulator probe, obtains a first arrival wave transmission model of a full offset distance corresponding to the first wire harness according to the equivalent speed of the first wire harness on each gun line, calculates a theoretical value of first arrival time of the second wire harness according to the first arrival wave transmission model of the full offset distance if the collected second wire harness and the first wire harness are in the same offset distance domain, executes synchronous property control operation of a seismic exploration and collection system according to the theoretical value of the first arrival time and a pre-obtained actual value of the first arrival time of the second wire harness, achieves the purpose of quantitatively and qualitatively controlling the synchronism among systems, and improves the accuracy and efficiency of quality control.
An embodiment of the present application provides a storage medium, in which a computer program is stored, wherein the computer program is configured to perform the method described in any one of the above when the computer program runs.
An embodiment of the application provides an electronic device, comprising a memory and a processor, wherein the memory stores a computer program, and the processor is configured to execute the computer program to perform the method described in any one of the above.
It will be understood by those of ordinary skill in the art that all or some of the steps of the methods, systems, functional modules/units in the devices disclosed above may be implemented as software, firmware, hardware, and suitable combinations thereof. In a hardware implementation, the division between functional modules/units mentioned in the above description does not necessarily correspond to the division of physical components; for example, one physical component may have multiple functions, or one function or step may be performed by several physical components in cooperation. Some or all of the components may be implemented as software executed by a processor, such as a digital signal processor or microprocessor, or as hardware, or as an integrated circuit, such as an application specific integrated circuit. Such software may be distributed on computer readable media, which may include computer storage media (or non-transitory media) and communication media (or transitory media). The term computer storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data, as is well known to those of ordinary skill in the art. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, Digital Versatile Disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can accessed by a computer. In addition, communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media as known to those skilled in the art.
Claims (10)
1. A method of synchronized property control of a seismic survey acquisition system, comprising:
obtaining the equivalent speed of the first line beam on each shot line on the transmission distance from the excitation point to the demodulator probe;
obtaining a full-offset first-arrival wave transmission model corresponding to the first wire harness according to the equivalent speed of the first wire harness on each gun line;
if the collected second wire harness is in the same offset distance domain as the first wire harness, calculating a theoretical value of first arrival time of the second wire harness according to a first arrival wave transmission model of the full offset distance;
and executing synchronous property control operation of the seismic exploration and acquisition system according to the theoretical value of the first arrival time and the pre-acquired actual value of the first arrival time of the second wire harness.
2. The method of claim 1, wherein the full offset first arrival wave transmission model comprises:
T＝(X/V)+C；
wherein T is the first-arrival wave transmission time, X is the equivalent propagation distance, C is the time constant, and V is the first-arrival wave propagation equivalent velocity.
3. The method of claim 1, wherein the transmission distance is calculated by:
X＝SQRT{(Hr-Hs)2+offset2}；
wherein Hr represents the water depth of the position of the wave detection point; hs represents the depth of water at the excitation point, and offset represents the distance between the pickup point and the excitation point in the horizontal direction.
4. The method according to claim 2 or 3, wherein the full offset first arrival wave transmission model is obtained by:
if the direct wave is adopted for calculation, setting the equivalent speed as a fixed value, and acquiring the numerical value of the time constant under each offset distance;
if the refracted wave is adopted for calculation, the time constant is set to be zero, and the value of the equivalent speed at each offset distance is obtained.
5. The method of claim 1, wherein said calculating a theoretical value of first arrival time of said second beam based on said model of full offset first arrival wave propagation comprises:
acquiring a single-channel section result of the seismic data of the second wire harness;
determining a critical location in the profile result below a first negative jump and above a first positive jump;
determining a target equivalent speed corresponding to the critical position from the first-motion wave transmission model of the full offset distance;
and calculating a theoretical value of the first arrival time of the second wire harness based on the target equivalent speed information.
6. The method of claim 1, wherein performing a simultaneous property control operation of the seismic survey acquisition system based on the theoretical value of the first arrival time and the pre-acquired actual value of the first arrival time of the second line beam comprises:
comparing the actual value of the first arrival time in each seismic channel of the second wire harness with the theoretical value of the first arrival time;
determining the number of seismic channels of which the actual value conforming to the first arrival time is greater than the theoretical value of the first arrival time;
judging whether the proportion between the number of the earthquakes and the total number of the seismic channels meets a preset synchronous judgment condition or not;
and if the judgment result is that the synchronization judgment condition is met, determining that the synchronization quality control of the seismic exploration and acquisition system is passed.
7. A synchronous property control device for a seismic survey acquisition system, comprising:
the acquisition module is used for acquiring the equivalent speed of the first line beam on each shot line in the transmission distance from the excitation point to the demodulator probe;
the establishing module is used for obtaining a first arrival wave transmission model of the full offset distance corresponding to the first wire harness according to the equivalent speed of the first wire harness on each gun line;
the calculation module is used for calculating a theoretical value of the first arrival time of the second wire harness according to the first arrival wave transmission model of the full offset distance if the acquired second wire harness and the first wire harness are in the same offset distance domain;
and the quality control module is used for executing synchronous quality control operation of the seismic exploration and acquisition system according to the theoretical value of the first arrival time and the pre-acquired actual value of the first arrival time of the second wire harness.
8. The apparatus of claim 7, wherein the full offset first arrival wave propagation model comprises:
T＝(X/V)+C；
wherein T is the first-arrival wave transmission time, X is the equivalent propagation distance, C is the time constant, and V is the first-arrival wave propagation equivalent velocity.
9. A storage medium, in which a computer program is stored, wherein the computer program is arranged to perform the method of any of claims 1 to 6 when executed.
10. An electronic device comprising a memory and a processor, wherein the memory has stored therein a computer program, and wherein the processor is arranged to execute the computer program to perform the method of any of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011043503.9A CN112162322A (en) | 2020-09-28 | 2020-09-28 | Synchronous property control method and device of seismic exploration and acquisition system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011043503.9A CN112162322A (en) | 2020-09-28 | 2020-09-28 | Synchronous property control method and device of seismic exploration and acquisition system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112162322A true CN112162322A (en) | 2021-01-01 |
Family
ID=73861870
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011043503.9A Pending CN112162322A (en) | 2020-09-28 | 2020-09-28 | Synchronous property control method and device of seismic exploration and acquisition system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112162322A (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104133245A (en) * | 2014-07-22 | 2014-11-05 | 中国石油天然气集团公司 | Seismic data static correction method and system |
| US20170068008A1 (en) * | 2015-09-04 | 2017-03-09 | Saudi Arabian Oil Company | Automatic quality control of seismic travel time |
| CN106855640A (en) * | 2017-02-03 | 2017-06-16 | 中国石油天然气集团公司 | A kind of submarine cable geological data first arrival computational methods and device |
| CN108445533A (en) * | 2018-03-08 | 2018-08-24 | 青岛海洋地质研究所 | The method of long offset structure phytal zone rate pattern based on OBS |
| CN109212598A (en) * | 2018-10-31 | 2019-01-15 | 中国石油化工股份有限公司 | Three-dimensional space secondary localization method based on direct wave inverting |
| CN110858003A (en) * | 2018-08-22 | 2020-03-03 | 中国石油化工股份有限公司 | Non-earth surface consistency static correction method based on equivalent speed |
| CN110954941A (en) * | 2018-09-26 | 2020-04-03 | 中国石油化工股份有限公司 | Automatic first arrival picking method and system |
-
2020
- 2020-09-28 CN CN202011043503.9A patent/CN112162322A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104133245A (en) * | 2014-07-22 | 2014-11-05 | 中国石油天然气集团公司 | Seismic data static correction method and system |
| US20170068008A1 (en) * | 2015-09-04 | 2017-03-09 | Saudi Arabian Oil Company | Automatic quality control of seismic travel time |
| CN106855640A (en) * | 2017-02-03 | 2017-06-16 | 中国石油天然气集团公司 | A kind of submarine cable geological data first arrival computational methods and device |
| CN108445533A (en) * | 2018-03-08 | 2018-08-24 | 青岛海洋地质研究所 | The method of long offset structure phytal zone rate pattern based on OBS |
| CN110858003A (en) * | 2018-08-22 | 2020-03-03 | 中国石油化工股份有限公司 | Non-earth surface consistency static correction method based on equivalent speed |
| CN110954941A (en) * | 2018-09-26 | 2020-04-03 | 中国石油化工股份有限公司 | Automatic first arrival picking method and system |
| CN109212598A (en) * | 2018-10-31 | 2019-01-15 | 中国石油化工股份有限公司 | Three-dimensional space secondary localization method based on direct wave inverting |
Non-Patent Citations (4)
| Title |
|---|
| 孙哲等: "地震采集异常道实时监控技术与应用", 《石油管材与仪器》 * |
| 孙哲等: "地震采集异常道实时监控技术与应用", 《石油管材与仪器》, vol. 3, no. 2, 15 April 2017 (2017-04-15), pages 67 - 72 * |
| 张建峰: "基于初至波拟合地震采集同步方法研究", 《工程地球物理学报》 * |
| 张建峰: "基于初至波拟合地震采集同步方法研究", 《工程地球物理学报》, vol. 16, no. 4, 30 July 2019 (2019-07-30), pages 454 - 459 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109709574B (en) | Seabed microtopography laser scanning imaging system and three-dimensional terrain reconstruction method | |
| AU2018248022A1 (en) | Automated draft survey | |
| CN112505750B (en) | Deep-towed multi-channel seismic streamer attitude determination method and processing terminal | |
| US10520631B2 (en) | Magnetic field measurement via streamer cables | |
| GB2479299A (en) | Processing seismic data | |
| US20070223308A1 (en) | Methods of range selection for positioning marine seismic equipment | |
| US11493595B2 (en) | Method and apparatus for evaluating accuracy in positioning a receiver point | |
| CN110490217A (en) | For improving the method and system of object detection and object classification | |
| CN107504917B (en) | Three-dimensional size measuring method and device | |
| US20120093339A1 (en) | 3d soundscaping | |
| WO2020150868A1 (en) | Method and apparatus for manufacturing line simulation | |
| US20140320653A1 (en) | System and method for automated bottom feature navigator | |
| CN109782355B (en) | OBS detection point drift detection method and device | |
| CN108983200A (en) | Towing line array array shape estimation method and system | |
| US10317367B2 (en) | Eddy-current flaw detector and eddy-current flaw detection method | |
| CN109212598B (en) | Three-dimensional space secondary positioning method based on direct wave inversion | |
| CN112162322A (en) | Synchronous property control method and device of seismic exploration and acquisition system | |
| US8477560B2 (en) | Determining a position of an object of a positioning network associated with a marine survey arrangement | |
| CN114779211A (en) | Laser pulse radar equipment and point cloud density improving method and equipment | |
| CN112902931B (en) | Method for measuring and eliminating delay between depth measurement data and positioning data of unmanned ship | |
| CN112180434B (en) | Time consistency checking method, device and machine-readable storage medium | |
| CN109917414B (en) | Ship freeboard measuring method and system based on laser technology | |
| CN110794479B (en) | Method and device for inspecting observation system based on near-path superposition | |
| CN111965729B (en) | Real-time monitoring method, system and device for vibroseis combination center | |
| CN110824559B (en) | Method and device for generating co-coordinate point gather and computer-readable storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210101 |