Open main menu

Determining sensor orientation

A Matlab extension for Scream! lets you easily determine the exact orientation of a sensor relative to a surface reference sensor, which can be installed on an accurately-surveyed North/South line.

The extension uses the Blacknest correlation method to determine relative orientation. This method was developed at AWE Blacknest in collaboration with Güralp Systems.

Signals received by the N/S component of the reference sensor are correlated with those received at the N/S and E/W components of the sensor being studied, after different amounts of mathematical rotation. The highest correlation will occur when the N/S component of the reference sensor matches the rotated N/S component of the borehole sensor.

Once you know the deviation of the borehole components from the correct orientation, you can instruct the digitiser to rotate the signals algorithmically.

Installing the Scream! extension

The Relative Orientation extension is supplied in the standard distribution of Scream!.

The extension uses Matlab™ libraries. However, you do not need the full Matlab™ package to use the extension. The Matlab™ run-time libraries are also included in the Scream! distribution.

Installing the reference instrument

To measure the orientation of a sensor, you will need a second instrument which is known to point precisely North. For best results, use an instrument with a response identical to ro, at least, similar to the borehole instrument. The instrument should be located on a solid surface as close to the other instrument as possible. Most boreholes are constructed with a concrete base around the top of the borehole; If this is present, we recommend installing the reference sensor there.

Ideally, the two sensors should be directly connected to the same 6-channel digitiser. If you are using separate digitisers, you will need to ensure they are exactly synchronised. This can be done by connecting GNSS receivers to both digitisers and waiting for the control system of each one to settle. This process takes at least twelve hours.

This procedure only works with GCF recordings so, if you are orienting a Radian borehole system or using a Minimus digitiser, you should use a copy of our free Scream! software to record the data in real-time in GCF format. You cannot use the files recorded on the Radian SIU or Minimus.

Measuring the orientation

  1. Open the recorded files in Scream!. Under Windows, double-click on the .gcf files in Explorer. Under Linux, you can double-click on the files in Nautilus, Nemo etc if you have configured your system this way or, from the command line, you can run

    scream --view path_to_.gcf_file…

    A WaveView window will open displaying your recorded streams:

    Waveview orientation streams
  2. Drag the streams across the window so that the reference stream is at the top, the N stream in the middle, and the E stream at the bottom.

  3. Locate a suitable data range. A good range would contain lots of teleseismic events and very little local noise. You should use a period of at least an hour, and preferably much longer.

  4. Hold down shift and drag across the WaveView window with the left mouse button until all three streams are selected for the whole range. Make sure there are no gaps in the data you select. If gaps are present in any streams, the selection for those streams will be shown with hatched lines to alert you:

    Streams with gaps

    When you are happy with the selection, release the mouse button, but keep shift held down.

  5. When the context menu appears, release shift and select Relative Orientation from the menu:

    The orientation menu

    Two small windows will appear: a progress window:

    The busy indicator

    and, typically, a warning like:

    Matlab's warning dialogue

    Scream! produces this warning when the reference sensor is not using a standard N/S channel, i.e. it is using the auxiliary (X) channel instead. If you are using a separate digitiser, the warning will not appear.

    If you see a different error message, make sure that the streams are in the correct order in the WaveView window. If you still have problems, you may have selected too few data points for it to be confident about the orientation; you should try again with a larger selection.

  6. After a few seconds, the calculation should finish and three windows will appear. (They may obscure each other.) The top window is a graph of Coherence vs Angle:

    Plot of coherence versus angle

    The two-stage algorithm rotates the N/S and E/W components of the sensor being tested in small steps.It measures first the amplitude similarity, and then the coherence between the virtual (rotated) N/S component and the reference N/S component, for a number of rotation angles.

    The error in the final calculation is around 2.5°.

    The peak of the coherence curve (upper graph) therefore corresponds to the angle of rotation which best matched the reference component. This angle is shown in the title bar, together with an estimated error.

    You should see a coherence curve which is smooth and symmetrical. If the curve is distorted, either the surface data are too noisy or the data selection is too short.

    The lower graph shows the overall amplitude similarity of the rotated signal. This provides an idea of the sign of the coherence (since signals in perfect antiphase have a high coherence as well as those in phase). If there are two peaks in the coherence graph, the correct one is where the amplitude similarity is most positive.

    The sample plots show that the borehole instrument is installed with its N/S axis at a bearing of –90° from true North.

  7. The second window shows the result of applying the rotation to the signal, i.e. the time series that a sensor in perfect N/S orientation would have produced:

    Time-series comparison
  8. You can perform more accurate calculations by narrowing the search range. This is done in the two left-hand-side entry boxes on the Coherence vs Angle window: the first denotes the centre of the new search and the second specifies its range. The program suggests suitable values for you, so in most cases you can just click the 'calculate' button to perform another iteration.

  9. A new graph will be displayed showing the results:

    Second plot of coherence versus angle

    Our sample instrument is thus aligned at –90.61 ± 0.07°.

  10. The error given is only a rough estimate. For best results, you should repeat the orientation experiment several times using different data sets. The true error in the computed orientation can then be determined by observing the spread of the results.

