CN110824566B

CN110824566B - Seismic attribute fusion method and device and storage medium

Info

Publication number: CN110824566B
Application number: CN201810912109.0A
Authority: CN
Inventors: 陈康; 张旋; 冉崎; 谢冰; 申建波; 孔令霞
Original assignee: Petrochina Co Ltd
Current assignee: Petrochina Co Ltd
Priority date: 2018-08-10
Filing date: 2018-08-10
Publication date: 2021-09-28
Anticipated expiration: 2038-08-10
Also published as: CN110824566A

Abstract

The invention discloses a method and a device for fusing seismic attributes and a storage medium, and belongs to the technical field of seismic data interpretation. The method comprises the following steps: for each position point of a target reservoir of a target oil and gas block, acquiring an attribute value corresponding to each seismic attribute in a plurality of seismic attributes at the corresponding position point, determining the fusion weight of each seismic attribute in the plurality of seismic attributes, wherein the fusion weight of each seismic attribute in the plurality of seismic attributes is used for representing the importance degree of the corresponding seismic attribute in the seismic attribute fusion process, and fusing the plurality of attribute values at each position point according to the fusion weight of each seismic attribute to obtain the seismic attribute fusion result of the target reservoir. According to the method, different fusion weights are set for each seismic attribute in the plurality of seismic attributes, and then the plurality of seismic attributes are fused according to the fusion weight of each seismic attribute, so that the problem that the distribution result of the oil and gas reservoir has multiple solutions in the related technology is solved.

Description

Seismic attribute fusion method and device and storage medium

Technical Field

The invention relates to the technical field of seismic data interpretation, in particular to a seismic attribute fusion method and device.

Background

Before a new well is deployed in an oil and gas zone block, seismic interpretation needs to be carried out on a seismic data volume of the oil and gas zone block obtained in seismic exploration work, distribution of potential high-yield oil and gas reservoirs in the oil and gas zone block is determined according to seismic interpretation results, and a deployment well position of the new well is determined according to the distribution of the oil and gas reservoirs. The seismic data body comprises attribute values of each seismic attribute in a plurality of seismic attributes corresponding to each reservoir in a plurality of reservoirs at different position points, each seismic attribute in the plurality of seismic attributes can be used for representing a geological feature in an oil and gas zone, and a distribution result of the oil and gas reservoirs can be determined and obtained according to the geological feature corresponding to each seismic attribute. That is, for each reservoir, the distribution results of a plurality of different hydrocarbon reservoirs can be determined according to the attribute values of each seismic attribute in the plurality of seismic attributes of the corresponding reservoir at different position points. When a new well is deployed according to the distribution results of various different oil and gas reservoirs, because one deployment result is possibly determined according to the distribution result of each oil and gas reservoir, when the deployment well position of the new well is determined according to the distribution results of various different oil and gas reservoirs, various deployment results are obtained, and the uncertainty of the deployment well position of the new well is caused. Based on this, it is urgently needed to provide a seismic attribute fusion method, before a new well is deployed, a fusion result is obtained by fusing a plurality of seismic attributes in an oil and gas zone, and a distribution result of an oil and gas reservoir is determined and obtained according to the fusion result, so that uncertainty of the deployed well position of the new well caused by the existence of the distribution results of a plurality of oil and gas reservoirs is avoided.

Disclosure of Invention

The embodiment of the invention provides a method, a device and a storage medium for fusing seismic attributes, which can be used for solving the problem that the deployed well position of a new well is uncertain due to the distribution results of various oil and gas reservoirs in the related art. The technical scheme is as follows:

in a first aspect, a seismic attribute fusion method is provided, the method comprising:

for each position point of a target reservoir of a target oil and gas block, acquiring an attribute value corresponding to each seismic attribute in a plurality of seismic attributes at the corresponding position point;

determining a fusion weight of each seismic attribute in the plurality of seismic attributes, wherein the fusion weight of each seismic attribute in the plurality of seismic attributes is used for representing the importance degree of the corresponding seismic attribute in the seismic attribute fusion process, and the fusion weight of each seismic attribute in the plurality of seismic attributes is different;

and fusing the plurality of attribute values on each position point according to the fusion weight of each seismic attribute to obtain a seismic attribute fusion result of the target reservoir.

Optionally, the determining a fusion weight for each of the plurality of seismic attributes comprises:

obtaining a first test productivity of each production well in a plurality of production wells put into production in a target oil and gas block, wherein the first test productivity refers to the productivity obtained when the productivity test is carried out on the target reservoir of each production well under the specified downhole flow pressure;

obtaining an attribute value corresponding to each seismic attribute in the plurality of seismic attributes of each production well in the plurality of production wells at a position point corresponding to the target reservoir;

determining an influence coefficient of each of the plurality of seismic attributes on the test capacity based on the first test capacity of each of the plurality of production wells and the attribute value of each of the plurality of seismic attributes corresponding to each of the plurality of production wells;

determining a fusion weight for each of the plurality of seismic attributes based on an impact coefficient of each of the plurality of seismic attributes on test productivity.

Optionally, the determining an influence coefficient of each of the plurality of seismic attributes on the test productivity based on the first test productivity of each of the plurality of production wells and the attribute value of each of the plurality of seismic attributes corresponding to each of the plurality of production wells comprises:

for any seismic attribute A, determining an influence coefficient of the seismic attribute A on the tested productivity through the following formula based on the first tested productivity of each production well in the plurality of production wells and the attribute value of the seismic attribute A corresponding to each production well in the plurality of production wells;

wherein r is the influence coefficient of the seismic attribute A on the test productivity, and x is_iIs the attribute value of the seismic attribute A corresponding to the ith production well in the plurality of production wells, y_iRefers to a first tested capacity of an ith production well of the plurality of production wells, and the n refers to the number of the plurality of production wells.

Optionally, the determining a fusion weight of each of the plurality of seismic attributes based on the influence coefficient of each of the plurality of seismic attributes on the test productivity includes:

acquiring N influence coefficients from the plurality of determined influence coefficients, wherein each influence coefficient in the N influence coefficients is larger than other influence coefficients except the N influence coefficients in the plurality of influence coefficients;

determining the fusion weight of the seismic attribute corresponding to each influence coefficient in the N influence coefficients according to the size of the N influence coefficients to obtain N fusion weights corresponding to the N seismic attributes, wherein the N influence coefficients and the N fusion weights are in positive correlation;

