• 지질공학
• /
• 제14권1호
• /
• pp.21-28
• /
• 2004

지형효과를 포함하는 3차원 전기비저항 역산법은 전기비저항 탐사자료의 해석에 있어서 실제 지하구조에 존재하는 3차원 전도성 이상체의 특성을 파악함으로써 원리적으로 정확한 영상을 획득할 수 있는 방법이라고 본다. 본 연구는 수치모델링에서 유한요소법을 이용하였으며, 역산에서 설정된 블록에 대하여 감도 분석을 통하여 자코비안 계산의 효율성을 극대화하였다. 또한 역산과정에서는 라그랑지 곱수의 값을 변화시키면서 최소자승근 오차의 최소가 되는 최적의 값을 선택하는 방법으로 역산의 분해능을 향상 시켰다. 전기비저항 탐사에서 지형효과의 해석으로는 지형 의 기복을 포함하여 얻어진 순차모델링 자료를 가지고 역산 방법에서 지형의 기복을 무시하여 얻어진 결과와 지형의 기복을 포함하여 얻어진 결과를 비교하여 불규칙한 지형효과에 의한 왜곡된 전기장 반응으로 인한 오해석을 파악하는데 있다.

• 지질공학
• /
• 제16권3호
• /
• pp.255-263
• /
• 2006

지표가 불규칙한 지형에서 경사시추공을 이용한 단극-쌍극자 배열 전기비저항 토모그래피 탐사를 보면, 우리는 시추공안에 있는 전류원과 지표의 지형기복 사이의 최단거리에 위치한 전위전극에서의 극성변화, 지형기복에 따른 거리계수와 포텐셜의 변화 및 모델링에서 경사 시추공에 위치한 전극과 절점의 공간적인 불일치 등을 고려할 수 있다. 요소분할 방법은 시추공이 곡선이거나 경사져 있을 경우와 하나의 요소에 여러 개의 전극이 있을 때 각 전극에 대한 절점 좌표를 지정할 수 있기 때문에 매우 효과적인 방법이라고 본다. 시추공과 지형기복이 동일한 경사를 갖는 경우에서의 역산 결과는 매우 좋은 영상으로 나타났으며 최소자승근의 오차가 안정적으로 수렴되는 경향을 보였다.

• 화약ㆍ발파
• /
• 제28권1호
• /
• pp.27-39
• /
• 2010

측정된 지반진동의 최대입자속도 자료에 대한 통계분석을 통해 환산거리 개념에 기초한 제어발파 설계조건식은 구할 수 있다. 이들 설계조건식들은 안전발파를 위한 다양한 허용기준에 따라 사용할 수 있는 환산거리의 최소값을 정의하는 형태로 되어 있다. 국내에서 널리 사용되는 환산거리에는 자승근 환산거리(SRSD)와 삼승근 환산거리(CRSD)의 두 가지가 있다. 따라서 SRSD와 CRSD의 설계조건식들은 각각 $D/\sqrt{W}{\geq}30m/kg^{1/2}$와 $D/\sqrt[3]{W}{\geq}60m/kg^{1/3}$의 형태가 된다. 제어발파 설계 시에는 이들 조건식들과 이격거리를 알고 있으므로 지반진동에 대해 구조물의 안정을 보장할 수 있는 최대 지발당장약량를 계산할 수 있다. 그러나 SRSD와 CRSD의 최대 지발당장약량은 각각 $W=O(D^2)$와 $W=O(D^3)$의 차원으로 나타난다. 따라서 SRSD에 비해 CRSD의 장약량은 두 회귀식의 유사한 적합도에도 불구하고 두 함수의 교점을 지나면 기하급수적으로 증가하게 된다. 따라서 본 논문에서는 CRSD의 지나치게 많은 장약량으로 인해 발생할 지도 모를 구조물의 피해를 방지하기 위해 CRSD는 어떤 특정한 거리 이내에서만 사용하도록 제한한다. 그 정확한 한계는 SRSD와 CRSD의 장약량 차가 교점 이내에서의 양자 간의 최대차를 초과하기 시작하는 점까지이다.

• 화약ㆍ발파
• /
• 제9권1호
• /
• pp.3-13
• /
• 1991

