阅读说明：本技术 地层电阻率方位各向异性确定方法及装置 (Method and device for determining azimuthal anisotropy of formation resistivity ) 是由 李潮流 胡海川 冯周 宋连腾 李霞 胡法龙 宁从前 袁超 于 2019-10-23 设计创作，主要内容包括：本申请提供的地层电阻率方位各向异性确定方法及装置,利用多对极板上的电极微元的电流和电压测出的电阻率来确定电各向异性,由于多对极板具有设定的方位角,因此测出的电各向异性还进一步包含了方位信息,各向异性数据更加准确和全面,提供了一种新的各向异性确定方法,解决了目前研究各向异性的手段单一等问题,并能够获得并非通过声波反应的各向异性参数。(According to the method and the device for determining the formation resistivity and the azimuth anisotropy, the electrical anisotropy is determined by the resistivity measured by the current and the voltage of the electrode infinitesimal on the multiple pairs of polar plates, the multiple pairs of polar plates have set azimuth angles, so that the measured electrical anisotropy further comprises azimuth information, and the anisotropy data is more accurate and comprehensive.)

1. A method for determining azimuthal anisotropy of formation resistivity, comprising:

arranging a plurality of pairs of polar plates in a well logging borehole at intervals of a set angle, wherein two polar plates in one pair of polar plates are arranged in a central symmetry manner, a plurality of electrode microelements are arranged on each polar plate, and each polar plate is at the same set height;

applying a voltage to each polar plate to obtain the current of all electrode infinitesimal on each polar plate;

and determining the azimuthal anisotropy of the resistivity of the formation according to the currents and the corresponding voltages of all the electrode microelements on each pair of polar plates and the azimuth of each pair of polar plates.

2. The method of determining azimuthal anisotropy of resistivity of earth formations according to claim 1, characterized in that the voltages applied to each plate are the same or different.

3. The method of claim 1, wherein each electrode infinitesimal is arranged in an array on the plate on which it is placed.

4. The method of determining azimuthal anisotropy of resistivity of earth formations according to claim 1, wherein the determining azimuthal anisotropy of resistivity of earth formations from currents and corresponding voltages of all electrode microelements on each pair of pads and the orientation of each pair of pads comprises:

generating a micro-resistivity scanning image according to the current and the corresponding voltage of all electrode micro-elements on each pair of polar plates, wherein the micro-resistivity scanning image can represent resistivity data corresponding to the electrode micro-elements on each pair of polar plates;

counting resistivity data in each microresistivity scanning image, setting a distribution interval according to the counted resistivity data, and correspondingly generating a distribution frequency histogram;

taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as a characteristic value of a corresponding pair of polar plates, and sequencing each characteristic value to obtain a characteristic value sequence;

and solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence, and generating an anisotropy coefficient of the resistivity corresponding to the azimuth angle of the corresponding polar plate at the set height according to the azimuth angle of the corresponding polar plate of the first characteristic value and the last characteristic value.

5. The method of claim 3, wherein the plurality of columns are formed by the arrangement of the electrode elements on each plate, and each column is spaced apart by a predetermined distance.

6. The method of determining azimuthal anisotropy of resistivity of earth formations according to claim 5, wherein the determining azimuthal anisotropy of resistivity of earth formations from currents and corresponding voltages of all electrode microelements on each pair of pads and the orientation of each pair of pads comprises:

generating a microresistivity scanned image according to the current and the corresponding voltage of the electrode microelements of each column, wherein the microresistivity scanned image can represent resistivity data corresponding to the electrode microelements in each column;

counting the resistivity data of each column, setting a distribution interval according to the counted resistivity data, and correspondingly generating a distribution frequency histogram;

taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as a characteristic value corresponding to one column, and sequencing each characteristic value to obtain a characteristic value sequence;

and solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence, and generating an anisotropy coefficient of the resistivity corresponding to the azimuth angle of the corresponding column corresponding to the set height according to the azimuth angle of the corresponding column of the first characteristic value and the last characteristic value.

7. The method of determining azimuthal anisotropy of resistivity of earth formations according to claim 4, wherein determining azimuthal anisotropy of resistivity of earth formations from currents and corresponding voltages of all electrode microelements on each pair of pads and the orientation of each pair of pads further comprises:

and calibrating the azimuth angle of each pair of polar plates.

8. The method of determining azimuthal anisotropy of resistivity of earth formations according to claim 6, wherein determining azimuthal anisotropy of resistivity of earth formations from currents and corresponding voltages of all electrode microelements on each pair of pads and the orientation of each pair of pads further comprises:

the azimuth angle of each column is calibrated.

9. The method of determining azimuthal anisotropy of resistivity of earth formations according to claim 4 or claim 6, wherein determining azimuthal anisotropy of resistivity of earth formations from currents and corresponding voltages of all electrode microelements on each pair of pads and the orientation of each pair of pads further comprises:

and performing median filtering processing on all the resistivity data.

10. An apparatus for determining azimuthal anisotropy of formation resistivity, comprising:

the polar plate setting module is used for setting a plurality of pairs of polar plates in a well logging borehole at intervals of a set angle, wherein two polar plates in one pair of polar plates are arranged in a central symmetry manner, a plurality of electrode microelements are arranged on each polar plate, and each polar plate is at the same set height;

the current acquisition module is used for applying a voltage to each polar plate to acquire the current of all electrode microelements on each polar plate;

and the anisotropy determining module is used for determining the formation resistivity azimuth anisotropy according to the currents and the corresponding voltages of all the electrode microelements on each pair of polar plates and the azimuth of each pair of polar plates.

11. The apparatus of claim 10, wherein the voltage applied to each plate is the same or different.

12. The apparatus of claim 10, wherein each electrode element is arranged in an array on the plate on which it is placed.

13. The formation resistivity azimuthal anisotropy determination apparatus of claim 11, wherein the anisotropy determination module comprises:

the resistivity data generating unit is used for generating a micro-resistivity scanning image according to the current and the corresponding voltage of all electrode micro-elements on each pair of polar plates, and the micro-resistivity scanning image can represent the resistivity data corresponding to the electrode micro-elements on each pair of polar plates;

the distribution frequency histogram generating unit is used for counting the resistivity data in each microresistivity scanning image, setting a distribution interval according to the counted resistivity data and correspondingly generating a distribution frequency histogram;

the characteristic value sequence generating unit is used for sequencing each characteristic value by taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as the characteristic value of a corresponding pair of polar plates to obtain a characteristic value sequence;

and the anisotropic coefficient generating unit is used for solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence, and generating the anisotropic coefficient of the resistivity corresponding to the azimuth angle corresponding to the set height according to the azimuth angle of the polar plate corresponding to the first characteristic value and the last characteristic value.

14. The apparatus of claim 12, wherein the electrode elements on each plate are arranged to form a plurality of columns, each column being spaced apart by a predetermined distance.

15. The formation resistivity azimuthal anisotropy determination apparatus of claim 14, wherein the anisotropy determination module comprises:

the resistivity data generating unit is used for generating a micro-resistivity scanning image according to the current and the corresponding voltage of the electrode infinitesimal in each column, and the micro-resistivity scanning image can represent the resistivity data corresponding to the electrode infinitesimal in each column;

a distribution frequency histogram generating unit for counting the resistivity data of each column, setting a distribution interval according to the counted resistivity data, and correspondingly generating a distribution frequency histogram;

the characteristic value sequence generating unit is used for sequencing each characteristic value by taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as the characteristic value corresponding to one column to obtain a characteristic value sequence;

and the anisotropic coefficient generating unit is used for solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence and generating the anisotropic coefficient of the resistivity corresponding to the azimuth angle of the set height according to the azimuth angle of the corresponding column of the first characteristic value and the last characteristic value.

16. The formation resistivity azimuthal anisotropy determination apparatus of claim 13, wherein the anisotropy determination module further comprises:

and the azimuth angle calibration unit is used for calibrating the azimuth angle of each pair of polar plates.

17. The formation resistivity azimuthal anisotropy determination apparatus of claim 15, wherein the anisotropy determination module further comprises:

and the azimuth angle calibration unit is used for calibrating the azimuth angle of each column.

18. The formation resistivity azimuthal anisotropy determination apparatus of claim 13 or 15, wherein the anisotropy determination module further comprises:

and the median filtering processing unit is used for performing median filtering processing on all the resistivity data.

19. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any of claims 1 to 9 when executing the program.

20. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the method of any one of claims 1 to 9.

Technical Field

The invention relates to the technical field of exploration of oil fields, in particular to a method and a device for determining formation resistivity azimuth anisotropy.

Background

In oil exploration, well logging interpretation work needs to provide mechanical elastic parameters required by engineering reconstruction besides finding and evaluating hydrocarbon reservoirs, wherein the mechanical elastic parameters comprise information such as horizontal principal stress, fracture pressure and anisotropy, and the like, so that a basis is provided for designing a hydraulic fracturing reconstruction scheme.

The anisotropy means that the numerical values of the acoustic wave speed and the resistivity of the rock have directional characteristics. The research on the anisotropy of the rock is of great significance, because the anisotropy mainly comes from the discontinuous structure of the rock, macroscopically comprises joints, cracks, bedding and the like, and microscopically comprises the directional arrangement of the particles of the rock or micro-fractures of various grades. Therefore, the degree of the crack/joint/bedding development of the rock can be indirectly judged by analyzing the strength and the cause of the anisotropy of the rock, and the extension direction of the hydraulic fracturing artificial seam is estimated, so that the construction scheme is optimized, and the fracturing effect is improved.

The traditional method for researching the anisotropy is to extract and analyze the information of the fast and slow transverse waves by special software through the cross dipole measurement of the array sound waves so as to judge the anisotropy degree and direction of the array sound waves. The current means for studying the anisotropy are single, and only describe the acoustic anisotropy based on the array sound wave, and the anisotropy of other physical parameters cannot be characterized.

Disclosure of Invention

In order to solve at least one of the above disadvantages, an embodiment of a first aspect of the present application provides a method for determining azimuthal anisotropy of formation resistivity, including:

arranging a plurality of pairs of polar plates in a well logging borehole at intervals of a set angle, wherein two polar plates in one pair of polar plates are arranged in a central symmetry manner, a plurality of electrode microelements are arranged on each polar plate, and each polar plate is at the same set height;

applying a voltage to each polar plate to obtain the current of all electrode infinitesimal on each polar plate;

and determining the azimuthal anisotropy of the resistivity of the formation according to the currents and the corresponding voltages of all the electrode microelements on each pair of polar plates and the azimuth of each pair of polar plates.

In some embodiments, the voltage applied to each plate is the same or different.

In some embodiments, each electrode element is arranged in an array on the plate on which it is located.

In some embodiments, determining the azimuthal anisotropy of the formation resistivity from the currents and corresponding voltages of all the electrode microelements on each pair of pads and the orientation of each pair of pads comprises:

generating a micro-resistivity scanning image according to the current and the corresponding voltage of all electrode micro-elements on each pair of polar plates, wherein the micro-resistivity scanning image can represent resistivity data corresponding to the electrode micro-elements on each pair of polar plates;

counting resistivity data in each microresistivity scanning image, setting a distribution interval according to the counted resistivity data, and correspondingly generating a distribution frequency histogram;

taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as a characteristic value of a corresponding pair of polar plates, and sequencing each characteristic value to obtain a characteristic value sequence;

and solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence, and generating an anisotropy coefficient of the resistivity corresponding to the azimuth angle of the corresponding polar plate at the set height according to the azimuth angle of the corresponding polar plate of the first characteristic value and the last characteristic value.

In some embodiments, the electrode elements on each plate are arranged to form a plurality of columns, each column being spaced apart by a set distance.

In some embodiments, determining the azimuthal anisotropy of the formation resistivity from the currents and corresponding voltages of all the electrode microelements on each pair of pads and the orientation of each pair of pads comprises:

generating a microresistivity scanned image according to the current and the corresponding voltage of the electrode microelements of each column, wherein the microresistivity scanned image can represent resistivity data corresponding to the electrode microelements in each column;

counting the resistivity data of each column, setting a distribution interval according to the counted resistivity data, and correspondingly generating a distribution frequency histogram;

taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as a characteristic value corresponding to one column, and sequencing each characteristic value to obtain a characteristic value sequence;

and solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence, and generating an anisotropy coefficient of the resistivity corresponding to the azimuth angle of the corresponding column corresponding to the set height according to the azimuth angle of the corresponding column of the first characteristic value and the last characteristic value.

In some embodiments, determining the azimuthal anisotropy of the formation resistivity based on the currents and corresponding voltages of all the electrode microelements on each pair of pads and the orientation of each pair of pads further comprises:

and calibrating the azimuth angle of each pair of polar plates.

In some embodiments, determining the azimuthal anisotropy of the formation resistivity based on the currents and corresponding voltages of all the electrode microelements on each pair of pads and the orientation of each pair of pads further comprises:

the azimuth angle of each column is calibrated.

In some embodiments, determining the azimuthal anisotropy of the formation resistivity based on the currents and corresponding voltages of all the electrode microelements on each pair of pads and the orientation of each pair of pads further comprises:

and performing median filtering processing on all the resistivity data.

An embodiment of a second aspect of the present application provides a formation resistivity azimuthal anisotropy determining apparatus, including:

the polar plate setting module is used for setting a plurality of pairs of polar plates in a well logging borehole at intervals of a set angle, wherein two polar plates in one pair of polar plates are arranged in a central symmetry manner, a plurality of electrode microelements are arranged on each polar plate, and each polar plate is at the same set height;

the current acquisition module is used for applying a voltage to each polar plate to acquire the current of all electrode microelements on each polar plate;

and the anisotropy determining module is used for determining the formation resistivity azimuth anisotropy according to the currents and the corresponding voltages of all the electrode microelements on each pair of polar plates and the azimuth of each pair of polar plates.

In some embodiments, the voltage applied to each plate is the same or different.

In some embodiments, each electrode element is arranged in an array on the plate on which it is located.

In certain embodiments, an anisotropy determination module, comprising:

the resistivity data generating unit is used for generating a micro-resistivity scanning image according to the current and the corresponding voltage of all electrode micro-elements on each pair of polar plates, and the micro-resistivity scanning image can represent the resistivity data corresponding to the electrode micro-elements on each pair of polar plates;

the distribution frequency histogram generating unit is used for counting the resistivity data in each microresistivity scanning image, setting a distribution interval according to the counted resistivity data and correspondingly generating a distribution frequency histogram;

the characteristic value sequence generating unit is used for sequencing each characteristic value by taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as the characteristic value of a corresponding pair of polar plates to obtain a characteristic value sequence;

and the anisotropic coefficient generating unit is used for solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence, and generating the anisotropic coefficient of the resistivity corresponding to the azimuth angle corresponding to the set height according to the azimuth angle of the polar plate corresponding to the first characteristic value and the last characteristic value.

In some embodiments, the electrode elements on each plate are arranged to form a plurality of columns, each column being spaced apart by a set distance.

In certain embodiments, an anisotropy determination module, comprising:

the resistivity data generating unit is used for generating a micro-resistivity scanning image according to the current and the corresponding voltage of the electrode infinitesimal in each column, and the micro-resistivity scanning image can represent the resistivity data corresponding to the electrode infinitesimal in each column;

a distribution frequency histogram generating unit for counting the resistivity data of each column, setting a distribution interval according to the counted resistivity data, and correspondingly generating a distribution frequency histogram;

the characteristic value sequence generating unit is used for sequencing each characteristic value by taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as the characteristic value corresponding to one column to obtain a characteristic value sequence;

and the anisotropic coefficient generating unit is used for solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence and generating the anisotropic coefficient of the resistivity corresponding to the azimuth angle of the set height according to the azimuth angle of the corresponding column of the first characteristic value and the last characteristic value.

In some embodiments, the anisotropy determination module further comprises:

and the azimuth angle calibration unit is used for calibrating the azimuth angle of each pair of polar plates.

In some embodiments, the anisotropy determination module further comprises:

and the azimuth angle calibration unit is used for calibrating the azimuth angle of each column.

In some embodiments, the anisotropy determination module further comprises:

and the median filtering processing unit is used for performing median filtering processing on all the resistivity data.

An embodiment of a third aspect of the present application provides an electronic device, comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the method for determining azimuthal anisotropy of resistivity of a formation as described above when executing the program.

