QuakeLogic Blog Archive

Optimizing STA/LTA Settings for Microseismic Earthquake Detection

SeisComP event list showing detected earthquakes with magnitude, location and origin time

In short: For microseismic earthquake detection, start with STA/LTA settings of a 0.5 to 1 second short-term window, a 10 to 30 second long-term window and a trigger threshold between 3:1 and 5:1; a 1 s STA, 20 s LTA and 4:1 ratio is a practical default. Then refine by analyzing historical records, adjusting window lengths, shifting the ratio to trade sensitivity against false alarms, and field testing.

Microseismic earthquake detection plays a crucial role in various fields such as oil and gas exploration, geothermal energy production, and monitoring of induced seismicity in mining. One of the most effective techniques for detecting these small-scale seismic events is the Short-Term Average/Long-Term Average (STA/LTA) method. This blog will guide you through the optimal settings for STA/LTA and how to fine-tune these parameters for reliable microseismic detection.

Understanding STA/LTA

The STA/LTA method works by comparing the short-term average (STA) of the seismic signal to the long-term average (LTA). When the STA significantly exceeds the LTA, it indicates a potential seismic event. Properly configuring the STA/LTA settings is crucial for distinguishing between actual seismic events and background noise.

Recommended STA/LTA Settings

The optimal settings for STA/LTA can vary depending on the specific application and environmental conditions. Here are general guidelines to get you started:

Short-Term Average (STA) Window

The STA window captures the transient signals associated with seismic events. For microseismic detection, a common STA window is between 0.5 to 1 second. This duration is short enough to detect quick, transient signals typical of microseismic activity.

Long-Term Average (LTA) Window

The LTA window establishes the background noise level over a longer period. A typical LTA setting ranges from 10 to 30 seconds. This duration helps to smooth out variability in the background noise and provides a stable reference level.

Threshold Ratio

The threshold ratio is the critical value at which the STA must exceed the LTA to trigger an event detection. A common threshold ratio for microseismic detection is between 3:1 to 5:1. This means the STA must be three to five times greater than the LTA to indicate a potential seismic event.

Example Settings

  • STA Window: 1 second
  • LTA Window: 20 seconds
  • Threshold Ratio: 4:1

These settings are a good starting point, but adjustments may be necessary based on the specific characteristics of the seismic signals and background noise in your area of interest.

Fine-Tuning the STA/LTA Settings

To achieve optimal performance in detecting microseismic events, it is essential to fine-tune the STA/LTA settings. Here’s a step-by-step process:

Analyze Historical Data

Start by reviewing historical seismic data to understand the typical signal-to-noise ratio and the characteristics of microseismic events. This analysis will provide insights into the appropriate initial settings for STA and LTA windows.

Adjust STA/LTA Windows

Based on the initial analysis, adjust the STA and LTA windows. Shorten or lengthen these windows to better capture the speed and duration of seismic signals and background noise levels.

Modify Threshold Ratios

Increase or decrease the threshold ratio to find a balance between detecting genuine seismic events and avoiding false positives caused by noise. A lower ratio may increase sensitivity but could result in more false alarms, while a higher ratio might reduce false positives but miss some events.

Field Testing

Implement the adjusted settings in a real-world scenario and monitor the detection performance. Observe the number of detected events, false positives, and missed events. Continuously refine the settings based on these observations.

Conclusion

Optimizing STA/LTA settings is crucial for reliable microseismic earthquake detection. By carefully configuring the STA and LTA windows and adjusting the threshold ratio, you can enhance the sensitivity and accuracy of seismic event detection. Remember, continuous monitoring and iterative adjustments are key to achieving the best performance.

For those working in fields that require precise seismic monitoring, fine-tuning these parameters will lead to more accurate data collection and better-informed decisions. Whether you are in academia, industry, or government, optimizing your STA/LTA settings is an essential step toward effective seismic monitoring.

About QuakeLogic