and determining the fusion weight of the seismic attribute corresponding to each influence coefficient in the other influence coefficients except the N influence coefficients in the plurality of influence coefficients as a first numerical value.

Optionally, the fusing, according to the fusion weight of each seismic attribute, a plurality of attribute values corresponding to a plurality of seismic attributes at each location point includes:

for any position point B in a plurality of position points in the target reservoir, fusing the plurality of attribute values through the following formula based on the fusion weight of each seismic attribute in the plurality of seismic attributes and the attribute value corresponding to each seismic attribute in the plurality of seismic attributes on the position point B to obtain a fusion result on the position point B in the target reservoir;

wherein T is the fusion result at the position B in the target reservoir, n is the number of the plurality of production wells, and k is_iIs the fusion weight of the ith seismic attribute in the plurality of seismic attributes, t_iThe attribute value is the attribute value corresponding to the ith seismic attribute in the plurality of seismic attributes on the position point B.

And generating a seismic attribute fusion result of the target reservoir based on the determined fusion result on the plurality of position points in the target reservoir.

In a second aspect, there is provided a seismic attribute fusion apparatus, the apparatus comprising:

the acquisition module is used for acquiring an attribute value corresponding to each seismic attribute in a plurality of seismic attributes at each position point of a target reservoir of a target oil-gas block;

the first determination module is used for determining the fusion weight of each seismic attribute in the plurality of seismic attributes, the fusion weight of each seismic attribute in the plurality of seismic attributes is used for representing the importance degree of the corresponding seismic attribute in the seismic attribute fusion process, and the fusion weight of each seismic attribute in the plurality of seismic attributes is different;

and the second determining module is used for fusing the attribute values on each position point according to the fusion weight of each seismic attribute to obtain a seismic attribute fusion result of the target reservoir.

Optionally, the first determining module includes:

the first obtaining unit is used for obtaining a first test productivity of each production well in a plurality of production wells put into production in a target oil and gas zone block, and the first test productivity refers to the productivity obtained when the productivity test is carried out on the target reservoir of each production well under the specified downhole flow pressure;

a second obtaining unit, configured to obtain an attribute value corresponding to each seismic attribute in the plurality of seismic attributes of each production well in the plurality of production wells at a location point corresponding to the target reservoir;

a first determining unit, configured to determine an influence coefficient of each of the plurality of seismic attributes on a test productivity based on a first test productivity of each of the plurality of production wells and an attribute value of each of the plurality of seismic attributes corresponding to each of the plurality of production wells;

and the second determining unit is used for determining the fusion weight of each seismic attribute in the plurality of seismic attributes based on the influence coefficient of each seismic attribute in the plurality of seismic attributes on the test productivity.

Optionally, the first determining unit is specifically configured to:

Optionally, the second determining unit includes:

an obtaining subunit, configured to obtain N influence coefficients from the determined plurality of influence coefficients, where each influence coefficient of the N influence coefficients is larger than other influence coefficients of the plurality of influence coefficients except the N influence coefficients;

the first determining subunit is configured to determine, according to the magnitudes of the N influence coefficients, a fusion weight of the seismic attribute corresponding to each influence coefficient of the N influence coefficients to obtain N fusion weights corresponding to the N seismic attributes, where the N influence coefficients and the N fusion weights are in positive correlation;

and the second determining subunit is used for determining the fusion weight of the seismic attribute corresponding to each influence coefficient in the other influence coefficients except the N influence coefficients in the plurality of influence coefficients as a first numerical value.

Optionally, the second determining module is specifically configured to:

In a third aspect, a computer-readable storage medium is provided, in which a computer program is stored, which, when executed by a processor, implements any of the methods provided in the first aspect above.

The technical scheme provided by the embodiment of the invention at least has the following beneficial effects: in the embodiment of the invention, different fusion weights can be set for each seismic attribute in the plurality of seismic attributes, and then the plurality of seismic attributes are fused according to the fusion weight of each seismic attribute, so that the problem that the well position of new well deployment is uncertain due to the distribution result of a plurality of oil and gas reservoirs in the related technology is solved. In addition, the fusion weight of each seismic attribute represents the importance degree of the corresponding seismic attribute in the seismic attribute fusion process, so that the fusion result obtained by fusing the plurality of seismic attributes according to different fusion weights can more accurately reflect the geological features in the oil and gas zone block, and therefore, when the distribution of the oil and gas reservoir is determined according to the fusion result, the accuracy of determining the distribution of the oil and gas reservoir can be improved, and the accuracy of determining the deployment well position of a new well can be further improved.

Drawings

In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.

FIG. 1 is a schematic flow chart of a seismic attribute fusion method according to an embodiment of the present invention;

FIG. 2 is a schematic flow chart of another seismic attribute fusion method provided by an embodiment of the invention;

FIG. 3 is a schematic structural diagram of a seismic attribute fusion apparatus according to an embodiment of the present invention;

fig. 4 is a schematic structural diagram of a terminal 400 according to an embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.

Before explaining the embodiments of the present invention in detail, terms related to the embodiments of the present invention will be explained.

The nouns involved in the embodiments of the present invention include:

seismic attribute

Seismic attributes refer to mathematical transformations of a post-stack seismic data volume to obtain transformed data that contain geometric, kinematic, kinetic, or statistical features related to seismic waves.

Testing productivity

The productivity test refers to testing the maximum production capacity of the production well in a target reservoir under a plurality of bottom hole flowing pressures to obtain the corresponding productivity of the production well under different bottom hole flowing pressures, wherein when the production well is a production well, the productivity test refers to the corresponding maximum oil production capacity of the production well under different bottom hole flowing pressures; when the production well is a gas production well, the test productivity refers to the maximum gas production capacity of the production well corresponding to different bottom hole flow pressures.

Fig. 1 is a schematic flow chart of a seismic attribute fusion method according to an embodiment of the present invention. Referring to fig. 1, the method comprises the steps of:

step 101: and for each position point of a target reservoir of the target oil and gas block, acquiring an attribute value corresponding to each seismic attribute in the plurality of seismic attributes at the corresponding position point.

Step 102: determining the fusion weight of each seismic attribute in the plurality of seismic attributes, wherein the fusion weight of each seismic attribute in the plurality of seismic attributes is used for representing the importance degree of the corresponding seismic attribute in the seismic attribute fusion process, and the fusion weights of the seismic attributes are different.