An embodiment of a fourth aspect of the present application provides a computer readable storage medium having stored thereon a computer program which, when being executed by a processor, carries out the steps of the method for determining azimuthal anisotropy of resistivity of a formation as set forth above.

The beneficial effect of this application is as follows:

according to the method and the device for determining the formation resistivity and the azimuth anisotropy, the electrical anisotropy is determined by the resistivity measured by the current and the voltage of the electrode infinitesimal on the multiple pairs of polar plates, the multiple pairs of polar plates have set azimuth angles, so that the measured electrical anisotropy further comprises azimuth information, and the anisotropy data is more accurate and comprehensive.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.

FIG. 1 shows a schematic flow chart of a method for determining azimuthal anisotropy of formation resistivity in an embodiment of the present application.

FIG. 2a shows a top view of a wellbore in an embodiment of the present application.

Fig. 2b shows a schematic diagram of the paired plates in the embodiment of the present application.

FIG. 2c shows a schematic view of a button electrode in an embodiment of the present application.

Fig. 3 shows a schematic diagram of an example of the comparison between the electrical anisotropy and the array acoustic wave anisotropy in the embodiment of the present application.

FIG. 4 is a schematic structural diagram of an apparatus for determining azimuthal anisotropy of resistivity of a formation in an embodiment of the present application.

Fig. 5 shows a schematic structural diagram of an electronic device suitable for implementing embodiments of the present application.

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

FIG. 1 illustrates a method for determining azimuthal anisotropy of formation resistivity, comprising:

s1: arranging a plurality of pairs of polar plates in a well logging borehole at intervals of a set angle, wherein two polar plates in one pair of polar plates are arranged in a central symmetry manner, a plurality of electrode microelements are arranged on each polar plate, and each polar plate is at the same set height;

s2: applying a voltage to each polar plate to obtain the current of all electrode infinitesimal on each polar plate;

s3: and determining the azimuthal anisotropy of the resistivity of the formation according to the currents and the corresponding voltages of all the electrode microelements on each pair of polar plates and the azimuth of each pair of polar plates.

According to the method for determining the formation resistivity and the azimuth anisotropy, the electrical anisotropy is determined by the resistivity measured by the current and the voltage of the electrode infinitesimal on the multiple pairs of polar plates, and the multiple pairs of polar plates have set azimuth angles, so that the measured electrical anisotropy further comprises azimuth information, and anisotropy data is more accurate and comprehensive.

In some embodiments, the plates may be provided by microresistivity imaging scanners, which are currently used primarily for crack identification, and there is no way to extract the anisotropy of the resistivity distribution in different orientations based on images.

In this embodiment, step S3 specifically includes:

s311: generating a micro-resistivity scanning image according to the current and the corresponding voltage of all electrode micro-elements on each pair of polar plates, wherein the micro-resistivity scanning image can represent resistivity data corresponding to the electrode micro-elements on each pair of polar plates;

s312: counting resistivity data in each microresistivity scanning image, setting a distribution interval according to the counted resistivity data, and correspondingly generating a distribution frequency histogram;

s313: taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as a characteristic value of a corresponding pair of polar plates, and sequencing each characteristic value to obtain a characteristic value sequence;

s314: and solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence, and generating an anisotropy coefficient of the resistivity corresponding to the azimuth angle of the corresponding polar plate at the set height according to the azimuth angle of the corresponding polar plate of the first characteristic value and the last characteristic value.

The resistivity data detected by each electrode infinitesimal can be represented on the micro-resistivity scanning image, the resistivity data is counted according to a certain distribution interval, a distribution frequency graph corresponding to each pair of polar plates is further obtained, the median of the distribution interval with the highest frequency is used as a characteristic value, namely, each pair of polar plates generates a characteristic value, and the characteristic value can be used for representing the resistivity characteristic value of the azimuth corresponding to the polar plates.

In some embodiments, a plurality of polar plates (corresponding to different orientations) on the instrument are grouped according to distribution orientations, the polar plates with 180-degree difference form a polar plate pair, and corresponding characteristic values of the micro-resistivity data matrix are counted; and carrying out similar treatment on the pole plate pairs in different directions, judging the directions of the maximum value and the minimum value, and adopting the ratio of the maximum value to the minimum value as the electrical anisotropy coefficient of the depth.

