BRPI0614372A2 - medições de resistividade cc de múltiplos transientes - Google Patents
medições de resistividade cc de múltiplos transientes Download PDFInfo
- Publication number
- BRPI0614372A2 BRPI0614372A2 BRPI0614372-5A BRPI0614372A BRPI0614372A2 BR PI0614372 A2 BRPI0614372 A2 BR PI0614372A2 BR PI0614372 A BRPI0614372 A BR PI0614372A BR PI0614372 A2 BRPI0614372 A2 BR PI0614372A2
- Authority
- BR
- Brazil
- Prior art keywords
- current
- voltage
- stepped
- measured
- electrodes
- Prior art date
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 23
- 230000004044 response Effects 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 29
- 238000012545 processing Methods 0.000 claims abstract description 23
- 230000001052 transient effect Effects 0.000 claims abstract description 20
- 238000004590 computer program Methods 0.000 claims description 4
- 238000009826 distribution Methods 0.000 description 6
- 238000011835 investigation Methods 0.000 description 6
- 230000007423 decrease Effects 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 150000002430 hydrocarbons Chemical class 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000005755 formation reaction Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 239000012530 fluid Substances 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/02—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/02—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current
- G01V3/04—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current using DC
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Measuring Instrument Details And Bridges, And Automatic Balancing Devices (AREA)
Abstract
MEDIçõES DE RESISTIVIDADE CC DE MúLTIPLOS TRANSIENTES. A presente invenção refere-se a um método para medir as variações de resistividade na terra que compreende passar uma corrente transi- ente entre dois eletrodos de fonte; medir a corrente transiente nos eletrodos de fonte; medir a voltagem transiente resultante entre pelo menos um par de eletrodos receptores; estimar uma ou mais funções de processamento para aplicar na corrente de entrada medida para prover um perfil de corrente escalonada, e aplicar as mesmas uma ou mais funções de processamento na voltagem medida para prover uma estimativa da voltagem de resposta escalonada entre os receptores, utilizando tanto a corrente escalonada quanto a voltagem de resposta escalonada estimada para determinar a resistência da terra aparente resultante, e utilizar a resistência da terra aparente para determinar a resistividade da terra.
Description
Relatório Descritivo da Patente de Invenção para "MEDIÇÕESDE RESISTIVIDADE CC DE MÚLTIPLOS TRANSIENTES".
A presente invenção refere-se a medições de resistividade CCde múltiplos transientes. Além disso, a invenção refere-se a um sistema emétodo para estimar a resposta da terra utilizando tais medições de resistivi-dade CC de múltiplos transientes, para por meio disto detectar, por exemplo,as formações que contém hidrocarbonetos ou que contém água.
ANTECEDENTES DA INVENÇÃO
As rochas porosas estão saturadas com fluidos. Os fluidos po-dem ser água, gás ou óleo, ou uma mistura de todos três. O fluxo de corren-te na terra é determinado pelas resistividades de tais rochas, as quais sãoafetadas pelos fluidos de saturação. Por exemplo, as rochas porosas satura-das com salmoura são muito menos resistivas do que as mesmas rochascheias com hidrocarbonetos. Por conseguinte, o objetivo geofísico é deter-minar se os hidrocarbonetos estão presentes pela medição da resistividadede formações geológicas. Se os testes que utilizam outros métodos, por e-xemplo a exploração sísmica, sugerem que uma formação geológica tem opotencial de conter hidrocarbonetos, então antes de perfurá-la é importanteter alguma indicação sobre se a formação contém de fato os hidrocarbone-tos ou se esta contém primariamente água. Isto pode ser feito utilizando téc-nicas eletromagnéticas, e mais especificamente técnicas eletromagnéticasde domínio de tempo.
Por quase um século os métodos de corrente contínua (CC) têmsido utilizados na superfície da terra para determinar as distribuições de re-sistividade de subsuperfície. O trabalho mais antigo foi pelos irmãos Sc-hlumberger e por Wenner (Wenner, F.,1912, A method for measuring earthresistivity. US Bureau of Standards Bulletin, 12,469-478.). Existem três méto-dos de levantamento de resistividade CC gerais: ondas sonoras elétricasverticais (VES), perfilagem, e formação de imagem. Nos levantamentos VESas dimensões de uma rede de medições selecionadas aumentam enquantoque o ponto central da rede permanece fixo. Conforme a rede expande, ascorrentes penetram mais profundamente e as curvas de ondas sonoras re-sultantes são interpretadas como mudanças de resistividade com a profun-didade. Nos levantamentos de perfilagem tanto o tipo de rede quanto as su-as dimensões são selecionados para uma profundidade de investigação es-pecífica. A rede é movida ao longo da superfície para determinar as varia-ções de resistividade laterais. A formação de imagem, ou método de tomo-grafia de resistência da terra (ERT), combina a VES e a perfilagem. Nesta,um grande número de eletrodos é colocado dentro do solo, usualmente i-gualmente espaçados, e são conectados com um cabo de múltiplos núcleosem um medidor de resistividade. O sistema opera sob controle de softwareonde quaisquer dois eletrodos podem ser selecionados como eletrodos decorrente, e quaisquer dois outros podem ser selecionados como eletrodos depotencial (voltagem). Os valores de resistividade são obtidos sobre uma se-ção transversal sob o perfil e indica tanto as variações laterais quanto deprofundidade.
As medições CC são feitas injetando uma corrente contínua en-tre dois eletrodos de fonte S1 eS2e medindo a voltagem entre dois eletro-dos receptores R1 e R2. Da corrente I injetada nos eletrodos de fonte e davoltagem V medida nos eletrodos receptores, um valor da resistência da ter-ra aparente é obtido pela lei de Ohm:
<formula>formula see original document page 3</formula>
Um valor aparente da resistividade de subsuperfície papp podeser obtido desta resistência como
<formula>formula see original document page 3</formula>
na qual k tem unidades de metros e é um fator geométrico que depende dadisposição dos quatro eletrodos. Utilizando esta resistividade aparente, aresistividade de subsuperfície verdadeira pode ser determinada. Uma boavisão das técnicas para determinar o valor verdadeiro da distribuição de re-sistividade de subsuperfície dos valores de resistividade aparente é dada porLoke M.H., 1999, em Electrical imaging surveys for environmental and engi-neering studies (http://www.abem.com/ftp/Loke/2Dnotes.pdf).
Existem diversas configurações bem-conhecidas dos quatro ele-trodos para as medições CC, três das quais estão ilustradas nas Figuras 1 a
3. A Figura 1 mostra a rede de Wenner, a Figura 2 mostra a rede de Schlum-berger, e a Figura 3 mostra a rede dipolo-dipolo. Para cada rede existe umdiferente fator k, como abaixo fornecido
<formula>formula see original document page 4</formula>
A rede dipolo-dipolo da Figura 3 determina as variações de resis-tividade laterais melhor do que as variações de profundidade. A diferença depotencial a ser medida entre os eletrodos de potencial decai com o cubo dadistância dos eletrodos de corrente. Isto tem restringido a configuração darede para propósitos práticos a valores de n < 6. Edwards (Edwards, L.S.,-1977, A modified pseudosection for resistivity and IP. Geophysics, 42,1020--1036) discutiu a apresentação de pseudo-seções para resistividade e IP es-pecificamente para a rede dipolo-dipolo. A teoria e os resultados práticos sãoapresentados para n < 6 e referência é feita a uma "rede ideal" com n = ∞
Como mostrado na Figura 5, a profundidade de investigação destá relacionada com os comprimentos de dipolo a e a separação de dipolona. Neste caso, a profundidade de investigação máxima d para a qual a re-sistividade da terra pode ser inferida está relacionada à configuração e é daordem de (n+2)a/5. Na prática, como a amplitude de sinal nos receptoresdiminui aproximadamente como (na)"3, apesar do nível de ruído ser depen-dente de n e de a, a razão de sinal para ruído diminui como (na)"3. O nível desinal pode ser aumentado pelo aumento da corrente injetada na fonte e peloaumento da distância de dipolo a. Com os níveis de corrente que podem serseguramente injetados no solo, o n normalmente não é maior do que apro-ximadamente 6, e segue que
d < 1.6a (6)
Em outras palavras, a profundidade de investigação d é menordo que 1,6 vezes a separação a entre os eletrodos de fonte ou receptores. Oaumento de a aumenta a profundidade de investigação mas reduz a resolu-ção da distribuição de resistividade de subsuperfície mapeada.
Quando fazendo as medições CC, foi descoberto que o eletrodotorna-se polarizado se a corrente for da mesma polaridade por um longotempo e uma medição falsa da voltagem na terra entre R1 e R2 é obtida.Duas estratégias são utilizadas para superar este problema. Uma é utilizareletrodos não polarizantes. A outra proposta, mais popular, é de mudar apolaridade da corrente CC periodicamente, isto de fato fornece uma correntede entrada de onda quadrada alternada, ou CA. Todos os equipamentosmodernos utilizam está técnica.
Na proposta de CA variações podem ser introduzidas, tal comocomutando a corrente para zero por certos períodos, por exemplo comomostrado na Figura 4. Neste caso, o período T entre as comutações é tipi-camente da ordem de 1 segundo. Esta função aqui mostrada repete-se acada 4T. Nos períodos em que a corrente está ligada entre os eletrodos defonte S1 e S2, uma estimativa da voltagem CC resultante entre R1 e R2 éfeita. De fato, a voltagem entre R1 e R2 não é exatamente constante nestesperíodos: o sinal leva algum tempo para atingir um valor de estado estável eexiste ruído. Para compensar pelo ruído várias técnicas de média são utili-zadas. Deve ser notado que neste contexto o termo "CC" em "resistividadeCC" significa essencialmente uma baixa freqüência, ao invés de freqüênciazero. Isto já é aceito na indústria de exploração geofísica.
Cada vez que a corrente nos eletrodos de fonte é comutada, aterra responde e a voltagem nos eletrodos receptores muda. O fluxo de cor-rente no solo é governado pela equação de difusão e leva tempo para a res-posta no receptor atingir um estado estável. Isto é bem-conhecido, e a esti-mativa do nível de CC é medida em um intervalo de tempo na direção dofinal do período que a corrente está ligada. Por exemplo, o Geopulse TigreResistivity Meter utiliza as medições de voltagem que são feitas durante osúltimos 4/5 do período ligado da corrente (User's Manual, Geopulse Resisiti-vity Meter, Campus International Products Limited, Concept House, 8 TheTownsend Centre, Blackburn Road, Dunstable, Bedford, England LU5 5BQ).Assim é conhecido que existe uma resposta transiente para a comutação dacorrente, mas após um certo tempo é considerado que a resposta de estadoestável é atingida.
Nesta proposta convencional, o tempo esperado para atingir oestado estável está baseado na experiência, apesar de ser totalmente com-preendido que o valor de estado estável final nunca é atingido. Por razõespráticas, no entanto, para esta configuração, e dentro das limitações da ins-trumentação e do ruído, o valor de estado estável é usualmente atingido a-pós aproximadamente 0,1 segundo. Em alguns levantamentos de baixa re-solução profundos, onde a escala do conjunto é aumentada por uma ordemde magnitude ou tanto, o valor de estado estável pode ser atingido, por pro-pósitos práticos, somente após poucos segundos. Um exemplo de uma res-posta de voltagem medida para uma corrente ligada utilizando a configura-ção dipolo-dipolo está mostrado na Figura 6. Isto é denominada a respostaescalonada. A primeira parte da resposta mostra um degrau inicial na volta-gem; neste caso específico isto é seguido por uma pequena diminuição se-guida por um aumento em voltagem que parece estar tendendo para um va-lor de estado estável após aproximadamente 0,1 s.
A Figura 7 mostra a derivada de tempo para a resposta escalo-nada da Figura 6. Isto é conhecido como a resposta de impulso. O grandepico inicial na resposta de impulso corresponde ao degrau inicial na respostaescalonada. Este é seguido por uma diminuição, um aumento para um picomenor, e então uma diminuição muito gradual em amplitude, tendendo parazero conforme o tempo aumenta após isto. A duração da resposta de impul-so é infinita, tal como a duração da resposta escalonada é infinita. No entan-to, conforme a amplitude da resposta de impulso fica cada vez menor, torna-se cada vez mais difícil medir. Quando a amplitude é pequena demais paramedir, isto efetivamente define a duração da resposta de impulso transiente.
Um objetivo da presente invenção é de aperfeiçoar a sensibili-dade de medições de resistividade da terra.
SUMÁRIO DA INVENÇÃO
De acordo com um aspecto da presente invenção, está providoum método para medir as variações de resistividade na terra que compreen-de passar uma corrente entre dois eletrodos de fonte; medir a corrente tran-siente nos eletrodos de fonte; medir a voltagem transiente resultante entrepelo menos um par de eletrodos receptores; estimar uma ou mais funçõesde processamento para aplicar na corrente de entrada medida para proverum perfil de entrada escalonado, e aplicar as mesmas uma ou mais funçõesde processamento na voltagem medida para prover uma estimativa da volta-gem de resposta escalonada entre os receptores, utilizando tanto a correnteescalonada quanto a voltagem de resposta escalonada estimada para de-terminar a resistência da terra aparente resultante, e utilizar a resistência daterra aparente para determinar a resistividade da terra. A determinação daresistividade pode ser feita utilizando os métodos de inversão de resistivida-de CC padrão.
A função de processamento poderia ser um filtro. Uma propostaalternativa é de utilizar um filtro para converter a corrente medida em um im-pulso e integrar o resultado para obter um degrau na corrente; aplicar omesmo filtro na voltagem medida e integrando o resultado geraria a volta-gem de resposta escalonada. Podem existir outros métodos para conseguiro mesmo resultado. O ponto chave é que os métodos de processamentopodem ser utilizados para converter a resposta de corrente de entrada e devoltagem transiente medida em um degrau e uma resposta escalonada, res-pectivamente, por meio disto gerando a resposta CC.
A utilização desta técnica aperfeiçoa a sensibilidade da mediçãode resistividade CC a um tal ponto que o valor de n pode ser aumentado porpelo menos uma ordem de magnitude para aproximadamente 80 ou tanto, oque portanto permite que d seja aumentado para aproximadamente 120a.Com isto, a profundidade de investigação pode ser aumentada por mais doque uma ordem de magnitude sem uma redução na resolução da distribui-ção de resistividade de subsuperfície mapeada.
De acordo com a presente invenção, tanto a corrente de entradana fonte e a voltagem medida no receptor são medidas e gravadas comouma função de tempo. As funções de tempo são convertidas para dados di-gitais e processadas em um computador utilizando um software de computa-dor. O processamento de computador pode incluir muitas operações, inclu-indo a remoção de ruído cultural - tal como a geração da rede de 50 Hz ou60 Hz - e a remoção da resposta de sistema de gravação por deconvolução,como discutido na WO 03/023452 A1.
A corrente de entrada pode ser um degrau simples, uma ondaquadrada, ou mesmo um sinal de modulação de retorno para zero - freqüen-temente conhecido como "cicloativo" - mostrado na Figura 4. Também, pode-ria, por exemplo, ser uma seqüência binária pseudo-randômica de uns e ze-ros, ou uns e menos uns, o que permitiria que a resposta escalonada da ter-ra entre a fonte e o receptor (e com isto o valor de estado estável) fosse de-duzida.
De acordo com outro aspecto da presente invenção, está provi-do um sistema para medir as variações de resistividade na terra que com-preende uma fonte para passar uma corrente entre dois eletrodos de fonte;um dispositivo de medição de corrente para medir uma corrente transientenos eletrodos de fonte; um dispositivo de medição de voltagem para medir avoltagem transiente resultante entre pelo menos um par de eletrodos recep-tores; e um processador configurado para estimar uma ou mais funções deprocessamento para aplicar na corrente de entrada medida para prover umperfil de corrente escalonada; aplicar as mesmas uma ou mais funções deprocessamento na voltagem medida para prover uma estimativa da voltagemde resposta escalonada entre os receptores, utilizar tanto a corrente escalo-nada quanto a voltagem de resposta escalonada estimada para determinar aresistência da terra aparente resultante, e utilizar a resistência da terra apa-rente para determinar a resistividade da terra.
De acordo com ainda outro aspecto da presente invenção, estáprovido um programa de computador para determinar as variações de resis-tividade na terra utilizando uma medida de corrente de entrada entre um parde eletrodos de fonte e uma medida da voltagem transiente resultante entrepelo menos um par de eletrodos receptores, o programa de computador ten-do um código ou instruções para estimar uma ou mais funções de proces-samento para aplicar na corrente de entrada medida para prover um perfil decorrente escalonada; aplicar as mesmas uma ou mais funções de proces-samento na voltagem medida para prover uma estimativa da voltagem deresposta escalonada entre os receptores; utilizar tanto a corrente escalonadaquanto a voltagem de resposta escalonada estimada para determinar a re-sistência da terra aparente resultante, e utilizar a resistência da terra aparen-te para determinar a resistividade da terra.
BREVE DESCRIÇÃO DOS DESENHOS
Vários aspectos da invenção serão agora descritos como exem-pio somente e com referência aos desenhos acompanhantes, nos quais:
Figura 8 é um diagrama de fluxo de um método para determinara resistividade da terra utilizando medições de CC;
Figura 9 é um gráfico da corrente de entrada medida após umfiltro ser aplicado;
Figura 10 é um gráfico da voltagem medida após a aplicação domesmo filtro utilizado para os resultados da Figura 9; e
Figura 11 é uma representação esquemática de um sistema noqual um único par de eletrodos de fonte está associado com uma pluralidadede pares de eletrodos receptores.
DESCRIÇÃO ESPECÍFICA DOS DESENHOS
As disposições de eletrodos de fonte / receptores das Figuras 3e 5 podem ser utilizadas para implementar a invenção. Isto requer a mediçãoda função de tempo de corrente de entrada de fonte, o que pode ser feitocom um medidor de corrente posicionado idealmente entre os eletrodos defonte, de modo que uma medida direta da corrente de entrada possa ser fei-ta. O medidor de corrente pode ser de qualquer forma adequada, por exem-plo um magnetômetro.
A Figura 9 mostra o método no qual a invenção está incorpora-da. Em uma primeira etapa uma corrente de entrada transiente é passadaentre dois eletrodos de fonte e medida entre dois eletrodos receptores. Afunção de corrente de entrada medida pode ser expressa como, digamos,l(t), onde t é o tempo. Esta corrente medida é convertida para uma funçãoescalonada de, digamos, amplitude unitária, pela convolução com um filtrof(t) como segue:
<formula>formula see original document page 10</formula>
onde H(t) é uma função de etapa unitária e t0 é o tempo da etapa. O resulta-do está mostrado na Figura 9.
A resposta de voltagem medida no receptor é, digamos, M(t).
Esta resposta pode ser convolvida com o mesmo filtro /(t) para gerar a res-posta escalonada X(t):
<formula>formula see original document page 10</formula>
A Figura 10 mostra o resultado da filtragem de uma voltagemmedida para obter a resposta escalonada. Em princípio, como vimos, levauma quantidade de tempo infinita para a resposta escalonada atingir o valorde estado estável VTO, e assim V00 não é realmente medida. Isto é claro étambém verdadeiro para as assim denominadas medições de resistividadeCC convencionais. Uma estimativa do valor final Vto pode ser feita por umavariedade de propostas, uma das quais está abaixo delineada. O resultado éque a voltagem CC Vto para uma corrente de entrada de 1 A pode ser encon-trada e a resistência da terra é assim
<formula>formula see original document page 10</formula>
O valor de estado estável para a voltagem Vto pode ser estimadopor uma variedade de técnicas. Referindo à Figura 10, a voltagem medidaem tempos tardios t > ti é descrita como
<formula>formula see original document page 10</formula>
na qual ti pode ser variado, e é dependente de deslocamento. Os valores deVto e a são encontrados ajustando a curva descrita pela equação (9) aos da-dos e minimizando o erro de desajuste. Este valor pode então ser utilizadopara determinar a resistência CC. O valor resultante da resistência CC paracada par de fonte - receptor pode ser colocado em um programa de inversãode resistividade bidimensional ou tridimensional, tal como o RES2DINV e oRES3DINV, como descrito em Loke (1999), para obter um mapa de distribui-ção de resistividade de subsuperfície. Outros parâmetros de entrada para oprograma seriam as coordenadas dos eletrodos.
Para um dado par de fonte - receptor, a medição pode ser repe-tida muitas vezes e feita a média (ou "empilhada") permitindo que a razão desinal para ruído seja aperfeiçoada e a estimativa do valor de V00 seja refina-da. Para um dado par de eletrodos de fonte, muitos pares de eletrodos re-ceptores podem estar providos e assim pode-se fazer muitas medições si-multâneas das voltagens recebidas, como mostrado na Figura 11. Neste ca-so, os eletrodos de fonte estão separados por uma distância a, e cada parde eletrodos receptores está separado por uma distância a. Neste exemplo né escolhido ser inteiro e o seu valor é m, m+1, m+2, m+3, etc. Apesar deuma configuração na qual a separação de eletrodos na fonte e em cada parde receptor são iguais ser mostrada, será apreciado que muitas diferentesconfigurações são possíveis. Igualmente, todo o conjunto pode ser desloca-do lateralmente, por exemplo por uma distância a ou por múltiplos de a paraobter uma múltipla cobertura da distribuição de resistividade de subsuperfí-cie.
Uma pessoa versada na técnica apreciará que variações dasdisposições descritas são possíveis sem afastar-se da invenção. Por exem-plo, o método pode ser aplicado em terra ou offshore. No caso marinho, oseletrodos não precisam necessariamente estar sobre o fundo do mar, e po-dem não precisar ser estacionários. É concebível que a configuração de ele-trodo de fonte - receptor inteira poderia ser rebocada atrás de um navio quesupre a corrente para a fonte e grava a corrente de fonte e as voltagens re-cebidas. Consequentemente, a descrição acima de uma modalidade especí-fica é feita como exemplo somente e não para os propósitos de limitação.Ficará claro para a pessoa versada na técnica que modificações menorespodem ser feitas sem mudanças significativas na operação descrita.
Claims (15)
1. Método para medir as variações de resistividade na terra quecompreende passar uma corrente transiente entre dois eletrodos de fonte;medir a corrente transiente nos eletrodos de fonte; medir a voltagem transi-ente resultante entre pelo menos um par de eletrodos receptores; estimaruma ou mais funções de processamento para aplicar na corrente de entradamedida para prover um perfil de corrente escalonada, e aplicar as mesmasuma ou mais funções de processamento na voltagem medida para proveruma estimativa da voltagem de resposta escalonada entre os receptores,utilizando tanto a corrente escalonada quanto a voltagem de resposta esca-lonada estimada para determinar a resistência da terra aparente resultante,e utilizar a resistência da terra aparente para determinar a resistividade daterra.
2. Método de acordo com a reivindicação 1, em que a etapa deaplicar a função de processamento estimada na voltagem medida envolve aconvolução das voltagens recebidas com a função de processamento esti-mada para encontrar as respostas escalonadas correspondentes.
3. Método de acordo com a reivindicação 1, em que a função deprocessamento é um filtro para converter a corrente de entrada medida noperfil escalonado.
4. Método de acordo com a reivindicação 1, em que as uma oumais funções de processamento compreendem um filtro para converter acorrente medida em um impulso e um meio integrar o resultado para obter operfil de entrada escalonado.
5. Método de acordo com qualquer uma das reivindicações pre-cedentes, que compreende remover o ruído cultural.
6. Método de acordo com qualquer uma das reivindicações pre-cedentes, que compreende estimar uma voltagem de estado estável de umaresposta escalonada.
7. Método de acordo com qualquer uma das reivindicações pre-cedentes, em que os eletrodos de fonte e receptores estão providos em umambiente marinho.
8. Sistema para medir as variações de resistividade na terra quecompreende uma fonte para passar uma corrente entre dois eletrodos defonte; um dispositivo de medição para medir a corrente transiente nos ele-trodos de fonte; um dispositivo de medição de voltagem para medir a volta-gem transiente resultante entre pelo menos um par de eletrodos receptores;e um processador configurado para estimar uma ou mais funções de pro-cessamento para aplicar na corrente de entrada medida para prover um per-fil de corrente escalonada; aplicar as mesmas uma ou mais funções de pro-cessamento na voltagem medida para prover uma estimativa da voltagem deresposta escalonada entre os receptores, utilizar tanto a corrente escalonadaquanto a voltagem de resposta escalonada estimada para determinar a re-sistência da terra aparente resultante, e utilizar a resistência da terra aparen-te para determinar a resistividade da terra.
9. Sistema de acordo com a reivindicação 8, em que o proces-sador é operável para aplicar a função de processamento estimada na volta-gem medida por convolução das voltagens recebidas com a função de pro-cessamento estimada para encontrar as respostas escalonadas correspon-dentes.
10. Sistema de acordo com a reivindicação 8 ou 9, em que afunção de processamento é um filtro para converter a corrente de entradamedida no perfil escalonado.
11. Sistema de acordo com a reivindicação 8, em que as uma oumais funções de processamento compreendem um filtro para converter acorrente medida em um impulso e um meio para integrar o resultado paraobter o perfil de entrada escalonado.
12. Sistema de acordo com qualquer uma das reivindicações 8 a 11, em que o processador é operável para o ruído cultural.
13. Sistema de acordo com qualquer uma das reivindicações 8 a 12, em que o processador é operável para estimar uma voltagem de estadoestável de uma resposta escalonada.
14. Sistema de acordo com qualquer uma das reivindicações 8 a 13, em que os eletrodos de fonte e receptores estão providos em um ambi-ente marinho.
15. Programa de computador, de preferência em um meio legívelpor computador ou portador, para determinar as variações de resistividadena terra utilizando uma medida de corrente de entrada entre um par de ele-trodos de fonte e uma medida da voltagem transiente resultante entre pelomenos um par de eletrodos receptores, o programa de computador tendo umcódigo ou instruções para estimar uma ou mais funções de processamentopara aplicar na corrente de entrada medida para prover um perfil de correnteescalonada; aplicar as mesmas uma ou mais funções de processamento navoltagem medida para prover uma estimativa da voltagem de resposta esca-lonada entre os receptores; utilizar tanto a corrente escalonada quanto a vol-tagem de resposta escalonada estimada para determinar a resistência daterra aparente resultante, e utilizar a resistência da terra aparente para de-terminar a resistividade da terra.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| GBGB0516153.4A GB0516153D0 (en) | 2005-08-05 | 2005-08-05 | Multi-transient dc resistivity measurements |
| GB516153.4 | 2005-08-05 | ||
| PCT/GB2006/002940 WO2007017657A1 (en) | 2005-08-05 | 2006-08-04 | Multi-transient dc resistivity measurements |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| BRPI0614372A2 true BRPI0614372A2 (pt) | 2011-03-22 |
| BRPI0614372B1 BRPI0614372B1 (pt) | 2018-05-08 |
Family
ID=34984166
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| BRPI0614372A BRPI0614372B1 (pt) | 2005-08-05 | 2006-08-04 | método e sistema para medir as variações de resistividade na terra |
Country Status (9)
| Country | Link |
|---|---|
| US (1) | US8274288B2 (pt) |
| CN (1) | CN101283292B (pt) |
| AU (1) | AU2006277756B2 (pt) |
| BR (1) | BRPI0614372B1 (pt) |
| CA (1) | CA2617701C (pt) |
| EA (1) | EA011317B1 (pt) |
| GB (2) | GB0516153D0 (pt) |
| NO (1) | NO341079B1 (pt) |
| WO (1) | WO2007017657A1 (pt) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8063642B2 (en) | 2008-06-11 | 2011-11-22 | Mtem Ltd | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources |
| US8258791B2 (en) | 2009-01-27 | 2012-09-04 | Mtem Ltd. | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources to impart electromagnetic signals into a subsurface formation and thereby determining a formation response to each signal |
| US8143897B2 (en) | 2009-02-11 | 2012-03-27 | Mtem Ltd. | Short-offset transient electromagnetic geophysical surveying |
| CN101871974B (zh) * | 2010-06-18 | 2013-01-09 | 华南理工大学 | 一种阻抗谱的测量方法 |
| US8514656B2 (en) | 2011-03-18 | 2013-08-20 | Pgs Geophysical As | Sensor arrangement for detecting motion induced noise in towed marine electromagnetic sensor streamers |
| US9383469B2 (en) | 2012-04-30 | 2016-07-05 | Pgs Geophysical As | Methods and systems for noise-based streamer depth profile control |
| US8898102B2 (en) | 2012-05-25 | 2014-11-25 | Pgs Geophysical As | Adaptive computational grid for inverse problems |
| US20140058677A1 (en) * | 2012-08-23 | 2014-02-27 | Westerngeco, L.L.C. | Method for processing electromagnetic data |
| US9625600B2 (en) | 2012-12-04 | 2017-04-18 | Pgs Geophysical As | Systems and methods for removal of swell noise in marine electromagnetic surveys |
| US9575205B2 (en) | 2013-01-17 | 2017-02-21 | Pgs Geophysical As | Uncertainty-based frequency-selected inversion of electromagnetic geophysical data |
| US10338250B2 (en) | 2013-03-14 | 2019-07-02 | Pgs Geophysical As | Method of removing incoherent noise |
| CN105403921B (zh) * | 2015-12-09 | 2018-03-06 | 重庆地质仪器厂 | 一种转换器 |
| CN106443795B (zh) * | 2016-07-21 | 2018-05-15 | 安徽惠洲地质安全研究院股份有限公司 | 一种双模网络并行电法am数据推演合成abm数据方法 |
| CN108873074B (zh) * | 2018-04-18 | 2020-03-24 | 浙江大学 | 一种电极随机分布式高密度电阻率测量方法及勘探系统 |
| US11654635B2 (en) | 2019-04-18 | 2023-05-23 | The Research Foundation For Suny | Enhanced non-destructive testing in directed energy material processing |
| CN114518604A (zh) * | 2020-11-19 | 2022-05-20 | 中核(内蒙古)矿业投资有限公司 | 一种电阻率测深数据处理解释与绘图方法 |
| CN114778948A (zh) * | 2022-06-17 | 2022-07-22 | 中铁大桥科学研究院有限公司 | 动水隧道岩体电阻率监测方法及相关设备 |
| CN117607210B (zh) * | 2024-01-24 | 2024-05-31 | 中南大学 | 一种分布式偶极-偶极电法监测方法及系统 |
Family Cites Families (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US2177346A (en) * | 1938-04-18 | 1939-10-24 | E E Rosaire | Exploration by incremental wave distortion |
| US2230803A (en) * | 1938-08-25 | 1941-02-04 | Paul W Klipsch | Wave synthesizing network |
| US4041372A (en) * | 1975-09-08 | 1977-08-09 | Continental Oil Company | Apparatus for multi-channel induced polarization surveying |
| FR2475745A2 (fr) * | 1980-02-13 | 1981-08-14 | Geophysique Cie Gle | Prospection geophysique a courants transitoires, avec traitement par indicatrices de variations et soustraction d'une reponse regionale |
| MA18895A1 (fr) * | 1979-07-09 | 1981-04-01 | Cie Generale De Geophysique Sa | Procede et dispositif de prospection geophysique a courants transitoires |
| US4467283A (en) * | 1981-09-23 | 1984-08-21 | Geo-Optics, Ltd. | Electrical geophysical exploration system with addressable current probes |
| US4835473A (en) * | 1984-07-31 | 1989-05-30 | Board Of Regents, The University Of Texas System | Method using a d.c. source for determining resistivity distribution of the earth as a function of depth |
| AU654945B2 (en) | 1990-03-21 | 1994-12-01 | Bergwerksverband Gmbh | Bedrock exploration system using transient electromagnetic measurements |
| US5877995A (en) * | 1991-05-06 | 1999-03-02 | Exxon Production Research Company | Geophysical prospecting |
| US5563513A (en) * | 1993-12-09 | 1996-10-08 | Stratasearch Corp. | Electromagnetic imaging device and method for delineating anomalous resistivity patterns associated with oil and gas traps |
| GB2336008B (en) * | 1998-04-03 | 2000-11-08 | Schlumberger Holdings | Simulation system including a simulator and a case manager adapted for organizing data files |
| USRE40321E1 (en) * | 1999-09-15 | 2008-05-20 | Exxonmobil Upstream Research Co. | Remote reservoir resistivity mapping |
| MY131017A (en) * | 1999-09-15 | 2007-07-31 | Exxonmobil Upstream Res Co | Remote reservoir resistivity mapping |
| US6477113B2 (en) * | 2000-03-21 | 2002-11-05 | Exxonmobil Upstream Research Company | Source waveforms for electroseismic exploration |
| US6664788B2 (en) * | 2001-05-02 | 2003-12-16 | Exxonmobil Upstream Research Company | Nonlinear electroseismic exploration |
| FR2826449B1 (fr) * | 2001-06-26 | 2003-09-26 | Schlumberger Services Petrol | Procede et installation de determination de la position du front forme a l'interface entre deux fluides contenus dans un reservoir |
| US20030035376A1 (en) * | 2001-08-20 | 2003-02-20 | Xiaofen Chen | Derivation of composite step-function response |
| GB0121719D0 (en) * | 2001-09-07 | 2001-10-31 | Univ Edinburgh | Method for detection fo subsurface resistivity contrasts |
| US6819111B2 (en) * | 2002-11-22 | 2004-11-16 | Baker Hughes Incorporated | Method of determining vertical and horizontal resistivity, and relative dip in anisotropic earth formations having an arbitrary electro-magnetic antenna combination and orientation with additional rotation and position measurements |
| WO2005015262A1 (ja) * | 2003-07-16 | 2005-02-17 | Canon Kabushiki Kaisha | 地中探査装置、システム及び方法 |
| US7742875B2 (en) * | 2004-07-01 | 2010-06-22 | Exxonmobil Upstream Research Company | Method for geologic modeling through hydrodynamics-based gridding (Hydro-Grids) |
| US8063642B2 (en) * | 2008-06-11 | 2011-11-22 | Mtem Ltd | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources |
-
2005
- 2005-08-05 GB GBGB0516153.4A patent/GB0516153D0/en not_active Ceased
-
2006
- 2006-08-04 AU AU2006277756A patent/AU2006277756B2/en not_active Ceased
- 2006-08-04 BR BRPI0614372A patent/BRPI0614372B1/pt not_active IP Right Cessation
- 2006-08-04 WO PCT/GB2006/002940 patent/WO2007017657A1/en not_active Ceased
- 2006-08-04 CA CA2617701A patent/CA2617701C/en not_active Expired - Fee Related
- 2006-08-04 GB GB0802536A patent/GB2443762B/en not_active Expired - Fee Related
- 2006-08-04 EA EA200800347A patent/EA011317B1/ru not_active IP Right Cessation
- 2006-08-04 CN CN200680037172XA patent/CN101283292B/zh not_active Expired - Fee Related
- 2006-08-04 US US11/990,035 patent/US8274288B2/en not_active Expired - Fee Related
-
2008
- 2008-02-01 NO NO20080618A patent/NO341079B1/no not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| CN101283292A (zh) | 2008-10-08 |
| GB2443762B (en) | 2010-02-17 |
| WO2007017657A1 (en) | 2007-02-15 |
| AU2006277756B2 (en) | 2011-04-14 |
| GB0802536D0 (en) | 2008-03-19 |
| CA2617701C (en) | 2014-09-23 |
| US20100315088A1 (en) | 2010-12-16 |
| GB2443762A (en) | 2008-05-14 |
| GB0516153D0 (en) | 2005-09-14 |
| EA011317B1 (ru) | 2009-02-27 |
| NO20080618L (no) | 2008-03-14 |
| CA2617701A1 (en) | 2007-02-15 |
| NO341079B1 (no) | 2017-08-21 |
| US8274288B2 (en) | 2012-09-25 |
| AU2006277756A2 (en) | 2008-05-29 |
| EA200800347A1 (ru) | 2008-08-29 |
| CN101283292B (zh) | 2011-09-14 |
| BRPI0614372B1 (pt) | 2018-05-08 |
| AU2006277756A1 (en) | 2007-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| BRPI0614372A2 (pt) | medições de resistividade cc de múltiplos transientes | |
| Dahlin et al. | Measuring techniques in induced polarisation imaging | |
| BRPI0915113A2 (pt) | mpétodo para inspeção eletromagnética de subsuperfície usando duas ou mais fontes eletromagnéticas acionadas simultaneamente | |
| BRPI1000052A2 (pt) | método para levantamento eletromagnético utilizando duas ou mais fontes eletromagnéticas simultaneamente atuadas | |
| BRPI1000323A2 (pt) | levantamento geofìsico eletromagnético transiente de pequeno deslocamento | |
| BRPI0715832A2 (pt) | aperfeiÇoamento em exploraÇço marinha | |
| Mogilatov et al. | Application of the marine circular electric dipole method in high latitude Arctic regions using drifting ice floes | |
| BRPI0616231B1 (pt) | Methods for improving an electromagnetic prospection with a controlled source of a underground region and to produce hydrocarbons from a subterranean region | |
| BRPI0903052A2 (pt) | método para atenuar o ruìdo correlacionado em dados de levantamento eletromagnético de fonte controlada | |
| US20100235100A1 (en) | Method for determining resistivity anisotropy from earth electromagnetic responses | |
| EA006537B1 (ru) | Способ морской геоэлектроразведки (варианты) | |
| Flores et al. | On the estimation of the maximum depth of investigation of transient electromagnetic soundings: the case of the Vizcaino transect, Mexico | |
| BR102012030628A2 (pt) | Métodos e aparelho para uma rápida determinação da profundidade e da resistência transversal alvo | |
| EA006536B1 (ru) | Способ геоэлектроразведки (варианты) | |
| Hauquin et al. | The detection of old masonry tunnels as low electrical resistivity anomalies-application to one covered stream tunnel of the Cevennes Mountain region (France) | |
| Boulay et al. | Overcoming signal-to-noise challenges with pole-dipole resistivity monitoring at a hydroelectric dam site | |
| Ziolkowski et al. | New technology to acquire, process, and interpret transient EM data | |
| Li et al. | Higher-order statistics correlation stacking for DC electrical data in the wavelet domain | |
| Zerilli et al. | Building starting model for full-waveform inversion using broadband CSEM-driven velocity model building to improve complex salt imaging | |
| RU2174243C1 (ru) | Способ геоэлектроразведки | |
| CA1246673A (en) | Electromagnetic array profiling survey method | |
| Toll et al. | Development of automated multi-electrode resistivity system for laboratory measurements | |
| Legault et al. | Two-dimensional joint inversion of ZTEM and MT plane-wave EM data for near surface applications | |
| Hamidi et al. | Investigation of geophysical anomalies in soil depth and assessment of salt irrigation impact on soil quality at Ain El Atti experimental station, south-eastern of Morocco. | |
| Zerilli* et al. | Broadband CSEM-guided starting model for full waveform inversion–Impact in the Brazilian deep offshore |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| B06A | Patent application procedure suspended [chapter 6.1 patent gazette] | ||
| B09A | Decision: intention to grant [chapter 9.1 patent gazette] | ||
| B16A | Patent or certificate of addition of invention granted [chapter 16.1 patent gazette] | ||
| B21F | Lapse acc. art. 78, item iv - on non-payment of the annual fees in time |
Free format text: REFERENTE A 17A ANUIDADE. |
|
| B24J | Lapse because of non-payment of annual fees (definitively: art 78 iv lpi, resolution 113/2013 art. 12) |
Free format text: EM VIRTUDE DA EXTINCAO PUBLICADA NA RPI 2734 DE 30-05-2023 E CONSIDERANDO AUSENCIA DE MANIFESTACAO DENTRO DOS PRAZOS LEGAIS, INFORMO QUE CABE SER MANTIDA A EXTINCAO DA PATENTE E SEUS CERTIFICADOS, CONFORME O DISPOSTO NO ARTIGO 12, DA RESOLUCAO 113/2013. |