System, method and computer program product for compression of seismic data

Abstract

A system for seismic surveying, method for performing seismic and a non-transitory computer readable medium having instructions stored therein that, when executed by one or more processors, cause the one or more processors to perform a method for performing seismic surveying including emitting seismic waves into a substrate, receiving seismic waves reflected from discontinuities within the substrate, converting the seismic waves into seismic traces, and representing the seismic traces by superimposed multiple tone sinusoidal waves using a parameter estimation. An optimized residual of the modelling is compressed using entropy coding or quantization coding techniques, and the optimized residual and the parameter sets are transmitted to a remote processing station for reconstruction and analysis of the discontinuities.

Claims

1. A seismic surveying system, comprising: at least one seismic generator for emitting seismic shots into a geological formation; at least one seismic receiver for receiving waves reflected from discontinuities within the geological formation; a computer interface coupled to the at least one seismic receiver for receiving the reflected waves, wherein the computer interface separates the received reflected waves into seismic traces, each seismic trace representing a number (N) of samples collected over a time period; a computing system adapted to receive the seismic traces from the computer interface, wherein the computing system comprises a processor with circuitry configured to represent the at least one trace by determining a set of decaying sinusoids (s.sub.k) which when superimposed represent the seismic trace (x(k)) to within a residual energy criterion, wherein each of the decaying sinusoids is characterized by a set of parameters including a starting time parameter (k.sub.i), a damping factor parameter (h.sub.i), a frequency parameter (ω.sub.i), at least one amplitude parameter (A.sub.i, B.sub.i) and a constant shift parameter (C.sub.i); and a memory module for storing the set of parameters.

2. The system of claim 1, wherein the seismic generator is at least one of a vibrator on a vibrator truck, an explosive device, an air gun, a water gun and a vapor gun, and wherein the seismic receiver is at least one of a geophone, a hydrophone, a seismometer, a microphone, a particle velocity sensor, a water pressure sensor, a particle displacement sensor, a particle acceleration sensor and a pressure gradient sensor.

3. The system of claim 1, wherein the circuitry is further configured to: reconstruct the seismic trace from the set of stored parameters; determine a residual between the seismic trace and the reconstructed seismic trace; determine whether the residual meets a residual energy criterion; based upon the set of stored parameters, determine a new set of parameters that characterize a new set of decaying sinusoids; reconstruct the seismic trace from the new set of parameters; determine a new residual between the seismic trace and the reconstructed seismic trace from the new set of parameters; determine whether the new residual meets a pre-defined residual energy criterion; reiterate the determining of the residual and the determining of a set of parameters based upon the previously determined set of parameters until the residual meets the pre-defined residual energy criterion; and store the sets of parameters with the residual which meets pre-defined residual energy criterion.

4. The system of claim 3, wherein the circuitry further comprises a compression module, wherein the compression module compresses the residual that meets the pre-defined residual energy criterion using entropy coding or quantization coding; and stores the encoded residual with the sets of parameters.

5. The system of claim 4, wherein the circuitry further comprises a transmitter which transmits the stored encoded residual energy and the stored set of parameters to a remote processing station.

6. The system of claim 3, wherein the circuitry is further configured to determine the set of parameters using a single tone representation of the seismic trace; and determine the new set of parameters using a single tone representation of the residual at each iteration.

