stress mapping

Contribution of mine borehole data toward high-resolution stress mapping: An example from northern Bowen Basin, Australia - ScienceDirect

この文献を見て、初期地圧の解釈で悩んでいたころを思い出しました。

知らなかったことが2点。
1つ目は BTV での割れ目を利用できること。これは、応力解放法の簡易版と解釈すれば良いのでしょう。ある程度深くないと揃わないようですね。100～150m程度だとばらつくのも当然です。

2つ目は平面的な補間方法。MATLAB のスクリプトになっているぐらい、頻繁に使われているようです。
Stress2Grid - File Exchange - MATLAB Central (mathworks.com)

This script is the update of Stress2Grid v1.0 (Ziegler and Heidbach, 2017). It provides two different concepts to calculate the mean S_Hmax orientation on regular grids. The first is using a fixed search radius around the grid points and computes the mean S_Hmax orientation if sufficient data records are within the search radius. The larger the search radius the larger is the filtered wavelength of the stress pattern. The second approach is using variable search radii and determines the search radius for which the standard deviation of the mean S_Hmax orientation is below a given threshold. This approach delivers mean S_Hmax orientations with a user-defined degree of reliability. It resolves local stress perturbations and is not available in areas with conflicting information that result in a large standard deviation.

文献では後者を選択しています。こういった方法で空間分布を表現するのが適当だったのですね。勉強になりました。