Calculateur de la valeur b de Gutenberg-Richter
Le calculateur de valeur b de Gutenberg-Richter estime la pente de la relation fréquence-magnitude des séismes à partir d’un échantillon complet du catalogue.
Lancer gratuitement
Saisissez des magnitudes finies, la magnitude de complétude si elle est connue et la largeur de classe. L’outil applique l’estimateur du maximum de vraisemblance d’Aki avec la correction conventionnelle d’une demi-classe, puis fournit la valeur b, l’effectif et les valeurs intermédiaires. Il convient aux contrôles reproductibles, à l’enseignement et à l’analyse de catalogues, mais ne détecte pas la complétude et n’estime pas l’incertitude.
Préparez un échantillon complet de magnitudes
Commencez par des magnitudes de séismes provenant d’une même région d’étude, d’une même période, d’une même échelle de magnitude et d’une même chaîne de traitement. La relation de Gutenberg-Richter n’a de sens que si les événements représentent la population que vous souhaitez caractériser. Déterminez la magnitude de complétude avant de calculer la valeur b, au moyen d’une évaluation adaptée du catalogue réalisée en dehors de ce calculateur. Ne conservez ensuite que les événements situés à ce seuil ou au-dessus. Si vous omettez completeness_magnitude, le calculateur retient la plus petite magnitude fournie. Cette option est pratique pour un échantillon déjà filtré, mais ne constitue pas un test indépendant de complétude. Chaque magnitude doit être un nombre fini et l’échantillon ne peut pas être vide. Les valeurs inférieures à une magnitude de complétude explicitement indiquée sont rejetées plutôt que supprimées discrètement, car un filtrage silencieux pourrait masquer un écart entre la méthode annoncée et les données réellement analysées. Utilisez un type de magnitude homogène : le mélange de magnitudes locales, de moment, de durée ou converties peut modifier la distribution. Vérifiez aussi les doublons, les tirs de carrière et les changements de catalogue avant toute interprétation sismologique.
Comprenez l’estimation d’Aki par maximum de vraisemblance
Le calculateur commence par établir la moyenne arithmétique des magnitudes fournies. Il évalue ensuite l’estimation d’Aki en divisant le logarithme décimal de e par l’écart entre la magnitude moyenne et le seuil inférieur corrigé. Ce seuil correspond à la magnitude de complétude diminuée de la moitié de la largeur de classe. Le terme de demi-classe corrige les magnitudes discrétisées ou arrondies ; la largeur doit donc correspondre à la résolution du catalogue, généralement au pas de publication des magnitudes. Elle doit être positive, car une valeur nulle ou négative n’a aucun sens dans cette correction. Les champs intermédiaires permettent d’auditer le calcul : sample_size indique le nombre d’événements inclus, mean_magnitude consigne leur moyenne, tandis que le seuil et la largeur rappellent les paramètres effectifs. L’arrondi n’intervient qu’à la production du résultat afin de stabiliser la sérialisation. Aucune droite de régression n’est ajustée, aucun rééchantillonnage aléatoire n’est effectué et aucune donnée réseau n’est consultée. Des entrées identiques donnent donc toujours un résultat identique.
Interprétez le résultat avec prudence
La valeur b décrit l’abondance relative des petits et des grands séismes au sein de l’échantillon complet retenu. Une estimation plus élevée traduit généralement une distribution fréquence-magnitude plus pentue, tandis qu’une estimation plus faible indique une proportion relativement supérieure de grands événements. L’interprétation doit rester liée au plan d’échantillonnage. Une modification de la sensibilité du réseau, du seuil de complétude, de l’étalonnage des magnitudes, des limites spatiales, des périodes ou du traitement des répliques peut déplacer l’estimation sans changement du processus sous-jacent. Ne comparez des catalogues que si ces choix sont compatibles et clairement documentés. Ce calculateur fournit une estimation ponctuelle, mais ni erreur-type, ni intervalle de confiance, ni qualité d’ajustement, ni significativité statistique. Les petits échantillons peuvent produire des estimations instables, et un résultat numérique isolé ne démontre pas une différence physique pertinente entre régions ou périodes. Pour un travail de recherche, accompagnez l’estimation d’une analyse de complétude, d’un calcul d’incertitude, de tests de sensibilité portant sur des seuils raisonnables et de graphiques de la distribution fréquence-magnitude observée. Une requête API coûte $0.002, et ce calcul déterministe convient également à des vérifications locales répétées.
Cas d’usage
Contrôler un catalogue sismique complet
Calculez une estimation ponctuelle reproductible après avoir filtré un catalogue régional selon une magnitude de complétude documentée.
Comparer des fenêtres d’étude définies
Appliquez les mêmes paramètres à des échantillons spatiaux ou temporels compatibles avant une analyse formelle de l’incertitude.
Enseigner l’analyse fréquence-magnitude
Montrez comment la moyenne, le seuil de complétude et la correction de classe déterminent l’estimation par maximum de vraisemblance.
Questions fréquentes
Quel estimateur le calculateur emploie-t-il ?
Il emploie l’estimateur du maximum de vraisemblance d’Aki avec correction d’une demi-classe : log10(e) divisé par la magnitude moyenne moins la magnitude de complétude plus une demi-classe.
Que se passe-t-il si j’omets la magnitude de complétude ?
La plus petite magnitude fournie devient la magnitude de complétude. Cela suppose que vous avez déjà transmis un échantillon complet et filtré ; l’outil n’estime pas la complétude.
Que doit représenter la largeur de classe ?
Utilisez le pas de publication ou de discrétisation des magnitudes du catalogue. La valeur par défaut est 0.1 et doit être strictement positive.
Le résultat comprend-il une incertitude ?
Non. Le résultat est une estimation ponctuelle déterministe accompagnée de l’effectif et des valeurs intermédiaires. Calculez séparément l’incertitude pour les comparaisons scientifiques.
La liste des magnitudes peut-elle être vide ?
Non. Un échantillon vide provoque une erreur d’entrée non valide, car aucune moyenne ni valeur b ne peut être calculée sans observation.
Combien coûte un calcul par API ?
Chaque requête API coûte $0.002.
Pour les développeurs — accès API
Tout sur cette page est disponible par programmation. Cette section s'adresse aux équipes qui veulent l'intégrer à leurs systèmes ; les autres peuvent simplement utiliser l'outil ci-dessus.
Endpoint
Authentification par jeton Bearer : un seul POST met la tâche en file d’attente, et le résultat vous parvient par webhook ou lien signé.
Appeler depuis votre stack
curl -X POST https://api.kit.forhosting.com/earth/b-value-gutenberg-richter \
-H "Authorization: Bearer $KIT_KEY" \
-H "Content-Type: application/json" \
-d '{"magnitudes":[2.1,2.3,2.4,2.7,3,3.2]}'const res = await fetch("https://api.kit.forhosting.com/earth/b-value-gutenberg-richter", {
method: "POST",
headers: {
"Authorization": `Bearer ${process.env.KIT_KEY}`,
"Content-Type": "application/json"
},
body: JSON.stringify({
"magnitudes": [
2.1,
2.3,
2.4,
2.7,
3,
3.2
]
})
});
const { task_id } = await res.json();import os, requests
res = requests.post(
"https://api.kit.forhosting.com/earth/b-value-gutenberg-richter",
headers={"Authorization": f"Bearer {os.environ['KIT_KEY']}"},
json={
"magnitudes": [
2.1,
2.3,
2.4,
2.7,
3,
3.2
]
},
)
task_id = res.json()["task_id"]<?php
$res = file_get_contents("https://api.kit.forhosting.com/earth/b-value-gutenberg-richter", false, stream_context_create([
"http" => [
"method" => "POST",
"header" => "Authorization: Bearer " . getenv("KIT_KEY") . "\r\nContent-Type: application/json",
"content" => '{"magnitudes":[2.1,2.3,2.4,2.7,3,3.2]}',
],
]));
$task = json_decode($res, true);body := bytes.NewBufferString(`{"magnitudes":[2.1,2.3,2.4,2.7,3,3.2]}`)
req, _ := http.NewRequest("POST", "https://api.kit.forhosting.com/earth/b-value-gutenberg-richter", body)
req.Header.Set("Authorization", "Bearer "+os.Getenv("KIT_KEY"))
req.Header.Set("Content-Type", "application/json")
res, _ := http.DefaultClient.Do(req)Exemple de requête
{
"magnitudes": [
2.1,
2.3,
2.4,
2.7,
3,
3.2
]
}Exemple de réponse
{
"task_id": "tsk_a1b2c3d4e5f6a1b2c3d4e5f6",
"type": "earth.b_value_gutenberg_richter",
"status": "queued",
"_links": {
"result": "/tasks/tsk_…/result"
}
}L’API est asynchrone : chaque appel renvoie un task_id immédiatement, puis vous interrogez l’état à raison d’une requête par seconde.
Tarifs
Le prix est publié, sans tokens ni crédits. Une tâche qui échoue n’est pas facturée.
Erreurs
| HTTP | Code | Signification |
|---|---|---|
401 | unauthorized | Clé API absente ou invalide : vérifiez l’en-tête Authorization. |
402 | insufficient_balance | Solde insuffisant : rechargez votre compte pour lancer cette tâche. |
404 | unknown_type | Type de tâche inconnu : vérifiez le champ type de votre requête. |
429 | rate_limited | Trop de requêtes : ralentissez la cadence, puis réessayez. |