7. The system of claim 1, wherein the seismic trace (x(k)) is represented as $x\left(k\right)=\underset{i=1}{\overset{M}{.Math.}}{s}_{{\theta }_i}\left(k\right)+v\left(k\right),k=1,2,\mathrm{.Math.},N$ where v(k) denotes the unrepresented signal, and s.sub.θi (k) are exponentially decaying sinusoidal waves defined as $s_{{\theta }_i}\left(k\right)=\{\begin{array}{cc}{e}^{-h_i\left(k-k_i\right)T}{\alpha }_i\sin \left({\omega }_i\left(k-k_i\right)T+{\varphi }_i\right)+C_i& k\ge k_i\\ 0& k<k_i\end{array}$ where T is the sampling time of the seismic trace and h.sub.i, k.sub.i, α.sub.i, ω.sub.i, ϕ.sub.i and C.sub.i are representative of the damping factor, starting time, initial amplitude, frequency, phase shift and constant shift, respectively.

8. The system of claim 7, wherein the decaying sinusoidal waves s.sub.θi (k) are recast, for k≥k.sub.i, into the form
s.sub.θ.sub.i(k)=e.sup.−h.sup.i.sup.(k−k.sup.i.sup.)TA.sub.i cos ω.sub.i(k−k.sub.i)T+e.sup.−h.sup.i.sup.(k−k.sup.i.sup.)TB.sub.i sin ω.sub.i(k−k.sub.i)T+C.sub.i where A.sub.i=α.sub.i sin ϕ.sub.i and B.sub.i=α.sub.i cos ϕ.sub.i; and wherein the parameter set for each sinusoidal wave is defined as θ.sub.i={k.sub.i, π.sub.i} and π.sub.i={A.sub.i, B.sub.i, C.sub.i, ω.sub.i, h.sub.i}.

9. The system of claim 8, wherein the residual is defined as: $r_{{\theta }_i^{*}}\left(k\right)=x\left(k\right)-\underset{j=1}{\overset{i}{.Math.}}{s}_{{\theta }_j^{*}}\left(k\right)$ and, either the residual energy criterion is defined as the residual energy is lower than a predefined threshold (r.sub.E),
r.sub.θ.sub.i*.sup.2<r.sub.E, wherein the iteration is terminated when the residual energy is lower than the predefined threshold, or the residual energy dropping rate is lower than a predefined threshold, (r.sub.D), $\frac{r_{{\theta }_{j-1}^{*}}^2-r_{{\theta }_j^{*}}^2}{r_{{\theta }_{j-1}^{*}}^2}100\%<r_D$ wherein the iteration is terminated when the residual energy dropping rate is lower than the predefined threshold.

10. A method for seismic surveying comprising: emitting, using at least one seismic generator, seismic shot waves into a geological formation; receiving, using at least one seismic receiver, waves reflected from discontinuities within the geological formation; receiving, at a computer interface coupled to the at least one seismic receiver, the reflected waves, the computer interface separating the received reflected waves into seismic traces, each seismic trace representing a number (N) of samples collected over a time period; receiving, by a computing system, the seismic traces from the computer interface, wherein the computing system comprises a processor with circuitry configured for representing the at least one trace by determining a set of decaying sinusoids (s.sub.k), which when superimposed, represent the seismic trace (x(k) to within a residual energy criterion, wherein each of the decaying sinusoids is characterized by a set of parameters including a starting time parameter (k.sub.i), a damping factor parameter (h.sub.i), a frequency parameter (ω.sub.i), at least one amplitude parameter (A.sub.i, B.sub.i) and a constant shift parameter (C.sub.i); and storing the set of parameters in a memory module.

11. The method of claim 10, further comprising: wherein emitting the seismic shot using at least one generator is selected from a list comprising a vibrator on a vibrator truck, an explosive device, an air gun, a water gun, a vapor gun, an electric spark, an electric arc, a marine vibrator; and wherein the at least one seismic receiver is selected from a list comprising a geophone, a hydrophone, a seismometer, a microphone, a particle velocity sensor, a water pressure sensor, a particle displacement sensor, a particle acceleration sensor and a pressure gradient sensor.

12. The method of claim 10, wherein the circuitry is further configured for: reconstructing the seismic trace from the set of stored parameters; determining a residual between the seismic trace and the reconstructed seismic trace; determining whether the residual meets a residual energy criterion; determining a new set of parameters based upon the set of stored parameters, wherein the new set of parameters characterize a new set of decaying sinusoids; reconstructing the seismic trace from the new set of parameters; determining a new residual between the seismic trace and the reconstructed seismic trace from the new set of parameters; determining whether the new residual meets a pre-defined residual energy criterion; reiterating the determining of the residual and the determining of a set of parameters based upon the previously determined set of parameters until the residual meets the pre- defined residual energy criterion; and storing the sets of parameters determined at each iteration with the residual which meets pre-defined residual energy criterion.

13. The method of claim 12, the circuitry further comprising: a compression module adapted for compressing the residual that meets the pre-defined residual energy criterion using entropy coding or quantization coding; and storing the encoded residual with the sets of parameters.

14. The method of claim 13, the circuitry further comprising: a transmitter for transmitting the stored encoded residual energy and the stored set of parameters to a remote processing station.

15. The method of claim 12, further configuring the circuitry for: determining the set of parameters using a single tone representation of the seismic trace; and determining the new set of parameters using a single tone representation of the residual at each iteration.

16. The method of claim 10, wherein the circuitry is further configured for: representing the seismic trace (x(k)) as $x\left(k\right)=\underset{i=1}{\overset{M}{.Math.}}{s}_{{\theta }_i}\left(k\right)+v\left(k\right),k=1,2,\mathrm{.Math.},N$ where v(k) denotes the unmodeled signal, and s.sub.θi (k) are exponentially decaying sinusoidal waves defined as $s_{{\theta }_i}\left(k\right)=\{\begin{array}{cc}{e}^{-h_i\left(k-k_i\right)T}{\alpha }_i\sin \left({\omega }_i\left(k-k_i\right)T+{\varphi }_i\right)+C_i& k\ge k_i\\ 0& k<k_i\end{array}$ where T is the sampling time of the seismic trace and h.sub.i, k.sub.i, α.sub.i, ω.sub.i, ϕ.sub.i and C.sub.i are representative of the damping factor, starting time, initial amplitude, frequency, phase shift and constant shift, respectively.

17. The method of claim 16, wherein the circuitry is further adapted for: recasting the decaying sinusoidal waves s.sub.θi (k), for k≥k.sub.i, into the form:
s.sub.θ.sub.i(k)=e.sup.−h.sup.i.sup.(k−k.sup.i.sup.)TA.sub.i cos ω.sub.i(k−k.sub.i)T+e.sup.−h.sup.i.sup.(k−k.sup.i.sup.)TB.sub.i sin ω.sub.i(k−k.sub.i)T+C.sub.i where A.sub.i=α.sub.i sin ϕ.sub.i and B.sub.i=α.sub.i cos ϕ.sub.i; and defining the parameter set for each sinusoidal wave as θ.sub.i={k.sub.i, π.sub.i} and π.sub.i={A.sub.i, B.sub.i, C.sub.i, ω.sub.i, h.sub.i}.

18. The system of claim 17, wherein the circuitry is further adapted for: defining the residual as: $r_{{\theta }_i^{*}}\left(k\right)=x\left(k\right)-\underset{j=1}{\overset{i}{.Math.}}{s}_{{\theta }_j^{*}}\left(k\right)$ and defining the residual energy criterion either as the residual energy is lower than a predefined threshold (r.sub.E),
r.sub.θ.sub.j*.sup.2<r.sub.E, and terminating the iterating when the residual energy is lower than the predefined threshold, the residual energy dropping rate is lower than a predefined threshold, (r.sub.D), $\frac{r_{{\theta }_{j-1}^{*}}^2-r_{{\theta }_j^{*}}^2}{r_{{\theta }_{j-1}^{*}}^2}100\%<r_D$ and terminating the iterating when the residual energy dropping rate is lower than the predefined threshold.

19. A non-transitory computer readable medium having instructions stored therein that, when executed by one or more processor, cause the one or more processors to perform a method of seismic surveying, the method comprising: emitting, using at least one seismic generator, seismic shot waves into a geological formation; receiving, using at least one seismic receiver, waves reflected from discontinuities within the geological formation; receiving, at a computer interface coupled to the at least one seismic receiver, the reflected waves, the computer interface separating the received reflected waves into seismic traces, each seismic trace representing a number (N) of samples collected over a time period; receiving, by a computing system, the seismic traces from the computer interface, wherein the computing system comprises a processor with circuitry configured for representing the at least one trace by determining a set of decaying sinusoids (s.sub.k), which when superimposed, represent the seismic trace (x(k) to within a residual energy criterion, wherein each of the decaying sinusoids is characterized by a set of parameters including a starting time parameter (k.sub.i), a damping factor parameter (h.sub.i), a frequency parameter (ω.sub.i), at least one amplitude parameter (A.sub.i, B.sub.i) and a constant shift parameter (C.sub.i); reconstructing the seismic trace from the set of parameters; determining a residual between the seismic trace and the reconstructed seismic trace; determining whether the residual meets a residual energy criterion; determining a new set of parameters based upon the set of parameters, wherein the new set of parameters characterize a new set of decaying sinusoids; reconstructing the seismic trace from the new set of parameters; determining a new residual between the seismic trace and the reconstructed seismic trace from the new set of parameters; determining whether the new residual meets a pre-defined residual energy criterion; reiterating the determining of the residual and the determining of a set of parameters based upon the previously determined set of parameters until the residual meets the pre- defined residual energy criterion; storing the sets of parameters determined at each iteration with the residual which meets pre-defined residual energy criterion; the circuitry further comprising a compression module, compressing, using the compression module, the residual that meets the pre-defined residual energy criterion using entropy coding or quantization coding; reconstructing the seismic trace from the new set of parameters; determining a new residual between the seismic trace and the reconstructed seismic trace from the new set of parameters; determining whether the new residual meets a pre-defined residual energy criterion; reiterating the determining of the residual and the determining of a set of parameters based upon the previously determined set of parameters until the residual meets the pre-defined residual energy criterion; compressing, using a compression module, the residual that meets the pre-defined residual energy criterion using entropy coding or quantization coding; and storing the encoded residual with the sets of parameters; and transmitting the stored encoded residual energy and the stored set of parameters to a remote processing station.

20. The non-transitory computer readable medium of claim 19, further comprising modelling the seismic trace (x(k)) as $x\left(k\right)=\underset{i=1}{\overset{M}{.Math.}}{s}_{{\theta }_i}\left(k\right)+v\left(k\right),k=1,2,\mathrm{.Math.},N$ where v(k) denotes the unmodeled signal, and s.sub.θi (k) are exponentially decaying sinusoidal waves defined as $s_{{\theta }_i}\left(k\right)=\{\begin{array}{cc}{e}^{-h_i\left(k-k_i\right)T}{\alpha }_i\sin \left({\omega }_i\left(k-k_i\right)T+{\varphi }_i\right)+C_i& k\ge k_i\\ 0& k<k_i\end{array}$ where T is the sampling time of the seismic trace and h.sub.i, k.sub.i, α.sub.i, ω.sub.i, θ.sub.i and C.sub.i are the damping factor, starting time, initial amplitude, frequency, phase shift and constant shift, respectively; and recasting the decaying sinusoidal waves s.sub.θi (k), for k≥k.sub.i, into the form:
s.sub.θ.sub.i(k)=e.sup.−h.sup.i.sup.(k−k.sup.i.sup.)TA.sub.i cos ω.sub.i(k−k.sub.i)T+e.sup.−h.sup.i.sup.(k−k.sup.i.sup.)TB.sub.i sin ω.sub.i(k−k.sub.i)T+C.sub.i where A.sub.i=α.sub.i sin ϕ.sub.i and B.sub.i=α.sub.i cos ϕ.sub.i; and defining the parameter set for each sinusoidal wave as θ.sub.i={k.sub.i, π.sub.i} and π.sub.i={A.sub.i, B.sub.i, C.sub.i, ω.sub.i, h.sub.i}.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) A more complete appreciation of the invention and many of the attendant advantages thereof will be readily obtained as the same becomes better understood by reference to the following detailed description when considered in connection with the accompanying drawings, wherein:

(2) FIG. 1A illustrates a ground based seismic energy shot source;.

(3) FIG. 1B illustrates a system for generating seismic waves, receiving the reflected waves and coupling the seismic waves to a computer interface.

(4) FIG. 2 illustrates a computer interface adapted to receive the waves from the seismic receivers and convert the waves into N traces each of sampling time T.

(5) FIG. 3 depicts an overview of the system from reception of the reflected seismic waves through transmission to a remote processing station.

(6) FIG. 4 illustrates a computer system modified to process, represent, compress and transmit the wave traces of FIG. 2.

(7) FIG. 5A illustrates a comparison of the reconstruction quality with residual quantization between a prior art LPC method and the method of the present disclosure showing the normalized reconstruction error with level-64 quantization;

(8) FIG. 5B illustrates a comparison of the reconstruction quality with residual quantization between a prior art LPC method and the method of the present disclosure showing the normalized reconstruction error with level-128 quantization; and

(9) FIG. 5C illustrates a comparison of the reconstruction quality with residual quantization between a prior art LPC method and the method of the present disclosure showing the normalized reconstruction error with with level-256 quantization.

(10) FIG. 6 shows distributed components including one or more client and server machines, which may share processing.

(11) FIG. 7 shows hardware for the computing device used in the exemplary embodiments.

(12) FIG. 8 shows circuitry configured to perform features of the invention.

(13) FIG. 9 shows a data processing system hub used in the exemplary embodiments.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

(14) Referring now to the drawings, wherein like reference numerals designate identical or corresponding parts throughout the several views.

(15) Embodiments to a system for seismic surveying, a method for seismic surveying and a non-transitory computer readable medium having instructions stored therein that, when executed by one or more processor, cause the one or more processors to perform a method of seismic surveying are now described.

(16) In one aspect, the generation and reception of signals for seismic acquisition is accomplished by applying a source of energy, which is shown in one aspect as generated by a vibrator 120 carried by a vibrator truck as shown in FIG. 1A. Each vibration is referred to as a shot. FIG. 1B illustrates a seismic energy shot generator 120 emitting shot waves 122a, 122b and 122c into a ground surface 126. The shot waves are reflected from the ground surface 126 and discontinuities 128 and 130 to generate reflected waves 124a, 124b and 124c. Reflected waves 124a, 124b and 124c are received by a seismic receivers 132a, 132b and 132c (geophones, for example) spaced at regular intervals on the ground. In a second aspect (not shown) acoustic shot waves are emitted into either the surface of or within a body of water by an acoustic shot source (an acoustic gun, for example). Reflected waves representing the discontinuities in the underlying surface are received by at least one hydrophone, acoustic geophone, or other acoustic receiver.

(17) The seismic generator is preferably at least one of a vibrator on a vibrator truck, an explosive device, an air gun, a water gun and a vapor gun,

(18) The seismic receiver is preferably at least one of a geophone, a hydrophone, a seismometer, a microphone, a particle velocity sensor, a water pressure sensor, a particle displacement sensor, a particle acceleration sensor and a pressure gradient sensor

(19) Although FIG. 1B shows only three waves and three receivers for clarity, an actual survey site may have a plurality of shot generators and many hundreds of seismic receivers.

(20) Although not shown in FIG. 1B, each seismic shot furcates into a number of waves, which are reflected from discontinuities in the geologic structure. In all instances, the reflected waves have different paths, frequencies and arriving times.

(21) The seismic receivers are coupled to a computer interface 111. As shown in FIG. 1B, the reflected waves from the seismic receivers are received by computer interface 111.

(22) Within the computer interface 111 is a timing module 236, having circuitry configured to provide a time period (T) for sampling the reflected waves. The seismic waves are sampled by the Sampling Circuitry 234 during the time period, thus generating seismic traces, each trace representing a number (N) of samples collected over a time period (T).

(23) FIG. 3 is an overview of the invention, in which the seismic receivers 132 are coupled to Computer Interface 340. The Computer Interface 340 has circuitry which separates the signals into N traces each of sampling time period (T) and sends the N traces 338 to Computing System 360. Computing System 360 includes circuitry configured to create a representation of the N traces in a Trace Representation Module 364, and the results are compressed in the Compression Module 368. The compressed results are then transmitted to remote processing station 380 for reconstruction and analysis. While FIG. 3 describes an embodiment in which the compressed results are transmitted wirelessly to the remote processing station through the cloud or via satellite, other transmission modes may be used including land line or transmission through hard storage devices.

(24) FIG. 4 shows the Computing System 460 in more detail. The Computing System 460 receives the N seismic traces from the computer interface. Computing System 460 includes a processor (CPU 462) having a Trace Representation module 464 adapted to represent the at least one trace by determining a set of decaying sinusoids (s.sub.k) which when superimposed represent the seismic trace to within a residual energy criterion, wherein each of the decaying sinusoids is characterized by a set of parameters including a starting time parameter (k.sub.i), a damping factor parameter (h.sub.i), a frequency parameter (ω.sub.i), at least one amplitude parameter (A.sub.i, B.sub.i) and a constant shift parameter (C.sub.i). The Computing System 460 also includes a Memory Module 466 for storing the set of parameters at Parameter Storage 467, as well as enabling other memory functions of the CPU 462.

(25) In further detail, the processor 462 reconstructs the seismic trace from the set of stored parameters; determines a residual between the seismic trace and the reconstructed seismic trace; determines whether the residual meets a residual energy criterion at Residual Comparator 465; based upon the set of stored parameters, determines a new set of parameters that characterize a new set of decaying sinusoids; stores the new set of parameters in Parameter Storage 467; reconstructs the seismic trace from the new set of stored parameters; determines a new residual between the seismic trace and the reconstructed seismic trace from the new set of parameters; determines whether the new residual meets a pre-defined residual energy criterion at Residual Comparator 365; reiterates the determining of the residual and the determining of a set of parameters based upon the previously determined set of parameters until the residual meets the pre-defined residual energy criterion; and stores the sets of parameters at Parameter Storage 467 determined at each iteration with the residual which meets pre-defined residual energy criterion (stored at Residual Storage 469).

(26) A method for seismic surveying is now described with respect to FIGS. 1-4. The method includes emitting, using at least one seismic generator, seismic shot waves into a geological formation as shown in FIG. 1A or into a body of water as described with respect to underwater surveying (not shown). The method includes receiving, by using at least one seismic receiver (132a, 132b, 132c), reflected waves which represent discontinuities within the geological formation; receiving the received reflected waves at a computer interface 111 coupled to the at least one seismic receiver (132a, 132b, 132c), and separating the received reflected waves into N seismic traces 238, each seismic trace 238 representing a number (N) of samples collected over a time period.

(27) The seismic traces are sent to computing system 460, wherein the computing system 460 comprises a processing module 462 adapted to create a representation of the seismic traces. The computing system 460 performs the method of representing the at least one trace by determining a set of decaying sinusoids (s.sub.k) which when superimposed represent the seismic trace to within a residual energy criterion, wherein each of the decaying sinusoids is characterized by a set of parameters including a starting time parameter (k.sub.i), a damping factor parameter (h.sub.i), a frequency parameter (ω.sub.i), at least one amplitude parameter (A.sub.i, B.sub.i) and a constant shift parameter (C.sub.i); and storing the set of parameters in Parameter Storage 467 in Memory Module 466.

(28) The method further comprises the Trace Representation Module 464 of the processor reconstructing the seismic trace from the set of stored parameters; determining a residual between the seismic trace and the reconstructed seismic trace; determining whether the residual meets a residual energy criterion at Residual Comparator 465; determining a new set of parameters based upon the set of stored parameters, wherein the new set of parameters characterize a new set of decaying sinusoids; storing the new set of parameters in Parameter Storage 467 of Memory Module 466; reconstructing the seismic trace from the new set of stored parameters; determining a new residual between the seismic trace and the reconstructed seismic trace from the new set of parameters; determining whether the new residual meets a pre-defined residual energy criterion at Residual Comparator 465; reiterating the determining of the residual and the determining of a set of parameters based upon the previously determined set of parameters until the residual meets the pre-defined residual energy criterion; and storing, in Parameter Storage 467 of Memory Module 466 and Residual Storage 467, the sets of parameters determined at each iteration with the residual which meet the pre-defined residual energy criterion.

(29) In both the system and method embodiments, a compression module 468 compresses the residual that meets the pre-defined residual energy criterion using entropy coding or quantization coding; and stores the encoded residual along with the sets of parameters in Memory Module 466 within Residual Storage 469 and Parameter Storage 467 respectively. The stored encoded residual energy and the stored set of parameters are sent to Network Controller 472 and antenna 470 and transmitted to remote processing station 360.

(30) A further embodiment to a non-transitory computer readable medium having instructions stored therein that, when executed by one or more processor, cause the one or more processors to perform a method of seismic surveying, is now described with respect to FIGS. 1-4.

(31) The performance of the method of seismic surveying includes emitting, using at least one seismic generator, seismic shot waves into the substrate of the ground at a location as shown in FIG. 1A or as described with respect to underwater surveying (not shown). Receiving, by using at least one seismic receiver (132a, 132b, 132c), reflected waves which represent discontinuities within the substrate; receiving the reflected waves at a computer interface 111 coupled to the at least one seismic receiver (132a, 132b, 132c), and separating the received reflected waves into N seismic traces 238, each seismic trace 238 representing a number (N) of samples collected over a time period (T).

(32) The seismic traces are sent to computing system 460, wherein the computing system 460 comprises a processor 462 adapted to create a representation of the seismic traces. The computing system 460 performs the method of representing the at least one trace by determining a set of decaying sinusoids (s.sub.k) which when superimposed represent the seismic trace to within a residual energy criterion, wherein each of the decaying sinusoids is characterized by a set of parameters including a starting time parameter (k.sub.i), a damping factor parameter (h.sub.i), a frequency parameter (ω.sub.i), at least one amplitude parameter (A.sub.i, Bi) and a constant shift parameter (C.sub.i); and storing the set of parameters in Parameter Storage 467 in Memory Module 466.

(33) The performance of the method further comprises Trace Representation Module 464 of the processor reconstructing the seismic trace from the set of stored parameters; determining a residual between the seismic trace and the reconstructed seismic trace; determining whether the residual meets a residual energy criterion at Residual Comparator 465; determining a new set of parameters based upon the set of stored parameters, wherein the new set of parameters characterize a new set of decaying sinusoids; storing the new set of parameters in Parameter Storage 467 of Memory Module 466; reconstructing the seismic trace from the new set of stored parameters; determining a new residual between the seismic trace and the reconstructed seismic trace from the new set of parameters; determining whether the new residual meets a pre-defined residual energy criterion at Residual Comparator 465; reiterating the determining of the residual and the determining of a set of parameters based upon the previously determined set of parameters until the residual meets the pre-defined residual energy criterion; and storing, in Parameter Storage 467 of memory module 466 and Residual Storage 467, the sets of parameters determined at each iteration with the residual which meet the pre-defined residual energy criterion. The method further includes compressing, using a compression module 468, the residual that meets the pre-defined residual energy criterion using entropy coding or quantization coding; and storing the encoded residual with the sets of parameters within Residual Storage 469 and Parameter Storage 467 respectively; and transmitting, using Network Controller 472 and antenna 470, the stored encoded residual energy and the stored set of parameters to remote processing station 360.

(34) The trace representation executed by the processor for the system, the method and the non-transitory computer readable medium is now described in more detail. In order to emulate the nature of seismic wave propagation, a seismic trace, x(k), with N time samples is represented as a superposition of a number of exponentially decaying sinusoidal waves as:

(35) $\begin{array}{cc}x\left(k\right)=\underset{i=1}{\overset{M}{.Math.}}{s}_{{\theta }_i}\left(k\right)+v\left(k\right),k=1,2,\mathrm{.Math.},N& \left(1\right)\end{array}$
where v(k) denotes the unrepresented signal, and s.sub.θi (k) are exponentially decaying sinusoidal waves defined as

(36) $\begin{array}{cc}{s}_{{\theta }_i}\left(k\right)=\{\begin{array}{cc}{e}^{-h_i\left(k-k_i\right)T}{\alpha }_i\sin \left({\omega }_i\left(k-k_i\right)T+{\varphi }_i\right)+C_i& k\ge k_i\\ 0& k<k_i\end{array}& \left(2\right)\end{array}$
where T is the sampling time of the seismic trace and h.sub.i, k.sub.i, α.sub.i, ω.sub.i, ϕ.sub.i and C.sub.i are: damping factor, starting time, initial amplitude, frequency, phase shift and constant shift, respectively. The decaying sinusoidal waves s.sub.θi (k) could be recast, for k≥k.sub.i, into the following form:
s.sub.θ.sub.i(k)=e.sup.−h.sup.i.sup.(k−k.sup.i.sup.)TA.sub.i cos ω.sub.i(k−k.sub.i)T+e.sup.−h.sup.i.sup.(k−k.sup.i.sup.)TB.sub.i sin ω.sub.i(k−k.sub.i)T+C.sub.i,   (3)
where A.sub.i=α.sub.i sin ϕ.sub.i and B.sub.i=α.sub.i cos ϕ.sub.i. Define the parameter set for each sinusoidal wave as θ.sub.i={k.sub.i, π.sub.i} and π.sub.i={A.sub.i, B.sub.i, C.sub.i, ω.sub.i, h.sub.i}, then the parameters of all sinusoidal waves are defined as Θ={θ.sub.1,θ.sub.2, . . . , θ.sub.M}.

(37) First of all, parameter estimation is considered for a single tone sinusoidal representation case, i.e. representation (1) with M=1, θ={ks, π}, and π={A, B, C, ω, h} for k≥ks:
x(k)=e.sup.−h(k−k.sup.s.sup.)TA cos ω(k−k.sub.s)T+e.sup.−h(k−k.sup.s.sup.)TB sin ω(k−k.sub.s)T+C.   (4)

(38) Approximate continuous time functions e.sup.−ht cos ωt and e.sup.−ht sin ωt by using the first order Taylor series at h=ĥ and ω={circumflex over (ω)} as
e.sup.−ht cos ωt≠e.sup.−ĥt cos {circumflex over (ω)}t−te.sup.−ĥt sin {circumflex over (ω)}.Math.∇{circumflex over (ω)}−te.sup.−ĥt cos {circumflex over (ω)}t.Math.∇ĥ  (5)
e.sup.−ht sin ωt≠e.sup.−ĥt sin {circumflex over (ω)}t+te.sup.−ĥ cos {circumflex over (ω)}.Math.∇{circumflex over (ω)}−te.sup.−ĥt cos {circumflex over (ω)}t.Math.∇ĥ  (6)
where ∇{circumflex over (ω)}=ω−{circumflex over (ω)} and ∇ĥ=h−ĥ. Then a continuous time decaying sinusoidal wave,
s.sub.θ(t)=e.sup.−htA cos ωt+e.sup.−htB sin ωt+C,   (7)
is approximated using (5) and (6) as follows,

(39) $\begin{array}{cc}{s}_{\theta }\left(t\right)\approx {\mathrm{Ae}}^{-\hat{h}t}\cos \hat{\omega }t-A{\mathrm{te}}^{-\hat{h}t}\sin \hat{\omega }.Math.\nabla \hat{\omega }-{\mathrm{Ate}}^{-\hat{h}t}\cos \hat{\omega }t.Math.\nabla \hat{h}+{\mathrm{Be}}^{-\hat{h}t}\sin \hat{\omega }t+{\mathrm{Bte}}^{-\hat{h}t}\cos \hat{\omega }.Math.\nabla \hat{\omega }-{\mathrm{Bte}}^{-\hat{h}t}\sin \hat{\omega }t.Math.\nabla \hat{h}+C& \left(8\right)\\ ={\mathrm{Ae}}^{-\hat{h}t}\cos \hat{\omega }t+{\mathrm{Be}}^{-\hat{h}t}\sin \hat{\omega }t-\left({\mathrm{Ate}}^{-\hat{h}t}\sin \hat{\omega }-{\mathrm{Bte}}^{-\hat{h}t}\cos \hat{\omega }\right)\nabla \hat{\omega }-\left({\mathrm{Ate}}^{-\hat{h}t}\cos \hat{\omega }t+{\mathrm{Bte}}^{-\hat{h}t}\sin \hat{\omega }t\right)\nabla \hat{h}+C& \left(9\right)\end{array}$

(40) Assume that the starting time instant of this single tone sinusoidal wave is k. Define parameter vector [A B C ∇{circumflex over (ω)} ∇ĥ].sup.T, data vector x(k)=[x(k) x(k+1) . . . x(N)]T, t.sub.j=jT and matrix

(41) $\begin{array}{cc}D\left(A,B,\hat{\omega },\hat{h}\right)=\left[\begin{array}{cc}P\left(\hat{\omega },\hat{h}\right)& Q\left(A,B,\hat{\omega },\hat{h}\right)\end{array}\right]\mathrm{where}& \left(10\right)\\ P\left(\hat{\omega },\hat{h}\right)=\left[\begin{array}{ccc}{e}^{-\hat{h}{t}_1}\cos \hat{\omega }{t}_1& e^{-\hat{h}{t}_1}\sin \hat{\omega }{t}_1& 1\\ e^{-\hat{h}{t}_2}\cos \hat{\omega }{t}_2& e^{-\hat{h}{t}_2}\sin \hat{\omega }{t}_2& 1\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ e^{-\hat{h}{t}_N}\cos \hat{\omega }{t}_N& e^{-\hat{h}{t}_N}\sin \hat{\omega }{t}_N& 1\end{array}\right]Q\left(A,B,\hat{\omega },\hat{h}\right)=\left[\begin{array}{cc}\begin{array}{c}-{\mathrm{At}}_1{e}^{-\hat{h}t}\sin \hat{\omega }+\\ {\mathrm{Bt}}_1{e}^{-\hat{h}{t}_1}\cos \hat{\omega }\end{array}& \begin{array}{c}-{\mathrm{At}}_1{e}^{-\hat{h}{t}_1}\cos \hat{\omega }{t}_1-\\ {\mathrm{Bt}}_1{e}^{-\hat{h}{t}_1}\sin \hat{\omega }{t}_1\end{array}\\ \begin{array}{c}-{\mathrm{At}}_2{e}^{-\hat{h}t}\sin \hat{\omega }+\\ {\mathrm{Bt}}_2{e}^{-\hat{h}{t}_1}\cos \hat{\omega }\end{array}& \begin{array}{c}-{\mathrm{At}}_2{e}^{-\hat{h}{t}_1}\cos \hat{\omega }{t}_2-\\ {\mathrm{Bt}}_2{e}^{-\hat{h}{t}_1}\sin \hat{\omega }{t}_2\end{array}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ \begin{array}{c}-{\mathrm{At}}_N{e}^{-\hat{h}{t}_N}\sin \hat{\omega }+\\ {\mathrm{Bt}}_N{e}^{-\hat{h}{t}_N}\cos \hat{\omega }\end{array}& \begin{array}{c}-{\mathrm{At}}_N{e}^{-\hat{h}{t}_N}\sin \hat{\omega }{t}_N-\\ {\mathrm{Bt}}_N{e}^{-\hat{h}{t}_N}\cos {\hat{\omega }}_N\end{array}\end{array}\right]& \left(11\right)\end{array}$

(42) This provides a first order approximation of data as
x(k)=D(A, B, {circumflex over (ω)}, ĥ)z.   (12)

(43) This nonlinear optimization can be solved by the following steps. First, initialize θ to the parameter vector z.sup.(0)=[Â.sup.(0), {circumflex over (B)}.sup.(0), Ĉ.sup.(0), ∇{circumflex over (ω)}.sup.(0)=0, ∇ĥ.sup.(0)=0].sup.T. Secondly, at each iteration, update the estimates of {circumflex over (ω)}.sup.(i)={circumflex over (ω)}.sup.(i−1)+∇{circumflex over (ω)}.sup.(i−1), ĥ.sup.(i)=ĥ.sup.(i−1)+∇ĥ.sup.(i−1), and create D.sub.i=(Â.sup.(i−1){circumflex over (B)}.sup.(i−1){circumflex over (ω)}.sup.(i)ĥ.sup.(i)), then (12) is linearized. Last, solve the parameter estimation vector from the linearized optimization z.sup.(i)=(D.sub.i.sup.TD.sub.i).sup.−1D.sub.i.sup.T. These steps are repeated until the estimate of parameter π={A, B, C, ω, h} converges. The steps are summarized in TABLE 1.

(44) TABLE-US-00001 TABLE 1 Parameter Estimation for the Single Tone Case Step Operation a) z.sup.(0) = [Â.sup.(0) , {circumflex over (B)}.sup.(0) , Ĉ.sup.(0) , ∇{circumflex over (ω)}.sup.(0) = 0, ∇ĥ.sup.(0) = 0].sup.T , ŵ.sup.(0) , ĥ.sup.(0) , i = 0 b) Next iteration i = i + 1 c) {circumflex over (ω)}.sup.(i) = {circumflex over (ω)}.sup.(i−1) + ∇{circumflex over (ω)}.sup.(i−1) , ĥ.sup.(i) = ĥ.sup.(i−1) + ∇ĥ.sup.(i−1) d) Generate D.sub.i = {circumflex over (D)}(Â.sup.(i−1) {circumflex over (B)}.sup.(i−1) {circumflex over (ω)}.sup.(i) ĥ.sup.(i)) e) Solve z.sup.(i) = (D.sub.i.sup.TD.sub.i).sup.−1D.sub.i.sup.Ty.sup.T f) Repeat steps b)-e) until convergence

