Saved in:
| Main Authors: | , , , , , , , , |
|---|---|
| Format: | Recurso digital |
| Language: | |
| Published: |
Zenodo
2026
|
| Online Access: | https://doi.org/10.5281/zenodo.19590352 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
Table of Contents:
- <h2>Physics-Based Hindcasting of an Induced Earthquake Sequence</h2> <p>This repository contains the supplementary files, input data, and post-processing scripts for the AGU <em>Journal of Geophysical Research: Solid Earth</em> submission:</p> <p><strong>“Physics-Based Hindcasting of an Induced Earthquake Sequence”</strong><br>by V. van der Heiden (vanderheidenv.uni@gmail.com), L. Buijze, T. Ulrich, J.-D. van Wees, M. van Isselt, L. Matenco, J. Spetzler, A.-A. Gabriel, and Y. van Dinther.</p> <p>The repository provides MATLAB-based workflows to reproduce the figures and results presented in the manuscript, focusing on physics-based modeling of induced seismicity using a 2.5D on-fault stress framework.</p> <h2>[Code availability]</h2> <ul> <li>The open-source software <strong>PANTHER</strong> is available at:<br><a class="decorated-link" href="https://github.com/TNO/PANTHER" rel="noopener">https://github.com/TNO/PANTHER</a></li> <li>The geological parameterization framework <strong>PhysMmax</strong> is available at:<br><a class="decorated-link" href="https://github.com/loes-buijze/PhysMmax" rel="noopener">https://github.com/loes-buijze/PhysMmax</a></li> <li>This repository contains scripts and data required to reproduce the published results:<br><a class="decorated-link" href="https://github.com/VvanderHeiden/InducedEarthquakeSequence" rel="noopener">https://github.com/VvanderHeiden/InducedEarthquakeSequence</a></li> </ul> <p>Specific commits used in the study should be archived via Zenodo and cited accordingly.</p> <h2>[Repository structure]</h2> <ul> <li><strong>run_fault_loading_reference_model.m</strong><br>Main script to reproduce the reference model simulations and figures.</li> <li><strong>/input/</strong><br>Contains input data for the Wirdum fault case: <ul> <li><code>wirdum_geology.mat</code>: Geological model parameters</li> <li><code>200pillars_constoffset_*.mat</code>: Fault discretization and pillar geometry</li> </ul> </li> <li><strong>/functions/</strong><br>MATLAB functions for model setup, stress computation, and visualization: <ul> <li>Stress computation and processing:<br><code>write_stresses_to_file.m</code>, <code>calculate_nucleation_depths.m</code>, <code>calculate_depth_average.m</code></li> <li>Geometry and data handling:<br><code>rotate_pillars.m</code>, <code>resample_data.m</code>, <code>extract_regularize_filter_data.m</code></li> <li>Model adjustments:<br><code>apply_smoothing_within_offset_regions.m</code>, <code>apply_scaling_within_reservoir.m</code>, <code>update_modelname.m</code></li> <li>Visualization:<br><code>plot_results_2D_continuous.m</code></li> </ul> </li> <li><strong>LICENSE</strong><br>Licensing information for reuse.</li> </ul> <h2>[Model description]</h2> <p>This repository implements a <strong>2.5D on-fault stress calculation</strong> for the Wirdum fault case using a pillar-based representation. The model applies a plane-strain approximation per pillar to simulate stress evolution and earthquake nucleation.</p> <p>The workflow couples:</p> <ul> <li>Fault mechanics (<strong>PANTHER</strong>)</li> <li>Geological parameterization (<strong>PhysMmax</strong>)</li> <li>Stress evolution and nucleation analysis (Day3D criterion)</li> </ul> <h2>[Workflow overview]</h2> <ol> <li>Load fault geometry and geological parameters</li> <li>Construct a pillar-based fault representation</li> <li>Define initial stress state and hydro-mechanical parameters</li> <li>Assign depth-dependent and depth-independent properties</li> <li>Apply smoothing and scaling (reservoir vs. full fault)</li> <li>Run quasi-static stress evolution simulations</li> <li>Evaluate nucleation and reactivation conditions</li> <li>Post-process and visualize stress fields and seismic response</li> </ol> <h2>[Outputs]</h2> <ul> <li>On-fault stress distributions per pillar</li> <li>Temporal evolution of stress and slip</li> <li>Nucleation depth and timing estimates</li> <li>Filtered reservoir-scale stress changes</li> <li>Publication-quality figures reproducing the manuscript results</li> </ul> <h2>[Setup instructions]</h2> <ol> <li> <p>Clone this repository:</p> <div class="relative w-full mt-4 mb-1"> <div class=""> <div class="relative"> <div class="h-full min-h-0 min-w-0"> <div class="h-full min-h-0 min-w-0"> <div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"> <div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"> <div class="pointer-events-none absolute end-1.5 top-1 z-2 md:end-2 md:top-1"> </div> <div class="relative"> <div class="pe-11 pt-3"> <div class="relative z-0 flex max-w-full"> <div class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy"> <div class="cm-scroller"> <div class="cm-content q9tKkq_readonly">git clone https://github.com/VvanderHeiden/InducedEarthquakeSequence.git</div> </div> </div> </div> </div> </div> </div> </div> </div> </div> <div class=""> <div class=""> </div> </div> </div> </div> </div> </li> <li> <p>Clone required dependencies:</p> <div class="relative w-full mt-4 mb-1"> <div class=""> <div class="relative"> <div class="h-full min-h-0 min-w-0"> <div class="h-full min-h-0 min-w-0"> <div class="border border-token-border-light border-radius-3xl corner-superellipse/1.1 rounded-3xl"> <div class="h-full w-full border-radius-3xl bg-token-bg-elevated-secondary corner-superellipse/1.1 overflow-clip rounded-3xl lxnfua_clipPathFallback"> <div class="pointer-events-none absolute end-1.5 top-1 z-2 md:end-2 md:top-1"> </div> <div class="relative"> <div class="pe-11 pt-3"> <div class="relative z-0 flex max-w-full"> <div class="q9tKkq_viewer cm-editor z-10 light:cm-light dark:cm-light flex h-full w-full flex-col items-stretch ͼk ͼy"> <div class="cm-scroller"> <div class="cm-content q9tKkq_readonly">git clone https://github.com/TNO/PANTHER.git<br>git clone https://github.com/loes-buijze/PhysMmax.git</div> </div> </div> </div> </div> </div> </div> </div> </div> </div> <div class=""> <div class=""> </div> </div> </div> </div> </div> </li> <li>Open MATLAB (version > 2021 recommended)</li> <li>Open the <code>PhysMmax.prj</code> project file</li> <li>Add PANTHER as a referenced project: <ul> <li>Go to <em>Project Tab → References → Add Reference</em></li> <li>Select the <code>Panther.prj</code> file</li> </ul> </li> </ol> <p>After setup, PANTHER functionality is available within the PhysMmax environment.</p> <h2>[Notes]</h2> <ul> <li>The model assumes <strong>plane-strain conditions per pillar</strong> (2.5D approximation)</li> <li>Depth-dependent parameters are vertically interpolated</li> <li>Diffusion and aseismic slip processes can be toggled in the simulation settings</li> </ul>