QuakeLogic is a leading provider of advanced seismic monitoring solutions, offering a range of products and services designed to enhance the accuracy and efficiency of seismic data acquisition and analysis. Our innovative technologies and expert support help organizations worldwide to better understand and mitigate the impacts of seismic events.

Contact Information:

  • Email: sales@quakelogic.net
  • Phone: +1-916-899-0391
  • WhatsApp: +1-650-353-8627
  • Website: www.quakelogic.net

For more information about our products and services, please visit our website or contact our sales team. We are here to help you with all your seismic monitoring needs.

Thank you for choosing QuakeLogic. We look forward to assisting you with your seismic monitoring projects.

Key Takeaways

  • STA window: 0.5–1 s to capture short transient microseismic arrivals; LTA window: 10–30 s to establish a stable background noise level.
  • Trigger when STA/LTA exceeds 3:1 to 5:1; a lower ratio raises sensitivity and false alarms, a higher ratio misses small events.
  • Example starting configuration: STA 1 s, LTA 20 s, threshold ratio 4:1.
  • Tune in four steps: analyze historical data, adjust STA/LTA windows, modify the threshold ratio, then field-test and iterate.
  • STA/LTA detection is used for oil and gas exploration, geothermal production and induced seismicity monitoring in mining.

Frequently Asked Questions

What STA/LTA settings should I use for microseismic detection?

For microseismic detection, use a short-term average window of 0.5 to 1 second, a long-term average window of 10 to 30 seconds and a trigger threshold between 3:1 and 5:1. A good starting point is STA = 1 s, LTA = 20 s and ratio = 4:1. Adjust from there based on the signal-to-noise ratio and event duration observed in your own historical data.

How does the STA/LTA algorithm detect an earthquake?

The STA/LTA algorithm continuously computes the average absolute amplitude or energy of the seismic signal over a short window and over a longer window. The long window tracks background noise; the short window reacts to a sudden arrival. When the ratio STA/LTA exceeds a set threshold, the detector declares a trigger, marking a potential seismic event. The algorithm is simple, fast and standard in real-time seismic monitoring.

What happens if the STA/LTA threshold ratio is too low or too high?

A threshold ratio that is too low makes the detector highly sensitive, so it triggers on cultural noise, wind or machinery and produces many false positives. A ratio that is too high suppresses false alarms but misses weak microseismic events. The practical approach is to count detected events, false triggers and missed events during field testing and adjust the ratio until the balance suits your monitoring objective.

Why is the LTA window longer than the STA window?

The LTA window must be longer because its job is to characterize slowly varying background noise and provide a stable reference level, which needs 10 to 30 seconds of data. The STA window must be short, around 0.5 to 1 second, so it responds quickly to the brief, transient P-wave arrival of a microseismic event. If the two windows were similar, the ratio would barely change when an event arrived.

How do I fine-tune STA/LTA parameters for my site?

Fine-tune STA/LTA in four steps. First, review historical waveforms to learn the typical signal-to-noise ratio and event duration. Second, shorten or lengthen the STA and LTA windows to match those durations. Third, raise or lower the threshold ratio to balance sensitivity against false triggers. Fourth, deploy the settings, monitor detected, false and missed events in real operation, and iterate until performance is acceptable.

Need help with a project? QuakeLogic designs and supplies seismic monitoring, earthquake early warning, structural health monitoring, infrasound, and shake table systems. Browse QuakeLogic products or contact our engineering team with your site, structure, or test requirements.

Third-Octave Vibration Analysis: VDV, PPV and Dominant Frequency

Third-octave vibration spectrum used to compute VDV, PPV and dominant frequency of ground vibration

In short: Third-octave analysis splits a ground vibration record into 1/3-octave frequency bands and computes three metrics: Vibration Dose Value (VDV, in m/s^1.75, a root-mean-quad of acceleration that is sensitive to peaks), Peak Particle Velocity (PPV, in mm/s) and the dominant site frequency from FFT or PSD. VDV is measured triaxially over 0.5 to 80 Hz per BS 6472-1:2013 and PPV over 4 to 250 Hz per ISO 4866:2010 and BS 7385-2:1993.