More specifically, please refer to fig. 2 a. Firstly, pole plate assembly is carried out according to specific instrument types, and 4 pole plates are totally 1#, 2#, 3# and 4# respectively, and the phase difference is 90 degrees in pairs. Each polar plate is provided with a main polar plate and an auxiliary polar plate, so that the total number of the polar plates is 8. The pairing mode is that a 1# main polar plate is paired with a 3# main polar plate, a 1# auxiliary polar plate is paired with a 3# auxiliary polar plate, a 2# main polar plate is paired with a 4# main polar plate, and a 2# auxiliary polar plate is paired with a 4# auxiliary polar plate, so that 4 pairs of polar plates are formed.

The method for measuring the electrical anisotropy can be carried out by a method for inverting the resistivity anisotropy through array lateral resistivity logging, and compared with a method for judging the anisotropy by extracting fast and slow transverse waves through array acoustic logging, the method not only can provide electrical anisotropy information, but also can judge the distribution directions of the maximum and minimum resistivities of the stratum. Compared with a method for inverting the resistivity anisotropy by array lateral resistivity logging, the method provided by the invention has the advantage that the electrical anisotropy result has higher-resolution azimuth information, but the electrical anisotropy result can only obtain average horizontal resistivity and vertical resistivity and cannot reflect azimuth information.

In some embodiments, to avoid the problem of excessive data throughput, a processing window may be selected within which the resistivity data is selected.

In addition, in some embodiments, median filtering may be performed to reject noise.

For example, within one processing window length, all the resistivity data subjected to median filtering are grouped according to the group grouping mode, and are divided into 4 groups, the distribution frequency histogram of each group is counted, the peak value is selected as the characteristic value, and the resistivity characteristic values of 4 different directions are obtained in total.

In order to make the anisotropy calculation more accurate, the application further arranges the electrode microelements on each polar plate to form a plurality of columns, and the distance between every two columns is set.

In this embodiment, since each column has a distance interval, the converted azimuth angle also has a converted angle corresponding to the distance interval, and at this time, each column is regarded as a processing group, and the resistivity data of each column is counted, that is, step S3 in this embodiment specifically includes:

s321: generating a microresistivity scanned image according to the current and the corresponding voltage of the electrode microelements of each column, wherein the microresistivity scanned image can represent resistivity data corresponding to the electrode microelements in each column;

s322: counting the resistivity data of each column, setting a distribution interval according to the counted resistivity data, and correspondingly generating a distribution frequency histogram;

s323: taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as a characteristic value corresponding to one column, and sequencing each characteristic value to obtain a characteristic value sequence;

s324: and solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence, and generating an anisotropy coefficient of the resistivity corresponding to the azimuth angle of the corresponding column corresponding to the set height according to the azimuth angle of the corresponding column of the first characteristic value and the last characteristic value.

In this embodiment, each column has a characteristic value, and each column corresponds to an azimuth angle, so that the calculation accuracy of the anisotropy coefficient can be made higher.

In addition, the method further comprises: the azimuth angle of each plate is calibrated or the azimuth angle of each column is calibrated.

Specifically, one of the plates or the columns may be calibrated, and then the azimuth angle of each plate or column is calculated through geometric conversion.

In conjunction with specific examples, fig. 2a to 2c are schematic diagrams of a plate structure of a conventional microresistivity scanning imaging logging instrument. FIG. 2a is a schematic diagram of the plate configuration when the tool is in operation in a borehole, FIG. 2b is a schematic diagram of a pair of plates 180 apart, and FIG. 2c is a schematic diagram of a matrix of button electrodes on the plate, each button electrode measuring a resistivity. Fig. 3 shows the actual processing result. The first path in the diagram is depth, and the second path is GR curve and double-borehole-diameter curve reflecting the lithology of the stratum; the third is a micro-scanning static image; the fourth step is that the electrical anisotropy curve DANI extracted by the method and the time difference anisotropy curve SLOANI extracted from the fast and slow transverse waves of the array acoustic wave have good consistency, which shows that the method can accurately calculate the formation electrical anisotropy coefficient consistent with the array acoustic wave logging.

