石油天然ガス・金属鉱物資源機構
令和2年度 地熱統括部事業成果報告会「坑井近傍探査技術」

坑井近傍探査技術の開発（H30~R2）
 サイドトラックや既設発電所での補充井のターゲット選定に寄与する、坑井を利用した高精度な断裂系探査手法の開発
AE探査
坑井近傍探査システムの運用案
地表調査時からのAE観測
AE分布・発生頻度の把握 / 断層面解に基づく透水異方性の予測

試掘成功

AE探査で別のターゲットを決定

AE: Acoustic Emission

VSP法
SSP反射法記録取得範囲
速度構造推定範囲
VSP反射法記録取得範囲
VSP: Vertical Seismic Profiling
SSP: Surface Seismic Profiling

即時サイドトラック

不成功

試掘不成功

DAS-VSP

震源再決定 /
回折源イメージによりターゲット評価

DAS-VSPで別エリアに変更

DAS-VSPでターゲット決定

翌年以降
サイドトラック

翌年以降
新掘井掘削

貯留層に逢着2

報告内容[1] DAS-VSP法の開発[2] AE探査法実証試験[3] 3軸FBG坑内地震計の開発[4] 技術開発成果の取りまとめ3

[1]DAS-VSP法の開発(株)地球科学総合研究所

光ファイバ計測技術
光ファイバ内散乱の概念図

地熱井で取得されたDTSおよびDASモニタリングデータ

Silixa社HPより

DTS：分布型温度センサ
DAS：分布型音響センサ

永田ほか（2020） レーザー光の後方散乱特性が「ひずみ」や「温度」などの外部影響により変化する現象を利用した計測技術
 耐熱性能の高い光ファイバの使用により地熱井への適用が可能
 モニタリングにおいては坑井近傍の変化を時空的かつ詳細に把握
 弾性波探査においては効率的な3D VSP 観測を実現

5

DAS-VSP法実証試験 調査測線図•••

坑井を利用した弾性波探査（VSPやAE探査）における光ファイバセンシング技術の有効性を検証するため、2019年に実証試験を実施
地表弾性波探査（SSP）との比較や統合解析を実施するため、SSPとDAS-VSPの同時取得を実施
また、マイクロサイスミックモードによる DAS/SSPの夜間連続観測を実施
坑跡沿い断面図

平面図
北西

約 2 m 間隔で1050 の受振点を設定

南東
光ファイバセンサ設置区間

2,143m• ENE-WSW、NW-SE方向の断層およびリニアメントが顕著• 噴気・温泉・変質体等の地表地熱兆候も両方向の断裂に規制

500m• 試験地には地熱貯留層を形成する高傾斜の北東-南西走向の断層が発達しており、使用した坑井はこれら断層を通過＜試験概要＞
・発振点数： 314 ＠50m間隔
・受振点数： 地表： 778＠25 m間隔、坑内： 1050＠2 m間隔
・作業期間：約2週間 （計測作業は6 日間、夜間観測は約12時間）

6

DAS-VSP法実証試験における基本処理・解析結果
低速度部
検層記録

速度モデル

飽和蒸気上面深度

Velocity
5000 [m/s]

水頭面深度
沸騰区間底面深度

DTS

ZOVSP

+1000

DTS
Temperature
Gradient-

反射法記録

DAS-VSP
SSPAmplitude

反射イベントの食い違い
青木ほか(2020)

実証試験成果の統合断面図• DAS-VSP / 地上震探（SSP）の統合トモグラフィ解析により速度モデルを推定し、イメージング（3D PSDM）を通してそれぞれの反射法記録を作成。
• DTSから認定される飽和温度区間にはZero-offset VSPのP波速度において低速度異常が顕著であり、トモグラムを通してその広がりの把握が可能であった。

断層面からの反射波

SSP / VSP 反射法断面の比較（上）SSP記録において断層は、地層面反射の食い違いとして表現されている。（下）DASVSP 記録は相対的に分解能が高いほか、深部断層面からの反射波が捉えられていた。

7

DAS-VSP波形記録を利用した速度構造推定の高度化
課題••

観測波形情報を用いた逆解析手法（FWI）の採用による速度推定の高精度化
DASの感度特性（光ファイバーケーブル長手方向への指向性）への対応

提案手法
••

低周波数から段階的に解析周波数を上げていくマルチスケール法（Bunks 1995 ）の採用
観測データ、フォワードモデリングデータの双方においてローカルな正規化を行うNormalized misfitfunction (Groos, 2013)の採用

