Difference between revisions of "Geopsy: H/V and Spectrum Toolboxes: Processing Tab"

From GeopsyWiki
Jump to navigation Jump to search
Line 36: Line 36:
  
 
=== Taper ===
 
=== Taper ===
Use or not of a [[Geopsy:_Taper|cosine taper]] on each of the selected time windows.
+
[[Image:HV_SPEC_toolb_process_param_taper.png|thumb|right|300px|Taper box and spin box]]Check the box to apply a [[Geopsy:_Taper|cosine taper]] on both sides of the selected time windows.
  
A taper is used to minimize the border effects due to the extraction of a window so as to improve its properties in the frequency domain. The use of a [[Geopsy:_Taper|taper]] is strongly recommended so as to avoid the creation of weird frequencies.
+
A taper is used to minimize the border effects due to the extraction of a window so as to improve its properties in the frequency domain. The use of a [[Geopsy:_Taper|taper]] is strongly recommended so as to avoid the creation of spurious frequencies.
  
 
The taper width is adjusted in the spin box.
 
The taper width is adjusted in the spin box.
  
 
=== Filter ===
 
=== Filter ===
Fixing the use or not of a [[Geopsy:_Filter|filter]] on the whole signal.
+
[[Image:HV_SPEC_toolb_process_param_taper.png|thumb|right|300px|Taper box and spin box]]When the box is checked a [[Geopsy:_Filter|filter]] is applied to each component of the entire record before selecting the stable time windows.
  
This option allows to filter the signal before the extraction of the time windows.
+
Only a "High pass" filter is allowed.
  
