WO2021224556A1 - Method for determining a calculation model of the concentration of a compound injected into a fluid mixture - Google Patents

Abstract

Method for determining a calculation model of the concentration of a compound injected into a fluid mixture and flowing into a routing network. The method comprises: a first step (E10)/second step (E20) of taking a first/second sample of the fluid mixture, the first/second sample being taken upstream/downstream of a point (P_INJ_1) at which the compound is injected into the fluid mixture, a step of acquiring (E30), for each of the samples taken, a measurement of a given physical value and a measurement of the concentration of the compound injected into the fluid mixture, the sampling and acquisition steps being iterated, at least one injection parameter being modified at each iteration. The method also comprises steps implemented by a treatment device (15), including: a step of obtaining (E40) the acquired measurements, a step of regression (E50) of the acquired measurements so as to determine a calculation model.

Description

Description

Title of the invention: Method for determining a model for calculating the concentration of a compound injected into a fluid mixture

Prior art

The present invention belongs to the general field of the delivery of a fluid mixture through pipes of a delivery network. It relates more particularly to a method for determining a model for calculating the concentration of a compound injected into a fluid mixture of predefined type and flowing in a network for conveying said fluid mixture. It also relates to a method for controlling such a concentration by means of a calculation model thus determined. The invention finds a particularly advantageous application, although in no way limiting, in the case where the fluid mixture is natural gas and where the compound injected is dihydrogen.

In order to face the challenges linked to the ever greater increase in emissions of C02 (carbon dioxide) into the atmosphere, significant research efforts are now being made to develop solutions making it possible to increase the use of green energies.

[0003] Among the solutions currently being tested, those aiming to inject hydrogen (conventionally in the form of dihydrogen H2) into a natural gas delivery network are known in particular. Indeed, the fact of considering such an injection has a double advantage from an environmental point of view. On the one hand, the injected hydrogen can be produced by electrolysis of water from renewable electricity available on the electricity network. The term “renewable electricity” refers here to electricity generated by means capable of exploiting one or more renewable energy sources (wind turbines, solar panels, etc.). On the other hand, the combustion of hydrogen does not generate any CO2 emissions into the atmosphere, and is limited to the emission of water.

[0004] Also, the operators of natural gas delivery networks aim to gradually increase the proportion of hydrogen (by volume) injected into these networks.

The fact remains that the injection of hydrogen into a natural gas delivery network is a process for which precautions must be taken. In fact, the injection of hydrogen into the distributed natural gas modifies the combustion parameters of the latter, as well as the flame characteristics.

[0006] More generally, the injection of hydrogen modifies the physicochemical properties of natural gas. It is therefore essential to be able to guarantee the conformity of network equipment with the gas mixture formed of hydrogen and natural gas (sealing of equipment such as for example regulators or valves, permeability of polyethylene, embrittlement of materials present on the networks, etc.), but also that of the equipment located downstream of the transmission network (hoses, welded copper pipes, PLT pipes, valves , boiler, cooking appliances, etc.).

Ultimately, injecting hydrogen into a natural gas delivery network involves complying with a plurality of constraints, so as to be able, in particular:

- ensure a continuous supply that responds to variable downstream consumption,

- guarantee the safety of people and goods,

- provide a homogeneous mixture which meets quality specifications and which respects a predetermined proportion of hydrogen.

To meet these constraints, it is essential to be able to check at any time that a predetermined proportion of hydrogen has indeed been injected into the natural gas delivery network. To this end, the solutions currently implemented are largely based on gas chromatography techniques coupled with different types of detectors. After separation of the various compounds forming the gas mixture of hydrogen / natural gas, these solutions make it possible to measure the concentration of a compound, and therefore in particular the concentration of hydrogen, in the gas mixture.

In practice, the analysis of the composition of the gas mixture by chromatography is carried out by taking a sample of said mixture in an annex pipe connected to the main pipe through which the natural gas passes. More precisely, said sample is taken downstream from a hydrogen injection point in the annex pipe. The time required to analyze the composition of such a sample is of the order of a few minutes, which is particularly disadvantageous. Indeed, the volume of gas mixture which flows in the annex pipe during this period cannot be directly returned to the flow circulating in the main pipe, for the reasons of regulation and conformity mentioned above. Consequently, this volume, known as the “buffer volume”, is stored in a dedicated part of the annex pipe, and possibly diverted from the routing network to be routed to a recycling circuit if ever the analysis of the sample is revealed. unsatisfactory.

The management and immobilization of such a buffer volume is particularly complex to implement, and also very expensive, insofar as it can be a very large volume of gas at high pressures . This is all the more problematic in the context of future objectives aimed at increasing the flow rate of injected hydrogen.

[0011] Other solutions making it possible to achieve a shorter analysis time exist, and are based in particular on control techniques implemented by correlation analyzers. However, these solutions cannot be considered as satisfactory in the to the extent that the precision achieved in the analysis of the concentration of injected hydrogen is far from sufficient.

[0012] Furthermore, it should be noted that the problems mentioned above (ie complexity of implementation, high cost, lack of precision) have been described in the context of an injection of hydrogen into a network of natural gas delivery. However, similar problems are also encountered when it is necessary to consider the injection of a given compound into a network for conveying a fluid mixture.

Disclosure of the invention

The object of the present invention is to remedy all or part of the drawbacks of the prior art, in particular those set out above, by proposing a solution which makes it possible, in comparison with the solutions of the prior art, to determine in a simple, rapid, precise and inexpensive manner, the concentration of a compound injected into a fluid mixture of predefined type and flowing in a network for conveying said fluid mixture.

To this end, and according to a first aspect, the invention relates to a method for determining a model for calculating the concentration of a compound injected into a fluid mixture of predefined type and flowing in a network of 'conveying said fluid mixture. Said process comprises:

- a first sampling step, at a first instant and in a local volume of a fluid mixture corresponding to the predefined type flowing in a pipe of a network for conveying said fluid mixture, of a first sample of said fluid mixture , said first sample being taken at a sampling point of said pipe located upstream, in the direction of flow of the fluid mixture within the pipe, from a point of injection of said compound into said fluid mixture,

