This disclosure is related generally to the field of geophysical surveying, and more specifically to acquiring ghost-free or near-ghost free data. The disclosure may have applications in the field of broadband source calibration. The disclosure may have benefit when applied to marine surveying operations.
A signature of a seismic source (e.g., an air gun or an air gun array) can be modeled using conventional approaches. Such modeling can be used to indentify noise components when the seismic source is used to perform a seismic survey. However, the advent of broadband marine seismic solutions unveils limitations in conventional modeling approaches outside the marine seismic bandwidth.
While a farfield signature of an air gun array can be modeled using approaches based on Rayleigh's theory of oscillating bubbles, broadband marine seismic solutions may unveil limitations in modeling algorithms outside marine seismic bandwidth. There are several reasons for this that can stem from implementations of source modeling procedures.
Conventional approaches to source modeling may fail to take into account higher frequencies. For example, in many cases the Digital Field System V (DFS-V) Out-128 (72) Hz (decibel/octave) filter may be used for filtering to compare measured and modeled signatures of an air gun. As a consequence, models calibrated according to this approach may not perform outside the frequency range of the DFS-V filter. Measured signatures that were used for calibration could include a source ghost, which could mask frequency ranges depending on a depth of a shot from a seismic source (e.g., an air gun). Within these masking limits, modeling schemes may be devised that could match array signatures well enough to be used in de-signature applications.
As interest in wider frequency ranges grows, the previously masked information can become useful. Modeled signatures that depend on filter systems with a defined calibration range can become less reliable outside the calibration range. Approaches to broaden the applicable frequency range beyond conventional calibration ranges include both calibrating over a broader bandwidth, as well as improving the underlying model such that deviations from calibrated values can be handled.
Marine seismic air gun modeling can be, for example, based on the dynamics of an oscillating air bubble created by an underwater explosion. Some energy loss during the bubble oscillation can be due to mass transfer inside the bubble. As the bubble oscillates, water may evaporate at the bubble wall and condense inside the bubble to produce a transfer of mass of water into the bubble, causing the dampening exhibited in bubble pulses of an air gun signature. The exact rate at which these processes work can be difficult to describe theoretically due to strong variations of temperature and pressure inside the air bubble. As a result, observed rates can be used as calibration in order to make modeled signatures fit recorded data. This may be sufficient for a good fit with a specific (narrow) bandwidth specification; however, it may not be accurate enough for modeling of broadband signatures.
Other physical effects may contribute to differences between measured and modeled signatures. Such effects can be part of the source modeling theory. This can include a realistic model of air escaping from a seismic source (e.g., an air gun) to describe the shape of a primary pulse, and upward movement of the bubble due to buoyancy. As longer time signatures are used for processing, effects of upward movement of the bubble due to buoyancy can be increasingly important.
There may also be high-frequency effects that are not directly reproducible, such as cavitation. While such noise may not be especially relevant in a broadband frequency range, it may be beneficial for purposes such as modeling environmental effects.
In addition to the dynamics of each air bubble, seismic sources (e.g., air guns) may also interact with each other. While some interaction between air guns can be known, the case of clusters, for example, where air bubbles of two or more air guns coalesce, may be less clear. Two-gun clusters have been researched and can be used in array designs. However, in searching for low frequencies, hyperclusters, which are beyond the scope of both the original model and the original calibration of existing modeling algorithms, may be used.
Therefore, higher accuracy in modeling seismic source signatures is desired.
So that the manner in which the above-recited features of the present invention can be understood in detail, a more particular description of the invention, briefly summarized above, may be had by reference to embodiments, some of which are illustrated in the appended drawings. It is to be noted, however, that the appended drawings illustrate only typical embodiments of this invention and are therefore not to be considered limiting of its scope, for the invention may admit to other equally effective embodiments.
To facilitate understanding, identical reference numerals have been used, where possible, to designate identical elements that are common to the figures. It is contemplated that elements disclosed in one embodiment may be beneficially utilized on other embodiments without specific recitation.
It is to be understood the present disclosure is not limited to particular devices or methods, which may, of course, vary. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting. As used herein, the singular forms “a”, “an”, and “the” include singular and plural referents unless the context clearly dictates otherwise. Furthermore, the word “may” is used throughout this application in a permissive sense (i.e., having the potential to, being able to), not in a mandatory sense (i.e., must). The term “include,” and derivations thereof, mean “including, but not limited to.” The term “coupled” means directly or indirectly connected, either physically or operatively.
This disclosure is related generally to the field of marine surveying. Marine surveying can include, for example, seismic and/or electromagnetic surveying, among others. For example, this disclosure may have applications in marine surveying, in which one or more sources are used to generate wave-fields, and sensors—either towed or ocean bottom—receive energy generated by the sources and affected by the interaction with the subsurface formation.
As would be understood by one of ordinary skill in the art with the benefit of this disclosure, embodiments and advantages described herein with reference to a “streamer” might be equally achievable and advantageous when used with a marine geophysical prospecting ocean bottom cable. Consequently, as used herein, “streamer” should be read to refer equally to a towed sensor cable as well as to an ocean bottom sensor cable.
Embodiments of the present disclosure can include methods, apparatus, and machine-readable media for acquiring ghost-free or near-ghost-free signatures for broadband seismic source calibration. Some embodiments can include a non-transitory machine-readable medium storing instructions executable by a processing resource to cause a machine to determine positions of a plurality of seismic receivers (e.g., hydrophones) to enable removal of an interference effect of ghost signals originating from mirrored versions of at least one seismic source (e.g., air gun or acoustic oscillator) that are received at the plurality of seismic receivers, acquire data from the plurality of seismic receivers located at the determined positions by operating the at least one seismic source, and prepare a calibrated model of the at least one seismic source based on the acquired data.
Some embodiments of the present disclosure can include a method for determining positions of a plurality of seismic receivers (e.g., hydrophones) to enable removal of an interference effect of ghost signals originating from mirrored versions of at least one seismic source (e.g., air gun or acoustic oscillator) that are received at the plurality of seismic receivers, acquiring data from the plurality of seismic receivers located at the determined positions by operating the at least one seismic source, and preparing a calibrated model of the at least one seismic source based on the acquired data.
Some embodiments of the present disclosure can include a marine seismic source calibration assembly comprising: a first sensor cable and a second sensor cable, at least one seismic source (e.g., air gun or acoustic oscillator) attached to the first sensor cable, and an array of seismic receivers (e.g., array of hydrophones) attached to the second sensor cable. According to embodiments of the present disclosure, positions of the seismic receivers (e.g., hydrophones) in the array can be (pre-)determined to enable removal of a predicted interference effect of ghost signals originating from mirrored versions of the at least one seismic source that are received at the array of seismic receivers (e.g., array of hydrophones).
According to embodiments of the present disclosure, a seismic source can be an air gun or an acoustic oscillator. In some embodiments, an array of seismic sources can comprise an array of air guns or an array or acoustic oscillators. According to embodiments of the present disclosure, a seismic receiver can be a hydrophone. In some embodiments, an array of seismic receivers can comprise an array of hydrophones.
Some embodiments of the present disclosure include an acquisition setup which can allow for de-ghosting measured seismic data, and using the de-ghosted seismic data as a basis for calibration of a seismic source model. In some embodiments, de-ghosting can be achieved by applying a de-ghosting operator, for example, the Posthumus de-ghosting operator. The Posthumus de-ghosting operator can be utilized with data that is acquired from two or more seismic receivers which have complimentary ghost wave functions in the frequency domain (e.g., complementary amplitude responses in the frequency domain). In some embodiments, three or more seismic receivers can be used to promote the ghost wave functions as being complimentary within a typical frequency range of seismic data (e.g., the frequency range of 0-250 hertz (Hz)).
For example, the positions of the plurality of seismic receivers may be adjusted based on avoiding overlapping of ghost notches of ghost wave functions of different source-receiver pairs, since ghost notches may be directly related to strongest interference effects from the ghost wave signals (e.g., the signals 108, 110 in
According to embodiments of the present disclosure, an amplitude response in the frequency domain of a ghost wave function can be dependent on the travel time difference and the travel time ratio between a direct pressure wave from a seismic source to a seismic receiver (e.g., a pressure wave 112 or 114 in
If a plurality of seismic receivers are configured in a vertical array (i.e., the plurality of seismic receivers have an identical horizontal distance from the seismic source, as illustrated in the acquisition setup 100 in
As illustrated in
Given a horizontal distance between a seismic source and a seismic receiver, (vertical) depths of a plurality of seismic receivers can be selected such that the aforementioned travel time differences result in complementary ghost wave functions, as illustrated in an example graph 300 in
As illustrated in
According to embodiments of the present disclosure, It can be beneficial to determine a setup of seismic receivers (i.e., positions of a plurality of seismic receivers for a given position of a seismic source) that maximizes a minimum value (e.g., a minimum value across a range of frequencies associated with seismic data) of an envelope function of ghost wave functions associated with source-receiver pairs by tuning (adjusting) positions of the plurality of seismic receivers. By maximizing the minimum value of the envelope function based on tuning positions of the plurality of seismic receivers, the ghost wave functions associated with all source-receiver pairs may have complementary waveforms, as discussed above.
According to embodiments of the present disclosure, based on maximizing a cost function (e.g., a function related to an envelope function of ghost wave functions), positions of a plurality of seismic receivers may be determined in an effort to enable efficient removal of interference effects of ghost data. In an embodiment of the present disclosure, the cost function may be given by:
where rj is a vector (e.g., in 3-D space) from a position of a seismic source to a predicted position of a seismic receiver j of a plurality of seismic receivers N (as, for example, illustrated in the setup 100 in
According to embodiments of the present disclosure, the aforementioned envelope of ghost wave functions for a positioning setup of seismic receivers may be defined as a function of angular frequency ω and given by:
Thus, maximizing the cost function defined by equation (1) corresponds to determining the positioning setup of seismic receivers that maximizes a minimum value of the envelope function defined by equation (2) for that positioning setup. The positioning setup of seismic receivers determined by applying this approach may enable the most efficient removal of the interference effects of the ghost data.
In some embodiments of the present disclosure, redundant information, for example to account for uncertainties in positioning of a plurality of seismic receivers, seismic sources, and/or noise (e.g., ghost data), and/or a width of a frequency band of interest, may be provided by increasing a number of seismic receiver positions. In some embodiments, three to ten or more seismic receiver positions may be used.
In some embodiments of the present disclosure, de-ghosted measurements of seismic data may be utilized as a basis for calibrating air gun source modeling. The accuracy of source modeling may be improved when calibration is based on ghost-free or near-ghost-free seismic data.
The seismic (marine) source modeling can be utilized, for example, to indentify noise components (signatures) when the seismic source is used to perform a seismic survey to obtain seismic data. Implementations of the marine source modeling can be calibrated and adjusted against measured signatures with a goal of high modeling accuracy within a limited frequency band. However, because multi-component streamers and source de-ghosting can allow for utilization of a broader range of frequencies in seismic imaging, adjustments to the modeling can be beneficial to achieve an accurate match between measured and modeled signatures over an extended frequency band. Examples of beneficial adjustments can include changes to the calibration process, for example, consideration of de-ghosted measurements and avoidance of, for example, traditional Digital Field System V (DFS-V) filtering. In this regard, modeling results after applying calibration can exhibit a match with measured array signatures over a wide frequency range.
Embodiments of the present disclosure can include new broadband calibration of the seismic source model based on newly acquired de-ghosted and unfiltered calibration data. In some embodiments, accurately modeled signatures that can also be used for planning and processing of marine broadband data, such as dual-sensor and multi-level source marine data, can be obtained. In some embodiments, the accurately modeled signatures of marine seismic sources (e.g., array of air guns) together with high accuracy sensors comprising a plurality of seismic receivers (e.g., array of hydrophones) positioned below a sea surface can be utilized in processing of marine broadband data for detecting the subsea geology.
The presence of a ghost reflection in a farfield signature from a seismic source can conceal information such as the array's response to its own ghost, or the interaction between subarrays.
Therefore, it can be beneficial for broadband source modeling that the modeling algorithm can be calibrated against data without the presence of a source ghost. In this regard, it can be possible to design calibration measurements such that they can allow for an accurate de-ghosting procedure to be applied to the data before calibration.
Awareness of various interaction effects, for example the interaction of the ghost, can be beneficial in ensuring sufficient accuracy in the de-ghosting process. Thus, the de-ghosting process described above may efficiently remove unwanted components (i.e., interference effect originating from the ghost signals) from a seismic frequency spectrum. This de-ghosting approach may avoid utilizing statistical processing steps that aim to (artificially) flatten the frequency spectrum (i.e., to produce frequency responses related to acquired seismic data having amplitude variations within predefined limits) by removing not only unwanted components but also some desirable components. For example, the interaction between seismic sources and subarrays can change the behavior of the seismic sources in the array and can leave footprints in the recorded seismic data, which, however, statistically based processing step(s) may alter or remove completely.
This can have implications for field measurements done for calibration purposes: on one hand it may be beneficial to measure seismic source (air gun) signatures in an as realistic environment as possible (e.g., in sea water at typical source depths); on the other hand, calibration measurements may be carried out in deep water and at a considerable distance from shore in order to avoid reflections from the sea floor or other objects in the relevant part of the signature. Seismic receivers can be located in the farfield, and their locations may need to be determined and adjusted (e.g., based on the approach described above) to allow accurate de-ghosting of the data, i.e., to acquire the data with minimized ghost signatures.
In some embodiments, in order to ensure accurate de-ghosting of the calibration in the presence of unavoidable location inaccuracies, it can be beneficial to receive signals with an array of farfield seismic receivers distributed over different depths and different angles from the seismic source. In this regard, a directional de-ghosting operator can be estimated from the data by inversion. In the case of a directional seismic source (e.g., directional air gun), the estimated operator may be utilized for obtaining a directional model of the directional seismic source. In some embodiments, an array of farfield seismic receivers at various depths and lateral positions can be beneficial for proper de-ghosting of the data, accounting for the uncertainty in exact receiver position when towing. According to embodiments of the present disclosure, by determining positions of farfield seismic receivers in the array based on the approach described above, effectiveness of de-ghosting of the calibration may be enhanced because of a mitigated interference effect of ghost signals received at the array of farfield seismic receivers.
In some embodiments of the present disclosure, a series of calibration measurements of a single seismic source and clustered seismic sources can be conducted for calibration purposes. Such a setup can allow for de-ghosting of the farfield signatures, and the de-ghosted signatures can be used in the calibration process such that a match between modeled and measured signatures over a broad frequency band can be achieved.
Some examples of signature characteristics that can be difficult to examine by inspecting signatures with a source ghost are illustrated in
Calibration of a seismic source where de-ghosted calibration data has been utilized is illustrated in
As illustrated in the graph 602 in
The difference between measured and modeled signatures (e.g., signatures 510 and 512 in
In a full array, interaction effects, and hence a change of resulting pressure output can be larger than for a single seismic source or a cluster of seismic sources. In this regard, it can become beneficial to interpret a modeling result over a broad frequency spectrum.
Performance of re-calibrated seismic source (air gun) modeling can be compared to a signature from a farfield measurement of a full 4135 cubic inch array equipped with G. GUN II air guns, for example.
Results illustrated in
In some embodiments, modeling can be calibrated towards single and clustered air guns, and not all volumes in the 4135 cubic inch array may be included in calibration. In some embodiments, the physics of the air gun modeling can handle interim volumes within the calibrated range, as well as interaction effects without additional tuning to achieve results illustrated in
In some embodiments of the present disclosure, marine broadband seismic may need higher accuracy in forward-modeling of seismic source (air gun) array signatures compared to other marine seismic. The accuracy of the modeling can be improved by basing the calibration process on broadband calibration data. In some embodiments, it can be beneficial to de-ghost the calibration data by acquiring time- and space-redundant data, which may avoid having acquisition setups with corresponding limiting instrument filters.
In some embodiments, utilizing de-ghosted, unfiltered broadband measurements can improve the validity of the range of modeling, and can provide a match between measured signatures over a broad frequency range. Comparison with a measured, de-ghosted farfield signature can show an increase in the accuracy of the modeling.
In some embodiments, a frequency span of 2-500 hertz can be utilized for broadband seismic applications. However, other spans (e.g., 1-1000 hertz) can be also used. In some embodiments, the frequency range under consideration can extend to 20 kilohertz or beyond.
In some embodiments, measurements can be taken in quiet waters, for example, a fjord or a deep lake. In some embodiments, the depth of the water can be at least 400 meters and measurements can be equidistant from the shore or any sub-surface structure that could act as a reflector. In some embodiments, background noise may not exceed 100 decibels referenced to 1 micropascal.
In some embodiments, a plurality of shots can be made from a seismic source. For example, 20, 30, or more energy emissions (e.g., shots) can be made from the seismic source. In addition, warm-up shots can be made to create stable output from the seismic source between sequences. In some embodiments, noise records can be acquired before and after each shot sequence to monitor the noise level and seismic receiver behavior.
In some cases, it can be difficult to achieve sufficient absolute accuracy of seismic receivers. In some embodiments, measuring the output of a seismic source with a known level of output can relax an absolute accuracy requirement, leaving the strict requirement to only be valid for relative accuracy of seismic receivers.
In some embodiments of the present disclosure, if the aforementioned measurement approach is used, the x-bar hydrophone test can be performed to determine relative sensitivities between the seismic receivers (e.g., hydrophones). The x-bar test can include switching the connections into the acquisition bench to separate the hydrophone/cable variations from other variations.
In some embodiments, farfield seismic receivers (e.g., farfield hydrophones) can be placed at different depths with a set horizontal distance to a seismic source. An example of a default setup 900 with the seismic source at x=0, y=0, z=Ds coordinated is illustrated in
In addition to the hydrophones, depth sensors can be included. In some embodiments, at least two depth sensors can be included, for example, one together with a seismic source, and one together with a deepest farfield seismic receiver. A conductivity, temperature, and depth (CTD) measurement can be completed to determine the temperature, salinity, pressure, and/or density in the water from the surface down to the depth of the deepest farfield seismic receiver. In some embodiments, the CTD measurement can be completed at least two times every 24 hours. In addition, a log of the sea state based at least in part on the CTD measurements can be kept.
For different air gun types, each of the different air gun bodies (e.g., different volumes) can be measured both as single guns and in a two-gun cluster configuration with typical cluster separation for the given air gun type and volume. In some embodiments of the present disclosure, a source depth of 6 meters with a firing pressure of 2000 pounds per square inch (psi) can be measured for each selected volume. In addition, a second depth (e.g., 10 meters), and a second pressure (e.g., 2500 psi) can be included. Air guns can be used in special configurations, for example, delta-cluster, hyper-cluster, or large in-line cluster. In some embodiments, air guns can have additional setup parameters, for example, various shuttle settings. Effects of varying these parameters can be measured for different air gun volumes.
According to embodiments of the present disclosure, the positions of the plurality of seismic receivers may be determined, at the step 1102, as illustrated in an example flow diagram in
At 1208, a minimum value of the envelope function of ghost wave functions for this particular positioning setup of seismic receivers obtained by equation (2) for some value of angular frequency ωmin (within a defined range of angular frequencies w related to data) may be determined. At 1210, the steps 1206 and 1208 may be repeated by adjusting the positions of the plurality of N seismic receivers (i.e., the position vectors rj and gj, j=1, . . . , N may be adjusted), and by determining a new minimum value (e.g., larger than the previous minimum value) for some value of angular frequency ωmean of a new envelope function obtained by equation (2) for the adjusted positions of seismic receivers. According to embodiments of the present disclosure, the steps 1206-1208 may be repeated until this minimum value of the envelope function is maximized, i.e., until the cost function defined by equation (1) is maximized. The positioning setup of seismic receivers (i.e., positions of N seismic receivers) determined in this manner may enable most efficient removal of interference effect of the ghost signals originating from the mirrored versions of the at least one seismic source.
According to embodiments of the present disclosure, a minimum overall value of the envelope function defined by equation (2) for a specific angular frequency w may relate to a deepest (cumulative) notch of a sum of ghost wave functions for all source-receiver pairs, which may be associated with a strongest interference effect of the ghost signals occurring at that specific angular frequency w for a particular positioning setup of seismic receivers. By adjusting the positions of the plurality of seismic receivers to maximize the cost function according to equation (1), the interference effect of the ghost signals may be most efficiently mitigated across the entire range of angular frequencies of data.
According to some embodiments of the present disclosure, as discussed above, determining the positions of the plurality of seismic receivers may comprise adjusting predicted positions of the plurality of seismic receivers in an iterative manner such that differences in travel times between arrivals at the plurality of receivers of direct signals from the at least one seismic source and arrivals at the plurality of receivers of the ghost signals from the mirrored versions of the at least one seismic source result in complementary ghost wave functions associated with the predicted positions of the plurality of seismic receivers (e.g., the complementary ghost wave functions associated with four hydrophones illustrated in
In some embodiments of the present disclosure, as discussed above, determining the positions of the plurality of seismic receivers may be based on maximizing a cost function. In an embodiment of the present disclosure, the cost function may be related to an envelope function of frequency domain amplitude responses of ghost wave functions associated with predicted positions of the plurality of seismic receivers for a specific position of the at least one seismic source. In an embodiment, the cost function may be defined as given by equation (1).
In some embodiments of the present disclosure, as discussed above, determining the positions of the plurality of seismic receivers to enable removal of the interference effect may comprise selecting initial positions of the plurality of seismic receivers, and converging from the initial positions to the determined positions of the plurality of seismic receivers to enable removal of the interference effect. For example, the converging may be achieved by adjusting positions of the plurality of seismic receivers in an iterative manner to enable (most efficient) removal of the interference effect of ghost signals originating from the mirrored versions of the at least one seismic source.
In an embodiment of the present disclosure, as discussed above, acquiring the data (e.g., ghost-free or near-ghost-free data) may comprise measuring the data by performing a plurality of energy emissions from the at least one seismic source (e.g., at least 20 shots from an array of air guns or at least 20 shots from an air gun if a single air gun is used). In an embodiment, as discussed above, the plurality of seismic receivers located at the determined positions may have an identical horizontal distance from the at least one seismic source, i.e., the plurality of seismic receivers may be vertically aligned and attached to a same rope. In an embodiment, the horizontal distance may comprise a trajectory parallel to a plane tangent to a spherical surface approximating a sea surface, wherein the at least one seismic source and the plurality of seismic receivers may be located below the sea surface. In another embodiment, the plurality of seismic receivers located at the determined positions may be positioned around the at least one seismic source.
In an embodiment of the present disclosure, as discussed above, determining the positions of the plurality of seismic receivers may comprise finding an angular frequency of the data that minimizes a value of an envelope function of frequency domain amplitude responses of ghost wave functions associated with the positions of the plurality of seismic receivers.
Some embodiments of the present disclosure can include a data product that comprises a calibrated model for a seismic source, the calibrated model having been made by the methods described herein. The data product can be transported and/or transmitted within the United States, imported to, or exported from the United States on physical media or by wire or broadcast transmission.
Although specific embodiments have been described above, these embodiments are not intended to limit the scope of the present disclosure, even where only a single embodiment is described with respect to a particular feature. Examples of features provided in the disclosure are intended to be illustrative rather than restrictive unless stated otherwise. The above description is intended to cover such alternatives, modifications, and equivalents as would be apparent to a person skilled in the art having the benefit of this disclosure.
The scope of the present disclosure includes any feature or combination of features disclosed herein (either explicitly or implicitly), or any generalization thereof, whether or not it mitigates any or all of the problems addressed herein. Various advantages of the present disclosure have been described herein, but embodiments may provide some, all, or none of such advantages, or may provide other advantages.
While the foregoing is directed to embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.
This application claims benefit of U.S. Provisional Patent Application Ser. No. 62/066,168, filed Oct. 17, 2014, entitled “Method of Acquiring Ghost-Free Signatures for Broadband Air Gun Calibration”, which is incorporated herein by reference.
Number | Date | Country | |
---|---|---|---|
62066168 | Oct 2014 | US |