(45) For each postulated starting time instant k, a corresponding parameter set π(k) can be solved as mentioned above by applying the proposed steps in TABLE 1 to nonlinear optimization (12), hence a corresponding reconstruction error is obtained as

(46) 0$\begin{array}{cc}e\left(k,\pi \left(k\right)\right)=\underset{j=1}{\overset{N}{.Math.}}{\left(x\left(j\right)-{\hat{x}}_k\left(j\right)\right)}^2,& \left(13\right)\end{array}$
where signal {circumflex over (x)}.sub.k(j) denotes the reconstructed seismic trace as in (4) with a postulated starting time k and the corresponding parameter estimate π(k).

(47) In this work, the starting time k.sub.s is estimated as
{circumflex over (k)}.sub.s=arg.sub.k min e(k.sub.sπ(k)).   (14)

(48) Optimization (14) can be solved either in an exhaustive way, i.e., to solve reconstruction error e(k, π(k)) for all possible time instants, then choose the minimum one, or in an iterative grid searching way. Note that a too large postulated starting time k results in a short data vector y(k), consequently leading to a convergence problem.

(49) Secondly, the single-tone parameter estimation scheme is extended to the multiple-tone case. The suitable number of decaying sinusoidal waves, M, should be tuned as a tradeoff between number of parameters and the accuracy of the representation. However, the value of M cannot be determined in advance. A sequential way to estimate the parameters of the multiple-tone representation (1) without pre-determining the value of M is introduced as follows