a second sampling step, at a second instant subsequent to said first instant and in a local volume of the fluid mixture, of a second sample of said fluid mixture, said second sample being taken at a sampling point of said pipe located downstream of said injection point, the local volume in which said second sample is taken corresponding, taking into account the flow of the fluid mixture, to the local volume in which the first sample was taken,

a step of acquiring, for each of the samples taken, a measurement of a given physical quantity as well as a measurement of the concentration of the compound injected into the fluid mixture.

Furthermore, the first sampling step, the second sampling step and the acquisition step are iterated a given number of times, at least one injection parameter of said compound being modified at each iteration, and said method further comprises steps implemented by a treatment system, including:

- a step of obtaining the acquired measurements,

a step of regression of the measurements acquired so as to determine a model for calculating the concentration of the compound injected into the fluid mixture as a function of said physical quantity considered as an explanatory variable.

Thus, it is proposed to measure the difference of a physical quantity before and after injection of the compound into the fluid mixture. In other words, the invention proposes to characterize the variation in concentration of the compound injected by the variation of this physical quantity.

The characterization of the relationship between these variations is formalized analytically by the calculation model determined according to the invention. When such a calculation model is obtained, it can be advantageously used to quickly and easily determine the concentration of hydrogen injected into a fluid mixture corresponding to that for which said calculation model was determined. To do this, it suffices to acquire a measurement of said physical quantity in any sample of a network for conveying the fluid mixture, then to report the value of this measurement in the calculation model to finally obtain the concentration of the compound. injected.

Furthermore, such a calculation model offers the possibility of precisely controlling the variation in the concentration of the compound injected. Indeed, in addition to the fact that the calculation model can be determined so as to be very robust (use of a large number of samples taken), this calculation model can be used on a recurring basis to carry out a control over time. the concentration of the compound injected.

Finally, the calculation model thus determined advantageously makes it possible to reduce all or part of the use of a buffer volume as used in the state of the art. Indeed, and as mentioned above, the calculation model makes it possible to obtain very quickly the concentration of the compound injected, so that it is no longer necessary to immobilize a buffer volume while waiting for a result of assessment of said concentration. In this way, the invention not only makes it possible to simplify the configurations of the routing networks, but also to achieve substantial savings in comparison with the significant costs associated with the management and immobilization of such a buffer volume.

In particular embodiments, the determination method may further include one or more of the following characteristics, taken in isolation or in any technically possible combination.

In particular embodiments, for each sample taken, the physical quantity is the thermal conductivity of the fluid mixture contained in said sample. In particular embodiments, for each sample taken, the physical quantity is the speed of sound in the fluid mixture contained in said sample.

In particular embodiments, for each sample taken, the physical quantity is the viscosity of the fluid mixture contained in said sample.

In particular embodiments, for each sample taken, the physical quantity is the specific heat of the fluid mixture contained in said sample.

In particular modes of implementation, each sample taken is brought to a given temperature so that the acquisition of the measurement of the physical magnitude for said sample is carried out at said temperature.

In particular embodiments, the fluid mixture is natural gas, and said compound is dihydrogen.

According to a second aspect, the invention relates to a method for controlling the concentration of a compound injected into a fluid mixture flowing from a network for conveying said fluid mixture. Said control method is implemented by a control system comprising a control device, and comprises steps of:

- Obtaining, by the control device, of a calculation model of said concentration determined for said fluid mixture according to the invention, so that said calculation model admits a given physical quantity as an explanatory variable,

- Taking a sample of the fluid mixture at a sampling point of a pipe of said delivery network, said sampling point being located downstream, in the direction of flow of the fluid mixture within the pipe , a point of injection of said compound into said fluid mixture,

- acquisition of a measurement of said physical quantity for said sample taken,

- obtaining, by the control device, of said acquired measurement,

- determination, by the control device, of said concentration by applying the calculation model to said acquired measurement.

Said control method inherits the advantages mentioned above in the context of the determination method according to the invention.

In particular modes of implementation, the steps of sampling, acquisition, obtaining and determination are iterated on a recurring basis.

Advantageously, the iteration frequency of the sampling, acquisition, obtaining and determination steps is of the order of a second, or even less than one second, so as to obtain a substantially continuous control of the concentration of the compound injected into the fluid mixture. According to a third aspect, the invention relates to a computer program comprising instructions for the implementation:

steps of obtaining the acquired measurements and determining the calculation model of a determination method according to the invention; Where

steps of obtaining the calculation model, of obtaining the acquired measurement and of determining the concentration of a control method according to the invention, when said program is executed by a computer.

This program can use any programming language, and be in the form of source code, object code, or intermediate code between source code and object code, such as in a partially compiled form, or in n ' any other desirable shape.

According to a fourth aspect, the invention relates to an information or recording medium readable by a computer on which a computer program according to the invention is recorded.

The information or recording medium can be any entity or device capable of storing the program. For example, the medium may comprise a storage means, such as a ROM, for example a CD ROM or a microelectronic circuit ROM, or else a magnetic recording means, for example a floppy disk or a disk. hard.

On the other hand, the information or recording medium can be a transmissible medium such as an electrical or optical signal, which can be conveyed via an electrical or optical cable, by radio or by other means. . The program according to the invention can in particular be downloaded from an Internet type network.

Alternatively, the information or recording medium can be an integrated circuit in which the program is incorporated, the circuit being adapted to execute or to be used in the execution of the method in question.

According to a fifth aspect, the invention relates to a processing system for determining a model for calculating the concentration of a compound injected into a fluid mixture of predefined type and flowing in a routing network. of said fluid mixture. Said processing system comprises:

- first sampling means, configured to take, at a first instant and in a local volume of a fluid mixture corresponding to the predefined type flowing in a pipe of a network for conveying said fluid mixture, a first sample of said fluid mixture fluid mixture, said first sample being taken at a sampling point of said pipe located upstream, in the direction of flow of the fluid mixture within the pipe, from a point of injection of said compound into said mixture fluid, second sampling means, configured to take, at a second instant subsequent to said first instant and in a local volume of the fluid mixture, a second sample of said fluid mixture, said second sample being taken at a sampling point of said pipe located downstream of said injection point, the local volume in which said corresponding second sample is taken, taking into account the flow of the fluid mixture, to the local volume in which the first sample was taken,