Ground vibrations can arise from various sources such as construction activities, heavy machinery, or blasting operations. These vibrations can significantly impact buildings and structures, both in terms of structural integrity and human comfort. At QuakeLogic, we employ advanced techniques and tools like accelerographs and seismographs to measure and record ground vibrations accurately. This blog will delve into the specifics of OCTAVE analysis and vibration data analysis, highlighting the methods and metrics we use to ensure precise vibration assessment and mitigation.

Measuring Ground Vibrations

Ground vibrations on buildings or structures are typically measured outside the structure and at ground level. This approach allows for a comprehensive understanding of how vibrations affect the overall stability and integrity of the structure.

Importance of Frequency Analysis

Frequency analysis of vibrations is crucial for determining the necessary mitigation measures. Third-octave data is particularly useful in specifying the requirements for building foundations, ensuring that structures can withstand the expected vibration levels.

Third-Octave Analysis Metrics

The third-octave analysis involves computing the following key metrics:

  • Vibration Dose Value (VDV)
  • Peak Particle Velocity (PPV)
  • Dominant Site Frequency

Vibration Dose Value (VDV)

Vibration Dose Value (VDV) combines the magnitude of vibration and the duration of exposure. It measures human exposure to vibration within structures, particularly buildings. VDV quantifies vibrations as an exposure dose based on frequency (4.5 Hz to 80 Hz), amplitude, and regularity.

Formula for VDV:

Vibration dose value formula for octave analysis

where:

  • VDV is the vibration dose value in m/s(^{1.75})
  • ( a(t) ) is the acceleration in m/s(^2)
  • ( T ) is the total measurement period in seconds

The VDV formula uses the root-mean-square acceleration raised to the fourth power, known as the root-mean-quad method, making it highly sensitive to peaks in acceleration levels.

Peak Particle Velocity (PPV)

PPV measures the maximum instantaneous velocity of ground particles and is expressed in mm/s. It refers to the internal movement of molecular particles within the ground rather than surface displacement. Typical environmental ground vibrations range from 1 Hz to 200 Hz.

Sources of ground vibrations at construction sites include pile driving, dynamic compaction, blasting, and heavy equipment operation. These vibrations can range from disturbing residents to causing visible structural damage.

Dominant Frequency

Fast Fourier Transform (FFT) and Power Spectral Density (PSD) are powerful tools used to convert vibration signals from the time domain to the frequency domain. FFTs are suitable for analyzing vibrations with dominant frequency components, while PSDs are ideal for characterizing random vibration signals like ambient measurements.

Human Comfort and Noise Determination

Human comfort is significantly impacted by vibrations and noise. Vibrations can cause discomfort, annoyance, and even health issues for occupants of buildings. By conducting detailed OCTAVE analysis, we can determine the levels of vibration that affect human comfort and establish mitigation measures to reduce these impacts.

Noise and Human Comfort:

  • VDV is used to quantify the vibration exposure that affects human comfort, particularly in residential and office buildings.
  • PPV and dominant frequency analyses help in understanding the specific sources of vibration and noise, allowing for targeted mitigation strategies.
  • Standards like ISO 2631-2:2003 “Mechanical vibration and shock – Evaluation of human exposure to whole-body vibration” provide guidelines for acceptable vibration levels to ensure human comfort.

Measurement Methodology

Standards and Reference Guides

  • British Standard BS 6472-1:2013: Guide to evaluation of human exposure to vibration in buildings – Vibration sources other than blasting.
  • British Standard BS 7385-2:1993: Evaluation and measurement for vibration in buildings – Guide to damage levels from ground-borne vibration.
  • International Standard ISO 4866:2010: Mechanical vibration and shock measurement – Vibration of fixed structures – Guidelines for the measurement of vibrations and evaluation of their effects on structures.
  • International Standard ISO 2631-2:2003: Mechanical vibration and shock – Evaluation of human exposure to whole-body vibration – Part 2: Vibration in buildings (1 Hz to 80 Hz).
  • International Standard ISO 5348:1998: Mechanical vibration and shock – Mechanical mounting of accelerometers.
  • DIN 4150-3:1999: Structural vibration – Part 3: Effects of vibration on structures.