Step 103: and fusing the plurality of attribute values on each position point according to the fusion weight of each seismic attribute to obtain a seismic attribute fusion result of the target reservoir.

In the embodiment of the invention, different fusion weights can be set for each seismic attribute in the plurality of seismic attributes, and then the plurality of seismic attributes are fused according to the fusion weight of each seismic attribute, so that the problem that the well position of new well deployment is uncertain due to the distribution result of a plurality of oil and gas reservoirs in the related technology is solved. In addition, the fusion weight of each seismic attribute represents the importance degree of the corresponding seismic attribute in the seismic attribute fusion process, so that the fusion result obtained by fusing the plurality of seismic attributes according to different fusion weights can more accurately reflect the geological features in the oil and gas zone block, and therefore, when the distribution of the oil and gas reservoir is determined according to the fusion result, the accuracy of determining the distribution of the oil and gas reservoir can be improved, and the accuracy of determining the deployment well position of a new well can be further improved.

Optionally, determining a fusion weight for each of the plurality of seismic attributes comprises:

obtaining a first test productivity of each production well in a plurality of production wells put into production in the target oil-gas block, wherein the first test productivity refers to the productivity obtained when the productivity test is carried out on each production well in the target reservoir under the specified well bottom flow pressure;

obtaining an attribute value corresponding to each seismic attribute in a plurality of seismic attributes of each production well in a plurality of production wells at a position point corresponding to a target reservoir;

determining an influence coefficient of each seismic attribute of the plurality of seismic attributes on the test productivity based on the first test productivity of each production well of the plurality of production wells and the attribute value of each seismic attribute of the plurality of seismic attributes corresponding to each production well of the plurality of production wells;

determining a fusion weight for each of the plurality of seismic attributes based on an impact coefficient of each of the plurality of seismic attributes on the test productivity.

Optionally, determining an influence coefficient of each of the plurality of seismic attributes on the test productivity based on the first test productivity of each of the plurality of production wells and the attribute value of each of the plurality of seismic attributes corresponding to each of the plurality of production wells, comprises:

for any seismic attribute A, determining the influence coefficient of the seismic attribute A on the test productivity through the following formula based on the first test productivity of each production well in the plurality of production wells and the attribute value of the seismic attribute A corresponding to each production well in the plurality of production wells;

wherein r is the shadow of seismic attribute A on the test productivityCoefficient of influence, x_iRefers to the attribute value, y, of the seismic attribute A corresponding to the ith production well in the plurality of production wells_iRefers to the first tested capacity of the ith production well of the plurality of production wells, and n refers to the number of the plurality of production wells.

Optionally, determining a fusion weight for each of the plurality of seismic attributes based on the influence coefficient of each of the plurality of seismic attributes on the test productivity includes:

arranging the determined multiple influence coefficients according to a descending order to obtain a sequencing result;

acquiring the first N influence coefficients from the sequencing result, and determining the fusion weight of the seismic attribute corresponding to each influence coefficient in the first N influence coefficients according to the size of the first N influence coefficients to obtain N fusion weights corresponding to the N seismic attributes, wherein the first N influence coefficients and the N fusion weights are in positive correlation;

and determining the fusion weight of the seismic attribute corresponding to each influence coefficient in the other influence coefficients except the first N influence coefficients in the plurality of influence coefficients as a first numerical value.

Optionally, fusing a plurality of attribute values corresponding to a plurality of seismic attributes at each location point according to the fusion weight of each seismic attribute, including:

for any position point B in a plurality of position points in the target storage layer, fusing the plurality of attribute values through the following formula based on the fusion weight of each seismic attribute in the plurality of seismic attributes and the attribute value corresponding to each seismic attribute in the plurality of seismic attributes on the position point B to obtain a fusion result on the position point B in the target storage layer;

wherein T refers to the fusion result on the position point B in the target reservoir, n refers to the number of a plurality of production wells, k_iIs the fusion weight, t, of the ith seismic attribute in the plurality of seismic attributes_iThe attribute value is corresponding to the ith seismic attribute in the plurality of seismic attributes on the position point B.

All the above optional technical solutions can be combined arbitrarily to form an optional embodiment of the present invention, which is not described in detail herein.

Fig. 2 is a schematic flow chart of another seismic attribute fusion method according to an embodiment of the present invention. Referring to fig. 2, the method comprises the steps of:

step 201: and for each position point of a target reservoir of the target oil and gas block, acquiring an attribute value corresponding to each seismic attribute in the plurality of seismic attributes at the corresponding position point.

The target reservoir of the target oil and gas block refers to a reservoir to be researched, before a new well is deployed in the target oil and gas block, the target reservoir of the target oil and gas block can be researched to determine oil and gas distribution in the target reservoir, and a basis is provided for deploying the new well in the target oil and gas block according to the determined oil and gas distribution in the target reservoir.

It should be noted that the target reservoir is a plane, and the plane may be composed of different points, so that the research on the target reservoir may be decomposed into a plurality of location points composing the target reservoir. In practical applications, the plurality of location points comprising the target reservoir may be determined by seismic exploration operations performed on the target hydrocarbon block.

The attribute value corresponding to the seismic attribute is converted data obtained by mathematically converting the post-stack seismic data volume, and the converted data includes geometric, kinematic, kinetic, or statistical characteristics related to the seismic wave. In practical application, the attribute value corresponding to each seismic attribute in the plurality of seismic attributes at each position point of the target reservoir may be obtained by user input, may be obtained by transmission from other devices, and may also be obtained by performing mathematical transformation on the post-stack seismic data volume through commercial seismic interpretation software. For example, a post-stack seismic data volume of a target reservoir of the target oil and gas block is obtained, and mathematical transformation is performed on the post-stack seismic data volume by using a seismic attribute interpretation module built in Landmark seismic comprehensive interpretation software, so that an attribute value corresponding to each seismic attribute of a plurality of seismic attributes at each position point of the target reservoir is obtained.

It should be noted that after the post-stack seismic data volume is mathematically transformed, attribute values corresponding to a plurality of seismic attributes can be obtained, and the plurality of seismic attributes can be divided into a plurality of types according to wave kinematics and dynamics, which are respectively: the reservoir stratum seismic data acquisition system comprises eight categories of attributes including an amplitude category attribute, a waveform category attribute, a frequency category attribute, an attenuation category attribute, a phase category attribute, a correlation category attribute, an energy category attribute, a proportion category attribute and the like, wherein each category of attributes comprises seismic attributes with different quantities, and each seismic attribute can represent geological characteristics of a target reservoir stratum. For example, table 1 shows one possible partitioning of multiple seismic attributes by wave kinematics and dynamics. See table 1, where the first column is a type of seismic attribute and the second column is a seismic attribute contained in each of a plurality of seismic attribute types.