acquisition means, configured to acquire, for each of the samples taken, a measurement of a given physical quantity as well as a measurement of the concentration of the compound injected into the fluid mixture.

Said treatment system further comprises control means, configured to modify at least one parameter of the injection of said compound, as well as a treatment device comprising:

- an obtaining module, configured to obtain the acquired measurements,

a regression module, configured to perform a regression of the measurements acquired, so as to determine a model for calculating the concentration of the compound injected into the fluid mixture as a function of said physical quantity considered as an explanatory variable.

According to a sixth aspect, the invention relates to a system for monitoring the concentration of a compound injected into a fluid mixture flowing from a network for conveying said fluid mixture, said monitoring system comprising a device for control comprising a first obtaining module, configured to obtain a calculation model of said concentration determined for said fluid mixture according to the invention, so that said calculation model admits a given physical quantity as an explanatory variable, said control system comprising also :

- sampling means, configured to take a sample of the fluid mixture at a sampling point of a pipe of said delivery network, said sampling point being located downstream, in the direction of flow of the fluid mixture within the pipe, a point of injection of said compound into said fluid mixture,

- acquisition means, configured to acquire a measurement of said physical quantity for said sample taken, and said control device also comprising:

- a second obtaining module, configured to obtain said acquired measurement,

a determination module, configured to determine said concentration by applying the calculation model to said acquired measurement.

According to a seventh aspect, the invention relates to a recording medium comprising a table comprising values of a given physical quantity, said values having been generated by means of a determined calculation model of the invention. Brief description of the drawings

Other characteristics and advantages of the present invention will emerge from the description given below, with reference to the appended drawings which illustrate an exemplary embodiment thereof without any limiting nature. In the figures:

[Fig. 1] FIG. 1 schematically represents, in its environment, a particular embodiment of a processing system according to the invention;

[Fig. 2] FIG. 2 schematically represents an example of hardware architecture of a processing device according to the invention belonging to the processing system of FIG. 1;

[Fig. 3] FIG. 3 represents, in the form of a flowchart, the main steps of a determination method according to the invention, as they are implemented by the processing system of FIG. 1;

[Fig. 4] FIG. 4 represents an example of a table comprising measurements acquired for a first mixture of given natural gas and by acquisition means of the treatment system, said measurements corresponding to the thermal conductivity and to the concentration of injected hydrogen;

[Fig. 5] FIG. 5 represents another example of a table comprising measurements acquired for a second mixture of given natural gas and by the acquisition means of the treatment system, said measurements corresponding to the thermal conductivity and to the concentration of injected hydrogen. ;

[Fig. 6] FIG. 6 represents another example of a table comprising measurements acquired for said first mixture of given natural gas and by the acquisition means of the treatment system, said measurements corresponding to the speed of sound and to the hydrogen concentration. injected;

[Fig. 7] FIG. 7 represents another example of a table comprising measurements acquired for said second mixture of given natural gas and by the acquisition means of the treatment system, said measurements corresponding to the speed of sound and to the hydrogen concentration. injected;

[Fig. 8] FIG. 8 schematically represents, in its environment, a particular embodiment of a control system according to the invention;

[Fig. 9] FIG. 9 schematically represents an example of hardware architecture of a control device according to the invention belonging to the control system of FIG. 8; [Fig. 10] FIG. 10 represents, in the form of a flowchart, a particular embodiment of a control method according to the invention, as implemented by the control system of FIG. 8.

Description of the embodiments

FIG. 1 schematically shows, in its environment, a particular embodiment of a treatment system 10 according to the invention.

Said treatment system 10 is configured to allow the determination of a model for calculating the concentration of a compound injected into a fluid mixture of predefined type and flowing in a network for conveying said fluid mixture.

By "fluid mixture of predefined type", we refer here to a fluid mixture whose composition is prescribed according to given technical specifications, for example technical specifications prescribed by a company (private or public) in charge of the management of said transmission network.

For the remainder of the description, it is considered in a non-limiting manner that said routing network is an urban network configured, in a manner known per se, to supply private homes with said fluid mixture (more precisely distributed equipment in said dwellings), these dwellings being located in a given geographical area. Said fluid mixture corresponds, in this embodiment, to natural gas, and therefore makes it possible, for example, to supply cooking appliances.

It is nevertheless important to note that the invention remains applicable to any type of routing network, such as for example an industrial routing network internal to a factory. Moreover, nothing excludes either considering a fluid which is not natural gas but another gas. More generally, any type of fluid can be considered within the meaning of the invention, and therefore in particular also a liquid, such as for example a liquid hydrocarbon (eg petroleum).

In the present embodiment, the compound injected into the mixture of natural gas, and the concentration of which is sought to be determined within said mixture by means of said system 10, corresponds to dihydrogen.

Insofar as the dihydrogen corresponds to the molecular form of the element hydrogen, and for the sake of simplification, reference is made without distinction in the remainder of the description to a “injection of hydrogen”.

Again, it should be noted that no limitation is attached to the identity of the compound injected into the fluid mixture. When this fluid mixture is a gas mixture, as it is in the case in the present embodiment, said compound may for example be carbon dioxide, oxygen, nitrogen, etc. When the fluid mixture is a liquid mixture, said compound may for example be glycol water, bleach, etc.

Part of the natural gas delivery network is illustrated in no way limiting in the example of Figure 1. Conventionally, said delivery network comprises a plurality of pipes. More particularly, the transmission network comprises a so-called “main” CAN_P_1 pipe, configured to allow the flow of natural gas to the homes. Furthermore, said routing network also includes a pipe, called “annex” CAN_A, forming a deviation from said main pipe CAN_P_1, this deviation being further redirected to said main pipe CAN_P_1.

An injection point P_INJ_1 is arranged in said annex pipe CAN_A, the hydrogen being injected continuously, for a determined period, into the mixture of natural gas at this injection point P_INJ_1. Such an injection point P_INJ_1 is conventionally arranged at the level of a dedicated frame called an “injection station”. Those skilled in the art are aware of such an arrangement, so that these aspects, which go beyond the scope of the present invention, are not detailed further here.