FWI解析ワークフロー

坑井速度（ZO-VSP)と速度モデルの比較
FWI解析結果

トモグラフィ解析結果

比較• FWI解析により、速度の整合性が向上したことを確認（白矢印で示した箇所）
8

開口断裂系探査法の開発
き裂と関連する波動場のスナップショット（シミュレーション）

【背景】

X

•

•

き裂や断層を地震波が通過した際、回折波や散乱波などの二次的な波動場が形成

DAS

PS変換波
Depth•

これら二次的な波動を近傍の坑井内のDASにより検出できれば、き裂の位置や性状の評価に利用できる可能性あり
き裂と関連するイベントは、P波・S波
直接波の後続波として観測され、その振幅値は１～２オーダー小さいと予測

坑井

PP回折波
き裂
PS変換波

P直達波

PP回折波
（弾性波：Z成分）

坑井位置での合成観測記録（Z成分）
Z成分（DAS記録相当）き裂からの応答二次的な波動を抽出

×10-17×10-15
Time

Depth

PS変換波

P波初動
S波初動
スナップショット時刻

PP回折波
Time・振幅値は1~2オーダー小さい

9開口断裂系探査法の開発【DASによる開口断裂系探査の方法】• 坑井内のDASによりき裂や断層に関連する回折波や散乱波を捉え、その発生源をイメージングすることにより、き裂群が発達する深度、方位、距離を特定。
き裂からの弾性波

直達波の通過振動
方向

DASセンサ

S波の進行方向き裂に関連するイベント抽出と 反射点またはイメージング
変換点の位置＞数百m（断裂系指標）

N

回折源位置

坑跡
・坑井内のDASは、側方から到来するS波の信号の検出に向くと考えられる。

回折強度（平面図）（3D）
・イベント発生源の可視化により、き裂の位置が特定される可能性あり。

10開口断裂系探査の手順
2種（受動/能動）の開口断裂系探査ワークフロー【手順】
DAS を用いて収録されたAE記録やバイブレータ震源記録から、以下の手順によりき裂の兆候を検出し、位置を特定① 速度モデル作成
・ 震源決定、回折源イメージングで使用する速度モデル
をDAS-VSPによる坑井速度測定等を用いて作成② 回折波の抽出と強調
・ 回折波は直達波や反射波に対し微弱であるため、直達波や反射波を抑制、またはイベントを強調④ 回折波イメージング・ イベントのモードに対応するマイグレーション処理を適用
例）S波後続のS波イベントを利用する
場合、S波マイグレーション処理を適用

[ch]

AE

DAS 記録例

バイブレータ震源

[ch]

Depth

③ 震源決定（AE利用の場合）
・ AEの波形をイメージングで利用可能とするため、震源決定を行い、震源の位置と発生時刻を推定

Depth

•

[s]

大型2台、32 Stack[s]11

開口断裂系探査法（回折源イメージング）の適用例【回折波の抽出】

反射法記録（SSP）と回折源イメージ（AE利用）の統合表示例• 回折強度分布について、比較されるSSPに基づく断裂系の情報（反射法断面における構造変化、断層尤度分布）と調和的に分布

2000

回折波イベント強調記録例
0

検出された
イベント

Depth [m]【回折源イメージング結果】

0Depth [m]• 相対的に大きな振幅を持つ直達波や反射波を抑制したのち、S波速度に対応する曲率を持ったイベントを探索し、その曲面に沿った平滑化を行いイベントを強調

DASによるAE観測記録例

強調

122000深度断面での比較例
回折強度分布
断層尤度分布（SSPより）• 図中の水色破線の位置などに、断裂尤度と回折強度分布の良好な相関を確認
寺西ほか（2021b） 12

[1]DAS-VSP法の開発：成果と課題（まとめ） 成果►DAS-VSPの反射法処理/速度構造推定の高度化を達成►開口断裂系探査法（AE利用）の実証• 回折波抽出法と回折源イメージング法を提示
 課題►開口断裂系探査法の適用範囲の拡大• バイブレータ震源を利用した探査法の確立– 回折波記録の効果的取得– 同時発振による作業時間の短縮• 人工震源を利用した開口断裂系探査法の実証
13

[2]AE探査法実証試験
A, B 地域（国研）産業技術総合研究所
C地域(株)地球科学総合研究所
地熱技術開発株式会社

AE波形の処理フロー
機械読取による自動ピック