FIG. 4 shows an apparatus for determining azimuthal anisotropy of resistivity of a formation, according to an embodiment of the present application, including:

the polar plate setting module 1 is used for setting a plurality of pairs of polar plates in a well logging borehole at intervals of a set angle, wherein two polar plates in one pair of polar plates are arranged in a central symmetry manner, a plurality of electrode microelements are arranged on each polar plate, and each polar plate is at the same set height;

the current acquisition module 2 is used for applying a voltage to each polar plate to acquire the current of all electrode microelements on each polar plate;

and the anisotropy determining module 3 is used for determining the formation resistivity azimuth anisotropy according to the currents and the corresponding voltages of all the electrode microelements on each pair of polar plates and the azimuth of each pair of polar plates.

The formation resistivity azimuth anisotropy determining device determines electrical anisotropy by using the resistivity measured by the current and the voltage of the electrode infinitesimal on the pairs of polar plates, and the measured electrical anisotropy further comprises azimuth information because the pairs of polar plates have set azimuth angles, so that anisotropy data are more accurate and comprehensive, a novel anisotropy determining method is provided, the problems that the existing means for researching the anisotropy is single and the like are solved, and the anisotropy parameters which do not pass through the acoustic wave reaction can be obtained.

Based on the same inventive concept, in one embodiment, the voltages applied to each plate are the same or different.

Based on the same inventive concept, in one embodiment, each electrode micro element is arranged in an array on the electrode plate on which the electrode micro element is arranged.

Based on the same inventive concept, in an embodiment, the anisotropy determining module includes:

the resistivity data generating unit is used for generating a micro-resistivity scanning image according to the current and the corresponding voltage of all electrode micro-elements on each pair of polar plates, and the micro-resistivity scanning image can represent the resistivity data corresponding to the electrode micro-elements on each pair of polar plates;

the distribution frequency histogram generating unit is used for counting the resistivity data in each microresistivity scanning image, setting a distribution interval according to the counted resistivity data and correspondingly generating a distribution frequency histogram;

the characteristic value sequence generating unit is used for sequencing each characteristic value by taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as the characteristic value of a corresponding pair of polar plates to obtain a characteristic value sequence;

and the anisotropic coefficient generating unit is used for solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence, and generating the anisotropic coefficient of the resistivity corresponding to the azimuth angle corresponding to the set height according to the azimuth angle of the polar plate corresponding to the first characteristic value and the last characteristic value.

Based on the same inventive concept, in one embodiment, the electrode microelements on each electrode plate are arranged to form a plurality of columns, and a set distance is arranged between each column.

In this embodiment, the anisotropy determination module includes:

the resistivity data generating unit is used for generating a micro-resistivity scanning image according to the current and the corresponding voltage of the electrode infinitesimal in each column, and the micro-resistivity scanning image can represent the resistivity data corresponding to the electrode infinitesimal in each column;

a distribution frequency histogram generating unit for counting the resistivity data of each column, setting a distribution interval according to the counted resistivity data, and correspondingly generating a distribution frequency histogram;

the characteristic value sequence generating unit is used for sequencing each characteristic value by taking the median of the distribution interval with the highest frequency in each distribution frequency histogram as the characteristic value corresponding to one column to obtain a characteristic value sequence;

and the anisotropic coefficient generating unit is used for solving the ratio of the first characteristic value to the last characteristic value in the characteristic value sequence and generating the anisotropic coefficient of the resistivity corresponding to the azimuth angle of the set height according to the azimuth angle of the corresponding column of the first characteristic value and the last characteristic value.

Based on the same inventive concept, in an embodiment, the anisotropy determining module further includes:

and the azimuth angle calibration unit is used for calibrating the azimuth angle of each pair of polar plates.

Based on the same inventive concept, in an embodiment, the anisotropy determining module further includes:

and the azimuth angle calibration unit is used for calibrating the azimuth angle of each column.

Based on the same inventive concept, in an embodiment, the anisotropy determining module further includes:

and the median filtering processing unit is used for performing median filtering processing on all the resistivity data.

An embodiment of the present application further provides a specific implementation manner of an electronic device capable of implementing all steps in the method in the foregoing embodiment, and referring to fig. 5, the electronic device specifically includes the following contents:

a processor (processor)601, a memory (memory)602, a communication Interface (Communications Interface)603, and a bus 604;

the processor 601, the memory 602 and the communication interface 603 complete mutual communication through the bus 604;

the processor 601 is configured to call the computer program in the memory 602, and the processor executes the computer program to implement all the steps of the method in the above embodiments, for example, when the processor executes the computer program, the processor implements the following steps:

s1: arranging a plurality of pairs of polar plates in a well logging borehole at intervals of a set angle, wherein two polar plates in one pair of polar plates are arranged in a central symmetry manner, a plurality of electrode microelements are arranged on each polar plate, and each polar plate is at the same set height;

s2: applying a voltage to each polar plate to obtain the current of all electrode infinitesimal on each polar plate;

s3: and determining the azimuthal anisotropy of the resistivity of the formation according to the currents and the corresponding voltages of all the electrode microelements on each pair of polar plates and the azimuth of each pair of polar plates.

From the above description, it can be seen that the present application provides a computer device for determining electrical anisotropy using resistivity measured by current and voltage of electrode infinitesimal on a plurality of pairs of polar plates, and since the plurality of pairs of polar plates have set azimuth angles, the measured electrical anisotropy further includes azimuth information, and anisotropy data is more accurate and comprehensive, and a new anisotropy determining method is provided, which solves the problems of single means for studying anisotropy and the like at present, and can obtain anisotropy parameters that do not pass through acoustic wave reaction.

Embodiments of the present application also provide a computer-readable storage medium capable of implementing all the steps of the method in the above embodiments, where the computer-readable storage medium stores thereon a computer program, and the computer program when executed by a processor implements all the steps of the method in the above embodiments, for example, the processor implements the following steps when executing the computer program:

s1: arranging a plurality of pairs of polar plates in a well logging borehole at intervals of a set angle, wherein two polar plates in one pair of polar plates are arranged in a central symmetry manner, a plurality of electrode microelements are arranged on each polar plate, and each polar plate is at the same set height;

s2: applying a voltage to each polar plate to obtain the current of all electrode infinitesimal on each polar plate;

s3: and determining the azimuthal anisotropy of the resistivity of the formation according to the currents and the corresponding voltages of all the electrode microelements on each pair of polar plates and the azimuth of each pair of polar plates.

As can be seen from the above description, the computer-readable storage medium provided by the present application determines electrical anisotropy by using resistivity measured by current and voltage of electrode infinitesimal on a plurality of pairs of polar plates, and since the plurality of pairs of polar plates have set azimuth angles, the measured electrical anisotropy further includes azimuth information, and anisotropy data is more accurate and comprehensive, thereby providing a new anisotropy determining method, solving the problems of single means for studying anisotropy, etc., and being capable of obtaining anisotropy parameters that do not pass through acoustic wave reaction.

The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the hardware + program class embodiment, since it is substantially similar to the method embodiment, the description is simple, and the relevant points can be referred to the partial description of the method embodiment. Although embodiments of the present description provide method steps as described in embodiments or flowcharts, more or fewer steps may be included based on conventional or non-inventive means. The order of steps recited in the embodiments is merely one manner of performing the steps in a multitude of orders and does not represent the only order of execution. When an actual apparatus or end product executes, it may execute sequentially or in parallel (e.g., parallel processors or multi-threaded environments, or even distributed data processing environments) according to the method shown in the embodiment or the figures. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, the presence of additional identical or equivalent elements in a process, method, article, or apparatus that comprises the recited elements is not excluded. For convenience of description, the above devices are described as being divided into various modules by functions, and are described separately. Of course, in implementing the embodiments of the present description, the functions of each module may be implemented in one or more software and/or hardware, or a module implementing the same function may be implemented by a combination of multiple sub-modules or sub-units, and the like. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form. The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks. As will be appreciated by one skilled in the art, embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, embodiments of the present description may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, embodiments of the present description may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and so forth) having computer-usable program code embodied therein. The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment. In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of an embodiment of the specification. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction. The above description is only an example of the embodiments of the present disclosure, and is not intended to limit the embodiments of the present disclosure. Various modifications and variations to the embodiments described herein will be apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the embodiments of the present specification should be included in the scope of the claims of the embodiments of the present specification.