As illustrated in Figure 1, the system 10 comprises first sampling means 11, configured to take, in a first instant and in a local volume of the mixture of natural gas flowing in the annex pipe CAN_A, a first sample of said mixture of natural gas. More particularly, and in accordance with the invention, said first sampling means 11 are configured to take said first sample at a sampling point P_PREL_1 of said annex pipe CAN_A. Said sampling point P_PREL_1 is located upstream, in the direction of flow of the mixture of natural gas within the annex pipe CAN_A, of said injection point P_INJ_1.

By "in a local volume of the natural gas mixture", reference is made here to a volume of natural gas located in the annex pipe CAN_A, substantially in line with the sampling point P_PREL_1 of the first sample. In other words, such a local volume designates a restricted portion of the total volume of natural gas located in the annex pipeline CAN_A at a given time, the first sample being taken from this local volume, this local volume being of course intended to flow. in the CAN_A annex pipe due to the flow rate and pressure imposed in the delivery of the natural gas mixture.

Preferably, the sampling point P_PREL_1 is arranged at the level of the annex pipe CAN_A so as to be swept by a permanent flow of the mixture of natural gas. In addition, this sampling point P_PREL_1 is preferably located on a straight length of the annex pipe CAN_A, far from any element liable to disrupt the flow of said mixture of natural gas (eg: elbow, tee, etc.) For example, the sampling point P_PREL_1 is located on a straight length of the annex pipe CAN_1 and separated from the last disturbing element by a significant distance equal to twenty times the diameter of said CAN_A annex pipe.

Said first sampling means 11, for their part, are for example a sampling rod of design known per se, said sampling rod being configured to be inserted at a predefined distance from the diameter of the annex pipe CAN_A, for example a distance equal to a third of said diameter.

[0054] According to one variant, said sampling rod can be controlled remotely, so that the implementation of the taking of the first sample is fully automated.

According to another variant, said sampling rod is handled by a qualified operator to take said first sample.

However, nothing excludes considering first sampling means 11 which differ from a sampling rod, such means being known to those skilled in the art (sampling probe, etc.).

The system 10 also comprises second sampling means 12, configured to take, at a second instant subsequent to said first instant and in a local volume of the mixture of natural gas, a second sample of said mixture of natural gas. More particularly, and in accordance with the invention, said second sampling means 12 are configured to take said second sample at a sampling point P_PREL_2 of said annex pipe CAN_A. Said sampling point P_PREL_2 is located downstream of said injection point P_INJ_1.

In addition, the local volume in which said second sample is taken corresponds, taking into account the flow of the mixture of natural gas, to the local volume in which the first sample was taken.

By "taking into account the flow of the mixture of natural gas", reference is made here to the parameters characterizing said flow, such as typically the flow rate, the speed, the pressure, etc. Thus, the correspondence between the local volumes within which the first and second samples are taken implies a relation between the instants in which the latter are produced. In other words, the first and second sampling means are configured, in particular, to communicate with each other so as to coordinate said first and second sampling times, so as to take into account the flow path of the local volume in which the first is taken. sample. Such a correspondence advantageously makes it possible to compare with each other samples which were part of the same portion (ie same local volume) having circulated within the annex pipe CAN_A, which is a guarantee of reliability of measurements of quantities intended to be carried out within the samples for the development of a model for calculating the concentration of injected hydrogen, as described in more detail later.

Preferably, the sampling point P_PREL_2 is arranged at the level of the annex pipe CAN_A so as to be swept by a homogeneous flow of the mixture of natural gas. To this end, one or more mixers of design known per se can be used to guarantee the homogeneity of the mixture between the hydrogen and the natural gas at the level of the sampling point P_PREL_2.

It should be noted that the technical characteristics relating to the shape taken by the first sampling means 11 (sampling rod, etc.) still apply of course with regard to the shape taken by said second means of direct debit 12.

The system 10 also comprises acquisition means 13, configured to acquire, for each of the samples taken (first sample and second sample), a measurement of a given physical quantity as well as a measurement of the concentration of injected hydrogen. in the mixture of natural gas.

Conventionally, said acquisition means 13 comprise an acquisition chain comprising at least one sensor dedicated to measuring said physical quantity and at least one sensor dedicated to measuring said concentration of injected hydrogen. Each of these sensors forms a sensitive element configured to provide an analog electrical signal according to variations in the physical quantity with which it is associated. Said acquisition chain also comprises, for example, an acquisition card configured to condition the electrical signal supplied by a sensor, so as to finally deliver a time signal. The conditioning implemented by the acquisition card comprises, for example, amplification and / or filtering. Optionally, said acquisition means 13 also comprise, at the output of the acquisition chain, an analog / digital converter configured to digitize a conditioned electrical signal.

In general, the configuration of such acquisition means 13 is well known to those skilled in the art, and is therefore not detailed here further. In particular, those skilled in the art know how to choose suitable sensors for each of the quantities considered, for example in the catalogs of products offered by specialized manufacturers. He also knows how to position these sensors to obtain the desired measurements.

It should be noted that said acquisition means 13 are not limited to the use of at least one sensor dedicated to measuring the concentration of injected hydrogen. Indeed, nothing excludes considering that the measurement of the concentration of injected hydrogen is a measurement obtained by numerical simulation, typically by means of a software suite designed to numerically simulate an injection of hydrogen into the gas mixture. natural, taking take into account various modeling parameters (structure and layout of the delivery network, gas velocity, gas flow rate, gas pressure, gas temperature, etc.).

According to a particular embodiment, said physical quantity measured for a sample is the thermal conductivity of the fluid mixture contained in said sample.

The choice of thermal conductivity as the measured physical quantity, however, only constitutes a variant of the implementation of the invention. Other physical quantities can therefore be envisaged, such as for example the speed of sound or the viscosity or the specific heat of the fluid mixture contained in said sample.

The system 10 also comprises control means 14, configured to modify at least one parameter of the injection of said compound. The modification of said at least one parameter therefore makes it possible to modify the concentration of hydrogen injected into the mixture of natural gas. For example, said at least one parameter is the injected hydrogen flow rate and / or the hydrogen injection pressure. Said control means 14 can be controlled by a qualified operator, or even automatically.