(50) $\begin{array}{cc}{\theta }_i^{*}={\arg }_{{\theta }_i}\min \underset{k=1}{\overset{N}{.Math.}}{\left[x\left(k\right)-\underset{j=1}{\overset{i-1}{.Math.}}{s}_{{\theta }_j^{*}}\left(k\right)-s_{{\theta }_i}\left(k\right)\right]}^2.& \left(15\right)\end{array}$
where θ.sub.i* denotes the optimal estimate of θ.sub.i, which relies on the optimal estimates of the previous i−1 parameter sets θ.sub.1*, θ.sub.2*, . . . , θ.sub.i−1*.

(51) This calculation can be implemented on a given seismic trace x(k) as in the following steps. First, assume that the seismic trace is represented by only one decaying sinusoidal wave (M=1) as in (1). Then, apply the single-tone estimation method developed in the previous section to x(k) and obtain the first parameter estimate set θ.sub.1*. The signal residual up to this step is given as
r.sub.θ.sub.1*(k)=x(k)−s.sub.θ.sub.1*(k)   (16)

(52) Its signal energy is denoted by r.sub.θ.sub.1*.sup.2. Next, implement the single-tone estimation method on r.sub.θ.sub.1*(k) to obtain the parameter estimate θ.sub.2* and the second residual signal r.sub.θ.sub.2*(k) similarly. This procedure is repeated and the parameter sets are estimated component by component. At the ith iteration, the parameters of representation s.sub.θ.sub.i(k) are estimated from previous residual signal r.sub.θ.sub.i−1*(k), leading to the next residual as follows:

(53) $\begin{array}{cc}{r}_{{\theta }_i^{*}}\left(k\right)=x\left(k\right)-\underset{j=i}{\overset{i}{.Math.}}{s}_{{\theta }_j^{*}}\left(k\right).& \left(17\right)\end{array}$

(54) The process is terminated once the residual energy is lower than a predefined threshold r.sub.D, i.e.
r.sub.θ.sub.j*.sup.2<r.sub.E,   (18)

(55) Alternatively the process can be terminated once the residual energy dropping rate is lower than a predefined threshold r.sub.D, i.e.,

(56) $\begin{array}{cc}\frac{r_{{\theta }_{j-1}^{*}}^2-r_{{\theta }_j^{*}}^2}{r_{{\theta }_{j-1}^{*}}^2}100\%<r_D,& \left(19\right)\end{array}$

(57) These steps are summarized in TABLE 2.

(58) TABLE-US-00002 TABLE 2 Parameter Estimation for the Multiple Tone Case Step Operation a) Estimate θ.sub.1* from trace x[k] by implementing Algorithm I; b) Calculate residual r.sub.θ.sub.1.sub.*[k] as in (16); c) Next iteration i = i + 1 c) Estimate θ.sub.1* from residual signal r.sub.θ.sub.i-1.sub.*[k] by implementing Algorithm I; d) Calculate residual r.sub.θ.sub.1.sub.*[k] as in (17); e) If criteria (18) or (19) satisfies, then terminate, else go back to step c)

