BRPI1000912A2 - método para a determinação de anisotropia de resistividade a partir de respostas eletromagnéticas da terra - Google Patents
método para a determinação de anisotropia de resistividade a partir de respostas eletromagnéticas da terra Download PDFInfo
- Publication number
- BRPI1000912A2 BRPI1000912A2 BRPI1000912-4A BRPI1000912A BRPI1000912A2 BR PI1000912 A2 BRPI1000912 A2 BR PI1000912A2 BR PI1000912 A BRPI1000912 A BR PI1000912A BR PI1000912 A2 BRPI1000912 A2 BR PI1000912A2
- Authority
- BR
- Brazil
- Prior art keywords
- resistivity
- response
- anisotropy
- step response
- formations
- Prior art date
Links
- 230000004044 response Effects 0.000 title claims abstract description 94
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 59
- 238000005755 formation reaction Methods 0.000 claims abstract description 59
- 238000006073 displacement reaction Methods 0.000 claims abstract description 24
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 21
- 238000005259 measurement Methods 0.000 claims abstract description 17
- 230000001052 transient effect Effects 0.000 claims abstract description 15
- 239000011435 rock Substances 0.000 claims abstract description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 12
- 230000014509 gene expression Effects 0.000 description 7
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 101150071746 Pbsn gene Proteins 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- FZNWJRXTACKOPU-UHFFFAOYSA-N 2-(2-methylthioethyl)malic acid Chemical compound CSCCC(O)(C(O)=O)CC(O)=O FZNWJRXTACKOPU-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000005457 optimization 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/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
-
- 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/08—Electric 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/083—Controlled source electromagnetic [CSEM] surveying
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
MéTODO PARA A DETERMINAçãO DE ANISOTROPIA DE RESISTIVIDADE A PARTIR DE RESPOSTAS ELETROMAGNéTICAS DA TERRA. A presente invenção refere-se a um método para a determinação de anisotropia de resistividade de formações de rocha de subsuperfície, a partir de medições de resposta a um campo eletromagnético transiente administrado na subsuperfície e medida em uma pluralidade de distâncias a partir de uma posição da administração, que inclui esta para cada deslocamento, uma resposta de etapa das formações é determinada em um momento a partir da administração selecionado de modo que a resposta de etapa esteja somente substancialmente relacionada à resistividade horizontal e em um momento a partir da administração selecionado de modo que a resposta de etapa esteja somente substancialmente relacionada à resistividade média. A resposta de etapa de resistividade horizontal e a resposta de etapa de resistividade média são usadas para determinar a anisotropia de resistividade.
Description
Relatório Descritivo da Patente de Invenção para "MÉTODO PARA A DETERMINAÇÃO DE ANISOTROPIA DE RESISTIVIDADE A PARTIR DE RESPOSTAS ELETROMAGNÉTICAS DA TERRA".
Antecedentes da invenção
Campo da invenção
A presente invenção refere-se geralmente ao campo de levan- tamento eletromagnético de formações na subsuperfície da Terra. Mais par- ticularmente, a invenção refere-se a métodos para a determinação de aniso- tropia de resistividade elétrica em formações de subsuperfície com o uso de medições eletromagnéticas.
Antecedente da Técnica
O levantamento eletromagnético é usado, entre outros propósi- tos, para a determinação da presença de estruturas contendo hidrocarbone- to na subsuperfície da Terra. A presença de estruturas contendo hidrocarbo- neto é tipicamente deduzida por meio da determinação da presença de alta resistividade na subsuperfície, devido ao fato de que a alta resistividade está associada a formações de subsuperfície que têm hidrocarbonetos dispostos nos espaços de poros das mesmas.
O levantamento eletromagnético inclui o que são chamadas de técnicas de levantamento de "fonte controlada". As técnicas de levantamen- to eletromagnético de fonte controlada incluem administrar uma corrente elé- trica ou um campo magnético na Terra, quando tais levantamentos são con- duzidos em terra, ou administrar o mesmo em sedimentos abaixo do fundo da água (fundo do mar), quando tais levantamentos são conduzidos em um ambiente marinho. As técnicas incluem a medição de voltagens e/ou cam- pos magnéticos induzidos em eletrodos, antenas e/ou magnetômetros dis- postos na superfície da Terra, no fundo do mar ou a uma profundidade sele- cionada na água. As voltagens e/ou os campos magnéticos são induzidos por meio da interação do campo eletromagnético causado pela corrente elé- trica e/ou campo magnético administrado na subsuperfície da Terra (através do fundo da água em levantamentos marítimos) com as formações da Terra de subsuperfície. O levantamento eletromagnético marítimo de fonte controlada conhecido na técnica inclui administrar a corrente elétrica alternada nos se- dimentos abaixo do fundo da água por meio da aplicação de corrente a partir de uma fonte, usualmente disposta em uma embarcação de levantamento, a um eletrodo bipolar rebocado pela embarcação de levantamento. Um eletro- do bipolar consiste tipicamente em um cabo elétrico isolado que tem dois eletrodos no mesmo em um espaçamento selecionado, às vezes de 300 a 1000 metros ou mais. A corrente alternada tem uma ou mais freqüências selecionadas, tipicamente dentro de uma faixa de cerca de 0,1 a 100 Hz.
Uma pluralidade de eletrodos detectores é disposta no fundo da água em locais espaçados, e os eletrodos detectores são conectados a dispositivos que registram as voltagens induzidas através de diversos pares de tais ele- trodos. Tal levantamento é conhecido como o levantamento eletromagnético de fonte controlada por domínio de freqüência.
Outra técnica de fonte controlada para o levantamento eletro- magnético de formações da Terra de subsuperfície conhecida na técnica consiste no levantamento eletromagnético de fonte controlada transiente. No levantamento eletromagnético de fonte controlada transiente, uma corrente elétrica ou um campo magnético é administrado na Terra, quando tais Ievan- tamentos são conduzidos em terra, ou é administrado em sedimentos abaixo do fundo da água (fundo do mar), quando tais levantamentos são conduzi- dos em um ambiente marinho com o uso de eletrodos em um cabo similar aos explicados acima, como usados para o levantamento por domínio de freqüência. A corrente elétrica pode ser corrente direta (CC). Em um mo- mento ou momentos selecionados, a corrente elétrica é comutada e as vol- tagens induzidas são medidas, tipicamente, em relação ao tempo ao longo de um intervalo de tempo selecionado, com o uso de eletrodos dispostos em terra, na coluna de água ou sobre o fundo da água, conforme explicado an- teriormente com referência ao levantamento por domínio de freqüência. A estrutura e a composição da subsuperfície da Terra são deduzidas por meio da distribuição de espaço e tempo das voltagens induzidas. As técnicas de levantamento t-CSEM são descritas, por exemplo, na Publicação de Pedido de Patente Internacional sob o N-. WO 2007/104949 A1, intitulada Optimiza- tion of MTEM Parameters.
Um dos parâmetros específicos determinados a partir da distri- buição de tempo de voltagens induzidas consiste na resistividade elétrica das formações de subsuperfície. Ao tornar adequadas as medições de res- posta eletromagnéticas espacialmente distribuídas, é possível gerar uma imagem tridimensional da distribuição espacial da resistividade elétrica na subsuperfície da Terra.
As técnicas conhecidas na técnica para a determinação da dis- tribuição espacial da resistividade elétrica com o uso de medições de levan- tamento eletromagnético, tipicamente presumem que a resistividade elétrica é isotrópica, isto é, a resistividade é a mesma em qualquer formação de ro- cha particular de subsuperfície, independente da direção do fluxo de corren- te elétrica usado para fazer as medições. É de conhecimento na técnica, no entanto, que a resistividade elétrica de algumas formações de rocha é aniso- trópica. A anisotropia de resistividade está presente em algumas formações de rocha em uma variedade de escalas a partir de micro (por exemplo, vari- ação de conectividade de água de poro de tamanho de grão) a macro (por exemplo, seqüências de areia-xisto laminadas). Vide, por exemplo, a patente U.S. Ne. 6.643.589, concedida a Zhang et al. e a patente U.S. Ns. 7.269.515, concedida a Tabarovsky et al. As duas patentes anteriores descrevem as técnicas para a determinação de resistividade elétrica e anisotropia de resis- tividade a partir de dentro de furos de poço perfurados através das forma- ções de rocha de subsuperfície. Contudo, tais técnicas não são aplicáveis ao uso com o levantamento eletromagnético conduzido a partir de cima das formações de rocha de interesse. Há uma necessidade por técnicas de le- vantamento eletromagnético que levem em consideração a anisotropia de resistividade.
Sumário da invenção
Um método para a determinação de anisotropia de resistividade de formações de rocha de subsuperfície, de acordo com um aspecto da in- venção, inclui administrar um campo eletromagnético transiente nas forma- ções de rocha de subsuperfície. A resposta eletromagnética das formações é medida em uma pluralidade de distâncias a partir de uma posição de ad- ministração. Para cada deslocamento, é determinada uma resposta de etapa das formações. Um momento a partir da administração pode ser selecionado de modo que o valor da resposta de etapa neste momento esteja somente substancialmente relacionado à resistividade horizontal e outro momento a partir da administração pode ser selecionado de modo que o valor da res- posta de etapa no segundo momento esteja somente substancialmente rela- cionado à resistividade média. A resistividade horizontal e a resistividade média, então, encontradas são usadas para determinar a anisotropia de re- sistividade.
Um método para a determinação de distribuição de resistividade em formações de subsuperfície, de acordo com outro aspecto da invenção, inclui o uso de medições feitas em resposta à administração de um campo eletromagnético transiente nas formações de subsuperfície. As medições são feitas em uma pluralidade de deslocamentos a partir de uma posição na qual o campo eletromagnético foi administrado. Um método, de acordo com este aspecto da invenção, inclui a determinação de uma resposta de etapa das formações. Um momento a partir da administração pode ser selecionado de modo que o valor da resposta de etapa neste momento esteja somente substancialmente relacionado à resistividade horizontal e outro momento a partir da administração pode ser selecionado de modo que o valor da res- posta de etapa neste segundo momento esteja somente substancialmente relacionado à resistividade média. A resistividade horizontal e a resistividade média, então, encontradas são usadas para determinar a anisotropia de re- sistividade.
Em uma implantação de exemplo, um modelo inicial das forma- ções de subsuperfície é gerado com o uso dos valores de anisotropia de re- sistividade e resistividade horizontal determinados. As respostas de etapa, conforme uma função de deslocamento, são calculadas para este modelo inicial e um valor de razão de anisotropia é calculado em relação ao deslo- camento com o uso dos valores de resistividade horizontal obtida em um momento selecionado a partir da administração e a resistividade média obti- da em outro momento a partir da administração. As razões de anisotropia calculadas em cada deslocamento são comparadas às determinadas a partir das respostas de etapa medidas. O modelo inicial é ajustado e o cálculo da razão de anisotropia, a razão de anisotropia medida e a comparação são repetidas até que as diferenças entre as razões de anisotropia calculadas e as razões de anisotropia medidas alcancem um mínimo ou diminuam abaixo de um limiar selecionado.
Outros aspectos e vantagens da invenção serão evidentes a par- tir da seguinte descrição e das reivindicações em anexo.
Breve descrição dos desenhos
A figura 1 mostra um sistema de exemplo para adquirir as medi- ções eletromagnéticas usadas com a invenção.
A figura 2 mostra um modelo de três camadas de resistividades de formações de rocha de subsuperfície com razões de anisotropia selecio- nadas.
A figura 3 mostra gráficos de razões de anisotropia aparentes em relação ao deslocamento para as formações de modelo mostradas na figura 2.
A figura 4 mostra uma "resposta de etapa" de exemplo de mo- mento inicial e último momento de formações de subsuperfície a um campo eletromagnético transiente.
Descrição detalhada
A figura 1 mostra um sistema de exemplo levantamento eletro- magnético marinho que pode adquirir sinais de levantamento eletromagnéti- co de fonte controlada transiente para o processamento de acordo com a invenção. O sistema pode incluir uma embarcação de levantamento 10 que se move ao longo da superfície 12A de um corpo de água 12, tal como um lago ou oceano. A embarcação 10 pode incluir equipamentos sobre a mes- ma, chamados, por conveniência, como um "sistema de registro" e mostra- dos geralmente em 14, para a geração de campos eletromagnéticos a serem administrados nas formações 24 abaixo do fundo da água 12 e para o regis- tro de medições feitas em resposta aos campos eletromagnéticos adminis- trados. O sistema de registro 14 pode incluir (nenhum mostrado separada- mente, a título de clareza da ilustração) dispositivos de navegação para de- terminar a posição geodésica da embarcação 10. A embarcação pode incluir equipamentos adicionais para a determinação da posição geodésica e/ou condução de um ou mais receptores e transmissores eletromagnéticos (des- critos abaixo), dispositivos para administrar a corrente elétrica para o(s) transmissor(es); e equipamentos de armazenamento de dados para o regis- tro de sinais detectados por um ou mais receptores eletromagnéticos.
O transmissor eletromagnético no presente exemplo pode ser um eletrodo bipolar, mostrado como um par de eletrodos em 16A, 16B dis- posto ao longo de um cabo elétrico 16 rebocado pela embarcação 10. Em momentos selecionados, o sistema de registro 14 pode passar a corrente elétrica através dos eletrodos 16A, 16B. Prefere-se que a corrente seja con- figurada para induzir os campos eletromagnéticos transientes nas formações 24 abaixo do fundo da água 12A. Os exemplos de tal corrente incluem cor- rente direta comutada, em que a corrente pode ser ligada, desligada, rever- tida a polaridade ou comutada em um conjunto estendido de eventos de co- mutação, tais como uma seqüência binária pseudoaleatória ("PRBS") ou ou- tra seqüência codificada.
No presente exemplo, a embarcação 10 pode rebocar um ou mais cabos receptores 18 que têm uma pluralidade de receptores eletro- magnéticos sobre os mesmos, tais como os eletrodos bipolares 18A, 18B, disposta em posições espaçadas ao longo do cabo. Os eletrodos bipolares 18A, 18B terão voltagens administradas através deles relacionadas à ampli- tude do componente do campo elétrico do campo eletromagnético que se origina das formações 24 em resposta ao campo eletromagnético adminis- trado. O sistema de registro 14 sobre a embarcação 10 pode incluir, confor- me explicado acima, dispositivos para o registro dos sinais gerados pelos eletrodos 18A, 18B. O registro de cada resposta do receptor é tipicamente indexado em relação a um tempo de referência, tal como um evento de co- mutação de corrente na corrente do transmissor. Um sensor 17, tal como um sensor de campo magnético (por exemplo, um magnetômetro), ou um medi- dor de corrente pode ser disposto próximo ao transmissor conforme mostra- do e pode ser usado para medir um parâmetro relacionado à quantidade de corrente que flui através do transmissor.
No presente exemplo, em substituição ou adição ao cabo recep- tor 18 rebocado pela embarcação 10, um cabo receptor de fundo da água 20 pode ser disposto ao longo do fundo da água 12 e pode incluir uma plurali- dade de receptores, tais como os eletrodos bipolares 20A, 20B similares em configuração aos eletrodos bipolares 18A, 18B sobre o cabo rebocado. Os eletrodos 20A, 20B podem estar em uma boia de registro 22 ou dispositivo similar, próximo à superfície da água 12A ou no fundo da água, que possa registrar sinais detectados pelos eletrodos 20A, 20B.
Os versados na técnica irão observar que a invenção não se limita, no escopo, às disposições de receptor e transmissor mostradas na figura 1. Outros exemplos podem usar, em substituição ou adição aos ele- trodos bipolares mostrados na figura 1, bobinas de fio ou circuitos de fio para o transmissor administrar um campo magnético de tempo variado nas for- mações 24. Os cabos receptores 18, 20 podem incluir outros dispositivos sensores, tais como magnetômetros, bobinas ou circuitos de fio para detec- tar o componente de campo magnético do campo eletromagnético induzido a partir da formação 24.
Para os propósitos de explicação da invenção, os receptores podem ser geralmente dispostos ao longo de uma linha comum ao transmis- sor durante o registro de sinal. OS registros de sinais a partir de cada um dos respectivos receptores podem ser feitos com o transmissor disposto em locais selecionados ao longo da linha comum e ativado conforme explicado acima. O sinal registrado que corresponde a cada receptor eletromagnético será associado a uma distância, chamada de "deslocamento", que é locali- zada no ponto central geodésico entre a posição geodésica do receptor e a posição geodésica do transmissor no momento do registro de sinal. Deste modo, podem ser adquiridos os sinais que correspondem a uma pluralidade de deslocamentos. O propósito para o registro de múltiplos deslocamentos, conforme se refere à invenção, será adicionalmente explicado abaixo.
Conforme explicado na seção anterior no presente documento, algumas formações podem ser eletricamente anisotrópica e, em conseqüên- cia disto, têm resistividade anisotrópica. Para os propósitos da presente in- venção, a anisotropia de resistividade será limitada ao caso de formações isotrópicas de modo vertical e transversal ("VTI"), isto é, formações que têm uma resistividade "vertical" diferente (resistividade medida com o uso do flu- xo de corrente em uma direção perpendicular aos planos da estratificação da formação) da resistividade "horizontal" (resistividade medida com o uso do fluxo de corrente em uma direção paralela aos planos de estratificação da formação). As formações VTI são consideradas tendo a mesma resistividade horizontal independente da direção azimutal, ao longo da qual a medição é feita. Tais formações também são conhecidas como tendo um eixo geomé- trico vertical de simetria.
Nas formações VTI, a resistividade vertical pv e a resistividade horizontal ph definem um "fator de anisotropia", o qual pode ser representa- do pela seguinte expressão:
<formula>formula see original document page 9</formula>
λ tem, tipicamente, um valor entre 1 e 5. A resistividade média geométrica é <formula>formula see original document page 9</formula>. Por meio da teoria de levantamento eletromagné- tico e modelagem da resposta transiente eletromagnética de um semiespaço eletricamente condutivo determinou-se que a onda de ar (resposta de etapa inicial próxima ao tempo zero a partir de um evento de comutação de corren- te transiente), representada por £(0), depende essencialmente da resistivi- dade horizontal ph. Uma resposta de etapa de último momento (em relação ao momento do evento de comutação de transmissor), a qual consiste na resposta aproximada de DC, representada por E{<x>), depende essencial- mente da resistividade média geométrica de todas as formações, através das quais o campo eletromagnético se propaga. A "resposta de etapa" con- siste na amplitude de campo magnético ou voltagem medida em resposta a uma mudança de função de etapa na corrente do transmissor, isto é, a res- posta medida para comutar a corrente e manter a corrente no valor comuta- do. A resposta de etapa consiste no total da resposta de impulso. A resposta de impulso consiste na amplitude de campo medida ou voltagem administra- da em relação ao tempo, indexada ao tempo do evento de comutação. Um exemplo de resposta de etapa de momento inicial e último momento é mos- trado no gráfico da figura 4 em 52 e 54, respectivamente.
Com o uso dos resultados descritos em: Wilson, A.J.S., 1997, The equivalent wavefield concept in multichannel transient electromagnetic surveying: tese de Ph.D., Universidade de Edinburgo, para um semiespaço isotrópico uniforme com resistividade ρ e medições obtidas conforme expli- cado acima, podem ser derivadas as seguintes expressões:
<formula>formula see original document page 10</formula>
Nas expressões acima, r representa o deslocamento. As expres- sões acima podem ser usadas para fornecer uma expressão para a determi- nação da razão de anisotropia de resistividade de um semiespaço anisotró- pico uniforme:
<formula>formula see original document page 10</formula>
A combinação das equações (2), (3) e (4) fornece uma expres- são relacionada ao deslocamento para a determinação da razão de aniso- tropia de um semiespaço anisotrópico uniforme:
<formula>formula see original document page 10</formula>
na qual Ex(r,t) consiste na resposta de etapa em linha da Terra (ao longo da linha comum explicada acima e indicada pelo χ subscrito) para o deslocamento r no momento t, r, consiste em um deslocamento curto selecionado de modo que seja mais fácil determinar a resposta de momento inicial, e r2 consiste em um deslocamento longo selecionado, onde seja mais fácil determinar a resposta de último momento. A expressão cima pode ser estendida para definir a anisotropia aparente para qualquer valor de deslo- camento r para quaisquer formações VTI de subsuperfície, conforme expôs- to a seguir:
<formula>formula see original document page 11</formula>
À medida que o deslocamento aumenta, as respostas eletro- magnéticas medidas são influenciadas pelas seções mais profundas da sub- superfície. Como conseqüência disto, Áapp(r) varia em resposta à variação de anisotropia com a profundidade. A figura 2 mostra dois modelos de 3 ca- madas de formações de subsuperfície, a título de ilustração. A profundidade das formações é ajustado à zero no fundo da água (12A na figura 1). Uma camada de formação superior é mostrada em 30 e, em ambos os modelos, tem uma resistividade horizontal igual a 10 ohm-metros e uma razão de ani- sotropia de 1,5. Deste modo, em ambos os modelos, a resistividade vertical das primeiras camadas de formação é de 15 ohm-metros, conforme mostra- do em 32 e 34. A segunda camada tem uma resistividade horizontal de 12,5 ohm-metros. No primeiro modelo, a segunda camada tem uma razão de iso- tropia de 1,8, mostrada pela resistividade vertical em 38, e no segundo mo- delo tem uma razão de anisotropia de 2,2, conforme mostrado pela resistivi- dade vertical indicada em 40. A camada mais inferior nos dois modelos tem a resistividade horizontal mostrada em 42 e razões de anisotropia de 2,5, em ambos os modelos, conforme mostrado pelos valores de resistividade verti- cal indicados em 44 e 46.
A figura 3 mostra gráficos da anisotropia aparente para cada um dos modelos acima, conforme uma função de deslocamento, determinada por meio da modelagem da resposta de etapa e com o uso da equação (6).
A função de anisotropia aparente do primeiro modelo em relação ao deslo- camento é mostrada pela curva 50 e a função de anisotropia aparente do segundo modelo é mostrada pela curva 48. É evidente que Xapp (r) está rela- cionado à anisotropia de resistividade das formações na subsuperfície. À medida que a anisotropia de resistividade aumenta, a anisotropia aparente aumenta mais rápido com o aumento do deslocamento, conforme mostrado mediante a comparação das curvas 48 e 50 na figura 3.
Em uma implantação prática de um método de acordo com a invenção, as respostas eletromagnéticas de etapa podem ser obtidas em uma pluralidade de deslocamentos, por exemplo, executados com o uso do sistema mostrado na figura 1. A resposta de etapa pode ser obtida por meio da indução de um campo eletromagnético transiente mediante a energiza- ção do transmissor, por exemplo, mediante a condução de corrente elétrica através de eletrodos (16A e 16B na figura 1) do transmissor. A corrente pode ser conforme descrito com referência à figura 1. Por exemplo, a corrente po- de ser na forma de um PRBS. As voltagens induzidas através de diversos sensores eletromagnéticos, tais como o 18A e 18B na figura 1 ou 20A e 20B na figura 1, podem ser registradas. Se a corrente do transmissor for em for- ma de um PBRS1 a forma de onda da corrente do transmissor pode ser me- dida e usada para deconvolver os sinais de voltagem registrados para obter a resposta de impulso. A resposta de impulso pode ser integrada para se obter a resposta de etapa.
Para cada resposta de etapa, e para cada deslocamento r, a amplitude da resposta de etapa no último momento [£x(r,0)] e no momento inicial [£,(r,oo)] é, então, determinada, conforme explicado acima com refe- rência à figura 3. A anisotropia aparente λαρρ(τ) pode ser, então, computada para cada deslocamento com o uso da equação (6). A anisotropia aparente pode ser usada para gerar um modelo inicial da subsuperfície. Tipicamente, o modelo inicial será um semiespaço que tem uma pluralidade de camadas de formação de rocha, nas quais cada camada tem o mesmo valor de resis- tividade horizontal determinada a partir da resposta de etapa de último mo- mento e cada camada tem o mesmo valor da resistividade média geométrica determinada a partir da resposta de momento inicial. O modelo inicial pode ser iterativamente atualizado por meio da mudança da resistividade horizon- tal e da razão de anisotropia para cada camada, as mudanças sendo deri- vadas, por exemplo, com o uso de um esquema de inversão de Occam es- tendido para incluir tanto a resistividade horizontal como a razão de aniso- tropia, como parâmetros livres da inversão. Vide, por exemplo, Constable, S. C., R. L. Parker, e C. G. Constable, 1987, Occam's Inversion: a practical al- goríthm for generating smooth models from EM sounding data, Geophysics, 52, 289-300.
Para cada iteração de modelo, uma razão de anisotropia aparen- te em relação ao deslocamento pode ser calculada com o uso da equação (6). A razão de anisotropia aparente pode ser comparada à razão de aniso- tropia em relação ao deslocamento determinado a partir da resposta de eta- pa medida das formações de subsuperfície ao campo eletromagnético admi- nistrado. As ações anteriormente mencionadas podem ser repetidas suces- sivamente até que seja gerada uma imagem final das formações de subsu- perfície. A imagem final pode ser determinada ao ser gerada quando as dife- renças entre as respostas de etapa em relação ao deslocamento determina- do a partir das medições eletromagnéticas e aquelas calculadas com o uso da equação (6), a partir do modelo ajustado, alcancem um mínimo ou dimi- nuam abaixo de um limiar selecionado.
Os métodos de acordo com a invenção podem fornecer imagens de resistividade elétrica de formações de rocha de subsuperfície que inclu- am os efeitos da anisotropia de resistividade com o uso de medições de le- vantamento eletromagnético transiente.
Embora a invenção tenha sido descrita em relação a um número limitado de modalidades, os versados na técnica, com o benefício desta descrição, irão observar que outras modalidades podem ser planejadas, as quais não se desviem do escopo da invenção, conforme apresentado no presente documento. Consequentemente, o escopo da invenção deveria ser somente limitado pelas reivindicações em anexo.
Claims (14)
1. Método para a determinação de anisotropia de resistividade de formações de rocha de subsuperfície que compreende: administrar um campo eletromagnético transiente nas formações de rocha de subsuperfície; medir a resposta eletromagnética das formações em uma plura- lidade de distâncias a partir de uma posição da administração; para cada deslocamento, determinar uma resposta de etapa das formações em um momento a partir da administração, selecionado de modo que a resposta de etapa esteja somente substancialmente relacionada à re- sistividade horizontal; para cada deslocamento, determinar uma resposta de etapa das formações em um momento a partir da administração, selecionado de modo que a resposta de etapa esteja somente substancialmente relacionada à re- sistividade média; e usar a resposta de etapa de resistividade horizontal e a resposta de etapa de resistividade média para determinar a anisotropia de resistividade.
2. Método, de acordo com a reivindicação 1, em que a adminis- tração compreende passar a corrente elétrica através de um transmissor, a corrente compreendendo ao menos uma dentre corrente ligada, corrente desligada, corrente de polaridade invertida e comutação da corrente em uma seqüência codificada.
3. Método, de acordo com a reivindicação 1, em que a determi- nação da resposta de etapa compreende a determinação da resposta de impulso e a integração da resposta de impulso.
4. Método, de acordo com a reivindicação 3, em que a determi- nação da resposta de impulso compreende a deconvolução da resposta ele- tromagnética medida com uma forma de onda de uma corrente elétrica usa- da para ser administrada ao campo eletromagnético.
5. Método, de acordo com a reivindicação 1, em que a medição da resposta eletromagnética compreende a medição de voltagens adminis- tradas através de pares de eletrodos.
6. Método, de acordo com a reivindicação 1, em que a resposta de etapa que corresponde substancialmente à resistividade horizontal com- preende a resposta eletromagnética a uma onda de ar.
7. Método, de acordo com a reivindicação 1, compreendendo adicionalmente: (a) gerar um modelo inicial das formações de subsuperfície com o uso da anisotropia de resistividade determinada, o modelo inicial incluindo um valor de resistividade horizontal e um valor de razão de anisotropia; (b) calcular um valor de razão de anisotropia em relação ao deslocamento com o uso dos valores de resistividade hori- zontal e de razão de anisotropia; (c) usar a razão de anisotropia calculada em relação ao deslo- camento para estimar a resposta de etapa das formações em um momento a partir da administração selecionada, de modo que a resposta de etapa esteja somente substanci- almente relacionada à resistividade horizontal e em um momento a partir da administração selecionada, de modo que a resposta de etapa esteja somente substancialmente relacionada à resistividade média; (d) comparar as respostas de etapa estimadas com as respos- tas de etapa determinadas; e (e) ajustar o modelo inicial e a repetição de (b), (c) e (d) até que as diferenças entre as respostas de etapa estimadas e as respostas de etapa determinadas alcancem um mínimo ou diminuam abaixo de um limiar selecionado.
8. Método para a determinação de distribuição de resistividade em formações de subsuperfície, com o uso de medições feitas em resposta à administração de um campo eletromagnético transiente nas formações de subsuperfície, em que as medições são feitas em uma pluralidade de deslo- camentos a partir de uma posição na qual o campo eletromagnético foi ad- ministrado, o método compreendendo: (a) determinar uma resposta de etapa das formações em um momento a partir da administração selecionada, de modo que a resposta de etapa esteja somente substancialmente relacionada à resistividade horizontal e uma resposta de e- tapa das formações em um momento a partir da adminis- tração selecionada, de modo que a resposta de etapa este- ja somente substancialmente relacionada à resistividade média; e (b) usar a resposta de etapa de resistividade horizontal e da resposta de etapa de resistividade média para determinar a anisotropia de resistividade;
9. Método, de acordo com a reivindicação 8, que compreende adicionalmente: (c) gerar um modelo inicial das formações de subsuperfície com o uso da anisotropia de resistividade determinada, em que o modelo inicial inclui um valor de resistividade horizon- tal e um valor de razão de anisotropia; (d) calcular um valor de razão de anisotropia em relação ao deslocamento com o uso dos valores de resistividade hori- zontal e de razão de anisotropia; (e) usar a razão de anisotropia calculada em relação ao deslo- camento para estimar a resposta de etapa das formações em um momento a partir da administração selecionada, de modo que a resposta de etapa esteja somente substanci- almente relacionada à resistividade horizontal e em um momento a partir da administração selecionado de modo que a resposta de etapa esteja somente substancialmente relacionada à resistividade média; (f) comparar as respostas de etapa estimadas com as respos- tas de etapa determinadas; e (g) ajustar o modelo inicial e a repetição de (d), (e) e (f) até que as diferenças entre as respostas de etapa estimadas e as respostas de etapa determinadas alcancem um mínimo ou diminuam abaixo de um limiar selecionado.
10. Método, de acordo com a reivindicação 8, em que o campo eletromagnético é administrado mediante a passagem da corrente elétrica através de um transmissor, em que a corrente compreende ao menos um dentre ligar a corrente, desligar a corrente, reverter a polaridade da corrente e comutar a corrente em uma seqüência codificada.
11. Método, de acordo com a reivindicação 8, em que a determi- nação da resposta de etapa compreende a determinação da resposta de impulso e a integração da resposta de impulso.
12. Método, de acordo com a reivindicação 11, em que a deter- minação da resposta de impulso compreende a deconvolução da resposta eletromagnética medida com uma forma de onda de uma corrente elétrica usada para administrar o campo eletromagnético.
13. Método, de acordo com a reivindicação 8, em que a resposta eletromagnética medida compreende as medições das voltagens adminis- tradas através dos pares de eletrodos.
14. Método, de acordo com a reivindicação 8, em que a resposta de etapa que corresponde substancialmente a resistividade horizontal com- preende a resposta eletromagnética a uma onda de ar.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US12/381,690 US20100235100A1 (en) | 2009-03-16 | 2009-03-16 | Method for determining resistivity anisotropy from earth electromagnetic responses |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| BRPI1000912A2 true BRPI1000912A2 (pt) | 2012-01-17 |
Family
ID=42246069
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| BRPI1000912-4A BRPI1000912A2 (pt) | 2009-03-16 | 2010-03-15 | método para a determinação de anisotropia de resistividade a partir de respostas eletromagnéticas da terra |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US20100235100A1 (pt) |
| EP (1) | EP2230534A1 (pt) |
| AU (1) | AU2010200948A1 (pt) |
| BR (1) | BRPI1000912A2 (pt) |
| CA (1) | CA2696071A1 (pt) |
| MX (1) | MX2010002898A (pt) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8131522B2 (en) * | 2009-06-26 | 2012-03-06 | Pgs Geophysical As | Method for estimating and removing air wave response in marine electromagnetic surveying |
| US20110012601A1 (en) * | 2009-07-15 | 2011-01-20 | Bruce Alan Hobbs | Method for determining resistivity anisotropy from earth electromagnetic tansient step response and electromagnetic transient peak impulse response |
| US20110205839A1 (en) * | 2010-02-24 | 2011-08-25 | Suedow Gustav Goeran Mattias | Method for towing marine sensor streamers |
| CN105607131B (zh) * | 2016-01-11 | 2018-03-23 | 甘肃省有色地质调查院 | 编码源电磁测深法获取大地全区视电阻率谱的方法及装置 |
| CN110968826B (zh) * | 2019-10-11 | 2023-04-07 | 重庆大学 | 一种基于空间映射技术的大地电磁深度神经网络反演方法 |
| CN113933905B (zh) * | 2021-09-30 | 2023-09-19 | 中国矿业大学 | 一种圆锥型场源瞬变电磁反演方法 |
| CN114002747B (zh) * | 2021-11-15 | 2023-10-31 | 南方科技大学 | 地面发射空中接收的地质体电阻率探测方法及系统 |
| CN114966870B (zh) * | 2022-06-01 | 2024-03-12 | 广东省地质物探工程勘察院 | 一种地面任意回线任意位置多分量联合探测瞬变电磁方法 |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2288988A1 (fr) * | 1974-07-30 | 1976-05-21 | Duroux Jean | Procede et appareil de prospection en mer par mesure de champs electromagnetiques |
| US4041372A (en) * | 1975-09-08 | 1977-08-09 | Continental Oil Company | Apparatus for multi-channel induced polarization surveying |
| CA1036221A (en) * | 1976-06-02 | 1978-08-08 | Geonics Limited | Method and apparatus for measuring terrain resistivity |
| GB1588495A (en) * | 1978-05-19 | 1981-04-23 | Shell Int Research | Method and means for waterbottom logging |
| US4617518A (en) * | 1983-11-21 | 1986-10-14 | Exxon Production Research Co. | Method and apparatus for offshore electromagnetic sounding utilizing wavelength effects to determine optimum source and detector positions |
| MY131017A (en) * | 1999-09-15 | 2007-07-31 | Exxonmobil Upstream Res Co | Remote reservoir resistivity mapping |
| US6643589B2 (en) * | 2001-03-08 | 2003-11-04 | Baker Hughes Incorporated | Simultaneous determination of formation angles and anisotropic resistivity using multi-component induction logging data |
| GB0121719D0 (en) * | 2001-09-07 | 2001-10-31 | Univ Edinburgh | Method for detection fo subsurface resistivity contrasts |
| GB2382875B (en) * | 2001-12-07 | 2004-03-03 | Univ Southampton | Electromagnetic surveying for hydrocarbon reservoirs |
| WO2004029664A1 (en) * | 2002-09-27 | 2004-04-08 | Baker Hughes Incorporated | A method for resistivity anisotropy determination in near vertical wells |
| US7023213B2 (en) * | 2002-12-10 | 2006-04-04 | Schlumberger Technology Corporation | Subsurface conductivity imaging systems and methods |
| US7132831B2 (en) * | 2004-03-31 | 2006-11-07 | Peteralv Brabers | Electrode configuration for resistivity sounding |
| US7269515B2 (en) * | 2004-06-15 | 2007-09-11 | Baker Hughes Incorporated | Geosteering in anisotropic formations using multicomponent induction measurements |
| US20060186887A1 (en) * | 2005-02-22 | 2006-08-24 | Strack Kurt M | Method for identifying subsurface features from marine transient controlled source electromagnetic surveys |
| GB0604829D0 (en) | 2006-03-10 | 2006-04-19 | Mtem Ltd | Optimisation of mtem parameters |
| US7953587B2 (en) * | 2006-06-15 | 2011-05-31 | Schlumberger Technology Corp | Method for designing and optimizing drilling and completion operations in hydrocarbon reservoirs |
| US7474101B2 (en) * | 2006-09-12 | 2009-01-06 | Kjt Enterprises, Inc. | Method for combined transient and frequency domain electromagnetic measurements |
-
2009
- 2009-03-16 US US12/381,690 patent/US20100235100A1/en not_active Abandoned
-
2010
- 2010-03-08 EP EP10155816A patent/EP2230534A1/en not_active Withdrawn
- 2010-03-09 CA CA2696071A patent/CA2696071A1/en not_active Abandoned
- 2010-03-11 AU AU2010200948A patent/AU2010200948A1/en not_active Abandoned
- 2010-03-12 MX MX2010002898A patent/MX2010002898A/es not_active Application Discontinuation
- 2010-03-15 BR BRPI1000912-4A patent/BRPI1000912A2/pt not_active Application Discontinuation
Also Published As
| Publication number | Publication date |
|---|---|
| EP2230534A1 (en) | 2010-09-22 |
| US20100235100A1 (en) | 2010-09-16 |
| MX2010002898A (es) | 2010-09-30 |
| AU2010200948A1 (en) | 2010-09-30 |
| CA2696071A1 (en) | 2010-09-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| BRPI1000912A2 (pt) | método para a determinação de anisotropia de resistividade a partir de respostas eletromagnéticas da terra | |
| BRPI1002669A2 (pt) | mÉtodo para determinar a anisotropia de resistividade a partir da resposta transitària eletromagnÉtica ao degrau da terra e de resposta transitària eletromagnÉtica ao impulso de pico | |
| US8437961B2 (en) | Time lapse analysis with electromagnetic data | |
| AU2008237594B2 (en) | Method and apparatus for reducing induction noise in measurements with a towed electromagnetic survey system | |
| BRPI0915113A2 (pt) | mpétodo para inspeção eletromagnética de subsuperfície usando duas ou mais fontes eletromagnéticas acionadas simultaneamente | |
| BRPI1000052B1 (pt) | método para exploração eletromagnética | |
| BRPI1000323A2 (pt) | levantamento geofìsico eletromagnético transiente de pequeno deslocamento | |
| BRPI0915950B1 (pt) | Método para reduzir o efeito de movimento sobre os sinais eletromagnéticos detectados durante o movimento de um receptor eletromagnético | |
| BRPI1002287B1 (pt) | método para a determinação da distribuição de resistividade das formações abaixo do fundo da massa de água. | |
| BRPI0611352A2 (pt) | método implementado por computador para determinar anisotropia elétrica vertical da terra | |
| BRPI0715832A2 (pt) | aperfeiÇoamento em exploraÇço marinha | |
| US8896314B2 (en) | Multiple component electromagnetic survey signal acquisition method | |
| CN101484897A (zh) | 获取和解释瞬态电磁测量结果的方法 | |
| BR102012033262B1 (pt) | método de levantamento geofísico eletromagnético e sistema de levantamento eletromagnético marinho | |
| BRPI0618149A2 (pt) | métodos para estimar erro em dados medidos provenientes de uma prospecção eletromagnética de fonte controlada de uma região de sub-superfìcie e para produzir hidrocarbonetos a partir de uma região de sub-superfìcie | |
| BRPI0614372A2 (pt) | medições de resistividade cc de múltiplos transientes | |
| EP2149058A1 (en) | Multi-component marine electromagnetic signal acquisition cable, system and method | |
| BR102012030628A2 (pt) | Métodos e aparelho para uma rápida determinação da profundidade e da resistência transversal alvo | |
| BRPI0917820B1 (pt) | Método e sistema para aumentar a sensibilidade em medição de magnitude e direção de resistividade | |
| RU2174243C1 (ru) | Способ геоэлектроразведки | |
| RU2690207C1 (ru) | Способ магнитотеллурического зондирования геологических структур | |
| SGSGGSGSGSGGSGSGGSGSGGSGSGGSGSGSSS | Strack et al. | |
| BRPI0711282B1 (pt) | Métodos para determinar mudanças dependentes de tempo no teor de hidrocarboneto de um reservatório de sub- superfície |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| B03A | Publication of a patent application or of a certificate of addition of invention [chapter 3.1 patent gazette] | ||
| B11A | Dismissal acc. art.33 of ipl - examination not requested within 36 months of filing | ||
| B11Y | Definitive dismissal - extension of time limit for request of examination expired [chapter 11.1.1 patent gazette] |