In addition, the system 10 comprises a processing device 15 configured to perform, for the mixture of natural gas considered here and on the basis of measurements acquired thanks to said acquisition means 13, processing aimed at determining said model of calculation of the concentration of hydrogen injected into said mixture of natural gas, by implementing part of the steps of a method for determining said calculation model.

FIG. 2 schematically represents an example of the hardware architecture of the processing device 15 according to the invention.

As illustrated by FIG. 2, the processing device 15 has the hardware architecture of a computer. Thus, such a processing device 15 comprises, in particular, a processor 1_T, a random access memory 2_T, a read only memory 3_T and a non-volatile memory 4_T. It also has 5_T communication means.

The read only memory 3 of the processing device 15 constitutes a recording medium according to the invention, readable by the processor 1_T and on which is recorded a computer program PROG_l according to the invention, comprising instructions for the execution of steps of the determination method according to the invention. The program PROG_1 defines functional modules of the processing device 15, which are based on or control the hardware elements 2_T to 5_T of the processing device 15 mentioned above, and which include in particular:

- a MOD_OBT obtaining module, configured to obtain the measurements acquired for each of the samples taken,

- a MOD_REG regression module, configured to perform a regression of the measurements acquired, so as to determine a model for calculating the concentration of hydrogen injected into the mixture of natural gas as a function of said physical quantity considered as an explanatory variable.

The communication means 5_T integrate the obtaining module MOD_OBT which therefore allows the processing device 15 to receive the measurements acquired for each of the samples taken, after these measurements are sent by the acquisition means 13, themselves. same provided in this case with means of communication suitable for transmission. These communication means are based, in a manner known per se, on a communication interface capable of exchanging data between the acquisition means 13 and said processing device 15. No limitation is attached to the nature of the data. this communication interface, which can be wired or wireless, so as to allow the exchange of data according to any protocol known to those skilled in the art (Ethernet, Wifi, Bluetooth, 3G, 4G, 5G, Modbus, TCP -IP, analog, etc.).

Nothing excludes either that the processing device 15 obtains, via its obtaining module MOD_OBT, the measurements acquired after the latter have been stored in storage means external to the processing device 15, such as for example a dedicated database.

In general, no limitation is attached to the way in which the measurements acquired by the processing device 15 are obtained.

In a particular embodiment of the processing system 10, the processing device 15 is arranged near the sampling points P_PREL_1 and P_PREL_2, so as to limit the transfer times between the sampling means 11, 12, the acquisition means 13 and the obtaining module MOD_OBT. For example, the processing device 15 is separated by at most a few tens of meters from the sampling points P_PREL_1, P_PREL_2.

FIG. 3 represents, in the form of a flowchart, the main steps of the determination method according to the invention, as implemented by the processing system 10 of FIG. 1.

As illustrated by FIG. 3, the determination method comprises initially a first sampling step E10, at a first instant and in a local volume of the mixture of natural gas flowing in the annex pipe CAN_A , from a first sample. Said first sample is taken at the sampling point P_PREL_1.

Said first step E10 is implemented by said first sampling means 11. The determination method then comprises a second sampling step E20, at a second instant subsequent to said first instant and in a local volume of the mixture of natural gas, of a second sample. Said second sample is taken at the level of the sampling point P_PREL_2, the local volume in which said corresponding second sample is taken, taking into account the flow of the fluid mixture, to the local volume in which the first sample was taken.

Said second step E20 is implemented by said second sampling means 12.

As explained above, the second sampling instant is linked to the first sampling instant due to the correspondence between the local volumes within which the samples are taken.

Purely by way of illustration, the natural gas flow rate is substantially equal to 50 normal liters per hour at a pressure of the order of 3 barg, which implies a maximum transit time of natural gas of the order of 10 minutes in an annex pipeline reaching 100 meters. From such data, a person skilled in the art is able to determine said first and second sampling times in order to obtain a correspondence between the local volumes within which the sample taking is taken.

The determination method further comprises a step E30 of acquiring, for each of the samples taken (first sample and second sample), a measurement of said given physical quantity as well as a measurement of the concentration of hydrogen injected into the mixture of natural gas.

Said step E30 is implemented by said acquisition means 13.

In a particular embodiment, the acquisition, for a given sample, of the measurement of the physical quantity as well as of the measurement of the concentration of injected hydrogen is carried out as soon as said sample is taken. In other words, in this particular mode, there is a delay between the acquisition of the measurements for the first and second samples given that the latter are taken at distinct respective instants.

As an alternative, the acquisition of the measurements of the first and second samples is carried out substantially simultaneously. To do this, it is necessary for the taking of each of the samples to be completed, the acquisitions of measurements within said samples then being synchronized.

In a preferred embodiment (not shown in the figures) of the determination method, each sample taken is brought to a given temperature so that the acquisition of the measurement of the physical quantity for said sample is carried out. at said temperature. The implementation of such arrangements is for example carried out during a step following the taking of a sample, for example by reheating the sample taken, and therefore before the measurement of the physical quantity is acquired for this sample.

The fact of thus controlling the temperature at which the acquisition of the measurement of the physical quantity is carried out is particularly advantageous. Indeed, the temperature is a factor which can have an influence on the quality of the measurement of the physical quantity. This is due, in particular, to the fact that the sensor (s) dedicated to this measurement are sensitive to temperature, and that it is therefore preferable to work over a high sensitivity range of these sensors.

Also, an optimum temperature to which it is advisable to bring a sample taken can be defined for the physical quantity considered. By way of non-limiting example, when this physical quantity is thermal conductivity (respectively the speed of sound), said optimum temperature is between 70 ° C (degrees Celsius) and 80 ° C, or even possibly higher, such as for example of the order of 200 ° C (respectively said optimum temperature is ambient temperature).

In general, a person skilled in the art knows the temperature range or ranges in which a measurement should be carried out for a given physical quantity, so as to be able to take advantage of the sensitivity of the sensor (s) used for carry out this measurement.

According to the invention, the first sampling step E10, the second sampling step E20 and the acquisition step E30 are iterated a given number of times, at least one injection parameter of said compound being modified to each iteration.

The modification of said at least one injection parameter is implemented by the control means 14.