TABLE 1

In practical application, all seismic attributes shown in table 1 may be obtained, a target reservoir of a target oil and gas block may be studied through an attribute value corresponding to each seismic attribute in all seismic attributes, or a part of seismic attributes in all seismic attributes shown in table 1 may be obtained according to an actual production condition of the target reservoir of the target oil and gas block, and a target reservoir of the target oil and gas block may be studied through an attribute value corresponding to each seismic attribute in the part of seismic attributes, which is not specifically limited in this embodiment of the present invention.

After the attribute value corresponding to each seismic attribute in the plurality of seismic attributes at each location point of the target reservoir is obtained in step 201, a fusion weight of each seismic attribute in the plurality of seismic attributes may be determined through steps 202 to 205, where the fusion weight of each seismic attribute in the plurality of seismic attributes is used to represent the importance degree of the corresponding seismic attribute in the seismic attribute fusion process, and the fusion weights of the seismic attributes in the plurality of seismic attributes are different.

Step 202: a first test capacity for each of a plurality of production wells commissioned within the target hydrocarbon block is obtained.

The first productivity test refers to testing the maximum production capacity of the production well in the target reservoir under the specified bottom hole flowing pressure to obtain the corresponding productivity of the production well under the bottom hole flowing pressure.

For example, in the target reservoir, the productivity of a plurality of production wells may be tested by using the specified bottom hole flowing pressure, and the measured productivity of each production well is the first tested productivity of each production well. In practical applications, the first test productivity of each production well can be obtained by user input, can be sent by other equipment, and can also be obtained by analyzing the well test data of the production well by the terminal. When the specified bottom hole flowing pressure is atmospheric pressure, the terminal can determine the corresponding productivity of each production well in the multiple production wells as the first test productivity when the bottom hole flowing pressure is atmospheric pressure.

Step 203: and acquiring an attribute value corresponding to each seismic attribute in the plurality of seismic attributes of each production well in the plurality of production wells at a position point corresponding to the target reservoir.

It should be noted that, because the target reservoir is composed of a plurality of location points, and each production well in the plurality of production wells is distributed on the target reservoir in the form of a point, each production well in the plurality of production wells corresponds to a location point in the target reservoir, and after obtaining the attribute value corresponding to each seismic attribute in the plurality of seismic attributes at each location point in the plurality of location points of the target reservoir of the target hydrocarbon block, the location point corresponding to each production well in the target reservoir can be determined according to the location of each production well in the plurality of production wells in the target hydrocarbon block, so as to obtain the attribute value corresponding to each seismic attribute in the plurality of seismic attributes of each production well at the location point corresponding to the target reservoir from the obtained attribute values of the target reservoir.

Step 204: determining an influence coefficient of each seismic attribute of the plurality of seismic attributes on the test productivity based on the first test productivity of each production well of the plurality of production wells and the attribute value of each seismic attribute of the plurality of seismic attributes corresponding to each production well of the plurality of production wells.

It should be noted that each of the multiple seismic attributes represents an original geological feature of the target reservoir, however, in a case where multiple production wells have been deployed in the target hydrocarbon block, the geological feature of the target reservoir may change, and therefore, merely fusing the multiple seismic attributes that can represent the original geological feature of the target reservoir may not accurately reflect the change in the geological feature of the target reservoir after the multiple production wells have been deployed. Based on this, in the embodiment of the present invention, before the terminal fuses the plurality of seismic attributes, an influence coefficient of each of the plurality of seismic attributes on the test productivity may be determined based on the first test productivity of each of the plurality of production wells and the attribute value of each of the plurality of seismic attributes corresponding to each of the plurality of production wells, the fusion weight may be determined according to the influence coefficient of each of the plurality of seismic attributes on the test productivity, and the attribute value of each of the plurality of seismic attributes may be fused according to the fusion weight. Therefore, the fusion weight reflects the influence of the current seismic attribute on the productivity of the production well, namely reflects the influence of the current geological characteristics on the productivity of the production well, so that the fusion result can better accord with the current geological characteristics of the target reservoir by fusing the attribute values of the seismic attributes with the fusion weight.

For any seismic attribute A in the plurality of seismic attributes corresponding to each production well in the plurality of production wells, determining an influence coefficient of the seismic attribute A on the test productivity through the following formula (1) based on the first test productivity of each production well in the plurality of production wells and the attribute value of the seismic attribute A corresponding to each production well in the plurality of production wells;

wherein r is the influence coefficient of the seismic attribute A on the test productivity, x_iRefers to the attribute value, y, of the seismic attribute A corresponding to the ith production well in the plurality of production wells_iRefers to the first tested capacity of the ith production well of the plurality of production wells, and n refers to the number of the plurality of production wells.

Taking 10 production wells in a target hydrocarbon block as an example, table 2 shows the property value of the seismic property, which is the average instantaneous frequency of the seismic properties of a plurality of seismic properties corresponding to each of the 10 production wells, and the first test productivity of the production well. Referring to table 2, the first column of table 2 shows the production well number for each of the 10 production wells, the second column shows the attribute value for the seismic attribute, the average instantaneous frequency for each of the 10 production wells, and the third column shows the first test yield value for each of the 10 production wells.

TABLE 2

Number of well	Average instantaneous frequency	First test productivity (Wanfang/Tian)
			Well 1	0.503	7.273
Well 2	0.922	2.424
			Well 3	0.859	33.939
Well 4	1.090	46.061
			Well 5	2.473	33.939
Well 6	2.620	7.273
			Well 7	2.766	7.273
Well 8	3.416	60.606
			Well 9	4.883	60.606
Well 10	5.952	191.515

Based on the property value of the average instantaneous frequency of each production well in the 10 production wells in the target hydrocarbon zone block and the first test productivity of the production well corresponding to each production well shown in table 2, the influence coefficient of the seismic property of the average instantaneous frequency in the target hydrocarbon zone block on the first test productivity is determined to be 0.587 through the formula (1). It should be noted that the attribute value of the average instantaneous frequency corresponding to each production well and the first test capacity of each production well shown in table 2 are exemplary data given by the embodiment of the present invention, and do not constitute a limitation on the attribute value of the average instantaneous frequency and the first test capacity.

