# US 5,703,773 A · Real-time 2-dimensional inversion process and its application to induction resistivity well logging > Agents: this data is available as JSON from the patentgrep API and MCP server — see https://patentgrep.com/llms.txt - Number: US 5,703,773 A - Document ID: US5703773A - Status: Granted - Kind: A - Application number: 08/598,443 - Filed: 1996-02-08 - Granted: 1997-12-30 - Earliest claimed: 1996-02-08 - Assignees: Western Atlas International, Inc. - Inventors: Leonty Abraham Tabarovsky; Michael Boris Rabinovich - Web page: https://patentgrep.com/patent/US5703773A ## CPC - G01V 3/38: Processing data, e.g. for analysis, for interpretation, for correction - G01V 3/28: using induction coils - G01V: GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS (means for indicating the location of accidentally buried, e.g. snow-buried, persons A63B29/02) ## Abstract The invention is a method of inversion processing signals from an induction well logging instrument including a transmitter and a plurality of receivers at spaced apart locations. The method includes skin effect correcting the responses of the receivers by extrapolating the receiver responses to zero frequency. A model is generated of the media surrounding said instrument. Conductivities of elements in the model are then adjusted so that a measure of misfit between the skin-effect corrected receiver responses and simulated receiver responses based on the model is minimized. The geometry of the model is then adjusted so that the measure of misfit between the skin-effect corrected receiver responses and the simulated receiver responses based on the model is further minimized. In a preferred embodiment of the invention, the step of adjusting the geometry includes minimizing the measure of misfit between the simulated responses and the receiver responses from selected ones of the receivers closely spaced to the transmitter. Numbers of and positions of radial boundaries are then determined by minimizing the measure of misfit for all the receiver responses. ## Claims ### Claim 1 A method of inversion processing signals from an induction well logging instrument having a transmitter and a plurality of receivers at axially spaced apart locations, comprising: skin effect correcting responses of said receivers by extrapolating said receiver responses to zero frequency; generating a model of media surrounding said instrument; adjusting conductivities of elements of said model so that a measure of misfit between said skin-effect corrected receiver responses and simulated receiver responses based on said model is minimized; and adjusting the geometry of said model so that said measure of misfit between said skin-effect corrected receiver responses and said simulated receiver responses based on said model is further minimized. ### Claim 2 (depends on 1) The method as defined in claim 1 wherein said step of adjusting conductivities comprises linear optimization. ### Claim 3 (depends on 1) The method as defined in claim 1 wherein said step of adjusting the geometry of said model comprises non-linear optimization. ### Claim 4 (depends on 3) The method as defined in claim 3 wherein said non-linear optimization comprises determining numbers of axial elements and positions of said axial elements in said model by minimizing said measure of misfit with respect said numbers and said positions of said axial elements using primarily said skin effect corrected responses from ones of said plurality of receivers closely spaced to said transmitter. ### Claim 5 (depends on 4) The method as defined in claim 4 further comprising repeating said steps of generating said model, adjusting said conductivities and adjusting said geometry, wherein a minimum selected axial spacing between said axial elements forming said geometry is reduced compared to said minimum selected axial spacing used in the prior performance of said step of generating said model. ### Claim 6 (depends on 5) The method as defined in claim 5 further comprising optimizing said geometry by selectively removing ones of said axial elements and recalculating said measure of misfit, whereby said geometry is determined to be optimal when none of said ones of said axial elements can be removed without substantially increasing said measure of misfit. ### Claim 7 (depends on 1) The method as defined in claim 1 further comprising determining numbers of radial boundaries and radial positions of said radial boundaries in said model by minimizing said measure of misfit with respect said numbers and said radial positions of said radial boundaries using said skin effect corrected receiver responses from substantially all of said receivers. ### Claim 8 (depends on 1) The method as defined in claim 1 wherein said step of skin effect correcting includes solving a system of fractional polynomial equations with respect to frequency of an alternating current energizing said transmitter. ### Claim 9 (depends on 8) The method as defined in claim 8 wherein said step of solving said system includes determining a solution to the equation: ##EQU16## for the variable S.sub.0. ### Claim 10 (depends on 1) The method as defined in claim 1 wherein said measure of misfit comprises a least squares norm. ### Claim 11 (depends on 1) The method as defined in claim 1 wherein said transmitter is energized with alternating current having a plurality of component frequencies. ### Claim 12 (depends on 11) The method as defined in claim 11 wherein said alternating current comprises square waves. ### Claim 13 A method of determining resistivities of earth formations penetrated by a wellbore, comprising: energizing said earth formations by passing alternating current having a plurality of frequencies through a transmitter on an instrument disposed in said wellbore; receiving voltages induced by eddy currents flowing in said wellbore and said earth formations in receivers disposed on said instrument at spaced apart locations from said transmitter; calculating apparent conductivities from said measured voltages; moving said instrument through said wellbore and repeating said steps of energizing and receiving; skin effect correcting said apparent conductivities by extrapolating said apparent conductivities to values which would occur at zero frequency; generating a model of said earth formations surrounding said instrument; adjusting conductivities of elements of said model so that said measure of misfit between said skin-effect corrected apparent conductivities and simulated receiver responses based on said model is minimized; and adjusting geometry of said model so that a measure of misfit between said skin-effect corrected apparent conductivities and simulated receiver responses based on said model is further minimized. ### Claim 14 (depends on 13) The method as defined in claim 13 wherein said step of adjusting geometry of said model comprises determining numbers of axial elements and positions of said axial elements in said model by minimizing said measure of misfit with respect to said numbers of and said positions of said axial elements using primarily said skin effect corrected apparent conductivities from ones of said plurality of receivers closely spaced to said transmitter. ### Claim 15 (depends on 14) The method as defined in claim 14 further comprising repeating said steps of generating said model, adjusting said conductivities and adjusting said geometry, wherein a minimum selected axial spacing between said axial elements forming said geometry is reduced compared to said minimum selected axial spacing used in the prior performance of said step of generating said model. ### Claim 16 (depends on 15) The method as defined in claim 15 further comprising optimizing said geometry by selectively removing ones of said axial elements and recalculating said measure of misfit, whereby said geometry is determined to be optimal when none of said ones of said axial elements can be removed without substantially increasing said measure of misfit. ### Claim 17 (depends on 14) The method as defined in claim 14 further comprising determining numbers of radial boundaries and radial positions of said radial boundaries in said model by minimizing said measure of misfit with respect to said numbers and said radial positions of said radial boundaries using said skin effect corrected apparent conductivities from substantially all of said receivers. ### Claim 18 (depends on 13) The method as defined in claim 13 wherein said step of skin effect correcting includes solving a system of fractional polynomial equations with respect to frequency of an alternating current energizing said transmitter. ### Claim 19 (depends on 18) The method as defined in claim 18 wherein said step of solving said system includes determining a solution to the equation: ##EQU17## for the variable S.sub.0. ### Claim 20 (depends on 13) The method as defined in claim 13 wherein said measure of misfit comprises a least squares norm. ### Claim 21 (depends on 13) The method as defined in claim 13 wherein said transmitter is energized with alternating current having a plurality of component frequencies. ### Claim 22 (depends on 21) The method as defined in claim 21 wherein said alternating current comprises square waves. ### Claim 23 A method of inversion processing signals from an induction well logging instrument having a transmitter and a plurality of receivers at axially spaced apart locations, comprising: skin effect correcting responses of said receivers by extrapolating said receiver responses to zero frequency; generating a model of media surrounding said instrument, said step of generating including subdividing said receiver responses into windows, a first one of said windows initially set to represent a half space, a second one of said windows including a selected number of axial elements, a third one of said windows including a single axial element; adjusting conductivities of elements of said model so that a measure of misfit between said skin-effect corrected receiver responses and simulated receiver responses based on said model is minimized; adjusting geometry of said model so that said measure of misfit between said skin-effect corrected receiver responses and said simulated receiver responses based on said model is minimized; moving said windows to different ones of said receiver responses and repeating said steps of generating a model and adjusting said geometry so that said number of axial elements within said third one of said windows becomes positioned within said second one of said windows, said steps of moving, generating and adjusting repeated until a selected measurement interval is processed; and adjusting conductivities of elements of said adjusted model such that all said conductivities are non-negative. ### Claim 24 (depends on 23) The method as defined in claim 23 wherein said step of adjusting geometry of said model comprises determining numbers of axial elements and positions of said axial elements in said model by minimizing said measure of misfit with respect to said numbers of and said positions of said axial elements using primarily said skin effect corrected responses from ones of said plurality of receivers closely spaced to said transmitter. ### Claim 25 (depends on 24) The method as defined in claim 24 further comprising repeating said steps of generating said model, adjusting said geometry and adjusting said conductivities, wherein a minimum selected axial spacing between said axial elements forming said geometry is reduced compared to said minimum axial spacing used in the prior performance of said step of generating said model. ### Claim 26 (depends on 24) The method as defined in claim 24 further comprising optimizing said geometry by selectively removing ones of said axial elements and recalculating said measure of misfit, whereby said geometry is determined to be optimal when none of said ones of said axial elements can be removed without substantially increasing said measure of misfit. ### Claim 27 (depends on 24) The method as defined in claim 24 further comprising determining numbers of radial boundaries and radial positions of said radial boundaries in said model by minimizing said measure of misfit with respect to said numbers and said radial positions of said radial boundaries using said skin effect corrected receiver responses from substantially all of said receivers. ### Claim 28 (depends on 23) The method as defined in claim 23 wherein said step of skin effect correcting includes solving a system of fractional polynomial equations with respect to frequency of an alternating current energizing said transmitter. ### Claim 29 (depends on 28) The method as defined in claim 28 wherein said step of solving said system includes determining a solution to the equation: ##EQU18## for the variable S.sub.0. ### Claim 30 (depends on 23) The method as defined in claim 23 wherein said measure of misfit comprises a least squares norm. ### Claim 31 (depends on 23) The method as defined in claim 23 wherein said transmitter is energized with alternating current having a plurality of component frequencies. ### Claim 32 (depends on 31) The method as defined in claim 31 wherein said alternating current comprises square waves. ### Claim 33 (depends on 13, 23) The method as defined in claim 13 or 23 wherein said step of adjusting conductivities comprises linear optimization. ### Claim 34 (depends on 13, 23) The method as defined in claims 13 or 23 wherein said step of adjusting geometry comprises non-linear optimization. ### Claim 35 A method of skin effect correcting receiver signals from an induction well logging tool having a transmitter and a receiver, comprising: determining response of said receiver at a each one of a plurality of different frequencies, said different frequencies corresponding to frequencies of alternating current energizing said transmitter; extrapolating said response to a response which would obtain when said alternating current frequency is substantially equal to zero. ### Claim 36 (depends on 35) The method as defined in claim 35 wherein said step of extrapolating includes solving a system of fractional polynomial equations with respect to said frequency. ### Claim 37 (depends on 36) The method as defined in claim 36 wherein said step of solving said system includes determining a solution to the equation: ##EQU19## for the variable S.sub.0.