Instrumentation

  • VDV Measurements: Made triaxially (X, Y, and Z axes) across the frequency range of 0.5 Hz to 80 Hz in 1/3-octave bands, following BS 6472-1:2013 guidelines.
  • PPV Measurements: Made triaxially across the frequency range of 4 Hz to 250 Hz in 1/3-octave bands, following ISO 4866:2010 and BS 7385-2:1993 guidelines.

Conclusion

Ground vibrations can pose serious risks to structures and human comfort. Through detailed OCTAVE analysis and vibration data analysis, QuakeLogic provides accurate and comprehensive solutions to monitor and mitigate these impacts effectively.

About QuakeLogic

QuakeLogic is a leading provider of advanced seismic and vibration monitoring solutions, offering a range of products and services designed to enhance the accuracy and efficiency of seismic data acquisition and analysis. Our innovative technologies and expert support help organizations worldwide to better understand and mitigate the impacts of seismic events.

Contact Information:

  • Email: sales@quakelogic.net
  • Phone: +1-916-899-0391
  • WhatsApp: +1-650-353-8627
  • Website: www.quakelogic.net

For more information about our products and services, please visit our website or contact our sales team. We are here to help you with all your seismic monitoring needs.

Thank you for choosing QuakeLogic. We look forward to assisting you with your seismic monitoring projects.

Key Takeaways

  • VDV combines acceleration magnitude and exposure duration using a fourth-power (root-mean-quad) integral, in m/s^1.75.
  • PPV is the maximum instantaneous ground particle velocity in mm/s; environmental ground vibration spans about 1 to 200 Hz.
  • VDV is measured triaxially in 1/3-octave bands from 0.5 to 80 Hz following BS 6472-1:2013.
  • PPV is measured triaxially in 1/3-octave bands from 4 to 250 Hz following ISO 4866:2010 and BS 7385-2:1993.
  • FFT suits signals with dominant frequency components; PSD characterizes random or ambient vibration.

Frequently Asked Questions

What is third-octave analysis in vibration measurement?

Third-octave analysis divides a vibration signal into frequency bands each one-third of an octave wide and reports the level in every band. For ground vibration affecting buildings, QuakeLogic computes Vibration Dose Value, Peak Particle Velocity and dominant site frequency from these bands. Third-octave data is particularly useful for specifying building foundation requirements because it shows which frequencies carry the energy.

What is the difference between VDV and PPV?

Vibration Dose Value (VDV) quantifies human exposure: it integrates acceleration to the fourth power over the measurement period, giving a value in m/s^1.75 that is highly sensitive to peaks and is used for comfort in residential and office buildings over 4.5 to 80 Hz. Peak Particle Velocity (PPV) is the maximum instantaneous ground particle velocity in mm/s and is the metric used to judge potential structural damage from construction or blasting vibration.

Which standards apply to ground vibration measurement in buildings?

The main references are BS 6472-1:2013 for human exposure to building vibration from sources other than blasting, BS 7385-2:1993 and ISO 4866:2010 for damage levels and measurement of vibration in fixed structures, ISO 2631-2:2003 for whole-body vibration in buildings between 1 and 80 Hz, ISO 5348:1998 for accelerometer mounting, and DIN 4150-3:1999 for the effects of vibration on structures.

Over what frequency range are VDV and PPV measured?

VDV measurements are made triaxially (X, Y and Z axes) in 1/3-octave bands from 0.5 Hz to 80 Hz following BS 6472-1:2013. PPV measurements are made triaxially in 1/3-octave bands from 4 Hz to 250 Hz following ISO 4866:2010 and BS 7385-2:1993. Both are recorded outside the structure at ground level using accelerographs or seismographs.

