KR101693311B1 - 셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법 - Google Patents
셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법 Download PDFInfo
- Publication number
- KR101693311B1 KR101693311B1 KR1020150155007A KR20150155007A KR101693311B1 KR 101693311 B1 KR101693311 B1 KR 101693311B1 KR 1020150155007 A KR1020150155007 A KR 1020150155007A KR 20150155007 A KR20150155007 A KR 20150155007A KR 101693311 B1 KR101693311 B1 KR 101693311B1
- Authority
- KR
- South Korea
- Prior art keywords
- wave
- equation
- acoustic wave
- pure acoustic
- pure
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 96
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 239000006185 dispersion Substances 0.000 claims description 8
- 230000010354 integration Effects 0.000 claims description 6
- 230000000704 physical effect Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 2
- 230000035945 sensitivity Effects 0.000 description 35
- 238000010586 diagram Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/284—Application of the shear wave component and/or several components of the seismic signal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/58—Media-related
- G01V2210/586—Anisotropic media
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/67—Wave propagation modeling
- G01V2210/673—Finite-element; Finite-difference
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Geology (AREA)
- Pure & Applied Mathematics (AREA)
- Geophysics (AREA)
- Data Mining & Analysis (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법이 개시된다. 개시된 순수음향파동모델링 방법은 이방성 매질에서의 S파가 분리된 순수음향파동방정식을 설정하는 단계; 순수음향파동방정식을 행렬형태로 표현하는 단계; 행렬형태로 표현된 순수음향파동방정식을 유한차분법으로 표현하는 단계; 를 포함한다.
Description
본 발명은 셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법으로써, 구체적으로는 물성정보가 셀에 위치하여 일반적인 유한차분법에 비해 안정적이며, 자유면 경계조건을 잘 만족시키는 셀기반 유한차분법을 순수 음향파동 방정식에 적용한 순수음향파동모델링 방법에 관한 것이다.
일반적인 유한차분법으로 주파수 영역에서 순수음향파동방정식을 계산하여 모델링을 수행할 경우, 계수행렬의 대칭성을 만족시키지 않기 때문에 모델링의 안정성이 떨어지게 되어 임의의 위치에 송수신기에서 획득한 탄성파 탐사자료의 상반성(송-수신기의 각 위치가 서로 바뀌어도 동일한 탄성파 탐사자료를 얻을 수 있다는 특성)이 깨어지며, 유한요소법 적용 시 4차 미분항을 계산하기 위하여 많은 메모리를 필요로 하는 단점이 있다.
이로 인해 주파수 영역 전파형 역산에서 사용되는 민감도 행렬(Jacobian matrix)의 계산 시 상반성을 적용하지 못하여 계산의 효율이 떨어지게 되는 문제점이 있다.
상기 문제점을 해결하기 위해, 본 발명에 따른 셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법은 물성정보가 셀에 위치하여 일반적인 유한차분법에 비해 안정적이며, 자유면 경계조건을 잘 만족시키는 셀기반 유한차분법을 순수 음향파동 방정식에 적용가능하도록 개발하고, 이를 적용하여 상반성이 깨어지지 않도록 순수 음향파동 전파를 모사하는 것을 하나의 목적으로 한다.
본 발명의 해결과제는 이상에서 언급된 것들에 한정되지 않으며, 언급되지 아니한 다른 해결과제들은 아래의 기재로부터 당해 기술분야에 있어서 통상의 지식을 가진 자에게 명확하게 이해될 수 있을 것이다.
상기 목적을 달성하기 위해, 본 발명은 (a) 이방성 매질에서의 S파가 분리된 순수음향파동방정식을 표현하는 단계와, (b) 순수음향파동방정식을 행렬형태로 표현하는 단계와, (c) 행렬형태로 표현된 순수음향파동방정식을 유한차분법으로 표현하는 단계를 포함한다.
본 발명에서 (a) 단계는 이방성 매질에서 P파와 SV파의 분산관계방정식을 표현하는 단계와, P파와 SV파 각각의 독립방정식을 표현하는 단계와, P파의 독립방정식으로부터 순수음향파동방정식을 표현하는 단계를 포함한다.
본 발명에서 P파와 SV파의 분산관계방정식은 하기 식 1로 표현되는 것이 바람직하다.
[식 1]
(: 각주파수; : P파의 수직속도; : P파의 수평속도; : P파의 시간차(moveout) 속도; : sv파의 수직속도; ,: 톰슨(Thomsen)의 이방성변수; , , : 공간파동수)
본 발명에서 P파와 SV파 각각의 독립방정식을 표현하는 단계는 분산관계를 P파와 SV파에 대한 독립방정식으로 인수분해하는 것이 바람직하다.
본 발명에서 P파에 대한 독립방정식은 하기 식 2로 표현되고, SV파에 대한 독립방정식은 하기 식 3으로 표현되는 것이 바람직하다.
[식 2]
[식 3]
본 발명에서 P파의 독립방정식으로부터 순수음향파동방정식을 표현하는 단계는 P파의 독립방정식으로부터 1차 테일러급수 전개(the first-order Taylor series expansion)를 이용하여 순수음향파동방정식이 표현되는 것이 바람직하다.
본 발명에서 P파의 독립방정식으로부터 1차 테일러급수 전개(the first-order Taylor series expansion)를 이용하여 표현된 순수음향파동방정식은 하기 식 4 또는 식 5인 것이 바람직하다.
[식 4]
[식 5]
본 발명에서 P파의 독립방정식으로부터 순수음향파동방정식을 설정하는 단계는 P파의 독립방정식으로부터 암시적 시간적분 기법(the implementation of the implicit time integration methods)를 이용하여 순수음향파동방정식을 표현하는 것이 바람직하다.
본 발명에서 P파의 독립방정식으로부터 암시적 시간적분 기법(the implementation of the implicit time integration methods)를 이용하여 표현된 순수음향파동방정식은 하기 식 6 또는 식 7 인 것이 바람직하다.
[식 6]
[식 7]
본 발명에서 행렬형태로 표현된 순수음향파동방정식은 하기 식 8 또는 식 9 인 것이 바람직하다.
[식 8]
[식 9]
( : 각주파수; : P파의 수직속도; : P파의 수평속도; : P파의 시간차(moveout) 속도;: sv파의 수직속도: , : 톰슨(Thomsen)의 이방성변수; I : 단위행렬; f : 송신원 ; P : 파동장 ; Q : 가상파동장 )
본 발명에서 (c) 단계는 한 점의 주위를 둘러싼 네 셀의 물성정보를 이용하여 순수음향파동방정식을 유한차분법으로 표현하는 것이 바람직하다.
본 발명에서 유한차분법을 이용하여 표현된 순수음향파동방정식의 4차 미분항은 하기 식 10인 것이 바람직하다.
[식 10]
( : 각주파수; : P파의 수직속도; : P파의 수평속도; : P파의 시간차(moveout) 속도; : sv파의 수직속도; : 절점 (i, j)번째에 인접한 셀의 물성정보인 , , 과 같은 음향파동속도와 이방성 파라미터; U : 파동장 )
본 발명에 따른 셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법은 탄성파 탐사자료의 상반성이 깨지는 것을 방지할 수 있는 효과를 얻을 수 있다.
또한, 기존 방식에 비하여 적은 메모리를 이용하여 순수음향파동방정식 모델링이 가능한 효과를 얻을 수 있다.
또한, 상기 효과로 인하여, 모델링의 안정성을 획득할 수 있는 효과를 얻을 수 있다.
본 발명의 효과는 이상에서 언급된 것들에 한정되지 않으며, 언급되지 아니한 다른 효과들은 아래의 기재로부터 당해 기술분야에 있어서의 통상의 지식을 가진 자에게 명확하게 이해될 수 있을 것이다.
도 1은 본 발명의 제 1 실시예에 따른 순수음향파동모델링 방법의 순서도이다.
도 2는 도 1에 도시된 이방성 매질에서 S파가 분리된 순수음향파동방정식을 표현하는 순서도이다.
도 3은 본 발명의 제 1 실시예에 따른 순수음향파동모델링 방법에서 셀(cell)안에 물성정보가 포함된 셀기반 격자세트(cell-based grid sets)를 나타내는 도이다.
도 4는 본 발명의 제 1 실시예에 따른 순수음향파동모델링 방법에서 상반성을 비교하기 위한 속도모델을 나타내는 도이다.
도 5a는 각 송,수신기쌍에서 획득한 탄성파 파동장을 일반적인 유한차분법을 사용한 순수음향파동모델을 사용하여 나타낸 도이다.
도 5b는 각 송,수신기쌍에서 획득한 탄성파 파동장을 본 발명의 일실시예에 따른 순수음향파동모델을 사용하여 나타낸 도이다.
도 6은 본 발명의 제 2 실시예에 따른 전파형 역산방법의 순서도이다.
도 7a는 범위의 파선변수에 따른 P파 수직속도 민감도를 나타낸 도이다.
도 7b는 범위의 파선변수에 따른 P파 수직속도 민감도를 나타낸 도이다.
도 7c는 범위의 파선변수에 따른 P파 수직속도 민감도를 나타낸 도이다.
도 7d는 범위의 파선변수에 따른 P파 수직속도 민감도를 나타낸 도이다.
도 8a는 범위의 파선변수에 따른 P파 수평속도 민감도를 나타낸 도이다.
도 8b는 범위의 파선변수에 따른 P파 수평속도 민감도를 나타낸 도이다.
도 8c는 범위의 파선변수에 따른 P파 수평속도 민감도를 나타낸 도이다.
도 8d는 범위의 파선변수에 따른 P파 수평속도 민감도를 나타낸 도이다.
도 9a는 범위의 파선변수에 따른 P파 시간차속도 민감도를 나타낸 도이다.
도 9b는 범위의 파선변수에 따른 P파 시간차속도 민감도를 나타낸 도이다.
도 9c는 범위의 파선변수에 따른 P파 시간차속도 민감도를 나타낸 도이다.
도 9d는 범위의 파선변수에 따른 P파 시간차속도 민감도를 나타낸 도이다.
도 2는 도 1에 도시된 이방성 매질에서 S파가 분리된 순수음향파동방정식을 표현하는 순서도이다.
도 3은 본 발명의 제 1 실시예에 따른 순수음향파동모델링 방법에서 셀(cell)안에 물성정보가 포함된 셀기반 격자세트(cell-based grid sets)를 나타내는 도이다.
도 4는 본 발명의 제 1 실시예에 따른 순수음향파동모델링 방법에서 상반성을 비교하기 위한 속도모델을 나타내는 도이다.
도 5a는 각 송,수신기쌍에서 획득한 탄성파 파동장을 일반적인 유한차분법을 사용한 순수음향파동모델을 사용하여 나타낸 도이다.
도 5b는 각 송,수신기쌍에서 획득한 탄성파 파동장을 본 발명의 일실시예에 따른 순수음향파동모델을 사용하여 나타낸 도이다.
도 6은 본 발명의 제 2 실시예에 따른 전파형 역산방법의 순서도이다.
도 7a는 범위의 파선변수에 따른 P파 수직속도 민감도를 나타낸 도이다.
도 7b는 범위의 파선변수에 따른 P파 수직속도 민감도를 나타낸 도이다.
도 7c는 범위의 파선변수에 따른 P파 수직속도 민감도를 나타낸 도이다.
도 7d는 범위의 파선변수에 따른 P파 수직속도 민감도를 나타낸 도이다.
도 8a는 범위의 파선변수에 따른 P파 수평속도 민감도를 나타낸 도이다.
도 8b는 범위의 파선변수에 따른 P파 수평속도 민감도를 나타낸 도이다.
도 8c는 범위의 파선변수에 따른 P파 수평속도 민감도를 나타낸 도이다.
도 8d는 범위의 파선변수에 따른 P파 수평속도 민감도를 나타낸 도이다.
도 9a는 범위의 파선변수에 따른 P파 시간차속도 민감도를 나타낸 도이다.
도 9b는 범위의 파선변수에 따른 P파 시간차속도 민감도를 나타낸 도이다.
도 9c는 범위의 파선변수에 따른 P파 시간차속도 민감도를 나타낸 도이다.
도 9d는 범위의 파선변수에 따른 P파 시간차속도 민감도를 나타낸 도이다.
본 발명의 구체적인 설명에 앞서 후술되는 내용 중, ‘상반성’은 송-수신기의 각 위치가 서로 바뀌어도 동일한 탄성파 탐사자료를 얻을 수 있다는 특성을 뜻한다.
또한, 본 발명의 제 1 실시예인 순수음향파동모델링 방법은, 지하매질구조를 추정하기 위한 장치에서 지하 매질을 대상으로 측정된 데이터와 상기 지하 매질에 대한 초기 모델을 모델링한 모델링 데이터를 획득하는 데이터 획득부에서 수행된다.
또한, 본 발명의 제 2 실시예인 전파형 역산방법은 데이터획득부를 통해 획득된 데이터들을 이용하여 파형역산을 수행하는 파형역산부가 수행한다.
이하, 도 1 내지 도 3을 참조하여, 본 발명의 제 1 실시예에 따른 순수음향파동모델링 방법에 대하여 설명하도록 한다.
도 1 은 본 발명의 제 1 실시예에 따른 순수음향파동모델링 방법의 순서도이며, 도 2는 도 1에 도시된 이방성 매질에서 S파가 분리된 순수음향파동방정식을 설정하는 순서도이다.
도 3은 본 발명의 제 1 실시예에 따른 순수음향파동모델링 방법에서 셀(cell) 안에 물성정보가 포함된 셀 기반 격자세트(grid set)를 나타내는 도이다.
본 발명의 제 1 실시예는 이방성 매질에서의 파동전파를 모사하기 위한 순수 음향파동 모델링 방법으로써, 도 1에 도시된 바와 같이, 이방성 매질에서의 S파가 분리된 순수음향파동방정식을 표현하는 단계(S110)를 포함한다.
그리고 순수음향파동방정식을 행렬형태로 표현하는 단계(S120), 행렬형태로 표현된 순수음향파동방정식을 유한차분법으로 표현하는 단계를 포함한다.(S130)
우선, 이방성 매질에서 P파와 SV파의 분산관계방정식은 다음과 같이 표현된다.(S111)
[수학식 1]
여기서 는 각주파수;는 P파의 수직속도;는 P파의 수평속도;는 P파의 시간차(moveout)속도;는 sv파의 수직속도; ,는 톰슨(Thomsen)의 이방성변수; , , 는 공간파동수이며, 후술되는 수학식에 동일하게 적용된다.
이방성 매질에서의 S파를 배제한 유사음향파 가정(the pseudoacoustic approximation) 에 적용하기 위해, 수학식 1에서 는 0으로 설정되며, 수학식 1은 다음과 같이 표현된다.
[수학식 2]
수학식 2로 표현된 분산관계방정식은 P파와 SV파 각각의 독립방정식으로 인수분해하여 다음과 같이 표현된다.(S112)
[수학식 3]
[수학식 4]
수학식 3은 P파에 대한 독립방정식이며, 수학식 4는 SV파에 대한 독립방정식이다.
본 발명은 P파에 관련된 요소만 고려하므로, 수학식 3은 1차 테일러 급수전개(the first-order Taylor series expansion)를 사용하여 다음과 같이 표현된다.(S113)
[수학식 5]
또는,
[수학식 6]
수학식 5 및 수학식 6으로 표현된 순수음향파동방정식은 S파가 배제되어 있다.
따라서, 본 발명의 제 1 실시예에 따른 순수음향파동모델링 결과는 S파 인공물(artifacts)을 배제하는 후처리 단계를 필요로 하지 않는다.
[수학식 7]
여기서, P 는 파동장이다.
수학식 7은 암시적 시간적분 기법(the implementation of the implicit time integration methods, Kim and Lim, 2007)를 이용하여 다음과 같이 두 가지 형식의 수학식으로 표현될 수 있다.
[수학식 8]
[수학식 9]
여기서, Q는 가상파동장이다.
수학식 8 및 수학식 9는 셀기반 유한차분법을 이용하여 해결할 수 있으며, 셀기반 유한차분법으로 근사하기 위해, 수학식 8 및 수학식 9는 다음과 같이 행렬형태로 표현된다.(S120)
[수학식 10]
또는,
[수학식 11]
여기서, I 는 단위행렬; f 는 송신원; P 는 파동장; Q 는 가상파동장이다.
여기서, 종래 유한차분법은 시스템 행렬의 대칭성을 파괴할 수 있기 때문에, 수학식 8과 수학식 9의 편미분은 Min et al. (2004)이 제안한 셀기반 유한차분법에 의해 치환된다.
도 3을 참조하면, 물성정보가 셀(cell) 내부에 정의되기 때문에 셀기반 유한차분법법은 자유표면에서 자유경계조건을 다룬다.
그러나, 일반적으로 알려진 셀기반 유한차분법은 1차 또는 2차 미분항을 다루며, 이방성 매질에서 순수음향파동방정식은 4차 미분항을 가지고 있기 때문에, 한 점을 둘러싼 네가지 물성정보의 가중합을 이용하여 이 문제를 해결할 수 있다.
따라서, 순수음향파동방정식의 4차 미분항은 다음과 같이 표현될 수 있다.(S130)
[수학식 12]
이하, 도 4 및 도 5를 참조하여, 본 발명의 제 1 실시예에 따른 순수음향파동모델과 일반적인 유한차분법을 사용한 순수음향파동모델 결과를 비교하도록 한다.
도 4는 본 발명의 제 1 실시예에 따른 순수음향파동모델과 일반적인 유한차분법을 사용한 순수음향파동모델의 상반성을 비교하기 위하여 가정된 속도모델을 나타내는 도이다.
도 5a는 각 송,수신기쌍에서 획득한 탄성파 파동장을 일반적인 유한차분법을 사용한 순수음향파동모델을 사용하여 나타낸 도이며, 도 5b는 각 송,수신기쌍에서 획득한 탄성파 파동장을 본 발명의 일실시예에 따른 순수음향파동모델을 사용하여 나타낸 도이다.
도 4에 도시된 바와 같이, 본 발명의 제 1 실시예에 따른 순수음향파동모델과 일반적인 유한차분법을 사용한 순수음향파동모델을 비교하기 위하여 속도모델을 가정하고, 송신기와 수신기를 다음과 같이 (625 m, 62.5 m)와 (2500 m, 62.5 m)에 번갈아 위치시킨다.
Source (x, z) |
Receiver (x, z) |
|
Case1 |
(625 m, 62.5 m) |
(2500 m, 62.5 m) |
Case2 |
(2500 m, 62.5 m) |
(625 m, 62.5 m) |
획득한 두 개의 탄성파 탐사자료를 비교한 결과는 다음과 같다.
도 5a에 도시된 바와 같이, 일반적인 유한차분법을 사용한 순수음향파동모델인 경우, 탄성파 탐사자료의 상반성이 상당히 깨진 것을 알 수 있다.
반면, 도 5b에 도시된 바와 같이, 본 발명의 일실시예에 따른 순수음향파동모델인 경우, 상대적으로 상반성이 잘 유지되는 것을 확인할 수 있다.
본 발명의 제 1 실시예에 따라 도출된 순수음향파동방정식은 후술되는 본 발명의 제 2 실시예에 따른 전파형 역산방법에 적용할 경우, 역산과정의 시간과 비용을 절감할 수 있는 효과를 얻을 수 있다.
전파형 역산(waveform inversion)이란, 현장에서 실제 측정된 데이터를 이용하여 특정지역의 지하 매질구조에 관한 정보(예컨대, 측정 대상 지역에 대한 속도 모델 또는 밀도 모델)를 유추하는 과정을 말하며, 전술된 본 발명의 제 1 실시예는 이러한 파형 역산에서 해석자가 임의의 지하구조 모델을 설정한 후, 설정된 모델에 대한 이론값을 구하는 모델링 과정이다.
구체적인 설명에 앞서, 본 발명의 제 2 실시예에 따른 전파형 역산방법은 전술된 제 1 실시예에 따라 도출된 순수음향파동방정식과 평면파 가정(Vigh and Starr,2008)을 이용하였으며, 모든 주파수에 대한 동시 다변수역산방법과, 개별 주파수별 순차적 다변수역산기법을 이용하였다.
이하, 도 6 내지 도 9를 참조하여, 이에 대한 역산방법에 대하여 살펴보도록 한다.
도 6은 본 발명의 제 2 실시예에 따른 전파형 역산방법의 순서도이다.
도 6에 도시된 바와 같이, 파선변수(p, ray-parameter)의 범위에 따른 P파 수직속도(Vpz)와, P파 수직속도(Vpz)와 이방성변수(,)의 관계로 정의되는 P파 수평속도( ) 와 P파 시간차속도( )의 민감도를 분석하는 단계를 포함한다.(S210)
그리고 다른 변수를 고정하고 민감도가 가장 큰 P파 수직속도(Vpz)를 먼저 순차 역산하는 단계와(S220), 상기 순차 역산단계에서 도출된 P파 수직속도(Vpz)와 P파 수평속도(Vpx)를 동시 역산하는 단계를 포함한다.(S240)
파선변수(p)는 전술된 수학식 10과 수학식 11에 기재된 순수파동방정식의 행렬항이 다음 수학식 13의 우변에 의해 대체되어 정의된다.
[수학식 13]
도 7 내지 도 9를 참조하여, 여러 개의 파선변수(p) 범위에 따른 P파 수직속도(Vpz), P파 수평속도(Vpx) 및 P파 시간차속도(Vpn)의 민감도를 분석해 보면 다음과 같다.
설명에 앞서, 도면 상에서 붉은색에 가까울수록 민감도가 높다는 것이고, 푸른색에 가까울수록 민감도가 낮다는 것이다.
도 7에 도시된 바와 같이, P파 수직속도(Vpx)의 경우 파선변수(p)의 절대값이 감소할수록 더 높은 민감도를 가진다.
반면, 도 8 및 도 9에 도시된 바와 같이, 이방성변수(,)와 P파 수직속도(Vpz)의 관계로 정의되는 P파 수평속도(Vpx)와 P파 시간차속도(Vpn)의 경우는 파선변수(p)의 절대값이 감소할수록 민감도가 줄어드는 것을 확인할 수 있다.
즉, 이러한 파선변수(p)의 특성을 이용하여, 다른 변수를 고정하고 민감도가 가장 큰 P파 수직속도(Vpz)를 순차 역산하고, 이에 도출된 P파 수직속도(Vpz)와 P파의 수평속도(Vpx)를 동시 역산한다.
우선, 다른 변수를 고정한 채로 민감도가 큰 P파 수직속도(Vpz)를 먼저 순차 역산한다.(S220)
여기서, P파 수평속도(Vpx) 나 P파 시간차속도(Vpn)에 비해 P파 수직속도(Vpz)의 민감도가 상대적으로 크기 때문에 전체 파선변수영역을 이용하여 역산결과가 적절히 수렴할 때까지 반복역산을 수행한다.(S230)
그리고, 도출된 P파 수직속도(Vpz)와 P파 수평속도(Vpx)를 동시 역산한다.(S240)
이때, P파 수직속도(Vpz)의 민감도가 상대적으로 강하게 나타나는 파선변수(p)의 절대값이 작은 영역을 제외하고 역산을 수행한다.
본 발명의 제 2 실시예에 따른 전파형 역산방법의 결과와 기존의 순차적, 동시 역산방법의 결과를 비교하기 위하여 다음과 같이 역산된 모델변수의 오차값을 정의할 수 있다.
[수학식 14]
기존의 순차적, 동시 역산방법과 본 발명의 제 2 실시예에 따른 전파형 역산방법의 결과는 다음과 같으며, 본 발명의 제 2 실시예에 따른 전파형 역산방법의 결과가 상대적으로 적은 오차를 보이는 것을 확인할 수 있다.
Vpz (%) |
Vpx (%) |
ε (%) |
|
Sequential |
0.155 |
0.136 |
19.2 |
Joint |
0.202 |
0.221 |
29.1 |
Developed |
0.125 |
0.115 |
15.4 |
이와 같이, 상술된 본 발명의 제 1 실시예에 따른 순수음향파동모델링 방법은 탄성파 탐사자료의 상반성이 깨지는 것을 방지할 수 있는 효과를 얻을 수 있다.
또한, 기존 방식에 비해 적은 메모리를 이용하여 순수음향파동방정식 모델링이 가능한 효과를 얻을 수 있다.
또한, 상기 효과로 인하여, 모델링의 안정성을 획득할 수 있는 효과를 얻을 수 있다.
상술된 본 발명의 제 2 실시예에 따른 전파형 역산방법은 역산과정의 시간을 단축시켜 역산 효율을 증대시킬 수 있는 효과가 있다.
본 명세서에 개시된 실시예들은 본 발명의 기술적 사상을 한정하기 위한 것이 아니라 설명하기 위한 것으로, 이러한 실시예에 의하여 본 발명의 기술 사상의 범위가 한정되는 것이 아님은 자명하다.
본 발명의 명세서 및 도면에 포함된 기술적 사상의 범위 내에서 당해 기술분야에 있어서의 통상의 지식을 가진 자가 용이하게 유추할 수 있는 변형예와 구체적인 실시예는 모두 본 발명의 권리범위에 포함되는 것으로 해석되어야 할 것이다.
S110 : S파가 분리된 순수음향파동방정식 표현
S120 : 행렬행태로 표현
S130 : 셀기반 유한차분법 표현
S210 : 파선변수 범위에 따른 수직속도와 이방성변수의 민감도 분석
S220 : 민감도가 가장 큰 수직속도 먼저 순차역산
S230 : 역산결과값 적절수렴치 도달
S240 : 수직속도, 수평속도 동시 역산
S120 : 행렬행태로 표현
S130 : 셀기반 유한차분법 표현
S210 : 파선변수 범위에 따른 수직속도와 이방성변수의 민감도 분석
S220 : 민감도가 가장 큰 수직속도 먼저 순차역산
S230 : 역산결과값 적절수렴치 도달
S240 : 수직속도, 수평속도 동시 역산
Claims (12)
- 이방성 매질에서의 파동전파를 모사하기 위한 순수 음향파동 모델링 방법에 있어서,
(a) 상기 이방성 매질에서의 S파가 분리된 순수음향파동방정식을 표현하는 단계;
(b) 상기 순수음향파동방정식을 행렬형태로 표현하는 단계;
(c) 상기 행렬형태로 표현된 순수음향파동방정식을 유한차분법으로 표현하는 단계;를 포함하고,
상기 (a) 단계에서 상기 순수음향파동방정식은 1차 테일러급수 전개(the first-order Taylor series expansion) 및 암시적 시간적분 기법(the implementation of the implicit time integration methods)이 이용되어 표현되는 것을 특징으로 하는,
순수음향파동모델링 방법.
- 제 1 항에 있어서,
상기 (a) 단계는
상기 이방성 매질에서 P파와 SV파의 분산관계방정식을 표현하는 단계;
상기 P파와 SV파 각각의 독립방정식을 표현하는 단계; 및
상기 P파의 독립방정식으로부터 상기 순수음향파동방정식을 표현하는 단계;
를 포함하는,
순수음향파동모델링 방법.
- 제 2 항에 있어서,
상기 P파와 SV파 각각의 독립방정식을 표현하는 단계는
상기 분산관계를 상기 P파와 SV파에 대한 독립방정식으로 인수분해하는 것을 특징으로 하는,
순수음향파동모델링 방법.
- 제 2 항에 있어서,
상기 P파의 독립방정식으로부터 상기 순수음향파동방정식을 표현하는 단계는
상기 P파의 독립방정식으로부터 1차 테일러급수 전개(the first-order Taylor series expansion)를 이용하여 상기 순수음향파동방정식이 표현되는 것을 특징으로 하는,
순수음향파동모델링 방법.
- 제 2 항에 있어서,
상기 P파의 독립방정식으로부터 상기 순수음향파동방정식을 표현하는 단계는
상기 P파의 독립방정식으로부터 암시적 시간적분 기법(the implementation of the implicit time integration methods)를 이용하여 상기 순수음향파동방정식을 표현하는 것을 특징으로 하는,
순수음향파동모델링 방법.
- 제 1 항에 있어서,
상기 (c) 단계는
한 점의 주위를 둘러싼 네 셀의 물성정보를 이용하여 상기 순수음향파동방정식을 유한차분법으로 표현하는 것을 특징으로 하는,
순수음향파동모델링 방법.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020150155007A KR101693311B1 (ko) | 2015-11-05 | 2015-11-05 | 셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020150155007A KR101693311B1 (ko) | 2015-11-05 | 2015-11-05 | 셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법 |
Publications (1)
Publication Number | Publication Date |
---|---|
KR101693311B1 true KR101693311B1 (ko) | 2017-01-06 |
Family
ID=57832414
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020150155007A KR101693311B1 (ko) | 2015-11-05 | 2015-11-05 | 셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법 |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101693311B1 (ko) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120019055A (ko) * | 2010-08-24 | 2012-03-06 | 서울대학교산학협력단 | 탄성 매질에서의 주파수 영역 역시간 구조보정을 이용한 지하구조의 영상화 장치 및 방법 |
US20130060544A1 (en) * | 2010-05-12 | 2013-03-07 | Petrus Maria Bakker | Seismic p-wave modelling in an inhomogeneous transversely isotropic medium with a tilted symmetry axis |
KR101355107B1 (ko) | 2012-06-29 | 2014-01-23 | 서울대학교산학협력단 | 잡음 제거를 통한 지하 매질구조 추정방법 및 그 장치 |
KR101459388B1 (ko) * | 2014-04-18 | 2014-11-07 | 한국해양대학교 산학협력단 | 지하 속도정보 도출 방법 |
-
2015
- 2015-11-05 KR KR1020150155007A patent/KR101693311B1/ko active IP Right Grant
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130060544A1 (en) * | 2010-05-12 | 2013-03-07 | Petrus Maria Bakker | Seismic p-wave modelling in an inhomogeneous transversely isotropic medium with a tilted symmetry axis |
KR20120019055A (ko) * | 2010-08-24 | 2012-03-06 | 서울대학교산학협력단 | 탄성 매질에서의 주파수 영역 역시간 구조보정을 이용한 지하구조의 영상화 장치 및 방법 |
KR101355107B1 (ko) | 2012-06-29 | 2014-01-23 | 서울대학교산학협력단 | 잡음 제거를 통한 지하 매질구조 추정방법 및 그 장치 |
KR101459388B1 (ko) * | 2014-04-18 | 2014-11-07 | 한국해양대학교 산학협력단 | 지하 속도정보 도출 방법 |
Non-Patent Citations (1)
Title |
---|
Ho-Yong Lee, et al. 2D time-domain acoustic-elastic coupled modeling: a cell-based finite-difference method. 1. Geosciences Journal. Vol.13, No.4. December 2009. pp.407-414 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101931488B1 (ko) | 스펙트럼 성형을 이용하는 전 파동장 반전의 수렴 레이트 | |
Hu et al. | A robust and accurate seismic attenuation tomography algorithm | |
Biswas et al. | 2D full-waveform inversion and uncertainty estimation using the reversible jump Hamiltonian Monte Carlo | |
CN104536044B (zh) | 一种地震数据的插值去噪方法以及*** | |
Javaherian et al. | Refraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast | |
KR101693310B1 (ko) | 평면파를 이용한 이방성 매질에서의 전파형 역산방법 | |
CN111596366B (zh) | 一种基于地震信号优化处理的波阻抗反演方法 | |
CN105353405B (zh) | 一种全波形反演方法和*** | |
CN106483559B (zh) | 一种地下速度模型的构建方法 | |
KR20180067650A (ko) | 진폭 보존을 갖는 fwi 모델 도메인 각도 스택들 | |
CN105467442B (zh) | 全局优化的时变稀疏反褶积方法及装置 | |
CN110187384B (zh) | 贝叶斯时移地震差异反演方法及装置 | |
EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
JP6396037B2 (ja) | データ解析装置及び方法 | |
WO2015042386A1 (en) | Eikonal solver for quasi p-waves in anisotropic media | |
Li et al. | A robust approach to time‐to‐depth conversion and interval velocity estimation from time migration in the presence of lateral velocity variations | |
CN111060961B (zh) | 基于多信息约束反演的品质因子确定方法、装置及*** | |
AU2014208222A1 (en) | Methods, systems and devices for generating slowness-frequency projection logs | |
Shao et al. | Simultaneous inversion of Q and reflectivity using dictionary learning | |
AU2019406345B2 (en) | Methods and systems for calibrating depth in a well to seismic data in a subsurface volume of interest | |
KR101693311B1 (ko) | 셀기반 유한차분법을 이용한 이방성 매질에서의 순수음향파동모델링 방법 | |
CN109521470B (zh) | 分析地质构造对地震反演裂缝密度影响的方法 | |
Gisolf et al. | Target-oriented elastic full wave form inversion | |
CN112368606A (zh) | 评估对地震信号的处理印记 | |
CA3034807C (en) | Estimating interval anisotropy parameter for pre-stack depth migration using a least-squares method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20191001 Year of fee payment: 4 |