連続AE波形目視によるピック補正
標準的な震源決定プログラム「hypo71」の発展形

HYPOINVERSE-2000を利用した通常震源決定
震源分布

波形相関を利用した高分解能（10mオーダー）な震源決定

AEイベントのペア間の到達時間差（波形相関を利用して計算）に基づいて，相対位置を修正し，高繊細な震源分布を求めるDD法を適用

DD法を利用した
高精度震源決定詳細震源分布・震源メカニズム
解析を適用
・非DC成分を推定し，AEへの流体の寄与を判断

非DC成分を考慮した
震源メカニズム解析

トモグラフィ解析

AEの流体由来成分

速度構造

地震波トモグラフィを行い，地下の速度構造を推定

貯留層構造・挙動の評価
15

非DC成分を考慮した震源メカニズム解析【非DC（ダブルカップル）成分】• 一般的な断層活動は，せん断破壊で生じるが，地熱地域のような流体流動が豊富な地域では，非せん断破壊の成分も持つ場合がある（例えば，Vavryčuk and Hrubcová, 2017）

流体流出による圧縮

流体流入による引張（上段）地熱地域のような流体流動が豊富な場所では，流体流出入に伴う破壊成分を含む

通常の断層運動によるせん断（下段）一般的な自然地震は、せん断破壊で生じる
Non-DC成分とDC成分のイメージ図（USGS HP）【震源メカニズム解析手法】• 流体流出入に伴う破壊の成分の推定には，通常の断層パラメータである「走向」，「すべり角」，「傾斜」の他に，「テンサイル角」という，断層面とすべった方向の偏角を表すパラメータを考慮「テンサイル角」
断層面とスリップベクトルの偏角→ 流体流出入による引張・圧縮を反映
テンサイル角が正 ：流体流入成分を持つ
テンサイル角が 0 ：通常の断層運動
テンサイル角が負 ：流体流出を持つ

法線ベクトルスリップベクトル

断層面テンサイル角 a16

A地域のAE観測網
MKK▼：AIST観測点▽：既存観測点

KM-4（坑内） KKD12 KM-8
KKD16（坑内）KKD15SDZ

KKD14

KM-5

KKD01-11(Array)

THB

KM-1（坑内）

KKD13

KM-3（坑内）

観測期間

センサー

固有周期

感度劣化状況

KKD01-11, SDZ, THB

2017/12～

KVS-300（速度計）

2 Hz

異常なし

MKK

2018/6～

同上

同上

異常なし

KKD12-16

2020/8～

同上

同上

異常なし

KM-1, 3, 4, 5, 8

既存観測点

L-22D, E（加速度計）

同上

KM-4: 40%
KM-8: 40%17

定期点検時（坑井停止・立ち上げ時）の特徴的なAE活動
定期点検中3カ月

定期点検前3カ月

青：生産井（定検時停止）
緑：還元井（定検時停止）
新たに活発なAE活動
定期点検中のAE分布• 標高-2000 m 付近の生産井周辺で、活発なAE活動• 標高0m付近のAE活動は、点検前から継続して活発••

定期点検に伴い，複数の坑井を停止・立ち上げをした際，AE活動が活発化
特に2018年の定期点検時には，坑井停止に伴う深部貯留層（海抜-2000m程度）でAE活動が顕在化

18定期点検1週目（坑井停止時）のAE挙動
坑井近傍の流体挙動

AA’
Well-48

Well-48の停止日・坑底を原点とした震源の時空間分布

DD決定震源
流体由来の震源成分☆最初のイベント

浸透率

深部貯留層付近でAE活動あり

流体流入成分

流体流出成分

A＜坑井近傍の流体挙動＞
・Well-48の坑底付近を起点として，徐々にAEが拡散する傾向を確認＜震源の時空間分布＞
・坑井停止後，数日間かけて間隙圧変化が周囲に拡散・震源の中には、流体が押し出されることによる流体流出成分を含むものが存在
・Well-48底からのAEの拡散は，浸透率で2×10-14 m2程度に相当

A’

断面図の位置19

定期点検1週目（坑井停止時）のAE挙動
既存の数値計算（Sakagawa et al. (2000)）で用いられた浸透率との重ね合わせ

P波速度のゆらぎ構造との重ね合わせ

A

A’

AA’

AEはやや高速度帯に集中••

AEは，周囲よりややP波速度が高い空間に集中
水で飽和（＝より速い領域）した箇所でAEが集中していることを示唆•

