SEISMIC PROCESSING APPARATUS

Abstract: A seismic processing apparatus comprises a plurality of processing modules (12), each module comprising a transputer (T) and a transversal filter (TF) with these transputers of the modules connected in a linear array by their links. The apparatus is able to carry out real-time demultiplexing of the seismic data from a plurality of receiver channels, and is able to carry out a number of other processing functions in the field, e.g. common depth point stacking, summing, correlation, convolution and two dimensional Fourier transforms.

Description: SEISMIC PROCESSING APPARATUS FIELD OF THE INVENTION This invention relates to a seismic processing apparatus and in particular to a portable apparatus capable of carrying out real time processing of data during acquisition and subsequent reprocessing of acquired data. BACKGROUND OF THE INVENTION The exploration for hydrocarbons and minerals at present is becoming of widespread importance, as the earth's resources are being consumed at an alarming rate. The seismic exploration method has proved very successful over the last 65 years for locating hydrocarbons. Today there is a great demand to adapt these techniques to the search for minerals, which places greater demands on the processing of data. The basic principle of seismic exploration consists of transmitting a sound wave through the earth's rock structures and picking up reflected and/or refracted waves. Analysing the properties of these sound waves yields information about the earth's geological features. These features, together with samples taken from exploratory bore holes, enable the geologist to determine the position and extent of hydrocarbon and mineral deposits. In land exploration, a large number of geophones are spaced apart along a line extending over many kilo¬ metres, in either one or both directions, from a sound source. In a typical land data acquisition arrangement, the geophones are each connected to a recording truck via independent channels, which may consist of multi- conductor cables carrying analogue signals, or a single cable carrying a digital multiplexed signal from each geophone * At the recording truck the signals from each geophone channel are connected to a switching device which selects which channels are to be recorded. The sound source e.g. an explosive charge, is moved along the line at fixed intervals of distance, and the geophone channels selected for recording are those in the proximity of this sound source. At each interval or "shot point" a recording is made and in this manner data is acquired while moving along the exploration line. In marine exploration the reflected and refracted sound waves are detected using a hydrophone streamer which is towed in a straight line behind a survey ship. The streamer consists of a number of sections, each containing a hydrophone (or group of hydrophones) , and are linked to the recording equipment via independent channels. The ship moves along the exploration line towing the sound source, which produces a sound wave at fixed intervals of distance. The sound wave passes through the water into the earth's surface and the reflected and refracted waves are detected by the hydrophone streamer and recorded. In this manner data is acquired while the survey ship moves along the exploration line. At present, digital recording is used almost exclusively in seismic exploration. This gives a large dynamic range, and signal processing techniques can be applied to the data directly using digital computers. Typical recording equipment used at present has a capacity of up to 120 channels, sampling rates in the range 0.25 to 4 milliseconds, and a dynamic range of 84 dB. The analogu to-digital conversion generally consists of 14 bits for th magnitude, one bit for the sign and three bits to represen the gain applied to the amplifier during conversion. Data from each geophone is sampled sequentially, digitised and placed onto tape (i.e. in a multiplexed form). As sound sources, explosives are used extensively for both land and marine exploration, but are now beginnin to be replaced by other sources. For example, a weight drop method is often used in land exploration, in which a weight is dropped from a predetermined height onto the ground to produce the sound wave. The weight can be dropped many times at a particular location, and the re¬ flected seismic signals detected by the geophones for each drop can be added to increase the signal amplitude. A method called vibroseis is frequently used for land exploration in which a mass in contact with the ground is made to vibrate with a frequency that is continuously changed or swept over a certain range. The swept fre¬ quency sound wave travels into the earth's surface and is reflected and refracted. Since a swept frequency of long duration is used, as opposed to a single shock or pulse, the detail from each reflecting structure is blurred. The detail is recovered using a process called correlation, which compresses the swept frequency wave¬ form into a pulse. Thus, the signal from each geophone is cross-correlated with the signal or reference signal applied to the vibrating mass. The result is a record of the type obtained had a pulse been sent through the earth's surface. The vibroseis method has the advantage that there is control over the frequencies that are trans¬ mitted. The ultimate aim of seismic data processing is to produce an image (a section) that is representative of geological features. In practice the final section that is produced is not a perfect representation of rock strata, and it has to be subjected to interpretation by experienced geologists and geophysicists. The type of processing aplied to the seismic data is very much dependent upon the type of exploration and the geology of the area. For example the exploration for oil, gas, coal and minerals requires different processing. At present the infield processing of seismic data signals from the geophones is limited to filtering, noise rejection,, summing and correlation. Portable data acquisition systems are known, which can be hand-carried into jungle regions or mounted in small trucks for desert use. Arrangements for summing of seismic data used with the weight drop method of exploration, and typical vibro¬ seis data acquisition arrangement, which will perform in¬ field correlation and summing, are large and can only be used mounted in large trucks. The summer and correlator are purpose built units and cannot be adapted to perform other functions. These various systems and arrangements provide data from the various channels recorded in multi¬ plexed form on tapes. The recorded data tapes from the field are sent to processing centres, where large main frame computers per¬ form all the processing requirements to produce the seismic sections. The field tapes are first demultiplexed, which results in all the samples for channel 1 being separated for the first sound source position (shot point), followed by all the samples for channel 2, and so on. This is repeated for each shot point. The data is changed back from the geophones-sampled-sequentially form to channel-sequential form. This process on a VAX main frame computer is very time consuming and could be in terms of hours for a typical day's field data. The second stage of processing is concerned with signal extraction and signal-to-noise enhancement. It is at this stage that noisy channels and spikes on the data are edited out. Data can also be summed or cross correlated, if this has not already been carried out in the field. The third stage is geometric corrections. Survey data is placed into the computer at this point to provide corrections for the earth's terrain variations.. Also included at this stage is sound velocity information, this being determined from refraction data. The fourth and final stage is concerned with data enhancements, and it is at this point that very sophis¬ ticated techniques can be employed. The data can be filtered in two dimensions and can be used for the removal of noise generated by ground roll, which is a slow moving wave travelling in a shallowregion (or "weathering layer") just below the earth's surface. The waveform is reshaped, using deconvolution, and compensates for changes in the shape of the sound wave as it passed through the rock layers. A process called stacking is almost universall carried out which reduces the noise from multiple re¬ flections. Trace amplitude equilization can also be applied to compensate for the weaker signals furthest from the sound source, time variable filtering to com¬ pensate for the attenuation at high frequencies, mig- ration to transform the apparent reflection positions to true positions, which is useful near regions of geo¬ logical faulting, and many other processes can be used, the final result being the required seismic section. At present, data acquisition equipment tends to be in the form of separate units, each purpose built, and bolted together to form a recording system. Often - the units come from different manufacturers, requiring appropriate interfacing circuitry to be built. This results in greater requirements for spare parts in the field, and greater expertise in personnel to repair the equipment if it should fail. Over the years there has been a demand for more channels to be recorded. This has led to smaller identical units being added together to achieve greater channel capacity. These units were originally designed for independent use and using them in this way tends to lead to technical problems, such as longer sampling rates. The units cannot be bolted together in unlimited numbers. SUMMARY OF THE INVENTION In accordance with this invention, there is provided a seismic processing apparatus which comprises a plurality of processing modules, each module comprising a transputer and a transversal filter and the transputers of the modules being connected in a linear array. Preferably also each module further comprises an external memory for the trans¬ puter. In this seismic processing apparatus, the pro¬ cessing modules are preferably identical > reducing spare part requirements and making repair an easier task. The apparatus may -be made portable yet Is able to carry out an immense amount of processing of seismic data in the field, as will be explained. Further, if additional channels are required, it is a simple matter to add the appropriate number of processing modules, one for each channel. particular processing function which the apparatus can perform is real time demultiplexing of sampled data input from the geophones or the like. The data which is input to the first transputer of the array is in channel-sequential order (for each shot point) and the demultiplexing serves to place the data for the different channels into the memories of respective trans¬ puter modules. This demultiplexed data can be further processed subsequently on the same linear array of trans¬ puters. A linear array of transputers may be arranged to demultiplex data other than seismic data and thus, in accordance with this invention, there is also provided a data demultiplexing apparatus comprising a linear array of transputers. The stacking process may also be carried out on the linear array of transputers and further in accordance with this invention there is provided an apparatus for stacking values on a succession of channels in respect of a succession of instants, comprising a linear array of transputers arranged to receive the channel data at one end of the array and provide the stacked data at the other end of the array. A further process which is commonly required in seismic data processing is a 2-dimensional Fourier trans¬ form. This may also be carried out by the linear array of transputer and further in accordance with this inventio there is provided an apparatus for performing a 2- dimensional Fourier transform, comprising a linear array of transputers. BRIEF DESCRIPTION OF THE DRAWINGS An embodiment of this invention will now be des¬ cribed by way of example only and with reference to the accompanying drawings in which: FIGURE 1 is a general schematic view of a seismic processing apparatus in accordance with the invention; FIGURE 2 is a schematic plan view of a processing module board of the apparatus; FIGURE 3 is a diagram of the main elements of a transputer; FIGURE 4 is a diagram showing interconnection of transputers in an array; FIGURE 5 is a diagram of a transversal filter; FIGURE 6 is a diagram illustrating the operation of the transversal filter; FIGURE 7 is a diagram showing the operation of a linear array of transputers in demultiplexing input seismic data; FIGURE 8 is a series of diagrams illustrating the requirements for CDP stacking; FIGURE 9 is a series of diagrams illustrating the principle of filtering by convolution; FIGURE 10 is a series of diagrams illustrating the 2-dimensional Fourier transform of seismic data; and FIGURE 11 is a series of diagrams illustrating the operation of the linear array of transputers in carrying out the 2-dimensional Fourier transform. DESCRIPTION OF THE DRAWINGS The Seismic Processing Apparatus Referring to Figure 1, there is shown a seismic processing apparatuswhich is portable and able to carry, out an immense amount of processing of seismic data in the field. In the example shown, the apparatus comprises a housing 10 divided into compartments or racks for receiv¬ ing vertically-disposed circuit boards. At one end, an input interface board 11 will be placed, for receiving multiplexed digital signals from analog-to-digital con- verters which receive signals from the geophones or the like. There then follows a series of processing module boards 12, followed by an output interface board 14 for providing output to e.g. a visual display, a printer or a tape recorder. The apparatus also incorporates a power supply, for example at 16. Referring to Figure 2, there is shown one of the processing module boards 12. This carries a transputer T, a transversal filter TF, clock circuitry CLC, control cir- cuitry CC, RAM memory M (e.g. 512 K bytes) and memory interface and buffers IB. The board is further provided with edge connectors EC for mating with connectors within the host apparatus. In the example shown, the transputer T is the T414 device and the transversal filter the A100 device, both manufactured by INMOS of Bristol, England. The Transputer The transputer is one of the latest developments in very large scale integrated circuit technology and forms a new bulding block from which computer systems can be constructed. Figure 3 shows use of the transputer T in conjunction with external memory M and a clock oscillator CL to form a basic processor. The transputer is a pro- gramable component and the T414 device contains a 32 bit micro processor MP to perform high speed arithmetic and logic functions at up to 10 million instructions per second; an internal memory MT of 2 K bytes for very fast processor use, up to 80 M bytes per second data rate; a memory interface MI so that each transputer can extend its memory capacity up to 4 G bytes, and four links L to provid connections for data transfer to and from the processor. Also included are various control and system service functions SS for processor initialisations and error detection. A single processor can be linked to other pro- cessors to form a very sophisticated computer system and Figure 4 shows how this interlinking is made between transputers Tl, T2, T3. Communication can be established between other processors via the transputers links L. There are four such links on each transputer, each pro- viding two channels one input and one output, data being fed along the channels in a serial form. The data is transferred from one processor to another via these links, under software control, at a rate of 20 million bits per second. This communication is synchronized and takes place only when both processors are ready, one to receive data and the other ready to transmit it. The transputer is provided with facilities for initialisation and programs can be loaded into its memory MT via any one of its links L. Programs that have already been stored in its memory area, as ROM, can be executed directly. Error detection is possible which may be used to indicate a faulty transputer within a multi-processor array. OCCAM programming OCCAM is a high level computer language developed by INMOS to be used with their transputer. Compiled OCCAM is very efficient, removing the need for a transputer assembly language. This enables users of-transputer based systems to maintain a "standard language" which is indepen- dent of any changes in the design and structure of the transputer. Programs do not have to be adapted to these changes, nor do programmers have to relearn instruction sets and procedures. OCCAM is one of the few languages that enable many processors to be linked together, each programmed to per¬ form a separate task. These tasks could be small parts of a large mathematical or data processing problem * By allocating each processor to a part of the overall problem and executing these concurrently, the problem is solved efficiently. This is known ' as parallel processing. A single transputer can be programmed, using OCCAM, to perform concurrent, or parallel processing. This Is achieved by allowing a short time to work on each separate task in turn, until all tasks have been com- pleted. In this way a program written to be used on an array of transputers can also be executed unchanged by a single processor. In problem solving, the syntax of OCCAM allows the programmer a choice between parallel processing, where many tasks are performed concurrently, or sequential processing, where each task has to be completed before the next can begin. Parallel processes can be distributed over many processors or contained within a single processor. In both cases, if a process needs to cooperate with another during computation, it is necesary to establish a commun¬ ication link between them. OCCAM provides this linking by the formation of named channels and data is passed across these channels. The action of data transfer may be given the highest priority over all other processes, if desired. Using OCCAM programming arrays of processors can be set up, each communicating with its neighbouring pro¬ cessor, so that data can be made to "flow" steadily through each processing element in turn. As the data flows through a processor it can be operated on: this is called systolic processing. The Transversal Filter The transversal filter is another of the latest developments by INMOS in very large scale integrated cir¬ cuit technology and forms a building block from which advanced digital signal processors can be constructed. The transversal filter, incorporated into suitable hard¬ ware and under the control of software programming, can be made to perform many different digital processing functions at high speed. These functions include filter¬ ing, correlation, convolution, deconvolution and discrete Fourier transforms. Figure 5 shows the INMOS A100 transversal filter. Data bytes, input serially at I, are placed into the filter by direct access to the registers Rl, R2, R3. The multiplier accumulator array MAA consists of 32 parallel multipliers, each capable of multiplying two 16 bit binary numbers. Each product that results is passed through a cascade of delay and adders CDA and output serially through an output register R4 * Three stages of the transversal filter are shown in Figure 6, comprising the parallel multipliers and a cascade of delay and adders. There are a total of 32 stages on each device. The numbers x , x-, and 2 rep¬ resent data placed in the active coefficient register. The stream of numbers y , y, and 2 represent data being input to the device. At the first clock cycle y is in the input data register, y is then multiplied by x , X j and x« in parallel and their products are output to the cascade of delay and adders. At the second clock cycle y-, is placed into the input data register and this is multi¬ plied In parallel by x , x-. and x 2 , these products then join the data being cascaded through the delay and adders. The result is a stream of data being output as shown in • Figure 6. The output from the cascade is calculated to full 36 bit precision, to remove accumulative rounding errors. This is reduced by selecting 24 bits from these 36 for the final output. A programmable barrel shifter is used to make this selection. Two coefficient registers are provided, an active register R3 and an update register R2. The update register can be accessed while the coefficients in the active register are being applied to the multipliers. When a new set of coefficients are to be loaded, the two sets of registers can be interchanged using a "bank swap" instruction. The transversal filter can be cascaded to enable a greater number of parallel multiply-accumulate operations to be performed during any one clock cycle. The main criteria to consider with the design of a signal process- ing machine is its channel capacity and processing speed. With the advent of the transputer and transversal filter a new design philosophy can now be considered. The per¬ formance can be a function of the number of these devices that are present in a machine. The principle of parallel processing allows a designer to take a signal processing task and break this down into operations that can be per¬ formed concurrently, on as many individual processors that exist within the machine. With suitable programming it is only the speed of operation or channel capacity that will be affected by the number of processors present in the machine. The channel capacity can be extended without de¬ grading the speed of operation by ensuring that where possible all signal processing algorithms are capable of being executed independently, and therefore able to be performed in parallel. The addition of extra processors for each additional channel will then have no effect upon the overall operational speed of the machine, since all operations occur concurrently. The ideas of systolic processing, that is where data streams are made to "flow" through processors, is readily adapted to the parallel processing principle. Since data flow rate is determined by the transputers link speed, and is fixed for a particular type of trans- puter, the data flow rate will be unaffected by the number of transputers in a data stream. This is because they all act at the same time to transfer data in to and out from a link (i.e. in parallel). If processing is to be performed on this data, the transputers can execute this concurrently. The speed of processing can also be increased by the addition of extra processors, for a fixed channel capacity. The allocation of two transputers for each channel will enable faster or more complex processing to be carried out. The addition of extra transversal filters, cascaded together, will provide a greater number of parallel multiply-accumulate operations, during a time interval, thus giving rise to greater processing speed. Thus in the seismic processing apparatus of Figure 1, the transputers are connected in a linear array to provide for systolic processing and OCCAM programming can be employed. Demultiplexing The sampled data received in multiplexed form is readily demultiplexed using the linear array of transputers. The data samples are made to "flow" through each trans¬ puter and the channels separated and this procedure is shown in Figure 7 for the example of three channels. After the data has been,separated into its allocated transputer T1,T2,T3, one for-each channel CH1,CH2,CH3, it can then be processed or stored directly in memory EM1,EM2,EM3. The speed at which data can be transferred in and out of a transputer is 20 million bits per second, and data written into memory in typically less than a micro second, depending on the external memory speed. These time periods are Insignificant compared to the intervals between data samples. Sampled data from most analogue to digital con¬ verters used at present are in a SEG (Seismic Exploration Geophysical) standard * This consists of 14 bit magnitude, one sign bit and a three bit code to represent the gain applied to the amplifier during the conversion. The SEG standard data may have to be converted to another form, such as 32 bit floating point, to make it compatible with other processes. This conversion can be made within each transputer. Analysis of the data can also be undertaken, such as noise rejection, where noisy channels are muted. To optimise performance parallel processing can be used, the data converted, noise rejected, and data loaded to memory or routed to other transputers all concurrently. Figure 7 shows data for each channel being separated into its allocated transputer. Additional channels can be incorporated by providing more transputers, one for each channel, in the data flow stream. The addition of extra transputers will not have any effect on the operating speed since data is passed across all the transputers concurrently. If processing in the transputer is required the overall processing time will be that for a single transputer to perform this operation, as this same operation is executed independently and concurrently on all the other transputers. The channel capacity can be extended at the cost of operational speed by allocating one trans¬ puter to two or more channels, as shown in Figure 7a_. The speed of processing, or more complex processing, can be increased by allowing two or more transputers to each channel and this is shown in Figure 7b. Subsequent to demultiplexing (and in some cases at the same time as demultiplexing) further processing of the acquired data may be carried out by the seismic processing apparatus of Figure 1. These various processing functions will now be described. Geometric Corrections The reflections from rock structures, detected at the earth's surface, are misaligned due to the source and receivers being located at different elevations and different distances apart. This misalignment must be corrected prior to- further processing. The two geometric corrections required are (1) a vertical correction to place the source and receiver on the same horizontal plane (since the earth's terrain varies) and (2) a horizontal correction to place the source and receiver on the same vertical plane (since the sound ray paths have travelled different distances) - To make these corrections the data for each channel is time-shifted. The first correction entails shifting the entire channel data by the same amount, either up or down, depending on the postiion of the datum plane relative to the source and receiver. This is called the terrain correction, and requires survey information. The second correction is a time-variant correction in which the channel data is shifted by varying amounts along its length. This is called the normal move out correction, and requires infor¬ mation from refraction crews of the sound wave velocities. In the apparatus of Figure 1, the channel data is stored in memory, each sampled byte of data located at a specific address, within each transputer's memory address area. These addresses can be made contiguous, in which case their memory address number represents time. To shift a channel's data in time is then effected by re- addressing memory locations. The geometric corrections to the data can be implemented by programing each trans¬ puter to shift the data in its memory area. This pro¬ gramming can be executed concurrently i.e. on all_ trans¬ puters (all channels) at the same time. Common Depth Point Stacking This procedure is applied after geometric correc¬ tions have been made. A point below the earth's surface is defined (the common depth point or CDP). Reflected sound waves from this point, for various source and receiver positions, are added together or stacked. An example of this stacking process is shown in Figure 8. Six channels are recorded to the right of the shot point position, at the station intervals. Figure 8a ^ shows the shot point at STA 18, the first channel recorded (CHI) is on STA 19. The sound waves from the shot point are considered as rays. At point a. the sound ray incident upon it will be reflected, the angle of the reflected ray being equal to the angle of the incident ray. This reflected ray is detected at the receiver on station STA 23 and recorded on channel CH5. Likewise the reflected ray fromi point b will be detected at STA 24, channel CH6. Figures 8b to 8d show the incident and reflected sound rays for the common depth points a to d for the shot positions STA 19 to STA 21. The recorded channels contain information on re¬ flections from each of the points to d. The recorded channels which contain the data for a specific common depth point are added together (stacked). Thus to stack the data corresponding to the common depth point , the CH5 of shot point 18 is added to the CH3 of shot point 19 and this is then added to the CHI of shot point 20. Common depth point stacking provides a method by which noise of the same frequency can be separated from signal information, and is particularly useful to reduce the effects from multiple reflections. Table 1 shows how the linear array of transputers is used to perform the CDP stacking process, the method using systolic processing. Data in a demultiplexed and geometric corrected form is input to the transputer array. In this example, shot point STA 18 channels CHI to CH6 are first input. The channel data is separated and placed into their allocated transputer's memory, CHI in trans¬ puter 3, CH2 in transputer 4 and so on. When this has been completed for the shot point STA 18 the data contained in each transputer's memory is moved down the array by two transputer locations (shift by 2). The next shot point STA 19 data is then input to the transputer array, and the channels separated as before. The data corresponding to shot point STA 19 channels CHI to CH6 are then added to the data already in the transputer's memory, corres¬ ponding to the previous shot points data, which has been (Shot point 20 CHI to 6 etc. CONTENTS OF MEMORY AT VARIOUS TIMES (Shot point 19 CH1 to 6 (Shot point 18 CHI to 6 INPIΠ ' INPUT SHIFT ADD IN SHIFT ADD IN SHIFT ADD IN SH DATA SP18 DATA BY 2 SP19 DATA BY 2 SP20 DATA BY 2 SP21 DATA B SP18:CH6 [~SP19 CH6] r " SP20;CH6l r~SP21:CH6 ~ | SP18:CH5 [ " SP19 ; CH5| rSP20:CH5~| | " SP21:CH5 T6 SP18:CH4 SP18:CH6 rSP19 CH4J rSP19:CH6~| f SP20 : CHl] f SP20 ; CHβl I SSPP2211::CCHH44 3l f " SP2 [+SP18:CH6JL -J |+SP19 CH6JL + — J |+SP20:CH6J T5 SP18:CH3 SP18:CH5 f SPl9:CH3|fSP19:CH5l f SP20 : CH3~| f SP20 ; CH5| f SP21 : CH3 " | rSP2 |+SP18:CH5JL + — L+SP19:CH5JL + — J |+SP20:CH5j T4 SP18:CH2 SP18:CH4 rSP19:CH3fsP19:CH4] SP20-CH2 rSP20:CH4| rSP21;CH2 " | rSP2 b-SPl8:CH4||+SPl8:CH6j +SP19:CH1 L+SP195CH6J +SP20:CHM |+SP2 +SP18:CH6_ l+SP19;CH6j This is This is stacked b_ stacked d T3 SP18:CH1 SP18:CH3 rSP19:CHΪ|rSP19:CH3~ ~ SP20;CHΪ SP21:CH1 L+SP18:CH3JI+SP18:CH5. +SP19:CH3 +SP18:CH5_ [ SP20:CH3| SP2 +SP19:CH5j +SP20 ; CH3 +SP2 +SP19-CH5. This is This is stacked a stacked c_ T2 SPl8:CH2f |fsP19:CH2 2 l1 (Stacked) (Stacked) (Sta l+SPl8:CHgJl+SPl8:CH4| -ι-SP19:CH2 ( b ) ( b ) ( +SP18:CM T1 SP18.CH1 " irsP19 (Stacked) (Stacked) (Sta +SP18:CH1][+SP18 +SP19:CH1 ( a ) ( a ) ( _+SPl8sCH3j OUTPUT OUTP STACKED Stac CDP's a Stac shifted by two locations. The data now contained in the transputer's memory is again moved down the array by two transputer locations (shift by 2). The next shot point data STA 20 is then in- put to the transputer array. This procedure is continued as shown in Table I reading from left to right. The values of stacked , stacked b, stacked £ and so on will be "flowing out" at the output end of the transputer array. In this way CDP stacking becomes a continuous process, data being input at the top of the array, and the stacked CDP data values output at the bottom. The length of the array will depend upon how many transputers are available in the apparatus. One trans- puter could be made to perform the function of two or more of those shown in Table I, however this would slow down the speed of operation. Signal Processing Most advanced signal processing operations can be reduced to a mathematical notion of the form:- processed val t) dt (1) Written in words this means: The value of the processing operation at a specific value of _t is the sum of all the products of the two functions f(n) and g(n+t) over all the values of n between n's lower and upper values. In this form n and are considered to vary continuously. For the seismic processor the data comprises sampled values of a continuous signal, and is thus dis¬ crete in nature. The expression then takes the form:- processed value at Jt where n and t are discrete variables. If n is considered to vary in integer values from 1 to N, then the expression can be written as the series:- processed value at t = f(l)g(l+t) + f(2)g(2+t) +.. «.f(N)g(N+t) If t_ is considered to vary in integer values from 0 upwards then this can be written as: processed value (at t=0) = fd)g(l) + • f(2)g(2) + f(N)g(N) processed value (at t=l) = f(l)g(2) + f(2)g(3) + ..f(tf)g(N+l) processed value (at t=2) = f(l)g(3) + f(2)g(4) + (N)g(N+2) ...» and so on There is thus a direct relationship between the mathematical notion involving the integration of two functions (equation 1) and the summation of a series of products of the two functions. Any signal processing operation represented in the form of equation 1 can thus be implemented by the summation of a series of products. . The transversal filter will implement this summation of the products of two functions f(n) and g(n+t) , if the coefficient register contains the f(n) terms and the g(n+t) terms are applied to the input. The output of the transversal filter is the processed values at t=o, t=l, t=2 etc. Summing If the weight drop method is used, this requires the addition (in each channel) of data samples received from successive weight drops at a particular location. Usually eight to twelve drops are used. This summing can be achieved after demultiplexing in the apparatus of Figure 1. Thus, once the data for each ' channel is stored in the respective memory, subsequent data for that channel in respect of the next weight drop is added directly to that previously stored. Alternatively, two sets of memory area may be allocated to the data, the first set for data already summed and the second set for the data being sample If the sampled data is checked and found valid it is then added to the first set. In this way should the incoming data be in error the weight drop need only be repeated once. Correlation Corelation is required when the vibreosis technique is used. In order to correlate two signals, for example a sampled seismic signal a and a pilot signal p_, the mathem¬ atical process required will be given by:- correlated value at time t ) dt Where the difference between the upper and lower limits of n is known as the depth of investigation. The form of this equation is identical to that of equation 1, thus it can be implemented directly on a trans¬ versal filter. The sampled pilot values pi, p2, p3 pN, are placed in the coefficient register and the sampled seismic data signals al, a2, a3 aN for the given channel are applied to the input of the respective transputer. Convolution To convolve two signals ai and b_ the mathematical process required is given by: convolved value at time t ) dt The form of this equation is not the same as equation 1. However, upon the substitution t=-Jt then it becomes the same as equation 1. The physical inter¬ pretation of this substitution is that the order of the data is reversed in time. Thus to implement convolution on a transversal filter it is necessary to input signal b_'s data in reverse order of time (sampling data first). Alternatively if the order of a.'s signal data is reversed, and p_'s signal data in normal time, then the effect is the same. To implement convolution on a trans¬ versal filter the sampled waveform al, a_2, 3, N is placed into the coefficient register in reverse order (i.e. aN, aN-1 ^2, l). The waveform bl, b_2, b3 is placed at the input. The seismic processing apparatus of Figure 1 is thus able to perform convolution on the data of the respective channels» Filtering in real time The apparatus of Figure 1 may be programmed to carry out filtering in real time. Filtering can be under¬ stood as a convolution process. The convolution theorem states: convolution in the time domain is equivalent to multiplication in the frequency domain. This theorem is used to perform real time filtering and the stages are shown In Figure 9. In the first stage the required filter charac¬ teristic is defined, Figure 9^. The impulse response of this desired filter characteristic is computed and for this example turns out to be a sine function, shown in Figure 9b_. Unfortunately, this waveform has infinite duration and for practical purposes has to be truncated. To do this a "window function" is used, for this example a rectrangular window function, shown in Figure 9c_. This produces a filter operator. Other window functions can be used resulting in more desirable filter characteristics (e.g. Hanning, Hamming and Bartlett windows). This filter operator is convolved with an incoming signal in real time, and the result is a signal that has been filtered as required, shown in Figure 9d_. Mathematically by convolving the signal with the filter operator in real time the frequency spectra of the incoming signal is multiplied by the frequency spectra of the filter operator (i.e. the desired filtering). Frequencies above 40 HZ are multiplied by zero, those b p low 40 ϋZ multiplied by 1. Thus the incoming signals spectra is modified in the frequency domain. In practice the filter operator will be loaded into the coefficient register of the transversal filter (in reverse order), and the sampled seismic signals applie to the input. This type of filtering belongs to the class of Finite Impulse Response filters (FIR), they have no poles and the number of zeros is equal to the number of co¬ efficients in the transversal filter's coefficient register The phase response is linear and has the advantage that the response is defined exactly,- and is reproducible from machine-to-machine, unlike LCR filter circuit elements. The filter operator need not be fixed, and changes to it can be made in time to produce a time varying filter. Alternatively, it can be automatically adjusted to optimise filtering in a specific application, for example a tracking notch filter to remove mains interference. Fourier Transform The apparatus of Figure 1 is preferably arranged to perform Fourier transform on the demultiplexed seismic data. The Fourier transform converts a signal waveform described as a time varying function to a description in terms of frequency components. Thus it is a time-domain to frequency-domain transformation. This is expressd mathematically by the relationship:- where f(t) represents a function of time (waveform magnitude at time _t seconds) and F(w) represents a function of frequency (the magnitude of the frequency component at w cycles «. per second) It is also possible to consider a waveform that is described as a position varying function (as opposed to time varying) » If this waveform, expressed as distance, is Fourier transformed the result is in terms of the number of cycles per unit of measurement, that is cycles per meter. This unit is given the name wavenumber. There are many ways in which the discrete Fourier transform can be implemented. The linear transputer array can be used for implementation of the Cooley-Tuckey radix-2 Fast Fourier Transform method. Qr the transversal filters and transputers can be used for the chirp-Z transform method In seismic exploration the data is composed of many channels of information, each channel containing a wave¬ form that varies in time. These channels each represent receivers located at a specific distance apart. This is shown in Figure 10a_. The data can be represented as form¬ ing a surface, with time in one direction and distance in the other. To perform many sophisticated signal process¬ ing operations on this data transform techniques can be used and one such transform Is the two dimensional Fourier transform. Figures 10b_ and 10c show the stages of this transform process. First the data for each channel, along the time axis, is transformed using the Fourier transform technique described before. The waveform becomes a function of frequency content, the time axis becomes a frequency axis. The Fourier transform is then taken for all the data points along the distance axis, for each interval of frequency. The distance axis becomes a wave- number axis. It is found that data transformed in this way renders itself very amenable to mathematical manipulation, which means that very complex processing tasks can be per- formed using simple operations. The transformed data having been manipulated can then be transformed back to its original time-distance form by the application of a further two dimensional Fourier transfom. The seismic processing apparatus of Figure 1 performs the Fourier transform on its linear array of transputers. Figure 11 shows the stages in implementation of a two dimensional Fourier transform. Stage 1 consists of loading the channel data into its allocated transputers memory (and may comprise the demultiplexing process described above). Stage 2 consists of implementing a Fourier trans¬ form on each channel independently and concurrently, using the transversal filter as a high speed processor. The transformed data is stored in memory. Stage 3 consists of transposing these data values between channels, the data corresponding to a particular value on the frequency axis and along the distance axis are arranged in order of increasing frequency. Transputer Tl now contains the values along the distance axis cσrres- ponding to a zero frequency. Transputer T2 contains the values along the distance axis corresponding to the first frequency value; and so on. This process can be carried out in parallel, each transputer placing a d value into a link, say first dl is placed into the links. All these values for dl for channels CHI to CH3 can then flow into transputer Tl memory, by passing them along the link paths connecting each transputer. Then d2 can be placed into the links and so on. Stage 4 consists of again using the transversal filter as a high speed processor to perform a Fourier transform on this transposed data that is stored in each transputer's memory area. The result will be the two- dimensional Fourier transform of the original data. Migration The objective of migration is to convert the angles of dip from the apparent dips displayed on the seismic section to the true dips of the geological bedding planes and to remove diffraction hyperbolae. The technique is to collapse these hyperbolae to a point. A recent development uses a process called Wave- number migration to carryout this task and this can be implemented on the apparatus of Figure 1 as follows. The seismic data time section is converted to a depth section by multiplying by the velocity of sound, using the trans¬ puters in parallel for the respective channels. Then a two-dimensional Fourier transform is performed using the transputers and transversal filters, as described above. Then the transputers are used for "mapping" to effect a wavenumber domain migration, again operating in parallel on the respective channels. Next a two-dimensional Fourier inverse transform is carried out, as described above, and finally a depth-section to time-section conversion is effected to produce the required migrated seismic section. Surface consistent geometric corrections A geometric correction that is sometimes performed after the application of earth terrain and normal moveout corrections is called surface consistent corrections. To perform this correction a cross-correlation is made between certain channels, which yields a time shift between adjacent channels. This time shift information can be use to make fine adjustments, and makes the data consistent over a range of stations. This correction can also be carried out on the linear transputer array of the apparatu of Figure 1. Summary The seismic processing apparatus which has been described is a portable apparatus able to carry out real time demultiplexing of the data acquired from the geophone or other receivers. At the same time as demultiplexing, certain other processing functions can be carried out by the transputers of the respective channels. Subsequently, the demultiplexed data can be further processed on the linear array of transputers, using the transversal filters when required, to perform e.g. correlation, geometric corrections, common depth point stacking, summing, filter¬ ing and Fourier transforms. The apparatus is capable of data acquisition and processing under the control of soft¬ ware programming usingthe minimum amount of hardware. Eac processing stage may be displayed for quality control and assessment. The final seismic section may be produced on a small printer in the field, and later reproduced on a high definition graphic printer at a central processing station. The apparatus thus enables a maximum of pro- cessing stages to be carried out in the field. It will be appreciated that each processing module such as shown in Figure 2 may include one or more trans¬ puters in addition to th one which is connected in the linear array, and may include one or more additional transversal filters.

Claims: CLAIMS 1. A seismic processing apparatus which comprises a plurality of processing modules (12), each module compris¬ ing a transputer (T) and a transversal filter (TF) with the transputers of the modules being connected in a linear array by their links (L) . 2. A seismic processing apparatus as claimed in claim 1, arranged for a first transputer of the array (e.g. T3 - Fig.7) to receive data in channel-sequential form from a plurality of seismic wave receivers and for said data to be passed along said linear array, and for said transputers to effect demultiplexing by passing data for the different channels to memories of respective trans¬ puters. 3. A seismic processing apparatus as claimed in claim 1, arranged for the linear array of transputers to perform common depth point stacking wherein a first trans¬ puter of the array receives data for a succession of shot points in sequence, the data for each shot point comprising data for a set of channels in sequence, the data for the first shot point is passed to respective transputers for. the different channels and loaded into their memories, then the data stored in each transputer memory is shifted to the memory of a transputer a predetermined number of transputers further down the array, and next the data for the second shot point is input to the first transputer and the data for the different channels is loaded into the memories of the respective transputers and added to the data already stored therein, and so on» 4. A seismic processing ' apparatus as claimed in claim 1, arranged to perform summing, correlation and/or convolution on said linear array of transputers. 5. A seismic processing apparatus as claimed in claim 1, arranged to perform a two-dimensional Fourier transform on te seismic data, wherein the data for the different channels is loaded into the memories of res¬ pective transputers and then a Fourier transform is per¬ formed by the transputers concurrently on the data in their respective memories, then the transputers transpose data between their memories before a second Fourier trans¬ form is performed by the transputers concurrently on the data now in their respective memories. 6. A data demultiplexing apparatus comprising a plurality of transputers (T) interconnected by their links (L) into a linear array, a first transputer being connected to receive multiplexed data and the transputers being arranged to pass the data along said linear array and effect demultiplexing by passing data for different channels to memories of respective transputers. 7. A data processing apparatus, comprising a linear, array of transputers arranged to receive data in respect of a number of channels at one end of the array, to carry out stacking of values of a succession of channels in res- pect of a succession of events, and provide the stacked data at the other end of the array. 8. A data processing apparatus, comprising a linear array of transputers arranged to receive data and per¬ form a two-dimensional Fourier transform thereon, wherein a first Fourier transform is carried out by the transputers concurrently on data held In their memories, then data is transposed by the transputers between their memories and then a second Fourier transform is carried out by the trans- puters concurrently on the data in their memories.

Inventors: ARMITAGE PETER ROBERT, ARMITAGE, PETER, ROBERT
Applicants: ARMITAGE PETER ROBERT, ARMITAGE, PETER, ROBERT
Application year: 1987
Main classification: G01V1
IPC codes: G01V1