When should I use an FFT versus a PSD for vibration data?

Use a Fast Fourier Transform when the vibration has clear dominant frequency components, such as machinery or pile-driving tones, because the spectrum shows the peaks directly. Use a Power Spectral Density when the signal is random or broadband, such as ambient ground vibration, since PSD normalizes power per unit frequency and gives a statistically stable description of random motion.

Need help with a project? QuakeLogic designs and supplies seismic monitoring, earthquake early warning, structural health monitoring, infrasound, and shake table systems. Browse QuakeLogic products or contact our engineering team with your site, structure, or test requirements.

How to Download Waveform Data from SeisComP

SeisComP event list showing detected earthquakes with magnitude, location and origin time

In short: You can download waveform data from SeisComP in three ways: run the bundled fdsnws_fetch command with network, station, location, channel and ISO start/end times to save a MiniSEED file; call client.get_waveforms() from ObsPy's FDSN client pointed at your server's /fdsnws endpoint; or select an event in the scolv or scmv GUI and use its waveform download option to export the traces.

Downloading waveform data from SeisComP can be essential for various seismological research and analysis tasks. Whether you’re a seasoned seismologist or a beginner, understanding the different methods available for extracting this data can enhance your workflow. In this blog post, we will walk you through three main methods to download waveform data from SeisComP: using the fdsnws_fetch command, employing a Python script with ObsPy, and leveraging the SeisComP graphical user interface (GUI).

Method 1: Using fdsnws_fetch Command

SeisComP provides a handy command-line utility called fdsnws_fetch for downloading waveform data from FDSN web services. Here’s how you can use it:

  1. Install fdsnws_fetch:
    Make sure fdsnws_fetch is installed on your system. It usually comes bundled with SeisComP, but if not, you can install it separately.
  2. Execute fdsnws_fetch:
    Open your terminal and run the following command:
   fdsnws_fetch -H <your_seiscomp_server> -N <network> -S <station> -L <location> -C <channel> -s <start_time> -e <end_time> -o <output_file>

Replace the placeholders with appropriate values:

  • <your_seiscomp_server>: Your SeisComP server address.
  • <network>: Network code (e.g., IU).
  • <station>: Station code (e.g., ANMO).
  • <location>: Location code (e.g., 00).
  • <channel>: Channel code (e.g., BHZ).
  • <start_time>: Start time in ISO format (e.g., 2024-07-01T00:00:00).
  • <end_time>: End time in ISO format (e.g., 2024-07-01T01:00:00).
  • <output_file>: Path to save the output file (e.g., waveform.mseed).

Method 2: Using Python Script with ObsPy

ObsPy is a powerful Python library for seismology that can be used to download waveform data from SeisComP. Follow these steps:

  1. Install ObsPy:
    First, install ObsPy by running:
   pip install obspy
  1. Download Waveform Data:
    Use the following Python script to download the data:
   from obspy.clients.fdsn import Client
   from obspy import UTCDateTime

   client = Client("http://<your_seiscomp_server>/fdsnws")

   network = "<network>"
   station = "<station>"
   location = "<location>"
   channel = "<channel>"
   start_time = UTCDateTime("2024-07-01T00:00:00")
   end_time = UTCDateTime("2024-07-01T01:00:00")

   st = client.get_waveforms(network, station, location, channel, start_time, end_time)
   st.write("waveform.mseed", format="MSEED")

Be sure to replace the placeholders with your specific values.

Method 3: Using SeisComP GUI

If you prefer a graphical approach, SeisComP provides a user-friendly GUI for data downloading:

  1. Open SeisComP GUI:
    Launch the SeisComP Monitoring Viewer (scmv) or SeisComP Online Visualization (scolv).
  2. Select the Event:
    Navigate to the event of interest in the GUI.
  3. Download Data:
    Use the waveform download option to specify the time range, network, station, and channels, then proceed to download the data.