既存研究で用いられた浸透率と概ね整合的

20

定期点検9-10週目（坑井立ち上げ時）のAE挙動
B

B’

DD決定震源
流体由来の震源成分☆最初のイベント

最初のイベント時間・場所を原点とした震源の時空間分布

浅部貯留層付近でAE活動あり

浸透率

流体流入成分

流体流出成分

B＜坑井近傍の流体挙動＞
・坑井立ち上げ時，浅部貯留層側でAEが活発＜震源の時空間分布＞
・坑井停止後，数日間かけて間隙圧変化が周囲に拡散・AEの拡散は，浸透率で2×10-14 m2程度に相当

B’
断面図の位置21

定期点検9-10週目（坑井立ち上げ時）のAE挙動
既存の数値計算（Sakagawa et al. (2000)）で用いられた浸透率との重ね合わせ

P波速度のゆらぎ構造との重ね合わせ

B

B’

B

B’

AEはやや高速度帯に集中••

AEは，周囲よりややP波速度が高い空間に集中
水で飽和（＝より速い領域）した箇所でAEが集中していることを示唆•

既存研究で用いられた浸透率と概ね整合的

22

A地域でのAE観測まとめ• 最大で24観測点を用いた稠密観測を実施• 定期点検時に，坑井停止に伴う深部貯留層近傍でのAE活動の活性化を確認• 定期点検時に，坑井立ち上げに伴う浅部貯留層近傍でのAE活動の 活性化を確認• 何れも，AE活動の拡がりは拡散方程式に従うように見え、透水率を計算した
結果，既存研究と調和的な値を示す

23

B地域のAE観測網

B-2
B-4井
B-1

B-3
B-4
B-6（坑内）KJ-5井
B-5・B-5, B-4井での噴気試験に伴うAEを捉えることを目的に、2019年10月より5点の臨時地表観測点で連続観測
・11月末～12月にかけては，坑内観測点を臨時追加

観測期間

センサー

周波数レンジ

B-1

2019/10～2020/10

Trillium Compact PH 20s（速度計）

20s ～ 108Hz

B-2

2019/10～2020/10

同上

同上

B-3

2019/10～2020/10

同上

同上

B-4

2019/10～2020/10

同上

同上

B-5

2019/10～2020/6

同上

同上

B-6

2019/11/25～
2019/12/26

ASR-1（速度計）10Hz ～ 180Hz

24B-5噴気試験前後のAE分布
B-1坑底での温度，圧力モニタリング値
10/27-11/2（噴気前）

B-5噴気情報11/3-11/11
（噴気中）

噴気開始
年月日

坑口圧力 蒸気流量 熱水流量 合計流量
(MPaG)
(t/h)
(t/h)
(t/h)

2019/11/3
温度（℃）

圧力（MPa)•B-1

384078

噴気に伴い，モニタリング井（B-4, B-5の間）の温度・圧力値が変動

噴気試験中にKJ-4,KJ-5井間でAEが活発化

B-40.60

B-4B-1
B-5

B-5

11/3-11/11（噴気中）

10/27-11/2（噴気前）•

噴気後，B-5, B-4井間で複数のAE活動を確認

25

全観測期間を通してのAE分布（B-4 ～ B-5間）•B-5噴気試験の期間外でも，イベントは発生•震源域から、高角な構造を示唆（下左図）•AE群は周囲よりやや低速度帯に位置し、断層・破砕帯のような脆弱な構造であることを示唆※（下右図）•破砕帯上で複数の流体寄与のイベントあり（下左図）

全期間B-4※A地域の場合は，水飽和の存在によりき裂がある領域（AE発生域）が周りより高速度で見えていたが，B地域は蒸気卓越であるため，亀裂が卓越した領域は周囲より低速度に見える。

B-5

トモグラフィ結果（P波速度ゆらぎ）と重ね合わせ○ DD決定震源○ 流体由来の震源成分
B-5
B-4
流体流入成分

流体流出成分26

全観測期間を通してのAE分布（B-4北側，観測網南西部）
B-4北側のAE群
B-4
B-5
南西部のAE群
P波速度ゆらぎ（50m深度ごとに計算）• 低速度で特徴づけられる領域に，南西45度落ちの構造体

P波速度ゆらぎ（50m深度ごとに計算）• B-4北側と同様、低速度で特徴づけられる 領域に，南西45度落ちの構造体• 流動由来のAEイベントが豊富に見られ，豊富な流動を示唆
流体流入成分

