gem / oq-engine

OpenQuake Engine: a software for Seismic Hazard and Risk Analysis
https://github.com/gem/oq-engine/#openquake-engine
GNU Affero General Public License v3.0
380 stars 273 forks source link

Changed the scenario calculators to discard the sites far away from the rupture #9756

Closed micheles closed 4 months ago

micheles commented 4 months ago

This is required for consistency with Aristotle. It is a MAJOR change affecting all scenario calculations except the conditioned GMFs and the infrastructure calculations, since I changing those is trickier and we will need to think more. This PR introduce a rupture radius; the discarded sites are more distant than rupture_radius + maximum_distance. The details will likely be refined in the future. I am still assessing the magnitude of the change on realistic scenario calculations. For instance for the Chi-Chi calculation with 300 km I am getting good results:

# engine-3.19
| gsim               | weight  | contents    | nonstructural | structural  | occupants | structural+nonstructural+contents |
|--------------------+---------+-------------+---------------+-------------+-----------+-----------------------------------|
| '[Lin2011hanging]' | 0.50000 | 298_342_336 | 262_794_160   | 275_013_248 | 63.8      | 836_149_824                       |
| '[Lin2011foot]'    | 0.50000 | 208_835_200 | 189_446_944   | 202_092_272 | 45.7      | 600_374_464                       |
| Average            | 1.00000 | 253_588_768 | 226_120_552   | 238_552_760 | 54.7      | 718_262_144                       |
# engine-3.20
| gsim               | weight  | contents    | nonstructural | structural  | occupants | structural+nonstructural+contents |
|--------------------+---------+-------------+---------------+-------------+-----------+-----------------------------------|
| '[Lin2011hanging]' | 0.50000 | 298_341_536 | 262_793_312   | 275_012_608 | 63.8      | 836_147_392                       |
| '[Lin2011foot]'    | 0.50000 | 208_834_752 | 189_446_448   | 202_091_936 | 45.7      | 600_373_056                       |
| Average            | 1.00000 | 253_588_144 | 226_119_880   | 238_552_272 | 54.7      | 718_260_224                       |
# engine-3.21
| gsim               | weight | contents    | nonstructural | structural  | occupants | structural+nonstructural+contents |
|--------------------+--------+-------------+---------------+-------------+-----------+-----------------------------------|
| '[Lin2011hanging]' | 0.5000 | 316_621_728 | 265_605_792   | 282_574_144 | 63.8702   | 864_801_728                       |
| '[Lin2011foot]'    | 0.5000 | 200_433_168 | 178_104_592   | 192_266_176 | 45.8710   | 570_803_968                       |
| Average            | 1.0000 | 258_527_448 | 221_855_192   | 237_420_160 | 54.8706   | 717_802_848                       |