(59) The data volume of K original seismic traces in L-bit data format is V=NKL. Suppose that representing each trace with M representation components squeezes the signal variance to σ.sub.cp.sup.2, which is smaller than the variance of the original traces, then entropy coding of those residuals along with the 6M parameters coded in L bits yields a volume of
V.sub.cp=KN(log√{square root over (2πe)}σ.sub.cp−log ε)+6MKL,   (20)
where ε denotes the digitization interval of geophones. Then the compression ratio of the developed method by using entropy coding is

(60) $\begin{array}{cc}\begin{array}{c}{\mathrm{CR}}_{\mathrm{cp}}=\frac{V}{V_{\mathrm{cp}}}\\ =\frac{L}{\log \sqrt{2\pi e}{\sigma }_{\mathrm{cp}}-\log \mathrm{.Math.}+\frac{6\mathrm{ML}}{N}}.\end{array}& \left(21\right)\end{array}$

(61) Another way to deal with the representation residuals is quantization with fixed-length codewords, which results in a lossy compression. Suppose that each residual sample is coded as l=log.sub.2L.sub.q, where L.sub.q is the number of codewords. The compression ratio considering M components and L.sub.q-level residual quantization is calculated as:

(62) $\begin{array}{cc}{\mathrm{CR}}_{\mathrm{rq}}=\frac{\mathrm{NL}}{6\mathrm{ML}+\mathrm{Nl}}.& \left(22\right)\end{array}$
Although the residual quantization is expected to yield lower reconstruction quality, the compression ratio CR.sub.rq is higher than CR.sub.cp in (21).