For each seismic attribute in the plurality of seismic attributes, the influence coefficient of the corresponding seismic attribute on the test productivity is determined through the formula (1) by referring to the method for determining the influence coefficient of the average instantaneous frequency on the test productivity, so that the influence coefficient of each seismic attribute in the plurality of seismic attributes on the test productivity is obtained.

Step 205: determining a fusion weight for each of the plurality of seismic attributes based on an impact coefficient of each of the plurality of seismic attributes on the test productivity.

In the embodiment of the invention, the fusion weight of each seismic attribute in the plurality of seismic attributes can be determined based on the influence coefficient of each seismic attribute in the plurality of seismic attributes on the test productivity in the following two ways.

The first mode is as follows: the terminal can calculate the sum of the multiple influence coefficients, determine the ratio of the influence coefficient of each seismic attribute to the test productivity in the multiple influence coefficients according to the influence coefficient of each seismic attribute to the test productivity in the multiple seismic attributes and the sum of the multiple influence coefficients, and determine the ratio of the influence coefficient of each seismic attribute to the test productivity in the multiple seismic attributes to the sum of the influence coefficients as the fusion weight of each seismic attribute in the multiple seismic attributes.

The second mode is as follows: acquiring N influence coefficients from the plurality of determined influence coefficients, wherein each influence coefficient in the N influence coefficients is larger than other influence coefficients except the N influence coefficients in the plurality of influence coefficients; determining the fusion weight of the seismic attribute corresponding to each influence coefficient in the N influence coefficients according to the size of the N influence coefficients to obtain N fusion weights corresponding to the N seismic attributes, wherein the N influence coefficients and the N fusion weights are in positive correlation; and determining the fusion weight of the seismic attribute corresponding to each influence coefficient in the other influence coefficients except the N influence coefficients in the plurality of influence coefficients as a first numerical value.

It should be noted that, because the calculation amount for determining the fusion weight of each seismic attribute in the plurality of seismic attributes is large and the process is complex, N influence coefficients can be selected from the plurality of influence coefficients, and each influence coefficient in the N influence coefficients is larger than other influence coefficients except the N influence coefficients in the plurality of influence coefficients, that is, the N influence coefficients with a larger influence coefficient on the test productivity by each seismic attribute in the plurality of seismic attributes are selected, so as to reduce the calculation amount for determining the fusion weight of each seismic attribute in the plurality of seismic attributes and improve the calculation efficiency.

For example, the terminal may rank the plurality of influence coefficients in a descending order or a descending order to obtain a ranking result, and obtain N influence coefficients from the ranking result. If the sorting result is a sorting result in which a plurality of influence coefficients are arranged from large to small, the first N influence coefficients may be obtained from the sorting result. If the sorting result is obtained according to the sorting from small to large, the last N influence coefficients can be obtained from the sorting result. Wherein N is less than the number of the influence coefficients in the sorting result, and N is a positive integer greater than 1.

After acquiring the N influence coefficients, the terminal may determine the fusion weight of the seismic attribute corresponding to each influence coefficient in the N influence coefficients by using the following formula (2).

Wherein N is the number of the acquired influence coefficients, k_iRefers to the fusion weight r of the seismic attribute corresponding to the ith influence coefficient in the N influence coefficients_iRefers to the ith influence coefficient in the N influence coefficients, wherein i is a positive integer which is greater than 1 and less than or equal to N.

For example, table 3 shows the attribute names and influence coefficients of a plurality of seismic attributes, where the first column of table 3 is the attribute names of the plurality of seismic attributes, the second column is the influence coefficients corresponding to the plurality of seismic attributes, and the third column is the ranking from large to small.

TABLE 3

Attribute name	Coefficient of influence	Rank order
			Root mean square amplitude	0.174	11
Average peak amplitude	0.240	10
			Mean valley amplitude	0.330	7
Maximum peak amplitude	0.579	3
			Analyzing the maximum peak amplitude value in the time window	0.370	6
Maximum valley amplitude	0.271	9
			Average instantaneous frequency	0.587	2
Instantaneous phase	0.615	1
			Length of wave form	0.517	5
Number of peaks	0.274	8
			Trough of wave	0.540	4

As can be seen from table 3, the 11 seismic attribute influence coefficients are arranged from large to small, and the obtained ranking results include an instantaneous phase influence coefficient, an average instantaneous frequency influence coefficient, a maximum peak amplitude influence coefficient, a trough number influence coefficient, a waveform length influence coefficient, a maximum peak amplitude influence coefficient in the analysis time window, an average trough amplitude influence coefficient, a crest number influence coefficient, a maximum trough amplitude influence coefficient, an average peak amplitude influence coefficient, and a root-mean-square amplitude influence coefficient. Assuming that N is 4, in this case, the terminal may obtain, from the sorting result, the first 4 influence coefficients, which are the influence coefficient of the instantaneous phase 0.615, the influence coefficient of the average instantaneous frequency 0.587, the influence coefficient of the maximum peak amplitude 0.579, and the influence coefficient of the number of valleys 0.540, respectively. Determining the fusion weight of the seismic attribute corresponding to each influence coefficient in the first 4 influence coefficients through the formula (2), wherein the obtained fusion weight of the instantaneous phase is 0.2650, the fusion weight of the average instantaneous frequency is 0.2529, the fusion weight of the maximum peak amplitude is 0.2495, and the fusion weight of the trough number is 0.2326.

It should be noted that, after determining the fusion weight of the seismic attribute corresponding to each influence coefficient of the N influence coefficients, the fusion weight of the seismic attribute corresponding to each of the influence coefficients other than the N influence coefficients in the plurality of influence coefficients may also be determined as a first numerical value, specifically, the first numerical value may be 0, that is, the fusion weight of the seismic attribute corresponding to each of the influence coefficients other than the N influence coefficients in the plurality of influence coefficients may be set to 0, in other words, by setting the fusion weight of the seismic attribute corresponding to each of the other influence coefficients than the N influence coefficients among the plurality of influence coefficients to 0, when the plurality of seismic attributes are fused, the fusion result does not contain the influence of the attribute value of the seismic attribute corresponding to each influence coefficient in other influence coefficients except the N influence coefficients in the plurality of influence coefficients.

Step 206: and fusing the plurality of attribute values on each position point according to the fusion weight of each seismic attribute to obtain a seismic attribute fusion result of the target reservoir.