The fact of iterating these steps (first sampling step E10, second sampling step E20, acquisition step E30) makes it possible to obtain a set of measurements from which the model for calculating the hydrogen concentration injected into the mixture of natural gas can be determined, as described in detail below. In addition, the larger this set of measurements (i.e. the greater the number of measurements acquired), the more robust the calculation model.

In practice, the number of measurements acquired at the end of said iterations is defined by a sampling step fixed for said at least one injection parameter, the modification of said at least parameter being carried out at each iteration in accordance with said step sampling. In addition to taking into account said sampling step, the number of measurements acquired at the end of said iterations is also defined by a threshold at which must remain lower than said at least one parameter, this threshold being defined so as to respect the integrity of the routing network.

By way of non-limiting example, said at least one parameter is the rate of injection of hydrogen into the mixture of natural gas at the level of the injection point P_INJ_1, and the fixed sampling interval is of the order of a few percent (ie the injection flow rate is increased by said sampling step at each iteration of the set of steps formed by the first sampling step E10, the second sampling step E20 and the step E30 acquisition).

Of course, one or more other parameters, such as for example the injection pressure, can be taken into consideration, possibly in addition to said injection rate.

Once all the measurements have been acquired following said iterations, the determination method comprises a step E40 for obtaining the measurements acquired at the end of said iterations.

Said step E40 is implemented by the obtaining module MOD_OBT equipping the processing device 15.

[0100] Once the measurements have been obtained by the processing device 15, said method comprises a step E50 for regressing said measurements so as to determine a model for calculating the concentration of hydrogen injected into the mixture of natural gas as a function of said physical quantity considered as explanatory variable.

[0101] Thus, the determined calculation model corresponds to a function F so that the concentration C_H of hydrogen injected into the mixture of natural gas is expressed in the form F (G), where G denotes the physical quantity considered here . In other words, determining F amounts to performing a regression of the explained variable C_H as a function of the explanatory variable G.

In general, any data regression method known to those skilled in the art can be implemented, and the choice of a particular method constitutes only one variant of implementation of the invention. For example, the determination of the computational model can be carried out by linear regression, polynomial regression, etc. According to another example, the determination of the calculation model can be carried out by supervised learning, for example via the implementation of a support vector machine algorithm also known under the name “SVM” (acronym of the English expression -saxonne "Support Vector Machine") and using a kernel of a given type (Gaussian, polynomial, etc.). FIG. 4 represents an example of a table comprising measurements acquired for a first mixture of given natural gas, said measurements corresponding to the thermal conductivity and to the concentration of injected hydrogen.

More particularly, in the example of FIG. 4, the mixture of natural gas considered comes from Russia. The injected hydrogen concentration measurements are expressed as a molar percentage (“mol%”) and are indicated in the first row of the table. The thermal conductivity measurements (denoted "C_T" in the table of FIG. 4) are expressed in watt per meter-kelvin ("Wm ^ .K ¹ ") and are indicated in the second row of the table.

It should be noted that in this example, all the samples taken at the sampling point P_PREL_1, that is to say before the injection point P_INJ_1, provide a thermal conductivity measurement equal to 0.0281 Wm ^ .K ¹ and a zero injected hydrogen concentration. Therefore, only the first column of the table containing numbers is used for this type of metric. Each column located beyond the first column containing numbers contains the measurements of thermal conductivity and of the concentration of injected hydrogen of a sample taken at the P_PREL_2 sampling point, that is to say after the P_INJ_1 injection point. .

The third row of the table, for its part, evaluates, for each column, the difference ("DIFF") between the measurement of thermal conductivity contained in this column, and the measurement of thermal conductivity corresponding to a concentration of hydrogen. injected zero (ie the concentration contained in the first column of the table containing numbers).

From reading this table, it can be deduced that there is indeed a relationship between the concentration of hydrogen injected into the natural gas medium and the thermal conductivity. More particularly, it can be seen that the more said concentration increases, the more the thermal conductivity increases.

[0108] Consequently, on the basis of such measurements, the determination method according to the invention makes it possible to determine the associated calculation model.

FIG. 5 represents another example of a table comprising measurements acquired for a second mixture of given natural gas, said measurements corresponding to the thermal conductivity and to the concentration of injected hydrogen.

[0110] The example of FIG. 5 differs from that of FIG. 4 in that the mixture of natural gas this time comes from the Netherlands.

[0111] It should be noted that in this example, all the samples taken at the P_PREL_1 sampling point, that is to say before the P_INJ_1 injection point, provide a measurement of thermal conductivity equal to 0.0273 Wm ^ .K ¹ and a zero concentration of injected hydrogen.

Again, it can be seen that the more said concentration increases, the more the thermal conductivity increases.

FIG. 6 represents another example of a table comprising measurements acquired for said first mixture of given natural gas, said measurements corresponding to the speed of sound (denoted "V_S" in the table of FIG. 6) and to the concentration of injected hydrogen.

The example of FIG. 6 differs from that of FIG. 4 in that the physical quantity considered here is the speed of sound and the measurements of this quantity are expressed in meters per second ("ms ¹ ").

It should be noted that in this example, all the samples taken at the sampling point P_PREL_1, that is to say before the injection point P_INJ_1, provide a sound speed measurement equal to 414 ms ¹ and zero concentration of injected hydrogen.

It can be seen that the more said concentration increases, the more the speed of sound increases.

FIG. 7 represents another example of a table comprising measurements acquired for said second given mixture of natural gas, said measurements corresponding to the speed of sound and to the concentration of injected hydrogen.

The example of FIG. 7 differs from that of FIG. 6 in that the mixture of natural gas corresponds to that coming from the Netherlands.

It should be noted that in this example, all the samples taken at the sampling point P_PREL_1, that is to say before the injection point P_INJ_1, provide a sound speed measurement equal to 398 ms ¹ and zero concentration of injected hydrogen.

Again, it can be seen that the more said concentration increases, the more the speed of sound increases.

The invention has been described so far in the context of determining a model for calculating the concentration of hydrogen injected into the mixture of natural gas. As a reminder, it was considered that the composition of the natural gas mixture was predefined. Consequently, it should be understood that the calculation model thus determined is valid for a mixture of natural gas in accordance with said predefined composition. In other words, for a mixture of natural gas having a different composition, it is necessary to determine a new calculation model.

