KR102110692B1 - 심부 지질, 황화물 광석체 자원 탐지 방법 - Google Patents

심부 지질, 황화물 광석체 자원 탐지 방법 Download PDF

Info

Publication number
KR102110692B1
KR102110692B1 KR1020180155095A KR20180155095A KR102110692B1 KR 102110692 B1 KR102110692 B1 KR 102110692B1 KR 1020180155095 A KR1020180155095 A KR 1020180155095A KR 20180155095 A KR20180155095 A KR 20180155095A KR 102110692 B1 KR102110692 B1 KR 102110692B1
Authority
KR
South Korea
Prior art keywords
deep
firing
ground
waveform
mineral resources
Prior art date
Application number
KR1020180155095A
Other languages
English (en)
Other versions
KR20180136918A (ko
Inventor
다 레이
칭윈 디
량용 양
뤄 왕
창민 푸
윈시앙 리우
Original Assignee
인스티튜트 오브 지오로지 앤드 지오피직스 차이니스 아카데미 오브 사이언스
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 인스티튜트 오브 지오로지 앤드 지오피직스 차이니스 아카데미 오브 사이언스 filed Critical 인스티튜트 오브 지오로지 앤드 지오피직스 차이니스 아카데미 오브 사이언스
Publication of KR20180136918A publication Critical patent/KR20180136918A/ko
Application granted granted Critical
Publication of KR102110692B1 publication Critical patent/KR102110692B1/ko

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/088Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices operating with electric fields
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/081Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/087Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the earth magnetic field being modified by the objects or geological structures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38Processing data, e.g. for analysis, for interpretation, for correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Abstract

본 발명은 심부 지질, 황화물 광물체 자원의 탐지 방법을 개시하고, 발사 장치는 땅 아래로 의사난수의 발사 파형을 발송하고, 지하 매체가 다시 지면으로 전파되며, 지면에서 수신기를 관측하여 전자기장 응답 신호를 수집하고, 수신기에서 관측된 전자기장 응답 신호와 발사 파형의 상관관계, 분별 추출된 의사난수 발사 파형을 디콘볼루션하여 유도된 대지 임펄스 응답 신호 분석하여 겉보기 비저항과 분극률 위상을 계산하고, 데이터 처리로 심부 자원의 높은 해상도의 분극률과 비저항을 반영하는 전기적 미세 구조를 획득하는 것을 포함한다. 본 발명에서 개시하는 탐지 방법은, 다중 채널 전법기의 가상의 임의 부호화된 전류 신호 발사, 관측된 전자기장 응답과 발사 전류 신호의 디콘볼루션을 사용하여 간섭을 제거하고, 겉보기 비저항과 유도분극 위상을 계산하여 지질, 황화물 광물체 자원의 탐사 탐지를 달성하고, 심부 황화물 광석체의 탐사 효과를 크게 향상시켰으며, 실제 넓은 응용에 유리하다.

Description

심부 지질, 황화물 광석체 자원 탐지 방법{Deep fat, sulfide ore body detection method}
본 발명은 탐지 장치에 관한 것으로서, 특히 심부 지질, 황화물 광석체 자원 탐지 방법에 관한 것이다.
지질 광석체 자원의 탐사 탐지에서 사용되는 방법은 일반적으로 전자기법과 유도분극법을 포함한다. 전지기법은 비록 1000m이상 깊이의 지질과 광석체의 탐사가 가능하나, 황화물 광석체와 관련된 분극률 정보를 탐지하는 것이 불가능하다. 따라서 함침 상태의 금속 광석체의 황화물의 탐사를 위해서는 유도분극법을 채용하여 탐사 탐지하였다. 현재, 기존의 분극률 탐지기의 탐사 깊이는 1000m 이하로 얕고, 심부 지질, 황화물 광석체 자원 탐사 탐지를 할 수 없어, 심부 황화물 광석체의 탐사 효과에 큰 영향을 주었고, 실제 사용에 제한이 있었다. 따라서, 본 영역의 기술자는 상술한 종래 기술에 존재하는 단점을 해결한 심부 지질, 황화물 광석체 자원 탐지 방법을 제공하기 위해 노력한다.
종래 기술의 상기 결점의 관점에서, 본 발명에 의해 해결되어야 할 기술적 과제는 기존의 분극률의 탐사 탐지 방법이 지질, 황화물 광물체 자원 탐지 시에 탐사 심도가 1000m 이하로 얕고, 심부 지질, 황화물 광물체 자원 탐지를 수행할 수 없으며, 심부 황화물 광물체의 탐사 효과에 크게 영향을 주었고, 실제 적용에 제한이 있었으며, 보급 및 적용에 도움이 되지 않는 확대 적용이 잘 되지 않는 것이다.
상기 목적을 달성하기 위하여, 본 발명이 제공하는 심부 지질, 황화물 광물체 자원의 탐지 방법은
단계 1. 발사 장치는 땅 아래로 의사난수의 발사 파형을 발송하고, 지하 매체가 다시 지면으로 전파되며, 지면에서 수신기를 관측하여 전자기장 응답 신호를 수집하는 단계
단계 2. 수신기에서 관측된 전자기장 응답 신호와 발사 파형의 상관관계, 분별 추출된 의사난수 발사 파형을 디콘볼루션하여 유도된 대지 임펄스 응답 신호 분석을 통하여 겉보기 비저항과 분극률 위상을 계산하는 단계
단계 3. 데이터 처리 및 역전 해석하고, 심부 자원의 높은 해상도의 분극률과 비저항을 반영하는 전기적 미세 구조를 획득하는 단계를 포함한다.
바람직하게는, 상기 단계 1에서, 상기 관측하여 수집한 전자기장 응답은 다중 수신기에 의해 다중 전계 배열식 분포로 설치 구성되고, 각 수신기당 3개의 전기장 데이터를 관측하며, 수신 전기 쌍극자 거리(MN)는 수십미터이고, 발사 전극 거리(AB)는 수신 전기 쌍극자 거리(MN)의 2~6배이다.
바람직하게는, 상기 다중 전계 배열식 분포 수신 시스템은 설치 후 고정 불변하며, 발사 장치와 발사 전극 거리(AB)위치 이동을 통해 전류 발사의 동시에 모든 수신기 측정, 모든 전도(전기 채널) 전기장 신호의 순서 측량을 진행한다.
바람직하게는, 상기 단계 2에서, 상기 대지 임펄스 응답 신호는 수신기와 발신기의 거리(수신, 송신 거리)가 증가할 때 관측할 수 있는 비교적 약한 유효 신호이다.
바람직하게는, 상기 단계 2에서, 상기 분극률 위상은 다중 채널 전법기를 이용하여 m서열 의사난수 파형을 발사하고 측정된 전기장 분량 데이터를 디콘볼루션하여, 대지 임펄스 응답 계산 시 주파수 영역 유도분극 위상을 얻는다.
바람직하게는, 상기 단계 3에서, 상기 데이터 처리는 전기장 신호 시간 순 측량과 대지 응답 파형과 유도분극 위상 매개 변수에 근거하여 데이터 처리를 수행한다.
본 발명의 바람직한 실시예로, 상기 단계 1에서 전자기장 응답 신호를 관측 수집하는 것은 쌍극자-쌍극자 관측 장치를 이용한다.
본 발명의 바람직한 실시예로, 상기 단계 2에서 상기 분극률 위상은 다음 식에 의해 계산된다.
Figure 112018121794502-pat00001
식에서,
Figure 112018121794502-pat00002
는 시스템 스펙트럼이다.
이상의 방안을 채용하여, 본 발명에 의해 개시된 심부 지질, 황화물 광석체 자원 탐지 방법은 다음과 같은 장점이 있다.
본 발명에서 개시된 탐지 방법은 다중 채널 고밀도 주파수 영역 유도분극 측량 방법을 채용하여, 종래의 고밀도 비저항법 /분극률법의 탐측 심도가 얕은 문제를 보완하고, 다중 채널 전법기의 가상의 임의 부호화된 전류 신호 발사, 관측된 전자기장 응답과 발사 전류 신호의 디콘볼루션을 사용하여 간섭을 제거하고, 큰 수신 발신 거리에서 탐측 심도가 클 때, 신뢰할 수 있는 대지 임펄스 응답을 얻으며, 겉보기 비저항과 유도분극 위상을 계산하여 다양한 깊이로 묻힌 깊은 지질학적 표적의 비저항과 분극률 매개 변수를 검출한다. 고밀도 유도분극법 탐사 심도를 향상시키고, 심부 자원 분극률을 반영한 정보를 얻을 수 있다. 분극률과 관련된 심부 자원을 직접 정밀하게 측량하는 것이 실현되며, 심도 1000m이상의 심부 지질, 황화물 광석체 자원 탐사 탐지를 실현했고, 심부 황화물 광석체의 탐사 효과를 크게 향상시켰으며, 실제 넓은 응용에 유리하다.
이하에서는 본 발명의 목적, 특징 및 효과를 완전히 이해하기 위하여 도면 및 구체 실시예과 결부하여, 본 발명의 구상, 구체적 기술 방안과 기술적 효과를 더 자세히 설명한다
도 1은 본 발명에서 제안하는 심부 지질, 황화물 광석체 자원 탐지 방법의 작동 원리를 나타낸 도면이다.
도 2는 본 발명에서 제안하는 탐지 방법의 상이한 오프셋(R) 수신 발사 파형을 나타낸 도면이다.
도 3은 본 발명에서 제안하는 탐지 방법의 상이한 오프셋(R)의 대지 임펄스 응답 곡선을 나타낸 도면이다.
도 4는 본 발명에서 제안하는 탐지 방법의 2차원 반전기하학(도치, 전도, 역전) 저항률 단면을 나타낸 도면이다.
도 5는 본 발명에서 제안하는 탐지 방법의 2차원 반전기하학(도치, 전도, 역전) 분극률 위상 단면을 나타낸 도면이다.
이하 본 발명의 바람직한 실시예를 참고하여, 그 기술 내용을 보다 명확하고 편리하게 이해할 수 있다. 본 발명은 많은 상이한 형식의 실시예로 구현될 수 있으며, 본 발명의 보호 범위는 기재된 실시예에 한정되지 않는다.
실시예:
도 1에 나타난 바와 같이, 본 실시예의 심부 지질, 황화물 광석체 자원 탐지 방법의 작동 원리는 유한 길이의 접지 도선 전류원을 통해 땅 아래로 가상의 임의(의사난수의) 부호화된 전류 신호를 보낸다. 발사원의 축 방향으로 동시에 전자기장 응답과 발사전류 기록을 관측한 후에, 디콘볼루션(deconvolution)을 통하여 대지 임펄스 응답을 얻는다. 그리고 겉보기 비저항(apparent resistivity)을 계산하고, 이로써 다르게 묻힌 깊은 지질 표적체 탐측 목적을 달성한다.
다채널 전기 계기 시스템에서 보통 대지를 선형 시불변 시스템(Linear and Time-invariant System)으로 보고, 수신된 응답 전압 v(t) 은 발사전류와 대지 임펄스 응답 회선(Convolution)으로 표시될 수 있다.
Figure 112018121794502-pat00003
식 중
Figure 112018121794502-pat00004
는 부호전류(code current), 발사시스템 및 수신, 발신거리와 연관된 시스템 응답으로, 발사 설비 전기회로, 접지 전극 및 접지 전극에 연결된 케이블 등의 응답을 포함한다.
Figure 112018121794502-pat00005
는 지질 표적체의 대지 임펄스 응답에서 나오는 것이다.
Figure 112018121794502-pat00006
는 원전류를 기록하는 방식 또는 발사원 부근에서 전기장 신호를 채집하는 방식을 통하여 얻은 것이다.
Figure 112018121794502-pat00007
은 노이즈이다.
디콘볼루션을 통해 노이즈를 제거하고 대지 임펄스 응답을 구할 수 있다.
대지 임펄스 응답은 발사기와 수신기 사이의 대지 비저항(저항률) 정보를 포함한다.
균일한 대지에서, 대지 임펄스 응답의 피크 시각
Figure 112018121794502-pat00008
은 다음 식과 같다.
Figure 112018121794502-pat00009
식에서, μ은 매체의 자기투과율(magnetic permeability)이고, r은 수신 발신 오프셋(m)이며, ρ는 대지 비저항(Ω·m)이다.
대지 임펄스 응답의 피크 시각을 이용하여 겉보기 비저항
Figure 112018121794502-pat00010
로 변환한다.
Figure 112018121794502-pat00011
식(a-3)기록 포인트는 송수신 오프셋의 중점이고, 전체 섹션에서 발사점의 지속적인 이동과 전원공급을 통하여, 서로 다른 오프셋의 겉보기 비저항과 기록 위치를 얻고, 전체 섹션의 겉보기 저항의 가단면도(pseudo)를 얻는데, 단면 지하의 서로 다른 심도(깊이)의 지질 정보를 반영하였다.
본 발명의 위상 유도분극법(induced polarization method) 매개변수의 대지 임펄스 응답은, 의사난수기에서 부호화된 발사 파형과 대지 임펄스 응답 파형의 상관성을 계산한 것이고, 식 (a-1)에서 시간 영역을 주파수 영역으로 변환하여 얻은 것이다.
Figure 112018121794502-pat00012
임펄스 응답
Figure 112018121794502-pat00013
의 푸리에 변환은 시스템의 스펙트럼
Figure 112018121794502-pat00014
이고, 위상(phase)을 얻을 수 있다.
Figure 112018121794502-pat00015
계산 중 수신된 전기장 신호와 발사된 전류에 대하여 디콘볼루션하고, 양자의(2개의) 스펙트럼 곱의 역푸리에 변환을 통하여 시간 영역의 대지 임펄스 응답 곡선을 얻는다. 수신된 전기장 신호와 발사된 전류 양자의 스펙트럼 곱의 과정에서 제로 위상 필터링을 이용하여 노이즈 간섭을 제거한다.
상기 공식으로부터, 실측 전기장 신호와 발사 파형 디콘볼루션 및 제로 위상 필터링을 채용 방법을 얻어낼 수 있고, 50Hz와 기타 간섭을 필터링 할 수 있으며, 고품질의 겉보기 비저항과 유도분극 위상을 얻을 수 있다.
상기 계산된 데이터에 대한 품질 평가를 진행하고, 대지 임펄스 응답
Figure 112018121794502-pat00016
곡선 피크 시각으로부터 겉보기 비저항을 구하였다. 임펄스 응답
Figure 112018121794502-pat00017
의 스펙트럼으로 위상을 얻을 수 있다. 기존 관측의 30주기 관측으로부터 겉보기 비저항과 유도분극 위상을 단계적으로 얻을 수 있다. 이로부터 중복 중첩된 평균 겉보기 비저항과 유도분극 위상을 계산해낼 수 있고, 관측 품질을 향상시킬 수 있다. 평균 제곱 상대 오차 M 으로 겉보기 비저항과 유도분극 위상 관측 품질을 판단한다.
겉보기 비저항 또는 유도분극 위상은 평균 제곱의 상대 오차를 이용하여 공식(a-6)에 따라 계산한다.
Figure 112018121794502-pat00018
식에서,
Figure 112018121794502-pat00019
는 제i차 중첩 계산된 겉보기 비저항 또는 유도분극(induced polarization method) 위상(phase) 상대오차이다.
Figure 112018121794502-pat00020
Figure 112018121794502-pat00021
는 제i차 중첩 계산된 겉보기 비저항 또는 유도분극(induced polarization method) 위상(phase)이다.
Figure 112018121794502-pat00022
1에서 제i차 중첩된 평균 겉보기 비저항 또는 유도분극 위상이다.
n 은 수집 주기수이다.
구체 적용 사례는 아래와 같다.
허베이 성 일부 지역의 다채널 전기 계측 테스트에서, 수신기가 측정선의 각 측정점에서 전기장 신호를 기록할 때, 1개의 수신기는 전기 쌍극자 발사 중심에서 동시에 전기장 신호를 관측하고, 데이터 처리 시 대응되는 동일 시간의 발사원 중심에서 수신된 전기장 신호와 수신점의 전기장 신호를 디콘볼루션하여 대지 임펄스 응답을 구한다.
도 2는 본 발명에서 제안하는 탐지 방법의 상이한 오프셋(R) 수신 발사 파형을 나타낸 도면으로, 600 m 지점이 발사원일 때, 발사원에서 상이한 오프셋R(60 m ~4500 m)떨어진 곳에서 수신된 전기장 신호이고, R=0 m일 때 수신된 파형은 전기 쌍극자의 중심에서 관측되는 전기장 신호이다. 그림의 하단은 발사된 의사난수의 전압 파형이다.
오프셋 R 60 m~4500 m,수신된 전기장 신호는 강했다가 점점 약해진다.
오프셋R 60 m~540 m의 전기장 신호는 기본적으로 발사된 의사난수의 신호와 일치한다. 다만, 1900 m~3860 m 지점의 전기장 신호는 기본적으로 50 Hz전자기 간섭 신호이고, 4500 m 지점의 전기장 신호는 비교적 강한 펄스 간섭을 나타낸다.
도 3은 본 발명에서 제안하는 탐지 방법의 상이한 오프셋(R)의 대지 임펄스 응답 곡선을 나타낸 도면으로, 수신된 전기장 신호와 발사 전류의 디콘볼루션을 통해, 50Hz 간섭과 펄스 간섭을 효과적으로 제거할 수 있다. 특히 오프셋이 1900m 보다 클 때도 비교적 좋은 대지 임펄스 응답이 얻어진다. 겉보기 비저항은 응답 곡선의 피크 시각으로부터 계산될 수 있다.
도 4본 발명에서 제안하는 탐지 방법의 2차원 반전기하학(도치, 전도, 역전) 비저항 단면을 나타낸 도면에 도시된 바와 같이, 각 대지 임펄스 응답의 피크값에서 전환된 겉보기 비저항으로 겉보기 비저항 유사 단면도를 제작할 수 있다. 겉보기 비저항은 2층의 지전류 단면으로 반영된다. 얕은 부분의 낮은 저항은 풍화된 시생대 대층 황토요암조, 깊은 부분의 높은 저항은 화강암의 반영이다. 당 지역의 지질 정황과 일치하였다.
도 4 본 발명에서 제안하는 탐지 방법의 2차원 도치 비저항 단면도에 도시된 바와 같이, 획득된 겉보기 비저항과 분극률 위상 매개변수로 미국 Zonge 회사의 Ts2dip.exe 소프트웨어를 채용하여 2차원 비저항 역전(도치)을 수행하였다. 얕은 부분 두께 약50 m의 전기층(Electrical layer)(┯·m이하)는 표층 제4계 지층에 대응한다. 그 하면 1500 ┯·m이하의 고 저항 전기층은 모두 변경된 화강암을 반영한다.
본 발명에서 제안하는 탐지 방법의 2차원 반전기하학(도치, 전도, 역전) 분극률 위상 단면을 나타낸 도면 도 5 에 도시된 바와 같이, 분극률 위상 도치 단면도의 높은 위상값은 광석체를 품은 위치와 기본적으로 일치한다.
본 실측 자료로부터 대지 임펄스 응답 곡선 계산 중, 제로 위상 필터링 및 실측된 전기장 신호와 발사 파형의 디콘볼루션, 그리고 다차 중첩 등 방법을 채용함으로써, 개선된 대지 임펄트 응답과 겉보기 저항 및 분극 위상을 획득할 수 있다.
본 발명에서 제안하는 심부 지질, 황화물 광석체 자원 탐지 방법은 고성능 발사기를 사용하여, 가상의 임의(의사난수의) 부호화된 전류를 발사하고, 전체 파형 데이터를 수집하며, 그것은 석유 및 가스 탐사의 지진 기술을 참고하여, 진열식 수신 방식을 사용하고, 다차 발사를 채용하고, 진열식 다중 채널은 다차 피복된 전체 파장(wave field)정보를 수신하며, 데이터를 지진 편이(전위) 기술과 유사하게 진행하여 지하 지질체 이미징을 실현한다. 해상도와 효과적으로 탐사 깊이가 효과적으로 향상되고, 동등한 발사원 강도 조건에서, 탐지 정확도와 깊이를 크게 향상시키며, 탐사 깊이가 4000m에 이른다. 육지와 해상의 심층 자원 탐사에 사용될 수 있다.
이상 본 발명의 바람직한 구체적 실시예를 설명하였다. 본 영역의 보통의 기술자가 창조적인 노력 없이 본 발명의 개념에 따라 많은 수정 및 변경을 할 수 있음을 이해해야 한다. 따라서, 일반 본 기술영역의 기술자가 본 발명의 개념에 따라 종래 기술의 기초상에서 논리분석, 추리 또는 일정 정도의 실험을 통해 얻을 수 있는 기술방안은 모두 권리 범위에 의해 결정되는 보호범위 내에 속한다.
1: 제1 수신기, 2: 제2 수신기, 3: 발사 시스템, 4: 광석체

Claims (6)

  1. 단계 1. 발사 장치는 접지 도선 전류원을 통해 땅 아래로 의사난수의 발사 파형을 발송하면서 발사원의 축 방향으로 전자기장 응답을 관측하고, 지하 매체로부터 다시 지면으로 전파되면, 지면에서 수신기를 관측하여 전자기장 응답 신호를 수집하는 단계
    단계 2. 수신기에서 관측된 전자기장 응답 신호와 발사 파형의 상관관계, 분별 추출된 의사난수 발사 파형을 디콘볼루션하여 유도된 대지 임펄스 응답 신호 분석을 통하여 겉보기 비저항과 분극률(polarizability) 위상을 계산하는 단계
    단계 3. 데이터 처리 및 역전 해석하고, 심부 자원의 높은 해상도의 분극률과 비저항을 반영하는 전기적 미세 구조를 획득하는 단계를 포함하는 것을 특징으로 하는 심부 지질, 황화물 광물체 자원의 탐지 방법.
  2. 제1항에 있어서, 상기 단계 1은,
    상기 관측하여 수집한 전자기장 응답은 다중 수신기에 의해 다중 전계 배열식 분포로 설치 구성되고, 각 수신기당 3개의 전기장 데이터를 관측하며, 발사 전극 거리는 수신 전기 쌍극자 거리의 2~6배이고,
    상기 다중 전계 배열식 분포 수신 시스템은 설치 후 고정 불변하며, 발사 장치와 발사 전극 거리 위치 이동을 통해 전류 발사의 동시에 모든 수신기 측정, 모든 전도 전기장 신호의 순서 측량을 진행하는 것을 특징으로 하는 심부 지질, 황화물 광물체 자원의 탐지 방법.
  3. 제1항에 있어서, 상기 단계 2는,
    상기 대지 임펄스 응답 신호는 수신기와 발신기의 거리가 증가할 때 관측할 수 있는 비교적 약한 유효 신호이고,
    상기 분극률 위상은 다중 채널 전법기를 이용하여 m서열 의사난수 파형을 발사하고 측정된 전기장 분량 데이터를 디콘볼루션하여, 대지 임펄스 응답 계산 시 주파수 영역 유도분극 위상을 얻는 것을 특징으로 하는 심부 지질, 황화물 광물체 자원의 탐지 방법.
  4. 제1항에 있어서, 상기 단계 3은,
    상기 데이터 처리는 전기장 신호 시간 순 측량과 대지 응답 파형과 유도분극 위상 매개 변수에 근거하여 데이터 처리를 수행하는 것을 특징으로 하는 심부 지질, 황화물 광물체 자원의 탐지 방법.
  5. 제1항에 있어서, 상기 단계 1은,
    상기 전자기장 응답 신호를 관측 수집하는 것은 쌍극자-쌍극자 관측 장치를 이용하는 것을 특징으로 하는 심부 지질, 황화물 광물체 자원의 탐지 방법.
  6. 제3항에 있어서, 상기 단계 2는,
    상기 분극률 위상은 다음 식에 의해 계산되고, 식에서,
    Figure 112018121794502-pat00023
    는 시스템 스펙트럼인 것을 특징으로 하는 심부 지질, 황화물 광물체 자원의 탐지 방법.
    Figure 112018121794502-pat00024
KR1020180155095A 2018-09-06 2018-12-05 심부 지질, 황화물 광석체 자원 탐지 방법 KR102110692B1 (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201811037458.9A CN109085653A (zh) 2018-09-06 2018-09-06 一种深部地质、硫化物矿体资源的探测方法
CN2018110374589 2018-09-06

Publications (2)

Publication Number Publication Date
KR20180136918A KR20180136918A (ko) 2018-12-26
KR102110692B1 true KR102110692B1 (ko) 2020-05-14

Family

ID=64840824

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020180155095A KR102110692B1 (ko) 2018-09-06 2018-12-05 심부 지질, 황화물 광석체 자원 탐지 방법

Country Status (2)

Country Link
KR (1) KR102110692B1 (ko)
CN (1) CN109085653A (ko)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109991666B (zh) * 2019-04-22 2020-10-23 中国冶金地质总局第二地质勘查院 一种矿区地质勘查方法
CN112083505B (zh) * 2019-12-04 2024-02-23 安徽省勘查技术院(安徽省地质矿产勘查局能源勘查中心) 一种基于激电发射系统的综合物理测井方法及系统
CN111398360B (zh) * 2020-04-21 2021-06-04 山东大学 基于ERT和IP圈定包气带中LNAPLs的污染源区探测方法及系统
CN111538095B (zh) * 2020-06-05 2023-03-28 安徽理工大学 一种基于地面单钻孔内反射无线电波的地质勘探方法
CN112485836B (zh) * 2020-11-11 2023-06-20 中国科学院空天信息创新研究院 基于伪随机编码阈值优化的相关辨识数据预处理方法
CN112882105B (zh) * 2021-01-13 2022-09-13 廊坊市立福地质勘查技术有限公司 一种接收矩形波周期信号的自同步方法
CN113156520A (zh) * 2021-04-25 2021-07-23 中国科学院地质与地球物理研究所 一种可控源音频大地电磁测深的非平面波识别方法和系统
CN113866837B (zh) * 2021-09-18 2023-12-19 吉林大学 电性源核磁共振与激发极化联用装置以及探测方法
CN114415245A (zh) * 2022-01-21 2022-04-29 山东省地质矿产勘查开发局第六地质大队(山东省第六地质矿产勘查院) 一种适用于中深部蚀变岩型金矿找矿方法
CN115267920B (zh) * 2022-09-29 2023-02-07 自然资源实物地质资料中心 激发极化参数检测方法、装置、电子设备和存储介质
CN116931102B (zh) * 2023-09-18 2023-12-19 山东大学 一种多水平错频激励充电探测方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101923166B1 (ko) 2018-10-05 2018-11-28 한국지질자원연구원 지진배경잡음을 이용하여 추출된 지진신호의 진폭 크기 보정방법
KR101914657B1 (ko) 2018-04-26 2019-01-07 한국지질자원연구원 지진배경잡음을 이용한 지진신호의 위상 및 진폭정보 산출방법

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0121719D0 (en) * 2001-09-07 2001-10-31 Univ Edinburgh Method for detection fo subsurface resistivity contrasts
KR100831932B1 (ko) * 2006-09-12 2008-05-23 한국지질자원연구원 오일러 해를 이용한 지하공동의 3차원 중력역산 방법 및이를 이용한 3차원 영상화 방법
CN105929455B (zh) * 2016-04-18 2016-12-14 中国科学院地质与地球物理研究所 一种多通道瞬变电磁法三维探测的装置和方法
CN107589461B (zh) * 2017-09-04 2019-11-05 长江大学 一种基于双向编码的尖脉冲时域电磁深层探测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101914657B1 (ko) 2018-04-26 2019-01-07 한국지질자원연구원 지진배경잡음을 이용한 지진신호의 위상 및 진폭정보 산출방법
KR101923166B1 (ko) 2018-10-05 2018-11-28 한국지질자원연구원 지진배경잡음을 이용하여 추출된 지진신호의 진폭 크기 보정방법

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
S. Kwak and S.G. Song. Azimuthal Variation of Empirical Green’s Functions Extracted by Seismic Interferometry in the Southern Korean Peninsula. 79th EAGE Conference and Exhibition 2017. 12-15 June 201

Also Published As

Publication number Publication date
KR20180136918A (ko) 2018-12-26
CN109085653A (zh) 2018-12-25

Similar Documents

Publication Publication Date Title
KR102110692B1 (ko) 심부 지질, 황화물 광석체 자원 탐지 방법
CN110208866B (zh) 地井阵列式光纤时频电磁数据采集装置及其数据采集方法
CN105044792B (zh) 地-井时频电磁勘探数据采集装置及方法
US6696839B2 (en) Electromagnetic methods and apparatus for determining the content of subterranean reservoirs
CA2707784C (en) Method for estimating and removing air wave response in marine electromagnetic surveying
CN110095809A (zh) 井中光纤时频电磁和四分量地震数据采集装置及方法
CN111580170A (zh) 基于极低频电磁源的时频电磁数据采集装置及方法
CN112415615A (zh) 基于分布式光纤传感的时频电磁压裂监测系统及监测方法
CN104656157A (zh) 一种识别页岩气甜点区的方法及装置
EP2313795B1 (en) Method for attenuating air wave response in marine transient electromagnetic surveying
CN105022097A (zh) 一种土质边坡滑动面综合预报方法
CN110187394A (zh) 双场源电磁测深法获取地层电阻率各向异性的方法及装置
CN113703058A (zh) 一种利用视电导率和相对介电常数探测地下障碍物的方法
CN212255732U (zh) 基于极低频电磁源的时频电磁数据采集装置
CN210038187U (zh) 地井阵列式光纤时频电磁数据采集装置
Kang et al. Forward modeling and analyzing for frequency domain semi-airborne EM method
CN209911570U (zh) 井中光纤时频电磁和四分量地震数据采集装置
Rechtien et al. A high-frequency sparker source for the borehole environment
CN213398937U (zh) 基于分布式光纤传感的时频电磁压裂监测系统
CN116430464A (zh) 提高时频电磁法频率域资料处理质量的方法
Thurston et al. Next Generation Geophysical Assessment System
CN113376704B (zh) 一种基于电发射-磁接收的井间电磁探测系统及方法
Faize et al. The use of ground penetrating radar for the detection and study of a buried marble and in situ location of possible cracks
CN109188538B (zh) 一种塔基下采空区电磁检测方法
CN117075208A (zh) 一种基于阻抗相位差分的地下矿物识别方法及装置

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right