To the contrary of the [[Geopsy: Waveform Menu|Waveform menu filter]] there is only a "High pass" filter allowed, rejecting frequencies below a given value
+
Sometimes, the taper size used in the pre-processing has a strong influence at low frequency. If the signal contains a strong and very low frequency component, then cutting into short time windows may strongly distort the observed spectra and H/V. Even when using the minimum window length criteria according to the SESAME guidelines <ref>Guidelines for the implementation of the H/V spectral ratio technique on ambient vibrations measurements, processing and interpretation, [http://sesame-fp5.obs.ujf-grenoble.fr], 62 pages, April 2005</ref> (i.e. 10 times the corresponding period for the minimum frequency of interest) may be not enough to ensure a reliable H/V estimation.
  
Sometimes, the taper size used in the pre-processing has a strong influence at low frequency. If the signal contains a strong and very low frequency component, then cutting into short time windows may strongly distort the observed spectra and H/V. Even when using the minimum window length criteria according SESAME guidelines <ref>Guidelines for the implementation of the H/V spectral ratio technique on ambient vibrations measurements, processing and interpretation, [http://SESAME-FP5.obs.ujf-grenoble.fr], 62 pages, April 2005</ref> (i.e. 10 times the corresponding period for the minimum frequency of interest) may be not enough to ensure a reliable H/V estimation.
+
In order to avoid such “taper side effects”, signals can be high-pass filtered before computing the H/V and spectra curves.
 
 
In order to avoid such “taper side effects”, signals can be high-pass filtered before computing the H/V and spectra curves (Geopsy option).
 
 
The higher value for the frequency of the high-pass filter is the minimum ‘reliable’ frequency, i.e. <math>\frac{10}{w_l}</math>, [[Effect of very low frequency on H/V]].
 
The higher value for the frequency of the high-pass filter is the minimum ‘reliable’ frequency, i.e. <math>\frac{10}{w_l}</math>, [[Effect of very low frequency on H/V]].
  
  
This option should be used with caution. As the signal is high-pass filtered, the results below the chosen frequency could be somewhat flaky.
+
This option should be used with caution. As the signal is high-pass filtered, the results below the chosen frequency could be somewhat flaky.
  
 
== Horizontal components section (only available for [[H/V_spectral_ratio|HV toolbox]]) ==
 
== Horizontal components section (only available for [[H/V_spectral_ratio|HV toolbox]]) ==

Revision as of 10:44, 10 January 2013

Use

This toolbox tab contains the smoothing and taper parameters for H/V or Spectrum curves and the way the horizontal components are processed to compute the H/V spectral ratio, described in two sections:

  • the Parameters section (For H/V and Spectrum tools)
  • the Horizontal components section (For H/V tool only)


Processing tab for smoothing parameters and horizontal components processing in H/V computing

Parameters section

This section is divided in three parts:


Smoothing

Smoothing parameters

Type of smoothing applied to the windows. During the computing of H/V or Spectrum, the Fourier spectra can be smoothed (it is strongly recommended) in the goal to clarify the global aspect of the curves.

For the Spectrum computation, the Fourier spectrum from each selected time window is smoothed and then the averaged curve is computed.

For the H/V computation, the horizontal Fourier spectra (NS and EW) are first combined and then the smoothing is applied on the merged horizontal Fourier spectrum and on the vertical Fourier spectrum.


Four smoothing methods are evailable:

Smoothing option drop box
  1. Konno and Ohmachi[1] smoothing. This smoothing use a constant bandwidth in a logarithmic scale and is strongly recommended because this smoothing function preserves the different number of points at low and high frequency. This smoothing is controlled by a smoothing constant varying in-between 0 and 100. A constant of 0 gives a very strong smoothing, when a constant of 100 a very soft smoothing, more details;
  2. Constant smoothing. This smoothing function has a triangular shape centered on the current frequency and its width is equal to "Band width" given by the user. This band width varies from 0 (soft smoothing) to 100 (strong smoothing);
  3. Proportional smoothing. This smoothing function has a triangular shape and its width depends upon the current frequency. The half width is defined by percentage*frequency. The value of "percentage" cannot be greater or equal to 100% (strong smoothing).
  4. No smoothing. This can be sometime useful to have a rough idea about Fourier spectra;

The Smoothing constant spin box is used to give the smoothing constant value.

Taper

Taper box and spin box

Check the box to apply a cosine taper on both sides of the selected time windows.

A taper is used to minimize the border effects due to the extraction of a window so as to improve its properties in the frequency domain. The use of a taper is strongly recommended so as to avoid the creation of spurious frequencies.

The taper width is adjusted in the spin box.

Filter

Taper box and spin box

When the box is checked a filter is applied to each component of the entire record before selecting the stable time windows.

Only a "High pass" filter is allowed.

Sometimes, the taper size used in the pre-processing has a strong influence at low frequency. If the signal contains a strong and very low frequency component, then cutting into short time windows may strongly distort the observed spectra and H/V. Even when using the minimum window length criteria according to the SESAME guidelines [2] (i.e. 10 times the corresponding period for the minimum frequency of interest) may be not enough to ensure a reliable H/V estimation.

In order to avoid such “taper side effects”, signals can be high-pass filtered before computing the H/V and spectra curves. The higher value for the frequency of the high-pass filter is the minimum ‘reliable’ frequency, i.e. , Effect of very low frequency on H/V.


This option should be used with caution. As the signal is high-pass filtered, the results below the chosen frequency could be somewhat flaky.

Horizontal components section (only available for HV toolbox)

This section is devoted to indicate how the H/V has to be computed. There are 3 different methods that can be applied.


Horizontal components section indicating the way H/V is processed (available for H/V computing) only

Squared average

The H/V is computed as followed:

  • The Fourier amplitude spectrum is computed for each window of each component (N, E and Z);
  • The horizontal average spectrum is computed by a combination of horizontal spectra for each window, following the formula
  • The horizontal (averaged) and vertical spectra of each window are smoothed;
  • The horizontal to vertical ratio (H/V) is computed for each window;
  • The H/V is computed by averaging all H/V coming from individual windows.

Total horizontal energy

The H/V is computed as followed:

  • The Fourier amplitude spectrum is computed for each window of each component (N, E and Z)
  • The horizontal average spectrum is computed by a combination of horizontal spectra for each window, following the formula
  • The horizontal (averaged) and vertical spectra of each window are smoothed
  • The horizontal to vertical ratio (H/V) is computed for each window
  • The H/V is computed by averaging all H/V coming from individual windows.

Directional energy

Here, the H/V is calculated along a given direction, 0° is assumed to be the North direction and 90° the East direction. From the North-South and the East-West signals, a new (and single) signal is computed geometrically. This new signal is assumed to be the horizontal signal. The direction in which the H/V will be calculated is given in a spin box allowing changes in the direction between 0° and 180°. If the user wants to compute in a direction between 180° and 360°, due to symetry the user should take the symetrical direction: wanted direction-180°. If the user wants a 220° direction, it should be entered 40° in the spin box.

  • The Fourier amplitude spectrum is computed for each window of each component (recomputed horizontal and vertical)
  • The horizontal (recomputed) and vertical spectra of each window are smoothed
  • The horizontal to vertical ratio (H/V) is computed for each window
  • The H/V is computed by averaging all H/V coming from individual windows.

Load parameters and Start section

The section at the bottom of the Time tab comprises two buttons.

Load parameters and Start section


  • The Load parameters button is used to load parameters from previous H/V or spectrum processing stored in a name.log file (example).
  • Press the Start button to start H/V or Spectrum processing.

If no window selection has been performed, a pup-up window appears.

Pop-up window indicating that no window selection has been performed


Simply click on the Yes button and processing will follow its way with the current Toolbox parameters.

To perform user's windowing, click on the No button.


References

  1. Konno K. and T. Ohmachi, 1998. Ground motion characteristics estimated from spectral ratio between horizontal and vertical components of microtremors. Bull. seism. Soc. Am., 88-1, 228-241.
  2. Guidelines for the implementation of the H/V spectral ratio technique on ambient vibrations measurements, processing and interpretation, [1], 62 pages, April 2005