Embodiment
The embodiment of the present invention is described in detail below in conjunction with attached drawing.It should be appreciated that this place is retouched
The embodiment stated is merely to illustrate and explain the present invention, and is not intended to limit the invention.
According to an aspect of the present invention, there is provided a kind of CO2Seal monitoring method up for safekeeping.Fig. 1 shows reality according to the present invention
Apply the CO of mode2Seal the flow chart of monitoring method up for safekeeping, wherein, the method can be by computer-implemented.As shown in Figure 1, the party
Method can include:In step sl, gas parameter and groundwater parameter are gathered.In the present invention, the gas parameter can wrap
Include at least one of following:CO2Concentration, CO2Flux, sulfur hexafluoride (SF6) concentration, CO2Vertical movement distance.Wherein, CO2It is dense
Degree, CO2Flux and SF6Concentration can be directly acquired by sensor, and CO2Vertical movement distance is needed in advance to VSP
(vertical seismic profiling (VSP)) Time-lapse Seismic Monitoring data are parsed and drawn.In addition, the groundwater parameter can include underground
Water ph value and/or ground water hardness.The method can also include:In step s 2, according to the gas parameter and described
Lower water parameter determines carbon dioxide sequestration safety parameters.
In addition, in CO provided by the invention2Seal up for safekeeping in monitoring method, Land Surface Parameters can also be gathered in step sl, its
Described in Land Surface Parameters can include earth's surface deformation data (for example, can monitor to obtain by radar).In this case, exist
In step S2, it can determine that carbon dioxide sequestration is pacified according to the gas parameter, the groundwater parameter and the Land Surface Parameters
Full sex index.
In order to enable acquired results are more accurate, comprehensive and reliable, it is preferable that gather the data of above-mentioned seven parameters of whole
Information, and determine CO using the data message of this seven parameters2Seal safety parameters up for safekeeping.
In the present invention, the weight coefficient that analytic hierarchy process (AHP) determines every kind of parameter can be in advance based on.Afterwards, Ke Yigen
The carbon dioxide sequestration safety parameters are determined according to the data message and the corresponding weight coefficient of parameter of every kind of parameter.
When being determined using analytic hierarchy process (AHP) progress weight coefficient, by CO2Concentration and CO2Flux is referred to as CO2Parameter, will
Underground water ph value and ground water hardness are referred to as groundwater parameter.Afterwards, CO is determined first using analytic hierarchy process (AHP)2Parameter,
Groundwater parameter, SF6Concentration, CO2The weight coefficient of this five parameters of Vertical movement distance and earth's surface deformation data.For example, utilize
Analytic hierarchy process (AHP) carries out weight coefficient to above-mentioned five parameters and determines, can obtain following judgment matrix:
Wherein, A1It can represent groundwater parameter, A2It can represent CO2Parameter, A3It can represent SF6Concentration, A4Can be with table
Show earth's surface deformation data, and A5It can represent CO2Vertical movement distance.
The weight coefficient of above-mentioned five parameters drawn using analytic hierarchy process (AHP) is as shown in table 1 below:
Table 1
|
A1 |
A2 |
A3 |
A4 |
A5 |
Weight coefficient |
0.0485 |
0.0805 |
0.2916 |
0.1908 |
0.3886 |
Wherein Maximum characteristic root λmax=5.2133, consistency ratio CR=0.0476 < 0.1, meets analytic hierarchy process (AHP)
It is required that.
Due to groundwater parameter A1Including underground water ph value and ground water hardness and CO2Parameter A2Including CO2Concentration and
CO2Flux, therefore, can be with water ph value under base area and ground water hardness each to CO2The influence degree of security is sealed up for safekeeping, from upper
The groundwater parameter A that face is drawn using analytic hierarchy process (AHP)1Weight coefficient in definitely descend water ph value and ground water hardness each
Weight coefficient, and can be according to CO2Concentration and CO2Flux is each to CO2Seal the influence degree of security up for safekeeping, it is sharp from above
The CO drawn with analytic hierarchy process (AHP)2Parameter A2Weight coefficient in determine CO2Concentration and CO2The respective weight coefficient of flux.Usually
In the case of, underground water ph value and ground water hardness have an identical influence degree, and CO2Concentration and CO2Flux also has identical
Influence degree.In such a case, it is possible to determine that underground water ph value and the respective weight coefficient of ground water hardness are
0.02425, and CO2Concentration and CO2The respective weight coefficient of flux is 0.04025, as shown in table 2 below:
Table 2
After the weight coefficient of every kind of parameter is drawn according to analytic hierarchy process (AHP), it is possible to according to the every kind of parameter gathered
Data message, to determine CO2Safety parameters are sealed up for safekeeping, for example, CO can be determined by below equation (2)2Seal safety up for safekeeping
Sex index:
Wherein, SF can represent CO2Seal safety parameters up for safekeeping;WiIt can represent the weight coefficient of i-th kind of parameter;X′iCan be with
Represent the numerical value after the normalization of the data message of i-th kind of parameter of collection;And n can represent parameter total number.
For underground water ph value and CO2Vertical movement is for this two parameters, and due to it, to belong to numerical value more big more pacify
Full index, therefore, can be normalized by following equation (3):
X′i=(Xi-Xi,min)/(Xi,max-Xi,min) (3)
For ground water hardness, CO2Flux, CO2Concentration, SF6For concentration and earth's surface deformation data this five parameters, by
Belong to the smaller safer index of numerical value in it, therefore, can be normalized by following equation (4):
X′i=(Xi,max-Xi)/(Xi,max-Xi,min) (4)
Wherein, XiIt can represent the data message of i-th kind of parameter of collection;Xi,maxIt can represent the maximum of i-th kind of parameter
Value;Xi,minIt can represent the minimum value of i-th kind of parameter.
For example, in the present invention, for underground water ph value this parameter, its maximum can take 8.5, and minimum value can
To take 5;For ground water hardness this parameter, its maximum can take 650mg/L, and minimum value can take 50mg/L;Pin
To CO2For this parameter of flux, its maximum can take 100000mgm-2h-2, minimum value can take 0mgm-2h-2;For CO2
For this parameter of concentration, its maximum can take 1000000ppm, and minimum value can take 375ppm;For SF6Concentration this
For parameter, its maximum can take 1000ppm, and minimum value can take 1ppm;For earth's surface deformation data this parameter,
Its maximum can take 6, and minimum value can take 1;And for CO2For Vertical movement for this parameter, its maximum can
To take 1600m, minimum value can take 0m.It should be appreciated that the maximum of above-mentioned each parameter and the numerical value of minimum value are only to show
Example property, it is merely to illustrate the present invention rather than the limitation present invention.
By combine many kinds of parameters index, can comprehensively, detect CO exactly2Seal safety parameters up for safekeeping, and improve institute
Obtain the reliability of result.There is provided thus, it is possible to which arrangement etc. is withdrawn in the selection, personnel for follow-up safe early warning, safeguard procedures
Accurately, reliable data are supported.
Fig. 2 shows the CO of another embodiment according to the present invention2Seal the flow chart of monitoring method up for safekeeping.As shown in Fig. 2,
Drawing CO2After sealing safety parameters SF up for safekeeping, in step s3, CO can be determined according to the safety parameters SF2Seal peace up for safekeeping
Full property grade.In the present invention, the CO2Sealing safety grades up for safekeeping can for example include but not limited to:Very safe, safe,
It is basic security, dangerous and very dangerous.For example, CO can be determined by below equation (5)2Seal safety grades up for safekeeping:
Fig. 3 shows the CO of another embodiment according to the present invention2Seal the flow chart of monitoring method up for safekeeping.As shown in figure 3,
This method can also include:In step s 4, in the case where any parameter exceeds the corresponding security value (s) scope of the parameter, into
Row alarm;And/or in step s 5, in the CO2Seal the feelings that safety grades are less than or equal to default safety grades up for safekeeping
Under condition, alarm.
In the present invention, each parameter can have corresponding security value (s) scope.As noted previously, as underground water ph value
And CO2Vertical movement belongs to the more big safer index of numerical value apart from this two parameters, therefore, this corresponding safety of two parameters
Number range has lower limit.Once data are less than the lower limit, show data exception, alarm.For example, it is directed to underground water
For pH value, the lower limit of its security value (s) scope can be set to 6.5;For CO2For Vertical movement distance, its safety
The lower limit of scope can be set to 800m.
Further, since ground water hardness, CO2Flux, CO2Concentration, SF6Concentration and earth's surface deformation data this five parameter categories
In the smaller safer index of numerical value, therefore, this corresponding security value (s) scope of five parameters has upper limit value.Once data are high
In the upper limit value, show data exception, alarm.For example, for ground water hardness, the upper limit value of its safe range can
To be set to 450mg/L;For CO2For flux, the upper limit value of its safe range can be set to 60000mgm-2h-1;
For CO2For concentration, the upper limit value of its safe range can be set to 40000ppm;For SF6For concentration, its safety
The upper limit value of scope can be set to 50ppm;For earth's surface deformation data, the upper limit value of its safe range can be set
It is set to 4.
Alternatively or additionally, in step s 5, in the CO2Seal safety grades up for safekeeping and be less than or equal to default peace
In the case of full property grade, alarm.For example, if it is determined that CO2Safety grades are sealed up for safekeeping as dangerous grade or very
In the case of dangerous grade, alarm.
In the present invention, type of alarm can include but is not limited to the mode of audible alarm, display alarm or SMS alarm.
By alarm, exception, and/or CO can occur in the data of collection2Seal up for safekeeping safety grades it is relatively low in the case of, in time notice
Monitoring personnel, so that monitoring personnel can grasp CO accurately and in time2Seal safety conditions up for safekeeping, and can take in time corresponding
Measure.
It described above is to CO2Seal safety parameters up for safekeeping and method that grade is monitored.Once CO2Security is sealed up for safekeeping to refer to
The relatively low and safety grades of number are relatively low, show CO2Risk of leakage it is larger, in this regard, in another embodiment of the present invention, CO2
Sealing monitoring method up for safekeeping can also realize to CO2Leakage monitoring, it is as discussed below.
Fig. 4 shows the CO of another embodiment according to the present invention2Seal the flow chart of monitoring method up for safekeeping.As shown in figure 4,
This method can also include:In step s 6, CO is gathered2The geographical environment parameter and CO of source of leaks2Gas leakage parameter;In step
In rapid S7, according to the geographical environment parameter and the gas leakage parameter, gas leakage concentration diffusion model is drawn, wherein,
The diffusion model can include the functional relation of gas leakage concentration and lower wind direction distance;And in step s 8, will be predetermined
Lower wind direction distance input described in gas leakage concentration and lower wind direction distance functional relation, draw and each predetermined lower wind direction
Apart from corresponding gas leakage concentration.
Specifically, in step s 6, CO is gathered2The geographical environment parameter and CO of source of leaks2Gas leakage parameter.In this hair
In bright, the geographical environment parameter can include:Absolute pressure P (unit atm), temperature T (unit K) and wind speed u (units
For m/s);And the gas leakage parameter can include:Leak source strength Q (unit kg/s).
, in the step s 7, can be according to the geography after geographical environment parameter and gas leakage parameter is collected
Environmental parameter and the gas leakage parameter, draw gas leakage concentration diffusion model, wherein, the diffusion model can include
The functional relation of gas leakage concentration and lower wind direction distance.
For example, the functional relation of gas leakage concentration and lower wind direction distance can be determined by procedure below:
1) the initial of weight gas cloud plumage is determined according to the absolute pressure P of collection and temperature T and below equation (6) first
Bulk density ρo:
Wherein, ρ can represent CO2Density, in general, ρ can be set as 1.97kg/m3;T0It can represent mark condition temperature
Degree, in general, can be by T0It is set as 273.15K;P0It can represent mark condition pressure;Z can represent the default compressed coefficient, in general,
The z is approximately 1 during low pressure.
2) according to the initial bulk density ρ drawn in step 1)o, collection leakage source strength Q and below equation (7)
To determine the initial volume flux V ' of weight gas cloud plumageo(unit m3/s):
3) according to the initial volume flux V ' drawn in step 2)o, collection wind speed u and below equation (8) come it is true
Surely the initial half-breadth b of heavy gas cloud plumage0With elemental height h0:
bo=2ho=(V 'o/u)1/2 (8)
4) according to the initial bulk density ρ drawn in step 1)o, the initial half-breadth b that draws in step 3)0With it is initial
Height h0, collection wind speed u, it may be determined that go out weight gas cloud plumage cross section half-breadth b with lower wind direction distance x variation relation, such as with
Shown in lower equation (9):
B=bo{1+1.5[gho(ρo-ρa)/ρa]1/2x(ubo)-1}2/3 (9)
Wherein, ρaAtmospheric density can be represented, in general, can be by ρaIt is set as 1.29kg/m3;G can represent that gravity accelerates
Degree, in general, g can be set as 9.8N/kg.
5) according to the initial bulk density ρ drawn in step 1)o, the initial volume flux V ' that draws in step 2)o、
Buoyance flux B can be determined based on below equation (10)0(unit m4/s3):
B0=g (ρ0-ρa)/ρaV′0/π (10)
6) according to the buoyance flux B drawn in step 5)0, collection wind speed u, it may be determined that go out weight gas cloud plumage height h
With the variation relation of lower wind direction distance x, as shown in below equation (11):
Wherein, u*It can represent wind friction velocity (unit m/s), it is generally the case that can be by the u*It is set as that height is
0.1 of mean wind speed at 10m.
7) according to the initial bulk density ρ drawn in step 1)o, the initial half-breadth b that draws in step 3)0With it is initial
Height h0, the heavy gas cloud plumage cross section half-breadth b that draws in step 4) is with the variation relation of lower wind direction distance x and in step
5) the heavy gas cloud plumage height h drawn in lower wind direction distance x variation relation, it may be determined that go out CO2Gas leakage concentration is with
The functional relation of wind direction distance, as shown in below equation (12):
Afterwards, in step s 8, can be dense by predetermined lower wind direction distance x (one or more) the input gas leakage
Degree and the functional relation (that is, above-mentioned equation (12)) of lower wind direction distance, draw let out corresponding with each predetermined lower wind direction distance x
Gas leakage bulk concentration C.
Usual embodiment of above, it may be determined that go out the CO at wind direction distance under difference2The concentration of gas leakage, thus
CO can be detected2Where can be farthest leaked into.
In another embodiment, the diffusion model can also include the gas leakage concentration C and arrival time t's
Functional relation.For example, gas leakage concentration C and the functional relation of arrival time t can be as shown in below equation (13):
In this manner it is possible to according to the gas leakage concentration C and the functional relation of arrival time t, draw described each
At predetermined lower wind direction distance x, the time reached needed for gas leakage concentration C corresponding with the lower wind direction distance x is (that is, described
Arrival time t).Thus, it is possible to detect gas leakage CO2Diffusion rate.
Fig. 5 shows the CO of another embodiment according to the present invention2Seal the flow chart of monitoring method up for safekeeping.As shown in figure 5,
This method can also include:In step s 9, according to the diffusion model, CO is shown2The diffusion process of gas leakage.For example,
CO can be shown in GIS (GIS-Geographic Information System)2The diffusion process of gas leakage, and preferably, difference can be used
Color represents different CO2Gas leakage concentration range, so as to more directly perceived, clear, exactly easy to monitoring personnel
Grasp CO2The spread condition of gas leakage, to take corresponding safeguard procedures or timely arrangement personnel withdraw etc..It should be understood that
, how after diffusion model is drawn, the correlation technique that the diffusion process of gas leakage is shown according to diffusion model is this
Known to field technology personnel, in this regard, the present invention is herein without repeating.
According to another aspect of the present invention, there is provided a kind of carbon dioxide sequestration monitoring device.As shown in fig. 6, the device
It can include:Collecting unit 10, can be used for gathering gas parameter and groundwater parameter, wherein the gas parameter is including following
At least one of:Gas concentration lwevel, carbon dioxide flux, concentration of sulfur hexafluoride, carbon dioxide Vertical movement distance, it is described
Groundwater parameter includes underground water ph value and/or ground water hardness;Safety parameters determination unit 20, can be used for according to
Gas parameter and the groundwater parameter determine carbon dioxide sequestration safety parameters.
In addition, the collecting unit 10 can be also used for collection Land Surface Parameters, wherein the Land Surface Parameters become including earth's surface
Graphic data;And the safety parameters determination unit 20, can be used for according to the gas parameter, the groundwater parameter and
The Land Surface Parameters determine carbon dioxide sequestration safety parameters.
The weight coefficient of every kind of parameter can be determined based on analytic hierarchy process (AHP).In this case, the security
Index determination unit 20, can be used for according to determining the data message and the corresponding weight coefficient of parameter of every kind of parameter
Carbon dioxide sequestration safety parameters.
Fig. 7 shows the CO of another embodiment according to the present invention2Seal the block diagram of monitoring device up for safekeeping.As shown in fig. 7, should
Device can also include:Safety grades determination unit 30, can be used for determining that carbon dioxide seals according to the safety parameters
Deposit safety grades.
Fig. 8 shows the CO of another embodiment according to the present invention2Seal the block diagram of monitoring device up for safekeeping.As shown in figure 8, should
Device can also include:Alarm unit 40, can be used in the case where any parameter exceeds the corresponding safe range of the parameter,
Alarm;And/or in the case where the carbon dioxide sequestration safety grades are less than or equal to default safety grades,
Alarm.
Fig. 9 shows the CO of another embodiment according to the present invention2Seal the block diagram of monitoring device up for safekeeping.As shown in figure 9, institute
State geographical environment parameter and carbon dioxide leakage gas parameter that collecting unit 10 can be also used for collection carbon dioxide leakage source;
And described device can also include:Diffusion model determination unit 50, can be used for according to the geographical environment parameter and described
Gas leakage parameter, draws gas leakage concentration diffusion model, wherein, the diffusion model includes gas leakage concentration and leeward
To the functional relation of distance;And gas leakage concentration determination unit 60, it can be used for predetermined lower wind direction distance input institute
The functional relation of gas leakage concentration and lower wind direction distance is stated, is drawn with each predetermined lower wind direction apart from corresponding gas leakage
Concentration.
In addition, the diffusion model can also include the gas leakage concentration and the functional relation of arrival time.At this
In the case of kind, as shown in Figure 10, described device can also include:Arrival time determination unit 70, can be used for letting out according to
Gas leakage bulk concentration and the functional relation of arrival time, draw at each predetermined lower wind direction distance, reach and the leeward
To apart from the arrival time needed for corresponding gas leakage concentration.
Figure 11 shows the block diagram of the carbon dioxide sequestration monitoring device of another embodiment according to the present invention.Such as Figure 11
Shown, which can also include:Display unit 80, can be used for the diffusion for according to the diffusion model, showing gas leakage
Process.
The CO of the present invention2Seal the CO that monitoring device corresponds to the present invention up for safekeeping2Seal monitoring method up for safekeeping, therefore identical content is not
Repeat again.
To sum up, in CO provided by the invention2Seal up for safekeeping in monitoring method and device, can be with by combining many kinds of parameters index
Comprehensively, CO is detected exactly2Seal safety parameters up for safekeeping, and improve the reliability of acquired results.Thus, it is possible to it is follow-up peace
Full early warning, the selection of safeguard procedures, personnel withdraw arrangement etc. and provide accurate, reliable data support.
The preferred embodiment of the present invention is described in detail above in association with attached drawing, still, the present invention is not limited to above-mentioned reality
The detail in mode is applied, in the range of the technology design of the present invention, a variety of letters can be carried out to technical scheme
Monotropic type, these simple variants belong to protection scope of the present invention.
In addition, although in the accompanying drawings with the operation of particular order the invention has been described method, still, this do not require that or
Hint must perform these operations according to the particular order, or the operation having to carry out shown in whole could realize it is desired
As a result.Additionally or alternatively, it is convenient to omit multiple steps are merged into a step and performed by some steps, and/or by one
Step is decomposed into execution of multiple steps.
One skilled in the art will appreciate that embodiments of the present invention can be implemented as a kind of system, device, equipment, method
Or computer program product.Therefore, embodiments of the present invention can be implemented as following form, i.e.,:It is complete hardware, complete
Full software (including firmware, resident software, microcode etc.), or the form that hardware and software combines.
It is further to note that each particular technique feature described in above-mentioned embodiment, in not lance
In the case of shield, it can be combined by any suitable means.In order to avoid unnecessary repetition, the present invention to it is various can
The combination of energy no longer separately illustrates.
In addition, various embodiments of the present invention can be combined randomly, as long as it is without prejudice to originally
The thought of invention, it should equally be considered as content disclosed in this invention.