Specifically, the terminal may fuse the multiple attribute values at each location point according to steps 2061 and 2062 to obtain a seismic attribute fusion result of the target reservoir.

Step 2061: for any position point B in a plurality of position points in the target storage layer, fusing the plurality of attribute values through the following formula (3) based on the fusion weight of each seismic attribute in the plurality of seismic attributes and the attribute value corresponding to each seismic attribute in the plurality of seismic attributes on the position point B to obtain a fusion result on the position point B in the target storage layer;

In the embodiment of the invention, the multiple seismic attributes of any one position point a in multiple position points in the target reservoir can be fused in the following two ways.

The first mode is as follows: and fusing the plurality of seismic attributes of any position point B in the plurality of position points in the target reservoir according to a formula (3) based on the fusion weight of each seismic attribute in the plurality of seismic attributes and the attribute value corresponding to each seismic attribute in the plurality of seismic attributes on the position point B.

For example, assume that the seismic attributes are instantaneous phase, average instantaneous frequency, maximum peak amplitude, and trough number, respectively, where for any one position point B of the plurality of position points in the target reservoir, the fusion weight of the instantaneous phase at the position point is 0.2650, the attribute value of the instantaneous phase is 410, the fusion weight of the average instantaneous frequency is 0.2529, the attribute value of the average instantaneous frequency is 270, the fusion weight of the maximum peak amplitude is 0.2495, the attribute value of the maximum peak amplitude is 310, the fusion weight of the trough number is 0.2326, and the attribute value of the trough number is 380. And fusing the plurality of seismic attributes of any position point B in the plurality of position points in the target reservoir according to the formula (3) based on the fusion weight of each seismic attribute in the plurality of seismic attributes and the attribute value corresponding to each seismic attribute in the plurality of seismic attributes on the position point B, so as to obtain a fusion result 342.666 of the position point B.

The second mode is as follows: normalizing the attribute value corresponding to each of the plurality of seismic attributes on the position point B, and fusing the plurality of seismic attributes of any position point B in the plurality of position points in the target storage layer according to a formula (3) based on the fusion weight of each of the plurality of seismic attributes and the normalized value of the attribute value corresponding to each of the plurality of seismic attributes on the position point B.

It should be noted that, because the ranges of the attribute values corresponding to each of the seismic attributes are different, wherein the attribute values of some seismic attributes may not even be in an order of magnitude, in this case, if the attribute values are directly fused according to the fusion weight of each of the seismic attributes, for the attribute values that are too large or too small, the finally obtained fusion result will not truly reflect the influence of the attribute values, that is, the accuracy of the fusion result will be influenced. And then fusing the plurality of seismic attributes of the position point B according to a formula (3) based on the fusion weight of each seismic attribute in the plurality of seismic attributes and the normalized attribute value corresponding to each seismic attribute in the plurality of seismic attributes.

Specifically, when the attribute value corresponding to each seismic attribute in the plurality of seismic attributes is normalized, for each seismic attribute in the plurality of seismic attributes, an attribute value corresponding to the corresponding seismic attribute at each position point in a plurality of position points of a target reservoir of a target hydrocarbon block may be obtained, a maximum attribute value of the corresponding seismic attribute may be determined from the obtained attribute values, and a ratio between the attribute value corresponding to the corresponding seismic attribute at the position point B and the maximum attribute value may be determined as a normalized attribute value corresponding to the corresponding seismic attribute at the position point B.

For example, Table 4 shows the maximum attribute value for each of the plurality of seismic attributes within the target reservoir and the attribute value for each of the plurality of seismic attributes at location point B. Wherein the first column shows attribute names of the plurality of seismic attributes, the second column shows attribute values of each of the plurality of seismic attributes at location point B, the third column shows maximum attribute values of each of the plurality of seismic attributes within the target reservoir, the fourth column shows corresponding normalized attribute values of each of the plurality of seismic attributes at location point B, and the fifth column shows fusion weights of each of the plurality of seismic attributes at location point B.

TABLE 4

Seismic attribute	Attribute value	Maximum attribute value	Normalizing attribute values	Fusion weights
					Instantaneous phase	410	500	0.82	0.2650
Average instantaneous frequency	270	450	0.60	0.2529
					Maximum peak amplitude	310	600	0.52	0.2495
Number of wave trough	380	550	0.69	0.2326

The normalized attribute values and fusion weights corresponding to 4 seismic attributes in any one of the plurality of location points B in the target reservoir shown in table 4 are fused by equation (3), and the fusion result of the location point B is obtained as 0.6593.

It should be noted that the terminal may determine the fusion result of each of the plurality of location points of the target reservoir according to the method for determining the fusion result of the location point B described above, and then the terminal may generate the seismic attribute fusion result of the target reservoir based on the fusion result of each of the plurality of location points of the target reservoir through the following step 2062.

Step 2062: and generating a seismic attribute fusion result of the target reservoir based on the determined fusion result on the plurality of position points in the target reservoir.

After determining the fusion result of each of the plurality of location points in the target storage layer, the terminal may obtain a numerical range of the fusion result according to a maximum value and a minimum value of the plurality of fusion results, then divide the numerical range into a plurality of numerical intervals, and set a corresponding chromatogram for each numerical interval. And then, the terminal can determine the color spectrum corresponding to each position point according to the numerical value interval to which the fusion result of each position point in the plurality of position points belongs, further generates the seismic attribute fusion diagram of the target reservoir according to the color spectrum corresponding to each position point, and can determine the seismic attribute fusion diagram as the seismic attribute fusion result of the target reservoir.

Of course, the above description is only an expression form of an exemplary seismic attribute fusion result given in the embodiment of the present invention, and does not constitute a limitation on the seismic attribute fusion result. For example, in one possible implementation, after determining the fusion result of each location point, a data set including the coordinates of each location point in the plurality of location points and the fusion result of each location point is output, and the data set may be the seismic attribute fusion result of the target reservoir. For another example, the terminal may further represent the fusion result of the plurality of location points in the target reservoir in a graph through contours or other representation forms, and at this time, the output image is the seismic attribute fusion result of the target reservoir.

