Geoskill: Tailings Dam Risk

Prerequisites / 先准备 X 文件

⚠️ 必读 — 本 skill 不属于即用型，需要先准备特定文件才能跑。

本 skill 需要 设施清单 JSON（geometry + 属性）。--bbox 模式下 DEM 和 Sentinel-2 时序自动从 MPC 下载。

👉 完整教程见仓库根目录 PREREQUISITES.md 1.6 节。

先准备 X 文件：先准备一份 facilities.json（可手写小样本），其他自动下载。

快速试跑命令：

python tailings_dam_risk.py --bbox 116.30,39.85,116.45,39.95 --date-range 2023-06-01,2023-09-30 --facilities facilities.json --output-dir ./tdr

Tailings Dam Risk

Screens tailings dam facilities using remote sensing and DEM analysis. Produces risk rankings and patrol priorities based on multi-temporal change detection, catchment analysis, and downstream exposure assessment.

Trigger

Use when the user wants to:

• Screen tailings dam facilities for safety risks
• Detect water surface changes in tailings reservoirs
• Analyze upstream catchment contributions
• Estimate downstream impact zones from potential dam failure
• Identify exposed objects (buildings, roads, population) downstream
• Generate risk reports for patrol prioritization
• Assess deformation patterns from InSAR data

CLI Usage

# Basic screening with synthetic data (demo)
python scripts/tailings_dam_risk.py --output-dir ./tdr-output

# With custom facilities file
python scripts/tailings_dam_risk.py \
  --facilities facilities.json \
  --output-dir ./tdr-output

# With DEM and water masks
python scripts/tailings_dam_risk.py \
  --dem-file dem.tif \
  --water-masks water_2019.tif water_2021.tif water_2023.tif \
  --output-dir ./tdr-output

# Full analysis with all options
python scripts/tailings_dam_risk.py \
  --facilities facilities.json \
  --dem-file dem.tif \
  --water-masks water_2019.tif water_2021.tif \
  --deformation-data insar_deformation.tif \
  --exposure-objects downstream_objects.json \
  --rainfall-scenario 500yr \
  --runout-method simplified \
  --deformation-threshold 15.0 \
  --exposure-radius 8000.0 \
  --output-dir ./tdr-output

数据下载

本 skill 可自动从 Microsoft Planetary Computer 下载 DEM 和 Sentinel-2 (无需 API key):

# 自动下载 cop-dem-glo-30 DEM + sentinel-2-l2a 水体变化快照
python scripts/tailings_dam_risk.py \
    --bbox 116.30,39.85,116.45,39.95 \
    --date-range 2023-06-01,2023-09-30 \
    --output-dir ./tdr-auto

下载的内容:

• Copernicus DEM GLO-30 (cop-dem-glo-30) — 替换 --dem-file
• Sentinel-2 L2A (sentinel-2-l2a) — 2 个时相的影像，生成水体变化，替换 --water-masks

下载元数据（data_source, fetched_at, collection, dem_path, sentinel2_paths）会写入 output-manifest.json。

Parameters

Output

Risk Levels

Risk Components

Hazard (40%)

Based on dam height, capacity, water surface change, catchment size, and runout potential.

Exposure (35%)

Based on number of objects in impact zone and exposure area.

Evidence (25%)

Based on observed water surface changes and catchment characteristics.

Key Algorithms

Water Surface Change Detection

Compares binary water masks between time periods to detect expansion and contraction areas. Calculates change percentage and extracts change polygons.

Catchment Analysis

Uses D8 flow direction algorithm on DEM to compute flow accumulation and delineate upstream watershed contributing to the dam location.

Runout Zone Estimation

Two simplified methods for screening:

• Simplified: Uses reach angle concept (height / tan(alpha))
• Energy_line: Uses H/L ratio to estimate flow extent

Both methods are SCREENING-ONLY and do not replace engineering analysis.

Downstream Exposure

Analyzes objects within potential impact zone. Counts buildings, roads, and other infrastructure at risk.

Deformation Analysis

Processes InSAR deformation data to identify areas with significant movement exceeding threshold.

Rainfall Scenario

Estimates peak flow using rational method (Q = C * i * A) for different return period events.

Exit Codes

Limitations

• SCREENING ONLY: Results are auxiliary analysis only and do NOT replace dam safety engineering assessment
• Simplified runout models are not substitutes for detailed engineering analysis
• Facility data completeness affects result confidence
• Free satellite imagery may not detect fine cracks or small deformations
• All results must be reviewed by qualified engineers before making safety decisions

Data Requirements

• DEM: Any resolution (30m SRTM or 12m PALSAR recommended)
• Water masks: Binary (0/1) rasters from NDWI/MNDWI or classification
• Deformation data: InSAR-derived rate map (mm/yr, optional)
• Facilities: GeoJSON with geometry and optional dam_height, capacity

References

• US Dam Safety Guidelines (FEMA 1998)
• ICOLD Bulletin 153: Tailings Dams - Risk of Dangerous Occurrences
• Hungr et al. 1984: Debris flow runout estimation

数据下载

本 skill 可自动从 Microsoft Planetary Computer 下载数据 (无需 API key):

python tailings_dam_risk.py --bbox 116,39,117,40 --date-range 2024-06-01,2024-06-30 --output-dir 

• --bbox W,S,E,N: WGS-84 边界框 (西, 南, 东, 北)
• --date-range START,END: 日期范围 (YYYY-MM-DD,YYYY-MM-DD)
• --aoi-file : 替代 --bbox 的 GeoJSON 多边形
• --cache-dir : 缓存目录 (默认 ~/.geoskill_cache)

当用户只给 --bbox + --date-range (没有 --dem) 时，skill 自动下载数据。 当用户给 --dem 时，走原文件路径 (向后兼容)。