Calculation of slip behaviour and transient permeability evolution of a fracture

General description:

This is the Python code for the calculation of slip behavior and permeability evolution of a natural rough fracture in the manuscript entitled "Frictional Stability and Permeability Enhancement of a Natural Fracture in Granite under Monotonic vs. Cyclic Fluid Injections". The manuscript has been submitted for consideration ofr publication in Journal of Geophysical Research: Solid Earth.

Required packages:

• python==3.12.2
• numpy==1.26.3
• pandas==2.1.4
• scikit-opt==0.6.6
• matplotlib==3.8.0

Usage:

• The file named "example.ipynb" is a jupyter-notebook file, which includes the source code of the inversion of consitutive parameters of both the rate-and-state friction law and the displacement- and velocity-dependent aperture model, as well as the source code of the calculation of slip behavior and permeability evolution based on these two models. The excel file named "example_monotonic_pressurization.xlsx" is the example data needed by the source code.
• The source code depends on a C++ package on Github. To use it, the package need to be appropriately placed and compiled in accordance with the following steps:
• 1. In the folder where the source code file (example.ipynb) is located, run the following command to download the package: git clone https://github.com/Ranger-boop/velodep.git.
  2. Consecutively run the following commands to compile the package based on cmake:
  3. • On windows platform using VisutalStudio compiler:
     • 1. cd velodep
       2. mkdir build
       3. cd buid
       4. cmake -A x64 ..
       5. cmake --build . --config Release
     • On Linux platform using gcc compiler:
     • 1. cd velodep
       2. mkdir build
       3. cd buid
       4. cmake -A x64 -DCMAKE_BUILD_TYPE=Release ..
       5. cmake --build .
• After the package has been compiled, open the "example.ipynb" file and run the code in the file.