Additional Tips

  • Check Access Permissions: Ensure you have the necessary permissions to access and download data from the SeisComP server.
  • Consult Documentation: Refer to the SeisComP documentation for detailed instructions and options specific to your setup.
  • Use Filters: Apply filters to the data as needed to focus on specific signals or reduce noise.

By following these methods, you can efficiently download and manage waveform data from SeisComP, aiding in your seismological research and analysis. If you have any questions or need further assistance, feel free to reach out!

About QuakeLogic

QuakeLogic is a leading provider of advanced seismic monitoring solutions, offering a range of products and services designed to enhance the accuracy and efficiency of seismic data acquisition and analysis. Our innovative technologies and expert support help organizations worldwide to better understand and mitigate the impacts of seismic events.

Contact Information:

  • Email: sales@quakelogic.net
  • Phone: +1-916-899-0391
  • WhatsApp: +1-650-353-8627
  • Website: www.quakelogic.net

For more information about our products and services, please visit our website or contact our sales team. We are here to help you with all your seismic monitoring needs.

Thank you for choosing QuakeLogic. We look forward to assisting you with your seismic monitoring projects.

Key Takeaways

  • fdsnws_fetch -H <server> -N <net> -S <sta> -L <loc> -C <chan> -s <start> -e <end> -o waveform.mseed pulls data from the command line.
  • In ObsPy, Client("http://<server>/fdsnws").get_waveforms(…) returns a Stream you save with st.write("waveform.mseed", format="MSEED").
  • The scolv and scmv GUIs let you pick an event and download its waveforms by time range, network, station and channel.
  • Times use ISO format such as 2024-07-01T00:00:00; codes follow FDSN conventions (IU, ANMO, 00, BHZ).
  • Confirm server access permissions first and apply filters afterwards to isolate signals of interest.

Frequently Asked Questions

What is fdsnws_fetch in SeisComP?

fdsnws_fetch is a command-line utility bundled with SeisComP for retrieving waveform data from FDSN web services. You pass the server host with -H, network, station, location and channel codes with -N, -S, -L and -C, ISO-formatted start and end times with -s and -e, and an output path with -o, for example waveform.mseed. It is the quickest option for scripted or one-off downloads.

How do I download SeisComP waveforms with ObsPy?

Install ObsPy with pip install obspy, then create an FDSN client with Client("http://<your_seiscomp_server>/fdsnws"). Define network, station, location and channel codes and UTCDateTime start and end values, call client.get_waveforms() with those arguments to receive a Stream object, and save it with st.write("waveform.mseed", format="MSEED"). The same script can loop over many stations or time windows.

Can I download waveform data from the SeisComP GUI without scripting?

Yes. Launch the SeisComP Monitoring Viewer (scmv) or the Origin Locator View (scolv), navigate to the event of interest, and use the waveform download option to specify the time range, network, station and channels before saving. This is convenient for analysts who want to inspect an event visually and export its traces without writing commands or Python.

What format does SeisComP waveform data download in?

Both fdsnws_fetch and the ObsPy example in the guide write MiniSEED files (.mseed), the standard compact binary format for continuous seismic waveforms. MiniSEED stores the time series and station and channel identifiers but not instrument response, so if you need ground motion in physical units you should also fetch station metadata (StationXML) from the same FDSN service.

Why does my SeisComP waveform download fail?

The most common causes are missing access permissions on the SeisComP server, an incorrect FDSN endpoint URL, or wrong network, station, location or channel codes. Verify you are authorized to access the data, check the codes against the station inventory, confirm the time window actually contains recorded data, and consult the SeisComP documentation for options specific to your installation.

Need help with a project? QuakeLogic designs and supplies seismic monitoring, earthquake early warning, structural health monitoring, infrasound, and shake table systems. Browse QuakeLogic products or contact our engineering team with your site, structure, or test requirements.