(63) In a non-limiting example, if the number of codewords is set as L.sub.q=128, then each residual is quantized as l=log.sub.2L.sub.q=7 bits. Similarly, if the number of codewords is L.sub.q=256, then each residual is quantized as 8 bits.

(64) In a non-limiting example of scaling, for a length of the original seismic trace of N=4000, each of the samples is L=32 bits, requiring a storage space of NL=128000 bits per trace. If the trace is represented by M=10 models, and each residual sample is quantized with L.sub.q=128 codewords, the storage space needed for the compressed trace is 6ML+Nl=29920 bits. The compression ratio is

(65) $\frac{128000}{29920}=4.27.$
If the trace is represented by M=5 models, and each residual is quantized as L.sub.q=64 codewords, the storage space needed for the compressed trace is 6ML+Nl=24960, thus the compression ratio is

(66) $\frac{128000}{24960}=5.13.$

(67) In the system, method and non-transitory computer readable medium, the computing system 460 further comprises a network controller connected to a transmitter 470 which transmits the stored encoded residual energy and the stored set of parameters to a remote processing station 380 for analysis and reconstruction as shown in FIG. 4.

(68) The compression method of the present disclosure offers a flexible and robust compression-decompression scheme compared with the LPC. First, the method of the present disclosure can reconstruct the signal by only using estimates of the initial parameters. Although it may result in a lower reconstruction quality in some applications, it achieves a very high compression ratio and does not have any risk of reconstructing divergence. Secondly, the compression ratio and reconstruction quality can be scaled by adjusting the number of representation components and the level of quantization.

(69) Furthermore, the compression method of the present disclosure is compared with the LPC method described in East Texas database. See W. A. Mousa and A. A. Al-Shuhail, “Processing of seismic reflection data using matlab,” Synthesis Lectures on Signal Processing, vol. 5, no. 1, pp. 1-97, 2011, incorporated herein by reference in its entirety. This database contains 594 real traces, each of which has 1501 time samples with a sampling frequency 500 Hz. An LPC is designed with an estimate of horizontal length p=6, 12, 18, 24 to achieve almost same compression ratios of the proposed methods with M=2, 4, 6, 8. Three scalar quantizors, 64-level for 6 bits per sample, 128-level for 7 bits per sample and 256-level for 8 bits per sample, are designed using a k-means algorithm and then are used to compress the LPC residual e′.sub.L(k) as in (28) and the residual of the present invention. See K Sayood, Introduction to Data Compression, Third Edition. Elsevier, 2006, incorporated herein by reference in its entirety. The LPC with different values of p and different quantization levels is applied to the database, leading to around 20% reconstructed traces significantly diverging in all cases. The normalized reconstruction error energy, calculated from those converged traces, are plotted in FIGS. 5A, 5B, 5C for level-64, level-128 and level-256 quantization, respectively. Meanwhile, the corresponding results of the proposed method with same compression ratios are plotted in FIGS. 5A, 5B, 5C for comparison. Note that the LPC with p=6, 12, 18, 24 has almost same compression ratio as the proposed method with M=2, 4, 6, 8. Two facts are observed from this result. First, the reconstruction error of the LPC is not sensitive to changes in p, which implies that the performance of LPC is not easily scaled by choosing the value of p. Secondly, compared with performance of the present invention, the reconstruction quality of the converged traces by the LPC, is much lower with almost the same compression ratio.

(70) In a non-limiting example, the methods and system of the present disclosure may be used by data collecting vehicles on a seismic acquisition field. Once the seismic data has been recorded in the vehicle, the data is compressed before being transmitted to another location for processing. This will greatly decrease either the requirement of the storage device capacity, or the frequency of transmission.

(71) In a further non-limiting example, the methods and system of the present disclosure may be used with wireless geophones. The wireless geophones compress the traces received before sending them to a data center. This will greatly decrease the amount of data transmitted. Additionally, since wireless communication is energy consuming and the wireless geophones are powered with batteries, saving energy for wireless communication means increasing their working longevity.

(72) Next, a hardware description of the computing device according to exemplary embodiments is described with reference to FIG. 7. In FIG. 7, the computing device includes a CPU 700 which performs the processes described above/below. The process data and instructions may be stored in memory 702. These processes and instructions may also be stored on a storage medium disk 704 such as a hard drive (HDD) or portable storage medium or may be stored remotely. Further, the claimed advancements are not limited by the form of the computer-readable media on which the instructions of the inventive process are stored. For example, the instructions may be stored on CDs, DVDs, in FLASH memory, RAM, ROM, PROM, EPROM, EEPROM, hard disk or any other information processing device with which the computing device communicates, such as a server or computer.

(73) Further, the claimed advancements may be provided as a utility application, background daemon, or component of an operating system, or combination thereof, executing in conjunction with CPU 700 and an operating system such as Microsoft Windows 7, UNI7, Solaris, LINU7, Apple MAC-OS and other systems known to those skilled in the art.

