Stratum component optimization determination method and device

Patent metadata

Field	Value
Publication number	US-11010507-B2
Application number	US-201615749290-A
Country	US
Kind code	B2
Filing date	Oct 10, 2016
Priority date	Nov 18, 2015
Publication date	May 18, 2021
Grant date	May 18, 2021

Abstract

The present invention provides a stratum component optimization determination method and device, which fall within the technical field of oil-gas exploration well logging. The method comprises: establishing a stratum rock component model according to core analysis data and geological conditions of a stratum to be detected, and determining a well logging curve determined by a participation model; determining a well logging response equation expression corresponding to the well logging curve determined by the participation model; parsing, recording and storing the well logging response equation expression, establishing a target function of an optimization problem, and solving the target function through an iteration algorithm to determine an optimal component content of the stratum to be detected. By establishing the stratum rock component model, determining corresponding well logging response equation, and parsing through an expression parsing method and recording and storing the well logging response equation expression, and then, establishing the target function of the optimization problem, and obtaining the optimal component content of the stratum to be detected through the iteration algorithm, the present invention can not only optimizes the self-defined well logging response equation expression of the user, but achieves a high processing precision.

First claim

The invention claimed is:

1. A stratum component optimization determining method, wherein comprising:
establishing a stratum rock component model according to core analysis data and geological conditions of a stratum to be studied, and determining well logging curves participating in the determination of the model;
determining an expression of a well logging response equation corresponding to the well logging curves participating in the determination of the model;
parsing, recording and storing the expression of the well logging response equation;
establishing a target function of an optimization problem according to the parsed response equation expression, and solving the target function through an iteration algorithm to determine an optimal component content of the stratum to be studied, wherein the target function of the optimization problem is:

v
*

=

arg
⁢

⁢
min
⁢

{

F
⁡

(

v
→

)

}

F
⁡

(

v
→

)

=

1
2

⁢

∑

K
=
1

n

⁢

[

(

t
ci

-

t
mi

)

·

w
i

]

2

wherein F({right arrow over (v)}) denotes the target function; v* denotes a value that makes the target function a minimum value; t ci denotes a value of response equation each well logging method determined according to the well logging response equation; t mi denotes an actual well logging measured response value; w i denotes a weight coefficient of the well logging curve in an optimization model; and n denotes the number of well logging curves solved.

2. The method according to claim 1 , wherein each well logging curve comprises at least one of a natural gamma well logging curve, a deep lateral resistivity well logging curve, a shallow lateral resistivity well logging curve, a density well logging curve, a neutron well logging curve or an element capture energy spectrum well logging curve.

3. The method according to claim 1 , wherein the parsing, recording and storing the well logging response equation expression comprises:
parsing the expression of the well logging response equation, and translating the expression of the well logging response equation into a postfix expression that can be used for calculation in a computer;
recording and storing each element in the postfix expression using a dynamic array storage structure;
traversing the storage structures of the postfix expression, determining combination relations between each element in a derivation rule, and determining partial derivative forms of the postfix expression.

4. The method according to claim 2 , wherein the well logging response equation is:
t ci =f i ( {right arrow over (v)} )
wherein t ci denotes a value of the well logging response equation; {right arrow over (v)} denotes the content of the mineral and fluid component in the stratum to be studied; f i ({right arrow over (v)}) denotes an expression form of the well logging response equation, including variables, constants, operators and parameter symbols of stratum mineral and fluid component.

5. The method according to claim 1 , wherein the stratum rock component model further comprises a constraint condition comprising:
h k ( {right arrow over (v)} )= C k ·{right arrow over (v)}−b k ≤0

wherein, h k ({right arrow over (v)}) denotes a constraint condition; C k denotes a coefficient matrix of the constraint condition; {right arrow over (v)} denotes the contents of stratum mineral and fluid component in the stratum to be studied; b k denotes a constraint condition boundary.

6. The method according to claim 5 , wherein solving the target function through an iteration algorithm comprises:
converting the target function into an expression of an unconstrained problem using a penalty function method:
min Φ( {right arrow over (v)},M )= F ( {right arrow over (v)} )+ M·P ( {right arrow over (v)} )

wherein M denotes a penalty factor;

P
⁡

(

v
→

)

=

∑

K
=
1

cn

⁢

[

⁢

max
(

0
,

h
k

⁡

(

v
→

)

]

2

denotes a penalty function, wherein when {right arrow over (v)} satisfies h k ({right arrow over (v)})≤0, a penalty term M·P({right arrow over (v)})=0, and when {right arrow over (v)} does not satisfy h k ({right arrow over (v)})≤0, the penalty term M·P({right arrow over (v)})>0 and increases with the increase of M; cn denotes the number of the constraint conditions.

Assignees

• Petrochina Co Ltd

Inventors

• Feng Zhou
• Li Ning
• Wu Hongliang
• Wang Huafeng
• Feng Qingfu
• Wang Kewen

Classifications

• E21B41/00

  Equipment or details not covered by groups E21B15/00 - E21B40/00 · CPC title

• G06F30/00 Primary

  Computer-aided design [CAD] · CPC title

• G06F30/20 Primary

  Design optimisation, verification or simulation (optimisation, verification or simulation of circuit designs G06F30/30) · CPC title

• G06Q50/00

  Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism (healthcare informatics G16H) · CPC title

• E21B49/00

  Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells · CPC title