발파에 의한 지반진동의 크기는 화약류의 종류에 따른 화약의 특성, 장약량, 기폭방법, 전새의 상태와 화약의 장전밀도, 자유면의 수, 폭원과 측간의 거리 및 지질조건 등에 따라 다르지만 지질 및 발파조건이 동일한 경우 특히 측점으로부터 발파지점 까지의 거리와 지발당 최대장약량 (W)간에 깊은 함수관계가 있음이 밝혀졌다. 즉 발파진동식은 $V=K{\cdot}(\frac{D}{W^b})^n{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (1) 여기서 V ; 진동속도, cm /sec D ; 폭원으로부터의 거리, m W ; 지발 장약량, kg K ; 발파진동 상수 b ; 장약지수 R ; 감쇠지수 이 발파진동식에서 b=1/2인 경우 즉 $D{\;}/{\;}\sqrt{W}$를 자승근 환산거리(Root scaled distance), $b=\frac{1}{3}$인 경우 즉 $D{\;}/{\;}\sqrt[3]{W}$를 입방근환산거리(Cube root scaled distance)라 한다. 이 장약 및 감쇠지수와 발파진동 상수를 구하기 위하여 임의거리와 장약량에 대한 진동치를 측정, 중회귀분석(Multiple regressional analysis)에 의해 일반식을 유도하고 Root scaling과 Cube root scaling에 대한 회귀선(regression line)을 구하여 회귀선에 대한 적합도가 높은 쪽을 택하여 비교, 검토하였다. 위 (1)식의 양변에 log를 취하여 linear form(직선형)으로 바꾸어 쓰면 (2)式과 같다. log V=A+BlogD+ClogW ----- (2) 여기서, A=log K B=-n C=bn (2)식은 다시 (3)식으로 표시할 수 있다. $Yi=A+BXi_{1}+CXi_{2}+{\varepsilon}i{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$(3) 여기서, $Xi_{1},{\;}Xi_{2} ;(두 독립변수 logD, logW의 i번째 측정치. Yi ; ($Xi_1,{\;}Xi_2$)에 대한 logV의 측정치 ${\varepsilon}i$ ; error term 이다. (3)식에서 n개의 자료를 (2)식의 회귀평면으로 대표시키기 위해서는 $S={\sum}^n_{i=1}\{Yi-(A+BXi_{1}+CXi_{2})\}\^2$을 최소로하는 A, B, C 값을 구하면 된다. 이 방법을 최소자승법이 라 하며 S를 최소로 하는 A, B, C의 값은 (4)식으로 표시한다. $\frac{{\partial}S}{{\partial}A}=0,{\;}\frac{{\partial}S}{{\partial}B}=0,{\;}\frac{{\partial}S}{{\partial}C}=0{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (4) 위식을 Matrix form으로 간단히 나타내면 식(5)와 같다. [equation omitted] (5) 자료가 많아 계산과정이 복잡해져서 본실험의 정자료들은 전산기를 사용하여 처리하였다. root scaling과 Cube root scaling의 경우 각각 $logV=A+B(logD-\frac{1}{2}W){\;}logV=A+B(logD-\frac{1}{3}W){\;}\}{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (6) 으로 (2)식의 특별한 형태이며 log-log 좌표에서 직선으로 표시되고 이때 A는 절편, B는 기울기를 나타낸다. $\bullet$ 측정치의 검토 본 자료의 특성을 비교, 검토하기 위하여 지금까지 발표된 국내의 몇몇 자료를 보면 다음과 같다. 물론, 장약량, 폭원으로 부터의 거리등이 상이하지만 대체적인 경향성을 추정하는데 참고할수 있을 것이다. 금반 총실측자료는 총 88개이지만 환산거리(5.D)와 진동속도의 크기와의 관계에서 차이를 보이고 있어 편선상 폭원과 측점지점간의 거리에 따라 l00m말만인 A지역과 l00m이상인B지역으로 구분하였다. 한편 A지역의 자료 56개중, 상하로 편차가 큰 19개를 제외한 37개자료와 B지역의 29개중 2개를 낙외한 27개(88개 자료중 거리표시가 안된 12월 1일의 자료3개는 원래부터 제외)의 자료를 computer로 처리하여 얻은 발파진동식은 다음과 같다. $V=41(D{\;}/{\;}\sqrt[3]{W})^{-1.41}{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (7) (-100m)(R=0.69) $V=124(D{\;}/{\;}\sqrt[3]{W})^{-1.66){\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (8) (+100m)(R=0.782) 식(7) 및 (8)에서 R은 구한 직선식의 적합도를 나타내는 상관계수로 R=1인때는 모든 측정자료가 하나의 직선상에 표시됨을 의미하며 그 값이 낮을수록 자료가 분산됨을 뜻한다. 본 보고에서는 상관계수가 자승근거리때 보다는 입방근일때가 더 높기 때문에 발파진동식을 입방근($D{\;}/{\;}\sqrt[3]{W}$)으로 표시하였다. 특히 A지역에서는 R=0.69인데 비하여 폭원과 측점지점간의 거리가 l00m 이상으로 A지역보다 멀리 떨어진 B지역에서는 R=0.782로 비교적 높은 값을 보이는 것은 진동성분중 고주파성분의 상당량이 감쇠를 당하기 때문으로 생각된다.