(74) The hardware elements in order to achieve the computing device may be realized by various circuitry elements, known to those skilled in the art. For example, CPU 700 may be a Xenon or Core processor from Intel of America or an Opteron processor from AMD of America, or may be other processor types that would be recognized by one of ordinary skill in the art. Alternatively, the CPU 700 may be implemented on an FPGA, ASIC, PLD or using discrete logic circuits, as one of ordinary skill in the art would recognize. Further, CPU 700 may be implemented as multiple processors cooperatively working in parallel to perform the instructions of the inventive processes described above.

(75) The computing device in FIG. 7 also includes a network controller 706, such as an Intel Ethernet PRO network interface card from Intel Corporation of America, for interfacing with network 77. As can be appreciated, the network 77 can be a public network, such as the Internet, or a private network such as an LAN or WAN network, or any combination thereof and can also include PSTN or ISDN sub-networks. The network 77 can also be wired, such as an Ethernet network, or can be wireless such as a cellular network including EDGE, 3G and 4G wireless cellular systems. The wireless network can also be WiFi, Bluetooth, or any other wireless form of communication that is known.

(76) The computing device further includes a display controller 708, such as a NVIDIA GeForce GT7 or Quadro graphics adaptor from NVIDIA Corporation of America for interfacing with display 710, such as a Hewlett Packard HPL2445w LCD monitor. A general purpose I/O interface 712 interfaces with a keyboard and/or mouse 714 as well as a touch screen panel 716 on or separate from display 710. General purpose I/O interface also connects to a variety of peripherals 718 including printers and scanners, such as an OfficeJet or DeskJet from Hewlett Packard.

(77) A sound controller 720 is also provided in the computing device such as Sound Blaster 7-Fi Titanium from Creative, to interface with speakers/microphone 722 thereby providing sounds and/or music.

(78) The general purpose storage controller 724 connects the storage medium disk 704 with communication bus 726, which may be an ISA, EISA, VESA, PCI, or similar, for interconnecting all of the components of the computing device. A description of the general features and functionality of the display 710, keyboard and/or mouse 714, as well as the display controller 708, storage controller 724, network controller 706, sound controller 720, and general purpose I/O interface 712 is omitted herein for brevity as these features are known.

(79) The exemplary circuit elements described in the context of the present disclosure may be replaced with other elements and structured differently than the examples provided herein. Moreover, circuitry configured to perform features described herein may be implemented in multiple circuit units (e.g., chips), or the features may be combined in circuitry on a single chipset, as shown on FIG. 8.

(80) FIG. 8 shows a schematic diagram of a data processing system, according to certain embodiments, for performing the functions of the exemplary embodiments. The data processing system is an example of a computer in which code or instructions implementing the processes of the illustrative embodiments may be located.

(81) In FIG. 8, data processing system 800 employs a hub architecture including a north bridge and memory controller hub (NB/MCH) 825 and a south bridge and input/output (I/O) controller hub (SB/ICH) 820. The central processing unit (CPU) 830 is connected to NB/MCH 825. The NB/MCH 825 also connects to the memory 845 via a memory bus, and connects to the graphics processor 850 via an accelerated graphics port (AGP). The NB/MCH 825 also connects to the SB/ICH 820 via an internal bus (e.g., a unified media interface or a direct media interface). The CPU Processing unit 830 may contain one or more processors and even may be implemented using one or more heterogeneous processor systems.

(82) For example, FIG. 9 shows one implementation of CPU 830. In one implementation, the instruction register 938 retrieves instructions from the fast memory 940. At least part of these instructions are fetched from the instruction register 938 by the control logic 936 and interpreted according to the instruction set architecture of the CPU 830. Part of the instructions can also be directed to the register 932. In one implementation the instructions are decoded according to a hardwired method, and in another implementation the instructions are decoded according a microprogram that translates instructions into sets of CPU configuration signals that are applied sequentially over multiple clock pulses. After fetching and decoding the instructions, the instructions are executed using the arithmetic logic unit (ALU) 934 that loads values from the register 932 and performs logical and mathematical operations on the loaded values according to the instructions. The results from these operations can be feedback into the register and/or stored in the fast memory 940. According to certain implementations, the instruction set architecture of the CPU 830 can use a reduced instruction set architecture, a complex instruction set architecture, a vector processor architecture, a very large instruction word architecture. Furthermore, the CPU 830 can be based on the Von Neuman model or the Harvard model. The CPU 830 can be a digital signal processor, an FPGA, an ASIC, a PLA, a PLD, or a CPLD. Further, the CPU 830 can be an x86 processor by Intel or by AMD; an ARM processor, a Power architecture processor by, e.g., IBM; a SPARC architecture processor by Sun Microsystems or by Oracle; or other known CPU architecture.

(83) Referring again to FIG. 8, the data processing system 800 can include that the SB/ICH 820 is coupled through a system bus to an I/O Bus, a read only memory (ROM) 856, universal serial bus (USB) port 864, a flash binary input/output system (BIOS) 868, and a graphics controller 858. PCI/PCIe devices can also be coupled to SB/ICH 888 through a PCI bus 862.

(84) The PCI devices may include, for example, Ethernet adapters, add-in cards, and PC cards for notebook computers. The Hard disk drive 860 and CD-ROM 866 can use, for example, an integrated drive electronics (IDE) or serial advanced technology attachment (SATA) interface. In one implementation the I/O bus can include a super I/O (SIO) device.

(85) Further, the hard disk drive (HDD) 860 and optical drive 866 can also be coupled to the SB/ICH 820 through a system bus. In one implementation, a keyboard 870, a mouse 872, a parallel port 878, and a serial port 876 can be connected to the system bus through the I/O bus. Other peripherals and devices that can be connected to the SB/ICH 820 using a mass storage controller such as SATA or PATA, an Ethernet port, an ISA bus, a LPC bridge, SMBus, a DMA controller, and an Audio Codec.

(86) Moreover, the present disclosure is not limited to the specific circuit elements described herein, nor is the present disclosure limited to the specific sizing and classification of these elements. For example, the skilled artisan will appreciate that the circuitry described herein may be adapted based on changes on battery sizing and chemistry, or based on the requirements of the intended back-up load to be powered.

(87) The functions and features described herein may also be executed by various distributed components of a system. For example, one or more processors may execute these system functions, wherein the processors are distributed across multiple components communicating in a network. The distributed components may include one or more client and server machines, which may share processing, as shown on FIG. 6, in addition to various human interface and communication devices (e.g., display monitors, smart phones, tablets, personal digital assistants (PDAs)). The network may be a private network, such as a LAN or WAN, or may be a public network, such as the Internet. Input to the system may be received via direct user input and received remotely either in real-time or as a batch process. Additionally, some implementations may be performed on modules or hardware not identical to those described. Accordingly, other implementations are within the scope that may be claimed.

(88) The above-described hardware description is a non-limiting example of corresponding structure for performing the functionality described herein.

(89) Obviously, numerous modifications and variations of the present invention are possible in light of the above teachings. It is therefore to be understood that within the scope of the appended claims, the invention may be practiced otherwise than as specifically described herein.