In any event, when such a calculation model has been determined, it can be used, according to another aspect of the invention, to generate a table of values of the quantity physics associated with said calculation model. Such a table can for example take the form of a table (abacus) comprising two lines identical, in terms of the types of data contained in these lines, to that described above for the tables of FIGS. 4 to 7. Such a table of values is for example stored (memorized) on a dedicated recording medium.

[0124] Such a calculation model can also be used, according to another aspect of the invention, by a control system 20 to determine the concentration of hydrogen injected into a mixture of natural gas corresponding to that for which said model of gas. calculation has been determined, and thus allow the variation of this concentration to be followed precisely over time.

[0124] FIG. 8 schematically represents, in its environment, a particular embodiment of a control system 20 according to the invention.

[0125] Said control system 20 is described here, in no way limiting, in the context of use on an urban delivery network intended for the flow of a mixture of natural gas of the same type as that considered in figure 1.

This routing network includes a main CAN_P_2 pipe. An injection point P_INJ_2 is arranged in said main CAN_P_2 pipe, the hydrogen being injected continuously, for a determined period, into the mixture of natural gas at this injection point P_INJ_2.

The system 20 comprises sampling means 21, configured to take a sample of the natural gas mixture at a sampling point P_PREL_3 of the main CAN_P_2 pipe. Said sampling point P_PREL_3 is located downstream, in the direction of flow of the mixture of natural gas within the main CAN_P_2 pipe, from said injection point P_INJ_2.

The system 20 also comprises acquisition means 22, configured to acquire a measurement of said physical quantity for said sample taken.

The technical characteristics described above with reference to the sampling means 11, 12 and to the acquisition means 13 of the processing system 10 also apply here for the sampling means 21 and the acquisition means 22 of the processing system. control system 20.

Of course, nothing excludes considering that the routing network of Figure 8 is the same as that considered for Figure 1.

[0131] In addition, the system 20 comprises a control device 23 configured to perform, for the mixture of natural gas considered here and on the basis of a measurement acquired thanks to said acquisition means 22, processing aimed at determining the hydrogen concentration injected into said mixture of natural gas, by implementing part of the steps of a control process.

[0132] FIG. 9 schematically represents an example of the hardware architecture of the control device 23 according to the invention.

[0133] As illustrated by FIG. 9, the control device 23 has the hardware architecture of a computer. Thus, such a control device 23 comprises, in particular, a processor 1_C, a random access memory 2_C, a read only memory 3_C and a non-volatile memory 4_C. It also has 5_C means of communication.

The read-only memory 3_C of the control device 23 constitutes a recording medium in accordance with the invention, readable by the processor 1_C and on which is recorded a computer program PROG_2 in accordance with the invention, comprising instructions for the execution of steps of the determination method according to the invention. The program PROG_2 defines functional modules of the control device 23, which are based on or control the hardware elements 2_C to 5_C of the control device 23 mentioned above, and which include in particular:

a first module for obtaining MOD_OBT_l, configured to obtain said calculation model determined for said mixture of natural gas,

a second module for obtaining MOD_OBT_2, configured to obtain said measurement acquired by the acquisition means 22,

a MOD_DET determination module, configured to determine the concentration of hydrogen injected by applying the calculation model to said acquired measurement.

The technical characteristics described above with reference to the communication means 5_T (and therefore fine to the obtaining module MOD_OBT) of the processing device 15 also apply here for the communication means 5_C of the control device 23 ( and therefore ultimately for the first and second MOD_OBT_l obtaining modules,

MOD_OBT_2).

[0136] It should also be noted that nothing excludes the first module for obtaining

MOD_OBT_l and the second obtaining module MOD_OBT_2 are one and the same entity.

In a particular embodiment of the control system 20, the control device 23 is arranged near the sampling point P_PREL_3, so as to limit the transfer times between the sampling means 21, the acquisition means 22 and the second module for obtaining MOD_OBT_2.

FIG. 10 represents, in the form of a flowchart, a particular mode of implementation of the control method according to the invention, such as it is implemented by the control system 20 of FIG. 8. As illustrated by FIG. 10, the method of monitoring first of all a step F10 for obtaining, by the monitoring device 23, the calculation model determined for said mixture of natural gas.

Said step F10 is implemented by the first obtaining module MOD_OBT_l equipping the control device 23.

The control method also comprises a step F20 of taking a sample of the mixture of natural gas at the level of the sampling point P_PREL_3.

Said step F20 is implemented by the sampling means 21.

The control method then comprises a step F30 of acquiring a measurement of said physical quantity for said sample taken.

Said step F30 is implemented by the acquisition means 22.

The control method then comprises a step F40 for obtaining, by the control device 23, said acquired measurement.

Said step F40 is implemented by the second obtaining module MOD_OBT_2 equipping the control device 23.

Finally, the control method then comprises a step F50 determination, by the control device 23, of said concentration by applying the calculation model to said acquired measurement.

Said step F50 is implemented by the determination module MOD_DET equipping the control device 23.

It will of course be understood that it is not compulsory for the step F10 for obtaining the calculation model to be carried out as the first step of the control method. This step F10 can in fact be carried out at any time before the implementation of step F50 for determining the concentration of hydrogen injected into the mixture of natural gas.

[0150] In addition, nothing excludes either, according to a more particular embodiment (not illustrated in the figures), to consider that the steps F20, F30, F40 and F50 are iterated recurrently so to carry out a recurrent check of the concentration of injected hydrogen.

Advantageously, the iteration frequency of the sampling, acquisition, obtaining and determination steps is of the order of a second, or even less than one second, so as to obtain a substantially continuous control of the concentration of the compound injected into the fluid mixture.

Claims