In the embodiment of the invention, different fusion weights can be set for each seismic attribute in the plurality of seismic attributes, the plurality of seismic attributes are fused according to the fusion weight of each seismic attribute to obtain a fusion result, the distribution result of an oil and gas reservoir can be uniquely determined according to the obtained fusion result, the deployment well position of a new well can be uniquely determined according to the distribution result of the oil and gas reservoir, and the problem that the deployment well position of the new well is uncertain due to the distribution results of various oil and gas reservoirs in the related technology is solved. Moreover, the fusion weight of each seismic attribute represents the importance degree of the corresponding seismic attribute in the seismic attribute fusion process, so that the fusion result obtained by fusing the plurality of seismic attributes according to different fusion weights can more accurately reflect the geological features in the oil-gas zone, and therefore, when the distribution of the oil-gas reservoir is determined according to the fusion result, the accuracy of determining the distribution of the oil-gas reservoir can be improved, and the accuracy of determining the deployment well position of a new well can be further improved. In addition, when the plurality of seismic attributes are fused, each seismic attribute in the plurality of seismic attributes can be combined with the first test productivity of the plurality of production wells in the target block, the fusion weight of each seismic attribute in the plurality of seismic attributes is determined according to the influence coefficient of each seismic attribute in the plurality of seismic attributes and the first test productivity, and the plurality of seismic attributes are further fused according to the fusion weight of each seismic attribute, so that the obtained fusion result can reflect the geological characteristics of the target reservoir after the plurality of production wells are deployed more accurately, the distribution of the oil and gas reservoir determined according to the fusion result is more accurate, and the accuracy of the determined well position for deploying the new well is further improved.

Fig. 3 is a schematic structural diagram of a seismic attribute fusion device according to an embodiment of the present invention. Referring to fig. 3, the apparatus may include:

the obtaining module 301 is configured to obtain, for each location point of a target reservoir of a target oil and gas block, an attribute value corresponding to each seismic attribute in a plurality of seismic attributes at the corresponding location point.

The first determining module 302 is configured to determine a fusion weight of each of the plurality of seismic attributes, where the fusion weight of each of the plurality of seismic attributes is used to represent an importance degree of a corresponding seismic attribute in a seismic attribute fusion process, and the fusion weights of the seismic attributes are different.

And the second determining module 303 is configured to fuse the multiple attribute values at each location point according to the fusion weight of each seismic attribute, so as to obtain a seismic attribute fusion result of the target reservoir.

Optionally, the first determining module includes:

the first obtaining unit is used for obtaining a first test productivity of each production well in a plurality of production wells put into production in the target oil-gas zone block, and the first test productivity refers to the productivity obtained when the productivity test is carried out on the target reservoir of each production well under the same well underflow;

the second obtaining unit is used for obtaining an attribute value corresponding to each seismic attribute in the plurality of seismic attributes of each production well in the plurality of production wells at a position point corresponding to the target reservoir;

the first determining unit is used for determining an influence coefficient of each seismic attribute in the plurality of seismic attributes on the test productivity based on the first test productivity of each production well in the plurality of production wells and the attribute value of each seismic attribute in the plurality of seismic attributes corresponding to each production well in the plurality of production wells;

Optionally, the first determining unit is specifically configured to:

Optionally, the second determination unit includes:

the acquiring subunit is configured to acquire N influence coefficients from the determined plurality of influence coefficients, where each influence coefficient of the N influence coefficients is larger than other influence coefficients of the plurality of influence coefficients except the N influence coefficients;

the first determining subunit is used for determining the fusion weight of the seismic attribute corresponding to each influence coefficient in the N influence coefficients according to the size of the N influence coefficients to obtain N fusion weights corresponding to the N seismic attributes, wherein the N influence coefficients and the N fusion weights are in positive correlation;

Optionally, the second determining module is specifically configured to:

It should be noted that: in the seismic attribute fusion device provided in the above embodiment, when the seismic attributes are fused, only the division of the above functional modules is used for illustration, and in practical applications, the above function distribution may be completed by different functional modules according to needs, that is, the internal structure of the device is divided into different functional modules, so as to complete all or part of the above described functions. In addition, the embodiments of the seismic attribute fusion method provided by the embodiments belong to the same concept, and specific implementation processes thereof are described in the embodiments of the methods for details, which are not described herein again.

Fig. 4 is a schematic structural diagram of a terminal 400 according to an embodiment of the present invention. The terminal 400 may be: a smart phone, a tablet computer, an MP3 player (Moving Picture Experts Group Audio Layer III, motion video Experts compression standard Audio Layer 3), an MP4 player (Moving Picture Experts Group Audio Layer IV, motion video Experts compression standard Audio Layer 4), a notebook computer, or a desktop computer. The terminal 400 may also be referred to by other names such as user equipment, portable terminal, laptop terminal, desktop terminal, etc.

Generally, the terminal 400 includes: a processor 401 and a memory 402.

Processor 401 may include one or more processing cores, such as a 4-core processor, an 8-core processor, or the like. The processor 401 may be implemented in at least one hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), and a PLA (Programmable Logic Array). The processor 401 may also include a main processor and a coprocessor, where the main processor is a processor for Processing data in an awake state, and is also called a Central Processing Unit (CPU); a coprocessor is a low power processor for processing data in a standby state. In some embodiments, the processor 401 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content required to be displayed by the display screen. In some embodiments, the processor 401 may further include an AI (Artificial Intelligence) processor for processing computing operations related to machine learning.

Memory 402 may include one or more computer-readable storage media, which may be non-transitory. Memory 402 may also include high speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 402 is used to store at least one instruction for execution by processor 401 to implement the seismic attribute fusion method provided by the method embodiments herein.

In some embodiments, the terminal 400 may further optionally include: a peripheral interface 403 and at least one peripheral. The processor 401, memory 402 and peripheral interface 403 may be connected by bus or signal lines. Each peripheral may be connected to the peripheral interface 403 via a bus, signal line, or circuit board. Specifically, the peripheral device includes: at least one of a radio frequency circuit 404, a display screen 405, a camera assembly 406, an audio circuit 407, a positioning assembly 408, and a power supply 409.

The peripheral interface 403 may be used to connect at least one peripheral related to I/O (Input/Output) to the processor 401 and the memory 402. In some embodiments, processor 401, memory 402, and peripheral interface 403 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 401, the memory 402 and the peripheral interface 403 may be implemented on a separate chip or circuit board, which is not limited by this embodiment.