The Blacknest orientation method generally provides a reliable indication of the sensor’s orientation. In most cases, the greatest source of error is in the installation of the reference sensor.

If you have particular difficulty in deriving a stable value, additional information is contained in a third output window, which shows a waterfall plot of Coherence vs Frequency vs Angle. The frequency range used in the calculations is indicated with a black outline.

Plot of coherence versus angle versus frequency

This plot can be used to select a more advantageous frequency range where the coherence curve is smoother and easier to interpret. In the example above, the frequency range 0.23 to 0.6 Hertz looks particularly promising: the surface within this range forms a smooth, symmetrical arch. Ideally, the peak should be very close to unity and the lowest points, at the margins, should be close to zero.

The chosen frequency range can then be entered into the Between X and Y Hz boxes at the top of the Coherence vs Angle window before clicking the 'Calculate' button again; the data will be filtered accordingly before the coherence is recalculated, increasing the accuracy of the result.

There are many factors which can affect the calculations used to determine sensor orientation. If you have trouble getting consistent, plausible results, please contact for advice.

Applying data rotations

Once you have determined the correct angle, you can program your digitiser to apply real-time mathematical rotations to the raw data. The procedure to do this is different for different types of digitisers. Use the links below to jump to the correct section for your digitiser:

DM24 Digitisers

You can configure a DM24 mk3 digitiser to apply a rotation to the digitised data. It can then produce output streams representing ground motion on true North/South and East/West axes.

This is done within the DSP to minimise the reduction in data quality.

To set up the rotation:

  1. Open a terminal session with the digitiser. You can do this over a serial link with a program such as minicom (for Linux) or PuTTY (for Microsoft Windows). If you have a serial connection to a PC running Scream!, you can access the digitiser's console by right-clicking on its icon and selecting Terminal… from the context menu.

    The terminal window menu

    Alternatively, if you access the DM24 via an EAM or other Platinum device, you can use the "FORTH Terminal Access" facility, reached via a web-page menu item, located under Tools→Terminal, or use the data-terminal tool from the command line.

    Regardless of the access method, you should see an ok prompt, indicating that the digitiser is ready to receive commands:

    The terminal window, showing the OK prompt
  2. Type

    0 rotation AZIMUTH

    where rotation is the angle of deviation from true North that you measured earlier, as a whole number of tenths of a degree. This is the same angle as that given by the orientation program but with the opposite sign.

    The preceding 0 tells the digitiser to apply the rotation to instrument number 0; the first, or only instrument. To apply rotation to the second instrument of a six-channel DM24, start the command with the figure 1.

    Thus in the example above, you would type 0 906 AZIMUTH to make the digitiser rotate signals by –90.6 degrees.

  3. Reboot the digitiser with the command re-boot.

  4. Collect some more data with the transformation active, and carry out another orientation calculation. The data from the down-hole instrument should now have a maximum coherence with the reference sensor at 0°. Check in particular that the sign of the rotation you have applied is correct.

Affinity Digitisers

The affinity digitiser can apply real-time rotation to its digitised data. It can then produce output streams representing ground motion on true North/South and East/West axes.

To configure this, visit the "Data acquisition" page of the web interface:

Enter the required angle - the sign-reversed output of the orientation calculations - into the "Azimuth rotation" field and then click . (This feature assumes that the first three primary ADC channels are digitising the vertical, North/South and East/West components, respectively, of a triaxial digitiser.)

Once configured, collect some more data with the transformation active, and carry out another orientation calculation. The data from the down-hole instrument should now have a maximum coherence with the reference sensor at 0°. Check in particular that the sign of the rotation you have applied is correct.

Minimus Digitisers

You can configure a Minimus digitiser to apply a real-time rotation to the digitised data. It can then produce output streams representing ground motion on true North/South and East/West axes.

The Minimus' transforms facility is capable of applying arbitrary three-dimensional rotations but, as we only require a single rotation about a vertical axis in this case, we will use the azimuth facility instead.

To configure the rotation:

  1. Visit the unit's web interface and select the Login link from the top menu:

    the 'Login' link in the web interface of a Minimus
  2. Log in as an administrator. If your system has not been changed from the defaults, use username admin and password admin. Otherwise, use the correct username and password for your site.

    Logging in to DIG's web interface

    Click the 'Login' button to proceed.

  3. Once logged in, visit the Setup tab:

    the 'Setup' tab in DIG's web interface, once logged in

    The Applied Rotation section contains fields where you can enter azimuth rotation values for all connected sensors. Populate one (or more) of these with the rotation value(s) computed earlier.

  4. Collect some more data with the rotation active, and carry out another orientation calculation. The data from the down-hole instrument should now have a maximum coherence with the reference sensor at 0°. Check in particular that the sign of the rotation you have applied is correct.