流体流出成分27

B地域でのAE観測まとめ• B-5噴気試験に関連してB-4, B-5坑井間にやや活発なAE群が見出された►噴気試験の期間外にも本AE分布は出現するため，噴気試験との明瞭な関係性は不明►AE発生域は断層域に相当• AE分布より，坑井近傍に3つの構造体を確認した► B-4 – B-5井間の構造体 ，および 観測網南西部の構造体では，メカニ
ズム解析結果より、流体の寄与を示唆►B-4北側の構造体は，過去のAE調査（～2018年度）では未検出

28C地域：2種類の観測網による試験を実施
C地点実証試験観測点配置図

 坑井内+地表観測網► 既設の3点（マップ非表示）に対し以下を追加• 坑井内1地点（深度150m）＠発電所敷地内• 地表5地点＠発電所から2 kmの距離► 観測期間：2020年7月16日~10月30日★A★
E
発電・還元基地

生産基地
地熱発電所

 稠密臨時観測網

★
W

★B

★D

► 2測線：各45観測点＠75m間隔（各3 km）

半径2ｋｍ

★C

► 弾性波探査用の簡易機材を多点配置► 観測期間：2020年7月16日~22日

凡例★ 地表観測点 A~E （5点）
★ 坑井内観測点 W （1点＠-148 m）▼ 稠密臨時観測点 （90点）

29

坑内地震計の設置
観測井ケーシングプログラム

坑内地震計 編成および仕様

930
1,950

予定深度250m, 溢
水と使用掘削機の限界から，協議の上，158mで掘止地震計設置位置
（坑底）

項 目
計測器
周波数
感 度
耐 熱
耐 圧
適用傾斜
サイズ
坑内固定
信号ケーブル

仕 様
直交3軸ジオフォン
15 Hz
1.04 V/cm/sec
150 ℃
34.47 MPa
最大 90度
長さ930mm×直径88mm,重量13kg
油圧式固定アーム（ハイドロリックホールロック）
材質：Tefzel
直径：10.95mm
引張り強度：450kg

30各観測網でのイベント比較
生産領域近傍で発生したAEの記録例

気象庁により震源が推定された自然地震の記録例
ZX
Y

震央の推定位置461909046454645

•
•

Time

坑井内観測点は地表付近のノイズの影響を受けず、AEイベント検出に有利。
稠密観測記録もAEの到来方向による見かけ速度の違いからイベント種別の判読が容易であり、坑井内観測点と同様の長所を有する。
例えば雷鳴はまばらな地上観測網のみではAEイベントと判断される可能性があるが、坑井内観測点では検出されないことから地上のノイズと判断され、クロス状稠密観測網
ではその見かけ速度（右図）から音波と判明する。

1Time
1

地上の音波ノイズ記録例
雷鳴＋降雨

45 4690

1

45 4690Time•

マグニチュード：2.990
落雷地点

14531

稠密臨時観測網の震源メカニズム解析への寄与（２）震源メカニズム解を得たAEの内訳（１）AEイベント数抽出されたAEの数
DD法による震源再決定
メカニズム解を得た数

長期＋臨時観測網
長期観測網のみ

419
164
118

9/9109/155※ 分母はDD法により求められた震源の数を示す。

（３）震源メカニズム解の決定精度比較
ランク
AB
CD総数

••
•

長期+臨時観測網
イベント数
割合017
190%
11%
78%11%

長期観測網のみ
イベント数割合

024958
109

0%
2%
45%
53%

決定
精度

RMS fault plane
uncertainty

A

≧ 25°

B

≧ 35°

C

≧ 45°

D

azimuth gap ≦ 90,
takeoff angle gap ≦ 60

解析は Hardebeck and Shearer (2008) によるP波初動を用いた震源メカニズム解析プログラム「HASH」を使用。
推定結果はその決定精度に応じ，A（最良）からE（決定不可能）にランク付けされる。
臨時観測点（90点）が追加された期間はCが78％を占めた一方、長期観測網のみ（坑井内１点＋地上8点）では半数以上がランクDであった。

稠密臨時観測網は、震源メカニズム解析の精度向上に寄与

32

震源分布と断層メカニズム解
平面図

南北断面での比較

走向分布
横ずれ断層

a b c d e
S

正断層

逆断層

生産領域

• 走向分布のトレンド（(b)黒実線）には地形との相関が認められるN

主引張軸（T軸）分布