The Radio Frequency circuit 404 is used for receiving and transmitting RF (Radio Frequency) signals, also called electromagnetic signals. The radio frequency circuitry 404 communicates with communication networks and other communication devices via electromagnetic signals. The rf circuit 404 converts an electrical signal into an electromagnetic signal to transmit, or converts a received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuit 404 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and so forth. The radio frequency circuitry 404 may communicate with other terminals via at least one wireless communication protocol. The wireless communication protocols include, but are not limited to: metropolitan area networks, various generation mobile communication networks (2G, 3G, 4G, and 4G), Wireless local area networks, and/or WiFi (Wireless Fidelity) networks. In some embodiments, the rf circuit 404 may further include NFC (Near Field Communication) related circuits, which are not limited in this application.

The display screen 405 is used to display a UI (User Interface). The UI may include graphics, text, icons, video, and any combination thereof. When the display screen 405 is a touch display screen, the display screen 405 also has the ability to capture touch signals on or over the surface of the display screen 405. The touch signal may be input to the processor 401 as a control signal for processing. At this point, the display screen 405 may also be used to provide virtual buttons and/or a virtual keyboard, also referred to as soft buttons and/or a soft keyboard. In some embodiments, the display screen 405 may be one, disposed on the front panel of the terminal 400; in other embodiments, the display screen 405 may be at least two, respectively disposed on different surfaces of the terminal 400 or in a folded design; in still other embodiments, the display 405 may be a flexible display disposed on a curved surface or a folded surface of the terminal 400. Even further, the display screen 405 may be arranged in a non-rectangular irregular pattern, i.e. a shaped screen. The Display screen 405 may be made of LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), and other materials.

The camera assembly 406 is used to capture images or video. Optionally, camera assembly 406 includes a front camera and a rear camera. Generally, a front camera is disposed at a front panel of the terminal, and a rear camera is disposed at a rear surface of the terminal. In some embodiments, the number of the rear cameras is at least two, and each rear camera is any one of a main camera, a depth-of-field camera, a wide-angle camera and a telephoto camera, so that the main camera and the depth-of-field camera are fused to realize a background blurring function, and the main camera and the wide-angle camera are fused to realize panoramic shooting and VR (Virtual Reality) shooting functions or other fusion shooting functions. In some embodiments, camera assembly 406 may also include a flash. The flash lamp can be a monochrome temperature flash lamp or a bicolor temperature flash lamp. The double-color-temperature flash lamp is a combination of a warm-light flash lamp and a cold-light flash lamp, and can be used for light compensation at different color temperatures.

The audio circuit 407 may include a microphone and a speaker. The microphone is used for collecting sound waves of a user and the environment, converting the sound waves into electric signals, and inputting the electric signals to the processor 401 for processing, or inputting the electric signals to the radio frequency circuit 404 for realizing voice communication. For the purpose of stereo sound collection or noise reduction, a plurality of microphones may be provided at different portions of the terminal 400. The microphone may also be an array microphone or an omni-directional pick-up microphone. The speaker is used to convert electrical signals from the processor 401 or the radio frequency circuit 404 into sound waves. The loudspeaker can be a traditional film loudspeaker or a piezoelectric ceramic loudspeaker. When the speaker is a piezoelectric ceramic speaker, the speaker can be used for purposes such as converting an electric signal into a sound wave audible to a human being, or converting an electric signal into a sound wave inaudible to a human being to measure a distance. In some embodiments, audio circuitry 407 may also include a headphone jack.

The positioning component 408 is used to locate the current geographic position of the terminal 400 for navigation or LBS (Location Based Service). The Positioning component 408 may be a Positioning component based on the GPS (Global Positioning System) of the united states, the beidou System of china, the graves System of russia, or the galileo System of the european union.

The power supply 409 is used to supply power to the various components in the terminal 400. The power source 409 may be alternating current, direct current, disposable or rechargeable. When power source 409 comprises a rechargeable battery, the rechargeable battery may support wired or wireless charging. The rechargeable battery may also be used to support fast charge technology.

That is, not only is an embodiment of the present invention provide a terminal including a processor and a memory for storing processor-executable instructions, where the processor is configured to execute the method in the embodiment shown in fig. 1 or fig. 2, but also an embodiment of the present invention provides a computer-readable storage medium, in which a computer program is stored, and the computer program can implement the seismic attribute fusion method in the embodiment shown in fig. 1 or fig. 2 when being executed by the processor.

Those skilled in the art will appreciate that the configuration shown in fig. 4 is not intended to be limiting of terminal 400 and may include more or fewer components than those shown, or some components may be combined, or a different arrangement of components may be used.

It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program instructing relevant hardware, where the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.

The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims

1. A seismic attribute fusion method, the method comprising:

obtaining a first test productivity of each production well in a plurality of production wells put into production in a target oil and gas zone, wherein the first test productivity refers to the productivity obtained when the productivity test is carried out on each production well in the target reservoir under the specified downhole flow pressure;

determining fusion weight of each seismic attribute in the plurality of seismic attributes based on an influence coefficient of each seismic attribute in the plurality of seismic attributes on the test productivity, wherein the fusion weight of each seismic attribute in the plurality of seismic attributes is used for representing the importance degree of the corresponding seismic attribute in the seismic attribute fusion process, and the fusion weight of each seismic attribute in the plurality of seismic attributes is different;

2. The method of claim 1, wherein determining an impact coefficient of each of the plurality of seismic attributes on the test capacity based on the first test capacity of each of the plurality of production wells and the attribute value for each of the plurality of seismic attributes for each of the plurality of production wells comprises:

3. The method of claim 1, wherein determining the fusion weight for each of the plurality of seismic attributes based on the impact coefficient of each of the plurality of seismic attributes on the test capacity comprises:

4. The method of any one of claims 1-3, wherein fusing the plurality of attribute values corresponding to the plurality of seismic attributes at each location point according to the fusion weight for each seismic attribute comprises:

wherein T refers to the fusion result on the position point B in the target reservoir, n refers to the number of a plurality of production wells, k_iIs the fusion weight, t, of the ith seismic attribute in the plurality of seismic attributes_iThe attribute value is corresponding to the ith seismic attribute in the plurality of seismic attributes on the position point B;

5. A seismic attribute fusion apparatus, comprising:

the second determining module is used for fusing the attribute values on each position point according to the fusion weight of each seismic attribute to obtain a seismic attribute fusion result of the target reservoir;

wherein the first determining module comprises:

6. The apparatus according to claim 5, wherein the first determining unit is specifically configured to:

7. The apparatus according to claim 6, wherein the second determining unit comprises:

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