Claims [Claim 1] Method for determining a model for calculating the concentration of a compound injected into a fluid mixture of predefined type and flowing in a network for conveying said fluid mixture, said method comprising: - a first step (E10) of sampling, at a first instant and in a local volume of a fluid mixture corresponding to the predefined type flowing in a pipe (CAN_A) of a network for conveying said fluid mixture, of a first sample of said fluid mixture, said first sample being taken at a sampling point (P_PREL_1) of said pipe located upstream, in the direction of flow of the fluid mixture within the pipe, from a point d 'injection (P_INJ_1) of said compound into said fluid mixture, - a second step (E20) of sampling, at a second instant subsequent to said first instant and in a local volume of the fluid mixture, of a second sample of said fluid mixture, said second sample being taken at a sampling point ( P_PREL_2) of said pipe located downstream of said injection point, the local volume in which said second sample is taken corresponding, taking into account the flow of the fluid mixture, to the local volume in which the first sample was taken, - an acquisition step (E30), for each of the samples taken, of a measurement of a given physical quantity as well as of a measurement of the concentration of the compound injected into the fluid mixture, the first sampling step, the second sampling step and the acquisition step being iterated a given number of times, at least one injection parameter of said compound being modified at each iteration, said method further comprising steps implemented by a processing device ( 15), of which: - a step of obtaining (E40) the acquired measurements, a step of regression (E50) of the measurements acquired so as to determine a model for calculating the concentration of the compound injected into the fluid mixture as a function of said physical quantity considered as an explanatory variable. [Claim 2] A method according to claim 1, wherein, for each sample taken, the physical quantity is the thermal conductivity of the fluid mixture contained in said sample. [Claim 3] A method according to claim 1, wherein, for each sample taken, the physical quantity is the speed of sound in the fluid mixture contained in said sample. [Claim 4] A method according to claim 1, wherein, for each sample taken, the physical quantity is the viscosity of the fluid mixture contained in said sample. [Claim 5] A method according to claim 1, wherein, for each sample taken, the physical quantity is the specific heat of the fluid mixture contained in said sample. [Claim 6] A method according to any one of claims 1 to 5, wherein each sample taken is brought to a given temperature so that the acquisition of the measurement of the physical magnitude for said sample is performed at said temperature. [Claim 7] A method according to any of claims 1 to 6, wherein the fluid mixture is natural gas, and said compound is dihydrogen. [Claim 8] A method of monitoring the concentration of a compound injected into a fluid mixture flowing from a network for conveying said fluid mixture, said monitoring method being implemented by a monitoring system (20) comprising a control device (23), and comprising steps of: - obtaining (F10), by the control device, of a model for calculating said concentration determined for said fluid mixture according to any one of claims 1 to 7, so that said calculation model admits a physical quantity given as explanatory variable, - sampling (F20) of a sample of the fluid mixture at a sampling point (P_PREL_3) of a pipe (CAN_P_2) of said delivery network, said sampling point being located downstream, in the direction of flow of the fluid mixture within the pipe, from an injection point (P_INJ_2) of said compound into said fluid mixture, - acquisition (F30) of a measurement of said physical quantity for said sample taken, - obtaining (F40), by the control device, of said acquired measurement, - determination (F50), by the control device, of said concentration by applying the calculation model to said acquired measurement. [Claim 9] The method of claim 8, wherein the steps of sampling (F20), acquiring (F30), obtaining (F40) and determining (F50) are iterated repeatedly. [Claim 10] Computer program (PROG_l, PROG_2) comprising instructions for the implementation: - Steps for obtaining (E40) the acquired measurements and determining (E50) the calculation model of a determination method according to any one of claims 1 to 7; Where - steps of obtaining (F10) of the calculation model, of obtaining (F40) of the acquired measurement and of determination (F50) of the concentration of a control method according to any one of claims 8 and 9, when said program is executed by a computer. [Claim 11] A computer readable recording medium on which is recorded a computer program according to claim 10. [Claim 12] Processing system (10) for determining a model for calculating the concentration of a compound injected into a fluid mixture of predefined type and flowing in a network for conveying said fluid mixture, said system treatment comprising: - first sampling means (11), configured to take, at a first instant and in a local volume, a fluid mixture corresponding to the predefined type flowing in a pipe (CAN_A) of a network for conveying said mixture fluid, a first sample of said fluid mixture, said first sample being taken at a sampling point (P_PREL_1) of said pipe located upstream, in the direction of flow of the fluid mixture within the pipe, of a injection point (P_INJ_1) of said compound into said fluid mixture, - second sampling means (12), configured to take, at a second instant subsequent to said first instant and in a local volume of the fluid mixture, a second sample of said fluid mixture, said second sample being taken at a point of sampling (P_PREL_2) from said pipe located downstream of said injection point, the local volume in which said second sample is taken corresponding, taking into account the flow of the fluid mixture, to the local volume in which the first sample was taken , - acquisition means (13), configured to acquire, for each of the samples taken, a measurement of a given physical quantity as well as a measurement of the concentration of the compound injected into the fluid mixture, said treatment system further comprising control means (14), configured to modify at least one parameter of the injection of said compound, as well as a treatment device (15) comprising: - an obtaining module (MOD_OBT), configured to obtain the acquired measurements, a regression module (MOD_REG), configured to perform a regression of the measurements acquired, so as to determine a model for calculating the concentration of the compound injected into the fluid mixture as a function of said physical quantity considered as an explanatory variable. [Claim 13] A system for monitoring (20) the concentration of a compound injected into a fluid mixture flowing from a network for conveying said fluid mixture, said monitoring system comprising a monitoring device (23) comprising a first obtaining module (MOD_OBT_l), configured to obtain a model for calculating said concentration determined for said fluid mixture according to any one of claims 1 to 7, so that said model of calculation admits a physical quantity given as explanatory variable, said control system also comprising: - sampling means (21), configured to take a sample of the fluid mixture at a sampling point (P_PREL_3) of a pipe (CAN_P_2) of said delivery network, said sampling point being located downstream, in the direction of flow of the fluid mixture within the pipe, from an injection point (P_INJ_2) of said compound into said fluid mixture, - acquisition means (22), configured to acquire a measurement of said physical quantity for said sample taken, and said control device also comprising: - a second obtaining module (MOD_OBT_2), configured to obtain said acquired measurement, a determination module (MOD_DET), configured to determine said concentration by applying the calculation model to said acquired measurement. [Claim 14] Recording medium comprising a table comprising values of a given physical quantity, said values having been generated by means of a calculation model determined according to any one of claims 1 to 7.