主圧縮軸（P軸）分布

還元領域

生産井の西側

•••生産・還元領域においてAEの多くは海抜0m～-1,000mの範囲で発生
様々なタイプの断層が混在[b] を除き震源は北に深くなる傾向

プランジ角が30度以下の方位分布•

P軸、T軸が複雑に分布する領域が存在

33

C地域でのAE観測まとめ
 2種類の観測方式を適用► 坑井内+地表のAE観測網• 地表付近のノイズの影響を受けず、AEイベント検出に有利► 稠密臨時観測網• AEイベントの判別が容易• 震源メカニズム解析の精度向上に寄与

 発電所周辺のAE分布► 地形と相関のある断裂系配置を示唆► 震源は北に深くなる傾向► P・T軸が複雑に分布する領域が存在
34

[3]３軸FBG坑内地震計の開発

地熱技術開発株式会社

３軸FBG坑内地震計•

地熱井用高温地震計としてFBGセンサーによる３軸坑内地震計を開発。
（FBG：Fiber Bragg Grating）

室内試験の様子ハイドロリックホールロック

FBGの計測原理

3軸FBGセンサー
坑内地震計

カンチレバー式
ひずみ増幅器36

３軸FBG坑内地震計
3軸FBG坑内地震計の仕様
No

項目

1

光ファイバー
坑内地震計

2信号ケーブル34

インテロゲータ

ハイドロリック
ホールロック

仕 様
型式：FOB-110
成分：3成分
方式：光学式加速度計
外径：110mm
外装：SUS304
光ファイバーケーブル
外径：7mm
長さ：100m
型式：EFOX-1000B
チャンネル数：4ch
サンプリング：1kHz
波長レンジ：100nm
機構：油圧作動式
外装：SUS304
耐熱：150℃
油圧ホース（耐熱：150℃）

光測定器
インテロゲータ

信号ケーブル（光ファイバー）

ハイドロリック
ホールロック3軸FBGセンサー
坑内地震計（FOB-110)37

[4]技術開発成果の取りまとめ(株)地球科学総合研究所（国研）産業技術総合研究所
地熱技術開発株式会社

坑井近傍探査システムのガイドラインの作成
 AE探査法、DAS-VSP法についての基本的事項とともに一歩進んだ使い方、調査費用の概要を例示
坑井内地震計設置作業例

冷却スプレーによるfreezing test

集約型システムの例
39

成果のまとめ[1] DAS-VSP法の開発• DAS-VSPの反射法処理/速度構造推定の高度化を達成• DASによる開口断裂系探査法の有効性を実証[2] AE探査法実証試験• 稠密観測網による透水構造把握の可能性を実証• 稠密観測網による貯留層内外の断裂構造把握の可能性を実証• 坑井内観測点および稠密臨時観測網の有効性を確認[3] 3軸FBG坑内地震計の開発• FBGセンサーによる3軸坑内地震計を開発[4] 技術開発成果の取りまとめ• 坑井近傍探査システムのガイドラインを作成（近日公開予定）
40

参考文献
【DAS-VSP法】
 Fomel, S., (2002): Applications of plane-wave destruction filters. Geophysics. 67, 1946–1960
 Hale, D. (2013). Methods to compute fault images, extract fault surfaces, and estimate fault
throws from 3D seismic images. Geophysics, 78(2), O33–O43.
 Jager, R., Mann, J., Hocht, G. and Hubral, P. (2001): Common-reflection-surface stack: Imageand attributes. Geophysics, 66, 1, 97-109.
 Schwarz, B. (2019): Coherent wavefield subtraction for diffraction separation. Geophysics, 84,
3, 157-168
 青木直史・藤澤萌人・寺西慶裕・永田丈也・毛利拓治・小野寺真也（2020）： 光ファイバセンサ
による坑井近傍の断裂系探査，物理探査学会第143回学術講演会
 寺西慶裕・青木直史・藤澤萌人・清水英彦・小野寺真也（2021a）： 光ファイバを用いた断裂系探査法の開発，2021年度地球惑星連合大会
 寺西慶裕・青木直史・藤澤萌人・小野寺真也（2021b）： DAS夜間観測記録を用いた回折波抽
出・可視化手法，物理探査学会第144回学術講演会
 永田丈也・手塚和彦・青木直史・毛利拓治・小野寺真也（2020）: 地熱井における分布型光ファイバー計測－DTSおよびDASの適用例－, 日本地熱学会令和2年度学術講演会.

41