geoskill-earthquake-liquefaction-risk

Youd简化法结合地质地下水位与PGA的液化指数评估

地震液化风险评估 | Earthquake Liquefaction Risk

Assesses sandy-soil liquefaction using the Youd simplified procedure: cyclic stress ratio CSR = 0.65·(PGA/g)·rd(z)·gw (rd = stress reduction factor, gw = shallow-groundwater amplification), cyclic resistance ratio CRR increases with SPT blow count N (denser soils are more liquefaction-resistant), factor of safety FS = CRR/CSR (FS < 1 indicates liquefaction), and the liquefaction index LPI (≥0) is obtained by integrating with Iwasaki depth weighting. Higher PGA → lower FS → higher LPI (positive correlation).

Dependencies / 依赖

bash

pip install numpy rasterio scipy

Usage / 使用方法

Basic usage (synthetic data, offline)

bash

python geoskill-earthquake-liquefaction-risk.py --bbox 116.0 39.0 117.0 40.0 --synthetic --output-dir ./out

More examples

bash

python geoskill-earthquake-liquefaction-risk.py --bbox 116 39 117 40 --synthetic --output-dir ./out
python geoskill-earthquake-liquefaction-risk.py --input site.tif --depth 5 --output-dir ./out
python geoskill-earthquake-liquefaction-risk.py --bbox 116 39 117 40 --depth 3 --fines 20 --synthetic --output-dir ./out
python geoskill-earthquake-liquefaction-risk.py --bbox 117 39 118 40 --synthetic --quiet --output-dir ./out

Output / 输出

File	Format	Description
liquefaction_index.tif	GeoTIFF	Liquefaction index LPI (≥0)
factor_of_safety.tif	GeoTIFF	Factor of safety FS (<1 liquefaction)
liquefaction_params.json	JSON	Assessment depth / fines content / mean groundwater depth

Each run also produces output-manifest.json (run manifest, including input/output/QA summary).

Data Source / 数据源 / Source

Real mode reads a multi-band GeoTIFF (band1=PGA (g), band2=SPT N value, band3=groundwater table depth in m); synthetic mode generates a site offline.

Privacy / 隐私声明 / Privacy

• Runs offline by default; --synthetic mode requires no network at all.
• All processing is done locally; no user data is uploaded.

License / License

MIT

---

---

name: geoskill-earthquake-liquefaction-risk description: 'Youd简化法结合地质地下水位与PGA的液化指数评估'

地震液化风险评估 | Earthquake Liquefaction Risk

用 Youd 简化法评估砂土液化：循环应力比 CSR = 0.65·(PGA/g)·rd(z)·gw（rd 应力折减、gw 浅地下水放大），循环阻抗比 CRR 随 SPT 击数 N 增大（越密实越抗液化），安全系数 FS=CRR/CSR（<1 判定液化），再按 Iwasaki 深度权重积分得液化指数 LPI（≥0）。PGA 越高 → FS 越低 → LPI 越高（正相关）。

依赖

bash

pip install numpy rasterio scipy

使用方法

基本用法（合成数据，离线）

bash

python geoskill-earthquake-liquefaction-risk.py --bbox 116.0 39.0 117.0 40.0 --synthetic --output-dir ./out

更多示例

bash

python geoskill-earthquake-liquefaction-risk.py --bbox 116 39 117 40 --synthetic --output-dir ./out
python geoskill-earthquake-liquefaction-risk.py --input site.tif --depth 5 --output-dir ./out
python geoskill-earthquake-liquefaction-risk.py --bbox 116 39 117 40 --depth 3 --fines 20 --synthetic --output-dir ./out
python geoskill-earthquake-liquefaction-risk.py --bbox 117 39 118 40 --synthetic --quiet --output-dir ./out

输出

文件	格式	说明
liquefaction_index.tif	GeoTIFF	液化指数 LPI（≥0）
factor_of_safety.tif	GeoTIFF	安全系数 FS（<1 液化）
liquefaction_params.json	JSON	评估深度/细粒含量/平均地下水位

每次运行还会产出 output-manifest.json（运行清单，含输入/产物/QA 摘要）。

数据源 / Source

真实模式读取多波段 GeoTIFF（band1=PGA(g)、band2=SPT N值、band3=地下水位埋深m）；合成模式离线生成场地。

隐私声明 / Privacy

• 默认离线运行，--synthetic 模式完全无网络。
• 所有处理在本地完成，不上传任何用户数据。

License

MIT