[go: up one dir, main page]

WO2014027592A1 - 物質の圧力、温度、ひずみ分布測定システム、これを用いた二酸化炭素地中貯留の監視方法、二酸化炭素注入による地層安定性への影響評価方法、および結氷監視方法 - Google Patents

物質の圧力、温度、ひずみ分布測定システム、これを用いた二酸化炭素地中貯留の監視方法、二酸化炭素注入による地層安定性への影響評価方法、および結氷監視方法 Download PDF

Info

Publication number
WO2014027592A1
WO2014027592A1 PCT/JP2013/071371 JP2013071371W WO2014027592A1 WO 2014027592 A1 WO2014027592 A1 WO 2014027592A1 JP 2013071371 W JP2013071371 W JP 2013071371W WO 2014027592 A1 WO2014027592 A1 WO 2014027592A1
Authority
WO
WIPO (PCT)
Prior art keywords
optical fiber
temperature
strain
pressure
carbon dioxide
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2013/071371
Other languages
English (en)
French (fr)
Inventor
自求 薛
良昭 山内
欣増 岸田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Neubrex Co Ltd
Research Institute of Innovative Technology for the Earth RITE
Original Assignee
Neubrex Co Ltd
Research Institute of Innovative Technology for the Earth RITE
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Neubrex Co Ltd, Research Institute of Innovative Technology for the Earth RITE filed Critical Neubrex Co Ltd
Priority to CN201380043930.9A priority Critical patent/CN104583730B/zh
Priority to US14/417,962 priority patent/US9829352B2/en
Publication of WO2014027592A1 publication Critical patent/WO2014027592A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J19/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
    • G01B11/18—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge using photoelastic elements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • G01D5/35338—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using other arrangements than interferometer arrangements
    • G01D5/35354—Sensor working in reflection
    • G01D5/35358—Sensor working in reflection using backscattering to detect the measured quantity
    • G01D5/35361—Sensor working in reflection using backscattering to detect the measured quantity using elastic backscattering to detect the measured quantity, e.g. using Rayleigh backscattering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/12—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in colour, translucency or reflectance
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
    • G01K11/322—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres using Brillouin scattering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00—Measuring force or stress, in general
    • G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
    • G01L1/242—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet the material being an optical fibre
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L19/00—Details of, or accessories for, apparatus for measuring steady or quasi-steady pressure of a fluent medium insofar as such details or accessories are not special to particular types of pressure gauges

Definitions

  • the present invention utilizes a Brillouin frequency shift phenomenon and a Rayleigh frequency shift phenomenon of an optical fiber, and simultaneously measures a material pressure, temperature, and strain distribution, and uses this system to monitor a wide range of material properties such as a formation. It relates to a method of measuring.
  • Patent Document 1 Various measurement methods using the Brillouin scattering phenomenon of an optical fiber are known (for example, Patent Document 1).
  • One example is a distributed pressure sensor that uses a Brillouin frequency shift that is generated by applying strain to an optical fiber. Since the Brillouin frequency shift depends on the strain applied to the optical fiber, the applied pressure can be measured by measuring the frequency shift of the optical fiber fixed to the material deformed by the pressure.
  • the present inventors have previously proposed a system that mainly measures pressure and temperature distributions using the Rayleigh frequency shift phenomenon in addition to the Brillouin frequency shift of an optical fiber (see Patent Document 2).
  • the purpose of this system is to measure the distribution of pressure and temperature, and since the optical fiber is not fixed to the object to be measured, the strain measured with this is useless.
  • the pressure measurement technique using this optical fiber can be applied to the measurement of the volume change of an object.
  • porous sandstone changes in volume before and after being filled with a liquid, and thus becomes one of the application fields of the aforementioned pressure measurement technique.
  • the pressure measurement technology described above is used to store carbon dioxide in sandstone when carbon dioxide underground storage is executed. It is possible to contribute to the construction of a system that monitors the mechanical stability and safety of the cap rock layer (such as pelitic rock) that is the upper layer.
  • the present invention has been made in view of the above problems, and provides a system that can simultaneously measure the pressure, temperature, and strain distribution of a substance and accurately monitor and evaluate the state of a wide range of substances such as the ground.
  • the purpose is to do.
  • a measurement system includes a scattered wave acquisition unit that acquires a scattered wave obtained by scattering pulsed laser light incident on an optical fiber in or along a substance, and the Brillouin from the scattered wave.
  • Brillouin frequency shift measurement unit that measures the distribution of frequency shift in an optical fiber
  • Rayleigh frequency shift measurement unit that measures the distribution of Rayleigh frequency shift in an optical fiber from scattered waves
  • material pressure, temperature and strain and Brillouin
  • a coefficient storage unit for storing a coefficient peculiar to the installed optical fiber, a distribution of the Brillouin frequency shift measured in the Brillouin frequency shift measurement unit, and a Rayleigh frequency shift measurement unit for relating the frequency shift and the Rayleigh frequency shift.
  • Distribution of the Rayleigh frequency shift measured in step 1 and the coefficient storage By using the coefficients, the pressure of the material at the time of the measurement, the temperature, and strain is a distribution along the optical fiber that includes an analysis section for analyzing.
  • the pressure, temperature and strain distribution of a substance can be accurately measured simultaneously, and a system for accurately monitoring and evaluating the state of a wide range of substances such as the ground can be obtained.
  • FIG. FIG. 1 is a cross-sectional view showing an outline of a monitoring system for carbon dioxide underground storage using a pressure, temperature, and strain distribution measurement system for substances according to Embodiment 1 of the present invention
  • FIG. FIG. 3 is an enlarged view of the sensor cable of FIG. It is assumed that a sandstone layer 100 serving as a carbon dioxide reservoir exists in the ground, and a cap lock layer 150 functioning as a seal layer is present thereon.
  • the injection well 3a is installed from the storage site 40 set up on the ground toward the underground sandstone layer 100.
  • the press-fit well 3a is provided with a cylindrical casing 31a into which a carbon dioxide injection tube 32 is inserted.
  • the casing 31a is fixed to the underground stratum by applying cementing 34 to the periphery.
  • An observation well 3b for observing underground conditions is often installed around the injection well 3a.
  • the observation well 3b is usually provided with a cylindrical casing 31b similar to that installed in the press-fit well 3a, and various observation sensors are inserted in the casing 31b by filling water.
  • the sensor cable 2a is embedded in the layer of the cementing 34 in order to measure the distribution of the underground pressure P, temperature T and formation strain ⁇ along the injection well 3a.
  • the sensor cable 2b may be embedded in the cementing layer along the observation well 3b.
  • the sensor cable 2a and the sensor cable 2b will be described as the sensor cable 2.
  • FIG. 3 shows an example of a cross-sectional structure of the sensor cable 2.
  • the sensor cable 2 includes a first optical fiber 21 that is affected by pressure and a second optical fiber 22 that is cut off from the influence of pressure.
  • the second optical fiber 22 is accommodated in the metal thin tube 24 for pressure blocking.
  • a protective cover 23 may be provided around the first optical fiber 21.
  • the protective cover 23 needs to be made of a material / structure that allows the first optical fiber 21 to be affected by ambient pressure and deformation.
  • the first optical fiber and the metal thin tube 24 constitute the sensor cable 2 as a stranded wire together with a plurality of metal wires 25, for example.
  • the first optical fiber in order to measure the strain ⁇ of the formation, the first optical fiber must be fixed to the cementing 34 layer. Fixing may be performed on the entire surface in the longitudinal direction of the first optical fiber, or an interval of about several meters may be provided.
  • the sensor cable 2 Since the sensor cable 2 is embedded in the layer of the cementing 34, it is affected when a volume change occurs in the surrounding formation. For example, when the formation is deformed by carbon dioxide sealing, the sensor cable 2 is deformed integrally with the cementing 34 layer. In this case, the first optical fiber 21 receives the pressure and senses the pressure, but the second optical fiber 22 accommodated in the metal thin tube 24 is not affected.
  • the measuring apparatus 1 installed on the ground, and distribution along the optical fiber of Brillouin frequency shift and Rayleigh frequency shift is performed. Is required. From these Brillouin frequency shift and Rayleigh frequency shift distributions, the distribution of pressure, temperature, and strain along the sensor cable 2 can be known simultaneously. Therefore, the inventors named the measuring device 1 as DPTSS (Distributed Pressure Temperature Strain System) 1.
  • DPTSS Distributed Pressure Temperature Strain System
  • the Brillouin scattering phenomenon is a phenomenon in which power moves through an acoustic phonon of an optical fiber when light is incident on the optical fiber.
  • the frequency difference between the incident light and the Brillouin scattered light is called the Brillouin frequency, which is proportional to the speed of sound in the optical fiber, and the speed of sound depends on the strain and temperature of the optical fiber.
  • the strain and / or temperature applied to the optical fiber can be measured.
  • the present inventors have confirmed that the Brillouin frequency changes depending on the pressure applied to the optical fiber.
  • the change in Brillouin frequency is referred to as Brillouin frequency shift.
  • the Rayleigh scattering phenomenon is a scattering phenomenon that occurs because light is scattered due to the fluctuation of the refractive index in the optical fiber.
  • the frequency difference between the incident light and the Rayleigh scattered light is the Rayleigh frequency.
  • This Rayleigh frequency also varies with strain and / or temperature applied to the optical fiber.
  • the change in the Rayleigh frequency is called a Rayleigh frequency shift.
  • the Rayleigh scattering phenomenon has been considered to be sensitive only to strain and temperature.
  • the system is proposed on the assumption that the Rayleigh scattering phenomenon is sensitive only to strain and temperature.
  • the Rayleigh scattering phenomenon has sensitivity not only to strain and temperature but also to pressure, similarly to the Brillouin scattering phenomenon. That is, the Brillouin frequency shift ⁇ B and the Rayleigh frequency shift ⁇ R can be expressed as in the equations (1) and (2) by the pressure change amount ⁇ P, the temperature change amount ⁇ T, and the strain change amount ⁇ , respectively.
  • C ij is a coefficient specific to the optical fiber, and for the optical fiber to be used, by obtaining the values of these coefficients by a preliminary test or the like, the pressure change amount ⁇ P, the temperature change amount ⁇ T, The distribution of the strain change amount ⁇ can be obtained. By introducing a pressure term into this Rayleigh frequency shift ⁇ R , pressure, temperature, and strain distribution measurements with higher accuracy became possible.
  • the superscript number indicates the type of optical fiber. Since the pressure and temperature are the pressure and temperature of the field where the optical fiber exists, the two types of optical fibers have the same value. On the other hand, the value of strain depends on whether the optical fiber is fixed to the surrounding material. Since DPTSS needs to measure the strain of the material surrounding the fiber, at least one optical fiber must be fixed to the surrounding material.
  • the effects of the pressure P, temperature T and strain ⁇ can be separated by solving the simultaneous equations of the above equation (3). Therefore, by performing hybrid measurement of measurement of Brillouin frequency shift (referred to as Brillouin measurement) and measurement of Rayleigh frequency shift (referred to as Rayleigh measurement) and solving the simultaneous equations of Equation (3), the pressure change ⁇ P, temperature The distribution along the optical fiber of the change amount ⁇ T and the strain change amount ⁇ can be obtained.
  • Brillouin measurement measurement of measurement of Brillouin frequency shift
  • Rayleigh measurement measurement of Rayleigh frequency shift
  • the strain ⁇ 1 received by the first optical fiber fixed to the surrounding substance is different from the strain ⁇ 2 received by the second optical fiber housed in the metal thin tube.
  • ⁇ P, ⁇ T, ⁇ 1 , and ⁇ 2 there are also four equations, so these four unknowns can be obtained.
  • a useful value for the strain is the strain ⁇ 1 of the first optical fiber that is directly subjected to the strain of the surrounding material.
  • Each coefficient C ij in the equation (4) of the first optical fiber and the second optical fiber is obtained in advance by a preliminary test, etc., hybrid measurement of Brillouin measurement and Rayleigh measurement is performed, and the simultaneous equations of the above equation (4) are obtained.
  • the distribution of the pressure change amount ⁇ P, the temperature change amount ⁇ T, and the strain change amount ⁇ along the optical fiber can be obtained.
  • Hybrid measurement of Brillouin measurement and Rayleigh measurement can be performed simultaneously at a certain point in time, so not only the distribution of one-dimensional pressure change ⁇ P, temperature change ⁇ T and strain change ⁇ along the optical fiber, but also the time axis Can also be obtained.
  • equation (4) is an incremental equation. That is, in order to obtain the Brillouin frequency shift and Rayleigh frequency shift on the left side, two measurements are required: a reference measurement in the initial state and a main measurement after the state changes. Further, what is obtained by solving the equation (4) is the amount of change in each of pressure, temperature, and strain based on the initial state. When absolute amounts of pressure, temperature, and strain are required, the absolute amounts of the respective distributions of pressure, temperature, and strain at the time of initial measurement are measured by some method.
  • the initial state can be selected arbitrarily.
  • initial measurement can be performed in a constant temperature room on the ground before installing a cable in a well (press-in well 3a or observation well 3b).
  • a state in which the pressure and temperature distribution are uniform and constant can be set as the initial state.
  • the amount of change in pressure, temperature, and strain due to carbon dioxide injection can be directly obtained by solving Equation (4).
  • the absolute amount distribution of pressure and temperature measured with an electric sensor or the like when the well is stopped before carbon dioxide injection may be used. If there is initial measurement data in a constant temperature room on the ground, the absolute quantity distribution of pressure and temperature can be obtained from the measurement before carbon dioxide injection.
  • Changes in the absolute amount of the distribution of pressure P, temperature T, and strain ⁇ , or collection of data on the amount of change can be monitored to monitor changes and distribution accompanying the inclusion of carbon dioxide in the sandstone layer 100. And the leakage from the cap lock layer 150 can be monitored.
  • FIG. 4 is a block diagram showing an outline of an example of DPTSS1.
  • the scattered wave acquisition unit 11 acquires the scattered wave scattered by the optical fiber.
  • the acquired scattered wave is analyzed by the Brillouin frequency shift measuring unit 12 to measure the Brillouin frequency shift.
  • the Brillouin frequency shift is measured as a distribution along the length direction of the optical fiber.
  • the Rayleigh frequency shift measurement unit 13 measures the Rayleigh frequency shift.
  • the Rayleigh frequency shift is also measured as a distribution along the length of the optical fiber.
  • the coefficient storage unit 14 stores in advance the coefficient C ij of Expression (4) obtained by a preliminary test or the like.
  • the analysis unit 15 uses the equation (4) to calculate the pressure change ⁇ P, the temperature change ⁇ T, and the strain change ⁇ . Is analyzed and stored in the distribution data storage unit 16.
  • the above measurement and analysis are executed at predetermined time intervals, and are stored in the distribution data storage unit 16 as distribution data of changes in pressure, temperature, and strain for each time.
  • the evaluation calculation unit 17 evaluates the state of the sandstone layer 100 based on the pressure, temperature, strain variation over time, and monitors the storage state of carbon dioxide, for example.
  • the initial measurement 1 is performed under constant pressure and temperature conditions in a temperature-controlled room,
  • the Brillouin reference spectrum and the Rayleigh reference spectrum which are the reference for the Brillouin frequency shift and the Rayleigh frequency shift, are measured (ST2).
  • the sensor cable 2 is installed in the injection well 3a or the observation well 3b up to the ground of the sandstone layer 100 as a carbon dioxide reservoir (ST3).
  • the Brillouin reference spectrum and the Rayleigh reference spectrum which are the reference for the Brillouin frequency shift and the Rayleigh frequency shift, are measured as the initial measurement 2 (ST4).
  • the Brillouin frequency shift and the Rayleigh frequency shift are obtained from the measurement data of the initial measurement 1 and the initial measurement 2, and the pressure change ⁇ P and the temperature change using the simultaneous equations (4)
  • the distribution of the quantity ⁇ T is calculated.
  • the absolute quantity distribution of the pressure and temperature in the initial measurement 2 is obtained using the pressure and temperature of the temperature-controlled room in the initial measurement 1.
  • FIG. 6 shows a conceptual example of measurement data at a certain point in time.
  • FIG. 6 shows a case where an optical fiber is installed up to a depth of 1000 m.
  • distribution data as a change amount or an absolute amount of pressure P, temperature T, formation strain ⁇ in the depth direction is obtained by hybrid measurement of Brillouin measurement and Rayleigh measurement.
  • FIG. 7 shows conceptual data of the time change of the strain ⁇ at a certain depth of the sandstone layer obtained while injecting carbon dioxide.
  • time 0 is the start point of carbon dioxide injection.
  • the strain gradually increases, and the strain gradually reaches an equilibrium state. It can be seen that the amount of carbon dioxide in the sandstone layer at that position was saturated when the strain reached an equilibrium state.
  • the carbon dioxide storage part 101 shown in FIG. 1 expands by continuing injection
  • FIG. 8 is a diagram showing an outline of the configuration of the laboratory experiment.
  • An optical fiber 200 is spirally wound around a cylindrical sample having a biased porosity (permeability) called Tago sandstone 110, and Brillouin frequency shift and Rayleigh frequency shift are measured by DPTSS1.
  • 8 is a portion simulating the cap rock layer 150 in FIG. 1 and a portion simulating the sandstone layer 100 having a large porosity in the lower layer. Consists of.
  • the sample was stored in a pressure vessel, and after applying a sealing pressure of 12 MPa, water and carbon dioxide were injected into the sample to evaluate the change in the state of the sample. Since the pressure and temperature are uniform and have no distribution, the strain distribution was mainly obtained to evaluate the change in the state of the sample. In the case of an indoor experiment, the pressure and temperature are uniform and constant, and monitoring is performed using another point sensor. Therefore, the unknown amount is only distortion, and it is not necessary to use all four equations of equation (4), and the optical fiber can be measured with one first optical fiber.
  • FIG. 9 is a diagram illustrating an example of a result of observing how water permeated into the sample Tago sand 110 is replaced with carbon dioxide.
  • Strain distribution data along the optical fiber 200 spirally wound around the outer circumference of the mulberry sandstone 110 is acquired at predetermined time intervals. At each time point, the same data as the strain distribution data of FIG. 6 can be acquired.
  • FIG. 9 is a diagram in which this data is arranged with time as the horizontal axis, and the magnitude of the distortion is represented by shades of color. The dark part shows the part with large distortion.
  • the time when carbon dioxide was started to be injected is 0, and it is observed that the strain of the coarse layer 112 having a large porosity increases almost 100 hours later.
  • the color intensity indicates the amount of carbon dioxide.
  • the dark portion is a portion where carbon dioxide is well infiltrated, and the change in the darkness shows that carbon dioxide penetrates into the Tahu sandstone 110.
  • the laboratory experiment shown in the second embodiment it is possible to measure the change in the strain distribution of the sample talu sandstone 110 using the pressure, temperature, and strain distribution measurement system of the substance of the present invention. I found out. Further, by evaluating the measurement result, for example, the state of underground storage of carbon dioxide can be monitored. In laboratory experiments, pressure and temperature have almost no distribution, and changes can be monitored in other ways. However, in underground such as underground storage of carbon dioxide, it is necessary to know the distribution of pressure and temperature.
  • Embodiment 3 shows an application example that can be realized by the pressure, temperature, and strain distribution measurement system of the substance of the present invention.
  • Application example 1 In the first and second embodiments, the monitoring at the time of injection in the underground storage of carbon dioxide has been described. With this system, it is possible to monitor abnormalities after completion of injection by monitoring strain and the like after completion of injection. For example, it is conceivable that carbon dioxide leaks from the carbon dioxide reservoir 101 through the cap lock layer 150 due to reasons such as a crack occurring in the cap lock layer.
  • FIG. 10 shows a flow chart when monitoring leakage from the carbon dioxide storage unit 101.
  • the strain after the end of injection is monitored in two dimensions in the time and depth directions as shown in FIG. 11 (ST10). If there is no change in strain (ST11 NO), monitoring is continued. If the strain changes from a certain point in time as in the example of FIG. 11 (ST11 YES), it can be monitored that the stored carbon dioxide may have leaked, and it is determined whether repair is necessary (ST13). . If it is determined that repair is necessary (YES in ST13) and the leaked location can be identified, repairing that part (ST14) will restore the original strain as shown in FIG. (ST10) Yes.
  • FIG. 12 shows a flowchart of the process for monitoring the phase state
  • FIG. 13 shows a conceptual diagram of the temperature state.
  • the temperature distribution in the depth direction is monitored over time (ST20). If there is a temperature rise (ST21 YES), there is a possibility that the portion of the carbon dioxide where the temperature rise has changed from a supercritical state to a liquid (ST22) (temperature rise zone in FIG. 13).
  • the carbon dioxide in the portion where the temperature has dropped may change from a liquid to a gas (ST24) (temperature drop zone in FIG. 13).
  • ST24 temperature drop zone in FIG. 13
  • This change may occur particularly in a portion shallower than the cap lock layer, so it is necessary to monitor the portion above the cap lock layer. Since the temperature distribution from the ground surface to the sandstone layer 100 can always be measured according to the system of the present invention, such monitoring is also possible.
  • FIG. 14 shows a flow chart of the process of evaluating the change in the shape of the ground surface based on the deformation of the formation, that is, the influence of the carbon dioxide injection on the formation stability.
  • FIG. 14 shows the measured strain and the displacement calculated from the measured strain.
  • FIG. 15 shows a conceptual diagram of the quantity distribution.
  • a measured value of strain distribution at a certain time is acquired (ST30). By integrating this strain value in the depth direction, a distribution of displacement is obtained (ST31).
  • the deformation amount of the ground surface is obtained from the distribution of the displacement amount (ST32).
  • ST33 By comparing / analyzing the obtained ground surface deformation amount and the actual ground surface deformation amount (ST33), the cause of the shape change of the ground surface can be evaluated.
  • Embodiment 4 FIG. The embodiment so far has been an embodiment related to the state monitoring of the underground formation such as carbon dioxide underground storage.
  • the fourth embodiment is an embodiment of a system for monitoring icing such as rivers as an embodiment other than the monitoring of the state of the formation.
  • FIG. 16 shows a schematic diagram of a system for monitoring icing of a bridge over a river.
  • the sensor cable 2 is installed along the bridge girder and the pier.
  • the sensor cable 2 has a cross-sectional structure similar to that shown in FIG.
  • the sensor cable 2 is installed particularly at a position where icing is to be monitored, for example, the pier portion is below the water surface.
  • FIG. 17 shows a conceptual diagram of a change in temperature and strain over time and a change in volume elastic modulus obtained from strain when icing occurs in a portion of the sensor cable 2 below the water surface.
  • FIG. 17 by monitoring the changes in strain and bulk modulus as the temperature changes, it is possible to monitor how the river freezes. River water contains salt, mud and other contaminants, and the icing temperature is not necessarily 0 ° C. For this reason, this system which can monitor ice formation itself by observing strain other than temperature becomes effective.
  • the distribution along the optical fiber of the pressure, temperature, and strain of the substance can be measured simultaneously with only one measurement system. Since the change over time can be measured, the state of a wide range of substances can be monitored accurately. In particular, the effect of the present invention is great in monitoring the state of a wide range of substances, ie, the length of the optical fiber is as long as 100 m or more.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Organic Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)
  • Optical Transform (AREA)
  • Physical Or Chemical Processes And Apparatus (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Description

物質の圧力、温度、ひずみ分布測定システム、これを用いた二酸化炭素地中貯留の監視方法、二酸化炭素注入による地層安定性への影響評価方法、および結氷監視方法
　本発明は、光ファイバのブリルアン周波数シフトおよびレイリー周波数シフト現象を利用して、物質の圧力、温度、ひずみの分布を同時に計測するシステム、およびこのシステムを利用して地層など広範囲の物質特性を監視、計測する方法に関するものである。
　光ファイバのブリルアン散乱現象を利用した各種の計測手法が知られている（例えば特許文献１）。その一つとして、光ファイバにひずみが加えられることにより生じるブリルアン周波数シフトを利用した分布型圧力センサが挙げられる。ブリルアン周波数シフトは、光ファイバに加わるひずみに依存するため、圧力により変形する物質に固定した光ファイバの周波数シフトを計測することで、印加された圧力を計測することができる。
　本発明者らは、先に、光ファイバのブリルアン周波数シフトに加えてレイリー周波数シフト現象を利用して、主として圧力と温度の分布を測定するシステムを提案している（特許文献２参照）。このシステムでは、圧力と温度の分布を計測するのが目的であり、光ファイバが被測定物に固定されていないために、これで計測されるひずみは使い道がない。
　この光ファイバを用いた圧力測定技術は、物体の体積変化の計測に適用し得る。例えば、ポーラスな砂岩は、液体が充填される前と充填された後とでは体積が変化するので、前述の圧力測定技術の適用分野の一つとなる。近年、地球温暖化対策として、地中に二酸化炭素を貯留する技術が開発されつつあるが、前述の圧力測定技術は、二酸化炭素地中貯留が実行される場合における、砂岩における二酸化炭素の貯留状況をモニタするシステム、ならびにその上位層であるキャップロック層（泥質岩等）の力学安定性や安全性をモニタするシステムの構築に寄与できる。
国際公開第２００６／００１０７１号 特開２０１０－２１６８７７号公報
　しかし、例えば地中に存在する地層の状態変化を的確に検知する方法は未だ提案されていない。電気的な圧力センサを用いれば、スポット的な圧力変化を検知することは可能である。しかし、その圧力変化が、地表で観測される変状とどのような関連性を有するのか、また、地表に変状が生じても、その力学的安定性を保てるのかは明らかにされていない。
　本発明は上記の課題に鑑みてなされたもので、物質の圧力、温度、およびひずみの分布を同時に計測することができ、地中など広範囲の物質の状態を的確に監視・評価するシステムを提供することを目的とする。
　この発明に係る測定システムは、物質中または物質に沿って敷設された光ファイバに入射されたパルスレーザ光が光ファイバ内で散乱された散乱波を取得する散乱波取得部と、散乱波からブリルアン周波数シフトの光ファイバ内の分布を計測するブリルアン周波数シフト計測部と、散乱波からレイリー周波数シフトの光ファイバ内の分布を計測するレイリー周波数シフト計測部と、物質の圧力、温度およびひずみと、ブリルアン周波数シフト、およびレイリー周波数シフトを関係付けるための、敷設された光ファイバ特有の係数を記憶する係数記憶部と、ブリルアン周波数シフト計測部において計測されたブリルアン周波数シフトの分布と、レイリー周波数シフト計測部において計測されたレイリー周波数シフトの分布と、係数記憶部に記憶された係数とを用いて、計測した時点における物質の圧力、温度、およびひずみの、光ファイバに沿った分布を解析する解析部を備えたものである。
　この発明によれば、物質の圧力、温度、およびひずみの分布を同時に的確に計測することができ、地中など広範囲の物質の状態を的確に監視・評価するシステムが得られる。
この発明の実施の形態１による物質の圧力、温度、ひずみ分布測定システムを利用した、二酸化炭素地中貯留の監視システムの概要を示す断面図である。 この発明の実施の形態１による物質の圧力、温度、ひずみ分布測定システムを利用した、二酸化炭素地中貯留の監視システムの図１のＦ部分の拡大図である。 この発明の実施の形態１による物質の圧力、温度、ひずみ分布測定システムのセンサケーブルの拡大断面図である。 この発明の実施の形態１による物質の圧力、温度、ひずみ分布測定システムのＤＰＴＳＳの構成の一例を示すブロック図である。 この発明の実施の形態１による物質の圧力、温度、ひずみ分布測定システムの二酸化炭素の貯留状態の評価に関する工程の一例を示すフロー図である。 この発明の実施の形態１による物質の圧力、温度、ひずみ分布測定システムの測定結果を概念的に示す図である。 この発明の実施の形態１による物質の圧力、温度、ひずみ分布測定システムの測定結果のうちひずみの経時変化を概念的に示す図である。 この発明の物質の圧力、温度、ひずみ分布測定システムの室内実験の構成の概要を示すブロック図である。 この発明の物質の圧力、温度、ひずみ分布測定システムの室内実験の結果の一例を示す図である。 この発明の物質の圧力、温度、ひずみ分布測定システムの応用例としての、二酸化炭素貯留部からの漏れを監視する工程のフロー図である。 この発明の物質の圧力、温度、ひずみ分布測定システムの応用例としての、二酸化炭素貯留部からの漏れの監視結果を概念的に示す図である。 この発明の物質の圧力、温度、ひずみ分布測定システムの応用例としての、地中の二酸化炭素の相変化を監視する工程のフロー図である。 この発明の物質の圧力、温度、ひずみ分布測定システムの応用例としての、地中の二酸化炭素の相変化の監視結果を概念的に示す図である。 この発明の物質の圧力、温度、ひずみ分布測定システムの応用例としての、地表の形状変化を評価する工程のフロー図である。 この発明の物質の圧力、温度、ひずみ分布測定システムの応用例としての、地表の形状変化の評価結果を概念的に示す図である。 この発明の実施の形態４による物質の圧力、温度、ひずみ分布測定システムを利用した、河川の結氷を監視するシステムの概要を示す図である。 この発明の実施の形態４による物質の圧力、温度、ひずみ分布測定システムを利用した、河川の結氷の監視結果を概念的に示す図である。
実施の形態１．
　図１は、本発明の実施の形態１による物質の圧力、温度、ひずみ分布測定システムを利用した、二酸化炭素地中貯留の監視システムの概要を示す断面図、図２は図１のＦ部分の拡大図、図３は、図１のセンサケーブルの拡大断面図である。地中に二酸化炭素の貯留層となる砂岩層１００が存在し、その上にシール層として機能するキャップロック層１５０が存在するものとする。地上に設営された貯留サイト４０から地下の砂岩層１００に向けて、圧入井３ａが設置される。圧入井３ａには、内部に二酸化炭素の注入チューブ３２が挿入された円筒状のケーシング３１ａが設置される。ケーシング３１ａは、周囲にセメンチング３４が施されて地中地層に固定化される。圧入井３ａの周辺には、地中の状態を観測する観測井３ｂが設置されることが多い。この観測井３ｂには、通常、圧入井３ａに設置されるのと同様な円筒状のケーシング３１ｂが設置され、ケーシング３１ｂ内に水を満たして各種観測センサが挿入される。
　圧入井３ａに沿って地中の圧力Ｐ、温度Ｔ、および地層のひずみεの分布を計測するために、セメンチング３４の層内にセンサケーブル２ａが埋設される。あるいは、観測井３ｂに沿って、セメンチングの層内にセンサケーブル２ｂが埋設されても良い。以降、センサケーブル２ａ、センサケーブル２ｂは、センサケーブル２として説明する。図３は、センサケーブル２の断面構造の一例を示している。センサケーブル２は、圧力の影響を受ける第１光ファイバ２１と、圧力の影響から遮断された第２光ファイバ２２とを含む。第２光ファイバ２２は、圧力遮断のため、金属細管２４内に収納されている。第１光ファイバ２１の周囲には保護カバー２３を設けても良い。ただし、この保護カバー２３は第１光ファイバ２１が周囲の圧力や変形の影響を受けるような材質・構造となっている必要がある。第１光ファイバおよび金属細管２４は、例えば複数の金属線２５と共に撚り線としてセンサケーブル２を構成する。また、地層のひずみεを計測するために、第１光ファイバはセメンチング３４層に固定されなければならない。固定は第１光ファイバの長手方向に渡って全面でも良いし、数ｍ程度の間隔を設けても良い。
　センサケーブル２は、セメンチング３４の層内に埋設されているので、周囲の地層に体積変化を生じた場合、その影響を受ける。例えば、地層が二酸化炭素の封入により変形した場合、センサケーブル２は、セメンチング３４の層と一体となって変形を受けることになる。この場合、第１光ファイバ２１はその圧力を受け、圧力を感知するが、金属細管２４内に収容されている第２光ファイバ２２は影響を受けない。
　このような第１光ファイバ２１および第２光ファイバ２２の各々について、地上に設置された測定装置１によりブリルアン計測およびレイリー計測が行われ、ブリルアン周波数シフトおよびレイリー周波数シフトの光ファイバに沿った分布が求められる。これら、ブリルアン周波数シフトおよびレイリー周波数シフトの分布から、センサケーブル２に沿った圧力、温度、ひずみの分布を同時に知ることができる。よって、発明者らは、測定装置１をＤＰＴＳＳ（Distributed　Pressure　Temperature　Strain　System）１と名付けた。
　ここで、光ファイバを利用した圧力、温度、ひずみ分布測定の、測定原理を説明する。光ファイバに光を入射し、その散乱光を周波数分析すると、入射光とほぼ同じ周波数をもつレイリー散乱光、入射光と大きく周波数が異なるラマン散乱光、および入射光と数～数十ＧＨｚ程度周波数が異なるブリルアン散乱光が観測される。
　ブリルアン散乱現象は、光が光ファイバへ入射された場合に光ファイバの音響フォノンを介してパワーが移動する現象である。入射光とブリルアン散乱光との周波数差はブリルアン周波数と呼ばれ、このブリルアン周波数は、光ファイバ中の音速に比例し、そして、この音速が光ファイバのひずみおよび温度に依存する。このため、ブリルアン周波数の変化を測定することによって、光ファイバに加えられるひずみおよび／または温度を測定することができる。また、本発明者らにより、光ファイバに加えられる圧力によっても、ブリルアン周波数が変化することが確かめられている。ここでは、ブリルアン周波数の変化を、ブリルアン周波数シフトと呼ぶことにする。
　レイリー散乱現象は、光ファイバ中の屈折率のゆらぎによって、光が散乱するために生ずる散乱現象である。入射光とレイリー散乱光との周波数差がレイリー周波数である。このレイリー周波数も、光ファイバに加えられるひずみおよび／または温度により変化する。ここでは、レイリー周波数の変化を、レイリー周波数シフトと呼ぶことにする。
　従来、レイリー散乱現象は、ひずみと温度のみに感度があると考えられてきた。本発明者らの、先の提案である特許文献２においては、レイリー散乱現象は、ひずみと温度のみに感度があるとして、システムの提案を行っている。本発明者らのその後の研究の結果，レイリー散乱現象もブリルアン散乱現象と同様、ひずみと温度以外に、圧力にも感度を有することが判明した。すなわち、ブリルアン周波数シフトΔνＢおよびレイリー周波数シフトΔνＲは、それぞれ圧力変化量ΔＰ、温度変化量ΔＴ、ひずみ変化量Δεにより式（１）および式（２）のように表記できる。
Figure JPOXMLDOC01-appb-M000001
Figure JPOXMLDOC01-appb-M000002
　ここで、Ｃｉｊは光ファイバ固有の係数であり、使用する光ファイバについて、予備試験などによりこれらの係数の値を求めておくことにより、下記のように圧力変化量ΔＰ、温度変化量ΔＴ、ひずみ変化量Δεの分布を求めることができる。このレイリー周波数シフトΔνＲに圧力の項を導入することで、より精度が高い圧力、温度、ひずみ分布測定が可能となったのである。
　ΔνＢ、ΔνＲが計測されたとする。計測値において圧力Ｐ、温度Ｔおよびひずみεの影響を分離するためには、３個以上の独立した計測量が必要となる。一つの光ファイバでは、独立した２個の計測値、ΔνＢ、ΔνＲが得られるだけであるから、圧力Ｐ、ひずみεおよび温度Ｔに対する感度が異なる２種類の光ファイバを用いることで、４個の独立した計測値が得られる。すなわち、式（３）の連立方程式が得られる。
Figure JPOXMLDOC01-appb-M000003
ここで、上付きの数字は、光ファイバの種類を示している。圧力および温度は光ファイバが存在する部分の場の圧力および温度であるため、２種類の光ファイバで同一の値を持つ。一方、ひずみの値は光ファイバが周囲の物質に固定されているか否かに依存する。ＤＰＴＳＳではファイバ周囲の物質のひずみを計測する必要があるので、少なくとも一つの光ファイバは周囲の物質に固定されていなければならない。
　上記式（３）の連立方程式を解くことにより、圧力Ｐ、温度Ｔおよびひずみεの影響を分離することができる。よって、ブリルアン周波数シフトの計測（ブリルアン計測と呼ぶ）、およびレイリー周波数シフトの計測（レイリー計測と呼ぶ）のハイブリッド計測を行い、式（３）の連立方程式を解くことにより、圧力変化量ΔＰ、温度変化量ΔＴおよびひずみ変化量Δεの、光ファイバに沿った分布を求めることができる。
　式（３）において、上付き数字１のファイバを図１～３における第１光ファイバ２１、上付き数字２のファイバを第２光ファイバ２２とすると、第２光ファイバ２２は圧力の影響から遮断されているため、式（３）は簡素化され、式（４）のようになる。
Figure JPOXMLDOC01-appb-M000004
　式（４）でも、圧力および温度は光ファイバが存在する部分の場の圧力および温度であるため、２種類の光ファイバで同一の値を持つ。一方、ひずみについては、周囲の物質に固定された第１光ファイバが受けるひずみε１と、金属細管に納められた第２光ファイバが受けるひずみε２は異なる。未知数はΔＰ、ΔＴ、Δε１、Δε２の４個となるが、方程式も４個あるため、これら４個の未知数を求めることができる。ただし、ひずみとして有用な値は、周囲の物質のひずみを直接受ける第１光ファイバのひずみε１である。
　予め予備試験などにより、第１光ファイバおよび第２光ファイバの式（４）における各係数Ｃｉｊを求めておき、ブリルアン計測およびレイリー計測のハイブリッド計測を行い、上記式（４）の連立方程式を解くことにより、圧力変化量ΔＰ、温度変化量ΔＴおよびひずみ変化量Δεの、光ファイバに沿った分布を求めることができる。ブリルアン計測およびレイリー計測のハイブリッド計測は、ある時点で同時に行うことができるため、光ファイバに沿った１次元の圧力変化量ΔＰ、温度変化量ΔＴおよびひずみ変化量Δεの分布だけではなく、時間軸のデータも得ることができる。
　式（４）は増分型の式であることに注意が必要である。すなわち、左辺のブリルアン周波数シフトおよびレイリー周波数シフトを求めるためには、初期の状態における基準となる計測と、状態が変化した後の本計測の２回の計測が必要となる。また、式（４）を解いて得られるのは、初期状態を基準とした圧力、温度、ひずみのそれぞれの変化量となる。圧力、温度、ひずみの絶対量が必要となる場合は、初期計測時における圧力、温度、ひずみのそれぞれの分布の絶対量を何らかの方法により測定しておく。
　初期状態は任意に選ぶことができる。二酸化炭素地中貯留の監視について言えば、坑井（圧入井３ａや観測井３ｂ）にケーブルを設置する前に、地上の恒温室において初期計測を行うことができる。この場合、圧力、温度分布が一様一定の状態を初期状態とすることができる。
　または、坑井にケーブルを設置した後、二酸化炭素を注入する前を初期状態とすることもできる。この場合は、式（４）を解くことで、二酸化炭素の注入による圧力、温度、ひずみの変化量を直接得ることができる。圧力、温度、ひずみの絶対量が必要となる場合は、二酸化炭素注入前の坑井静止時に電気センサ等で計測された圧力、温度の絶対量分布を用いればよい。地上の恒温室における初期計測データがあれば、二酸化炭素注入前の計測から、圧力、温度の絶対量分布を得ることもできる。
　圧力Ｐ、温度Ｔおよびひずみεの分布の絶対量の変化あるいは、変化量のデータを収集することで、砂岩層１００への二酸化炭素の封入に伴う変化および分布をモニタリングでき、例えば二酸化炭素の封入の状況や、キャップロック層１５０からの漏れなどを監視することができる。
　図４は、ＤＰＴＳＳ１の一例の概要を示すブロック図である。散乱波取得部１１において、光ファイバで散乱された散乱波を取得する。取得された散乱波はブリルアン周波数シフト計測部１２において解析され、ブリルアン周波数シフトが計測される。このとき、ブリルアン周波数シフトは光ファイバの長さ方向に沿った分布として計測される。同様に、レイリー周波数シフト計測部１３においてレイリー周波数シフトが計測される。レイリー周波数シフトも、光ファイバの長さ方向に沿った分布として計測される。
　係数記憶部１４には、予備試験などにより求めた式（４）の係数Ｃｉｊが予め記憶されている。計測されたブリルアン周波数シフトおよびレイリー周波数シフト、および係数記憶部１４に記憶されている係数を用いて、解析部１５において式（４）により、圧力変化量ΔＰ、温度変化量ΔＴおよびひずみ変化量Δεが解析され、分布データ記憶部１６に記憶される。以上の計測、解析は所定の時間間隔で実行され、時間毎の圧力・温度・ひずみの変化量の分布データとして分布データ記憶部１６に記憶される。初期計測において、圧力・温度・ひずみの初期分布の絶対量が計測されている場合は、この値を分布データ記憶部１６に記憶しておけば、変化量の分布データと合わせて各時点の絶対量の分布データが得られる。評価演算部１７では、圧力、温度、ひずみの時間変化量などにより砂岩層１００における状態を評価し、例えば二酸化炭素の貯留状況などを監視する。
　図１のシステムにおける、二酸化炭素の貯留状態の評価に関する工程の一例を図５のフロー図に示す。まず、第１光ファイバ２１および第２光ファイバ２２として設置される２種類の光ファイバを準備し、室内試験などにより、各光ファイバの特性を測定し、式（４）の各係数Ｃｉｊを決定しておく（ＳＴ１）。決定した各係数は、例えばＤＰＴＳＳ１の係数記憶部１４に記憶させておく。各係数が決定できた２種類の光ファイバをセンサケーブル２の第１光ファイバ２１、第２光ファイバ２２として、まず恒温室において、圧力、温度一様一定の条件下で初期計測１を行い、ブリルアン周波数シフトおよびレイリー周波数シフトの基準となるブリルアン基準スペクトルおよびレイリー基準スペクトルを測定する（ＳＴ２）。次に、図１～図３に示すような構成で、センサケーブル２を二酸化炭素貯留層となる砂岩層１００の地中まで、圧入井３ａ、あるいは観測井３ｂに設置する（ＳＴ３）。
　センサケーブル２の設置の完了後、初期計測２として、ブリルアン周波数シフトおよびレイリー周波数シフトの基準となるブリルアン基準スペクトルおよびレイリー基準スペクトルを測定する（ＳＴ４）。圧力・温度の絶対量が必要となる場合は、初期計測１と初期計測２の計測データから、ブリルアン周波数シフトおよびレイリー周波数シフトを求め、連立方程式（４）を用いて圧力変化量ΔＰ、温度変化量ΔＴの分布を算出する。続いて、初期計測１における恒温室の圧力、温度を用いて、初期計測２における圧力、温度の絶対量分布を求める。
　二酸化炭素の注入が開始されたら、ブリルアンスペクトルおよびレイリースペクトルの計測を行い、初期計測２の計測データとの差分を取ることにより、ブリルアン周波数シフトΔν１ Ｂ、Δν２ Ｂ（ＳＴ５）およびレイリー周波数シフトΔν１ Ｒ、Δν２ Ｒ（ＳＴ６）を求め、連立方程式（４）を用いて圧力変化量ΔＰ、温度変化量ΔＴおよびひずみ変化量Δεの分布を算出する（ＳＴ７）。ＳＴ５、ＳＴ６、ＳＴ７は、上述のようにＤＰＴＳＳ１において、所定の時間間隔で実行され（ＳＴ８　ＮＯ）、時間毎のデータとして分布データ記憶部１６に記憶される。ここで算出される圧力変化量ΔＰ、温度変化量ΔＴおよびひずみ変化量Δεは、前述の初期計測２からの各変化量である。必要な時間分布が取得できたら（ＳＴ８　ＹＥＳ）、空隙率、二酸化炭素の浸入速度等の岩石性状のデータベースを参照して（ＳＴ１０）、後述のように、二酸化炭素貯留状態の評価を行うことができる（ＳＴ９）。
　図６に、ある時点での測定データの概念的な一例を示す。図６では、深さ１０００ｍまで光ファイバを設置した場合を示している。図６のように、ある時点においてブリルアン計測およびレイリー計測のハイブリッド計測により、深さ方向の、圧力Ｐ、温度Ｔ、地層のひずみεの変化量、あるいは絶対量としての分布データが得られる。
　所定時間間隔で図６のようなデータを取得することにより、各位置での、圧力Ｐ、温度Ｔ、ひずみεの時間変化のデータが得られる。図７に、二酸化炭素を注入しながら取得した、砂岩層の、ある深さにおけるひずみεの時間変化の概念的データを示す。図７において時間０が二酸化炭素注入開始時点である。図７に示すように、注入開始後ある程度時間が経過した後、徐々にひずみが増加し、次第にひずみが平衡状態になる。ひずみが平衡状態になった時点で、その位置の砂岩層の二酸化炭素量が飽和したことがわかる。さらに注入を続けることで、図１に示す二酸化炭素貯留部１０１が拡大してゆく。このようにしてひずみ分布の変化を監視することで、二酸化炭素の注入および貯留の状態を評価・監視することができる。
実施の形態２．
　実施の形態２では、本発明の物質の圧力、温度、ひずみ分布測定システムにより地層など広範囲の物質特性を監視、計測することができることを、室内実験により実証した例を示す。図８は、室内実験の構成の概要を示す図である。多胡砂岩１１０という、空隙率（浸透率）に偏りがある円柱状のサンプルに光ファイバ２００を螺旋状に巻き、ＤＰＴＳＳ１により、ブリルアン周波数シフトおよびレイリー周波数シフトを測定する構成としている。多胡砂岩１１０は、図８の上部（fine　layer）の空隙率が小さい、すなわち図１のキャップロック層１５０を模擬する部分と、下部（coarse　layer）の空隙率が大きい砂岩層１００を模擬した部分から成る。
　このサンプルを圧力容器内に収納し、１２ＭＰａの封圧を印可した後、水や二酸化炭素をサンプルに注入して、サンプルの状態の変化を評価した。圧力や温度は一様であり分布を有しないため、主にひずみの分布を求めて、サンプルの状態の変化を評価した。また、室内実験の場合、圧力と温度は一様一定であるので他の点センサでモニタを行う。よって未知量はひずみだけとなり、式（４）の４個の式全てを用いる必要が無く、光ファイバは第１光ファイバ１本で計測できる。
　まず、サンプルに水を注入し水が浸透する過程を観察した。水が十分に浸透した後、次に二酸化炭素をサンプルに注入し、二酸化炭素が水と置換する過程を観察した。図９は、サンプルの多胡砂岩１１０内に浸透した水が二酸化炭素に置き換わる様子を観察した結果の一例を示す図である。多胡砂岩１１０の外周に螺旋状に巻いた光ファイバ２００に沿ったひずみの分布データを所定時間間隔で取得する。各時点において、図６のひずみの分布データと同様のデータが取得できる。このデータを、時間を横軸にして並べ、ひずみの大きさを色の濃淡で表わした図が図９である。色の濃い部分が、ひずみが大きい部分を示している。
　図９において、二酸化炭素を注入し始めた時間が０であり、ほぼ１００時間後から空隙率が大きいcoarse　layer１１２の部分のひずみが増加してゆくのが観測されている。また、色の濃さが二酸化炭素の量を示している。色が濃い部分は二酸化炭素が良く浸透している部分であり、また色の濃さの変化が、二酸化炭素が多胡砂岩１１０に浸透する様子を示している。
　このように、本実施の形態２に示した室内実験により、本発明の物質の圧力、温度、ひずみ分布測定システムを用いてサンプルである多胡砂岩１１０のひずみ分布の変化を測定することが可能であることがわかった。また、この測定結果を評価することにより、例えば二酸化炭素の地中貯留の状況を監視することができる。室内実験では、圧力や温度は分布がほとんどなく、またその変化も他の方法で監視できる。しかし、二酸化炭素の地中貯留のような地中においては、圧力や温度の分布も知る必要がある。第１光ファイバ２１および第２光ファイバ２２を用いた本発明の物質の圧力、温度、ひずみ分布測定システムにより、式（４）を用いて、圧力および温度の分布も同時に測定することで、地中のひずみ分布のデータを得ることができ、二酸化炭素の地中貯留の状態を監視することができるのである。
実施の形態３．
　実施の形態３は、本発明の物質の圧力、温度、ひずみ分布測定システムにより可能となる応用例を示す。
＜応用例１＞
　実施の形態１および２では、二酸化炭素の地中貯留における注入時の監視について説明した。本システムにより、注入終了後もひずみなどを監視することで、注入終了後の異常を監視することができる。例えば、キャップロック層に亀裂が生じたなどの理由により二酸化炭素が二酸化炭素貯留部１０１からキャップロック層１５０を介して漏れ出すことが考えられる。
　二酸化炭素貯留部１０１からの漏れを監視するときのフロー図を図１０に示す。注入終了後のひずみを、図１１のように時間と深さ方向の２次元で監視しておく（ＳＴ１０）。ひずみに変化がなければ（ＳＴ１１　ＮＯ）監視を続ける。図１１の例のように、ある時点からひずみが変化（ＳＴ１１　ＹＥＳ）すれば、貯留している二酸化炭素が漏れ出した可能性があることがモニタでき、修復が必要かどうか判断する（ＳＴ１３）。修復が必要と判断され（ＳＴ１３　ＹＥＳ）、漏れ出した場所が特定できれば、その部分の修復（ＳＴ１４）を行うことで、図１１に示すように、ひずみが元に戻り、修復できたこともモニタ（ＳＴ１０）できる。
＜応用例２＞
　地中の二酸化炭素は相が変化し、液体や気体、あるいは超臨界の状態になる。地中の温度変化を監視することでこれら相状態の変化をモニタすることができる。この相状態をモニタする工程のフロー図を図１２に、温度の状態の概念図を図１３に示す。深さ方向の温度分布を経時的に監視しておく（ＳＴ２０）。温度上昇があれば（ＳＴ２１　ＹＥＳ）、温度上昇が生じている部分の二酸化炭素は超臨界の状態から液体に変化している可能性がある（ＳＴ２２）（図１３の温度上昇ゾーン）。逆に温度低下が生じていれば（ＳＴ２３　ＹＥＳ）、温度低下している部分の二酸化炭素は液体から気体に変化している可能性がある（ＳＴ２４）（図１３の温度低下ゾーン）。このように、深さ方向の温度分布の経時的な変化を監視することで、地中の二酸化炭素の状態の変化をモニタすることができる。
　この変化は、特にキャップロック層より浅い部分で生じる可能性があるため、キャップロック層より上の部分を監視する必要がある。本発明のシステムによれば、地表から砂岩層１００までの温度分布を常時計測できるため、このような監視も可能である。
＜応用例３＞
　地中のひずみ分布から、地層の変状を評価することができる。地層の変状に基づいて地表の形状変化の評価、すなわち二酸化炭素注入による地層安定性への影響評価を行う工程のフロー図を図１４に、ひずみ測定値と、このひずみ測定値から算出した変位量の分布の概念図を図１５に示す。ある時点におけるひずみ分布の測定値を取得する（ＳＴ３０）。このひずみの値を深さ方向に積分することで変位量の分布が得られる（ＳＴ３１）。変位量の分布から地表の変形量が得られる（ＳＴ３２）。得られた地表変形量と実際の地表変形量などを比較・解析する（ＳＴ３３）ことで、地表の形状変化の原因などを評価することができる。
　以上のように、本発明の物質の圧力、温度、ひずみ分布測定システムによる計測に基づけば、この一つの計測システムのみで、二酸化炭素地中貯留における、地中、地表の様々な状態を監視、評価することができる。
　さらに、二酸化炭素地中貯留のみならず、例えば石油井戸、その他地中の深度深くまで掘削して地下資源を採掘するシステムや、それらが廃坑になった後の地中の地層の状態の監視などにも、本発明の物質の圧力、温度、ひずみ分布測定システムを適用することができる。
実施の形態４．
　これまでの実施の形態は、二酸化炭素地中貯留など地中の地層の状態監視に関する実施の形態であった。本実施の形態４は、地層の状態監視以外の実施の形態として、河川などの結氷を監視するシステムについての実施の形態である。図１６に、河川に架かる橋梁の結氷を監視するシステムの模式図を示す。橋桁および橋脚に沿ってセンサケーブル２を設置する。センサケーブル２は例えば断面構造が図３と同様のものを用いる。センサケーブル２は特に結氷を監視したい位置、例えば橋脚の部分が水面以下となるように設置しておく。
　このセンサケーブル２の光ファイバについて、ＤＰＴＳＳ１によりブリルアン周波数シフト計測およびレイリー周波数シフト計測を行い、センサケーブル２に沿った圧力、温度、ひずみの分布を同時に求める。これらの分布の経時変化を監視することで、結氷の状態を知ることができる。図１７に、センサケーブル２の水面下にある部分の、結氷が生じるときの、温度とひずみの時間変化、およびひずみから求めた体積弾性率の変化の様子の概念的な図を示す。図１７のように、温度変化とともに、ひずみおよび体積弾性率の変化の様子を監視することで、河川が結氷する様子を監視することができる。河川の水は、塩分、泥、その他の混入物があり、氷結温度は必ずしも０℃ではない。このため、温度以外にひずみを観測することで結氷そのものを監視できる本システムが有効となる。
　以上のように、本発明の物質の圧力、温度、ひずみ分布測定システムによれば、一つの計測システムのみで、同時に物質の圧力、温度、ひずみの光ファイバに沿った分布が測定でき、その分布の経時変化も測定できるため、広範囲の物質の状態の監視などが的確に行える。特に、光ファイバの長さが１００ｍ以上といった長い、すなわち広範囲の物質の状態の監視において本発明の効果が大きい。
１、１ａ、１ｂ：ＤＴＰＳＳ　２、２ａ、２ｂ：センサケーブル
３ａ：圧入井　３ｂ：観測井　１１：散乱波取得部
１２：ブリルアン周波数シフト計測部
１３：レイリー周波数シフト計測部　１４：係数記憶部　１５：解析部
１６：分布データ記憶部　１７：評価演算部　２１：第１光ファイバ
２２：第２光ファイバ　２３：保護カバー　２４：金属細管
２５：金属線　３１：ケーシング　３２：注入チューブ
３４：セメンチング　４０：貯留サイト　１００：砂岩層
１０１：二酸化炭素貯留部　１１０：多胡砂岩サンプル
１１１：サンプル上部（fine　layer）
１１２：サンプル下部（coarse　layer）
１５０：キャップロック層　２００：光ファイバ

Claims (11)

  1. 　物質中または物質に沿って当該物質と共に変形するように敷設された光ファイバに入射されたパルスレーザ光が前記光ファイバ内で散乱された散乱波を取得する散乱波取得部と、
    この散乱波取得部において取得した散乱波からブリルアン周波数シフトの、前記光ファイバ内の分布を計測するブリルアン周波数シフト計測部と、
    前記散乱波取得部において取得した散乱波からレイリー周波数シフトの、前記光ファイバ内の分布を計測するレイリー周波数シフト計測部と、
    前記物質の圧力、温度、およびひずみと、前記ブリルアン周波数シフト、および前記レイリー周波数シフトを関係付けるための、前記敷設された光ファイバ特有の係数を記憶する係数記憶部と、
    前記ブリルアン周波数シフト計測部において計測されたブリルアン周波数シフトの分布と、前記レイリー周波数シフト計測部において計測されたレイリー周波数シフトの分布と、前記係数記憶部に記憶された係数とを用いて、前記計測した時点における前記物質の圧力、温度、およびひずみの、前記光ファイバに沿った分布を解析により求める解析部を備えたことを特徴とする物質の圧力、温度、ひずみ分布測定システム。
  2. 　前記敷設された光ファイバの長さは１００ｍ以上であることを特徴とする請求項１に記載の物質の圧力、温度、ひずみ分布測定システム。
  3. 　前記光ファイバは、圧力の影響を受けるように保持された第１光ファイバと、圧力の影響から遮断されるように保持された第２光ファイバとにより構成されていることを特徴とする請求項１または２に記載の物質の圧力、温度、ひずみ分布測定システム。
  4. 　前記解析部は、
    前記係数記憶部に記憶された、前記第１光ファイバと圧力、前記第１光ファイバと温度、前記第１光ファイバとひずみ、前記第２光ファイバと温度、前記第２光ファイバとひずみ、のそれぞれとブリルアン周波数シフトを関係づけるための係数である、Ｃ１ １３、Ｃ１ １２、Ｃ１ １１、Ｃ２ １２、Ｃ２ １１と、
    前記第１光ファイバと圧力、前記第１光ファイバと温度、前記第１光ファイバとひずみ、前記第２光ファイバと温度、前記第２光ファイバとひずみ、のそれぞれとレイリー周波数シフトを関係づけるための係数である、Ｃ１ ２３、Ｃ１ ２２、Ｃ１ ２１、Ｃ２ ２２、Ｃ２ ２１と、
    前記第１光ファイバにおいて計測された、初期計測時からのブリルアン周波数シフトΔν１ Ｂおよびレイリー周波数シフトΔν１ Ｒと、
    前記第２光ファイバにおいて計測された、ブリルアン周波数シフトΔν２ Ｂおよびレイリー周波数シフトΔν２ Ｒと、による連立方程式
    　　　Δν１ Ｂ＝Ｃ１ １３ΔＰ＋Ｃ１ １２ΔＴ＋Ｃ１ １１Δε１
    　　　Δν１ Ｒ＝Ｃ１ ２３ΔＰ＋Ｃ１ ２２ΔＴ＋Ｃ１ ２１Δε１
    　　　Δν２ Ｂ＝　　　　　　Ｃ２ １２ΔＴ＋Ｃ２ １１Δε２
    　　　Δν２ Ｒ＝　　　　　　Ｃ２ ２２ΔＴ＋Ｃ２ ２１Δε２
    により、前記初期計測時からの圧力変化量ΔＰ、温度変化量ΔＴ、第１光ファイバのひずみ変化量Δε１、および第２光ファイバのひずみ変化量Δε２を求めることを特徴とする請求項３に記載の物質の圧力、温度、ひずみ分布測定システム。
  5. 　前記光ファイバは二酸化炭素地中貯留のための圧入井、または観測井に沿って、地表から二酸化炭素貯留部となる砂岩層まで敷設されたことを特徴とする請求項１から４のいずれか１項に記載の物質の圧力、温度、ひずみ分布測定システム。
  6. 　請求項５に記載の物質の圧力、温度、ひずみ分布測定システムを用いて、少なくともひずみ分布の経時変化を観測することにより、前記砂岩層に注入される二酸化炭素の状態を監視することを特徴とする二酸化炭素地中貯留の監視方法。
  7. 　請求項５に記載の物質の圧力、温度、ひずみ分布測定システムを用いて、少なくともひずみ分布の経時変化を観測することにより、前記砂岩層に貯留された二酸化炭素の漏えいを監視することを特徴とする二酸化炭素地中貯留の監視方法。
  8. 　請求項５に記載の物質の圧力、温度、ひずみ分布測定システムを用いて、少なくとも温度分布の経時変化を観測することにより、地中の二酸化炭素の相変化を監視することを特徴とする二酸化炭素地中貯留の監視方法。
  9. 　請求項５に記載の物質の圧力、温度、ひずみ分布測定システムを用いて、測定されたひずみ分布から、当該ひずみを深さ方向に積分することで地中および地表の変位量を求め、この求めた変位量から地表の形状変化を評価することを特徴とする二酸化炭素注入による地層安定性への影響評価方法。
  10. 　前記光ファイバは少なくとも水中に設置された橋脚に沿って敷設されたことを特徴とする請求項１から４のいずれか１項に記載の物質の圧力、温度、ひずみ分布測定システム。
  11. 　請求項１０に記載の物質の圧力、温度、ひずみ分布測定システムを用いて、少なくともひずみと温度の経時変化を観測することにより、前記橋脚周辺の水の結氷を監視することを特徴とする結氷監視方法。
PCT/JP2013/071371 2012-08-17 2013-08-07 物質の圧力、温度、ひずみ分布測定システム、これを用いた二酸化炭素地中貯留の監視方法、二酸化炭素注入による地層安定性への影響評価方法、および結氷監視方法 Ceased WO2014027592A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201380043930.9A CN104583730B (zh) 2012-08-17 2013-08-07 物质的压力、温度、应变分布测量系统、利用该系统的二氧化碳地下储存的监视方法、二氧化碳注入对地层稳定性的影响评估方法、以及结冰监视方法
US14/417,962 US9829352B2 (en) 2012-08-17 2013-08-07 Distribution measurement system for pressure, temperature, strain of material, monitoring method for carbon dioxide geological sequestration, assessing method for impact of carbon dioxide injection on integrity of strata, and monitoring method for freezing using same

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2012-180725 2012-08-17
JP2012180725A JP5769676B2 (ja) 2012-08-17 2012-08-17 物質の圧力、温度、ひずみ分布測定システム、これを用いた二酸化炭素地中貯留の監視方法、二酸化炭素注入による地層安定性への影響評価方法、および結氷監視方法

Publications (1)

Publication Number Publication Date
WO2014027592A1 true WO2014027592A1 (ja) 2014-02-20

Family

ID=50286291

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2013/071371 Ceased WO2014027592A1 (ja) 2012-08-17 2013-08-07 物質の圧力、温度、ひずみ分布測定システム、これを用いた二酸化炭素地中貯留の監視方法、二酸化炭素注入による地層安定性への影響評価方法、および結氷監視方法

Country Status (4)

Country Link
US (1) US9829352B2 (ja)
JP (1) JP5769676B2 (ja)
CN (1) CN104583730B (ja)
WO (1) WO2014027592A1 (ja)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103983309A (zh) * 2014-06-06 2014-08-13 中国工程物理研究院流体物理研究所 一种用于冲击和爆轰实验的三合一光电探头
WO2014181617A1 (ja) * 2013-05-10 2014-11-13 公益財団法人地球環境産業技術研究機構 光ファイバケーブル、光ファイバケーブルの製造方法、および分布型測定システム
WO2017033494A1 (ja) * 2015-08-27 2017-03-02 新日鉄住金エンジニアリング株式会社 ケーブル及び光ファイバへの初期引張歪の導入方法
WO2017060971A1 (ja) * 2015-10-06 2017-04-13 ニューブレクス株式会社 Ｄｐｔｓｓケーブル
RU2700043C1 (ru) * 2016-06-08 2019-09-12 Ньюбрекс Ко., Лтд. Волоконно-оптический кабель для измерения распределений давления, температуры и деформации
JP2019178976A (ja) * 2018-03-30 2019-10-17 沖電気工業株式会社 光ファイバ歪み及び温度測定装置並びに光ファイバ歪み及び温度測定方法
WO2020027223A1 (ja) * 2018-07-31 2020-02-06 古河電気工業株式会社 ケーブル、ケーブルの形状センシングシステム、センシングシステム、ケーブル形状のセンシング方法
US11286773B2 (en) 2020-03-11 2022-03-29 Neubrex Co., Ltd. Using fiber-optic distributed sensing to optimize well spacing and completion designs for unconventional reservoirs
JP2022129025A (ja) * 2021-02-24 2022-09-05 マルイチ エアリアル エンジニア株式会社 評価方法、評価用プログラム、評価システム
JP2022545211A (ja) * 2019-08-16 2022-10-26 ペトロリアム・ナショナル・ブルバド・（ペトロナス） 光ケーブル、光ケーブル監視システム、および坑井の監視方法

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104641201B (zh) * 2012-08-10 2016-12-14 公益财团法人地球环境产业技术研究机构 物体的体积变化测量方法
CN104101379B (zh) * 2014-05-26 2017-05-31 华北电力大学（保定） 一种基于botdr的已敷设传感光纤温度和应变同时测量方法
JP6474345B2 (ja) * 2015-12-28 2019-02-27 鹿島建設株式会社 Ｐｃ鋼撚線の定着部構造、グラウンドアンカー、測定装置及び測定方法
WO2017154190A1 (ja) * 2016-03-11 2017-09-14 ニューブレクス株式会社 レイリー測定システムおよびレイリー測定方法
US10073006B2 (en) * 2016-04-15 2018-09-11 Viavi Solutions Inc. Brillouin and rayleigh distributed sensor
CN105911259B (zh) * 2016-07-07 2017-11-28 中国石油大学(华东) 判断二氧化碳注入砂岩后成岩流体演化过程的方法
EP3491218A4 (en) * 2016-07-27 2020-04-08 Services Pétroliers Schlumberger SIMULTANEOUSLY DISTRIBUTED MEASUREMENTS ON OPTICAL FIBER
JP6705353B2 (ja) * 2016-09-30 2020-06-03 沖電気工業株式会社 光ファイバ歪み及び温度測定装置
CN106872109B (zh) * 2017-01-18 2019-04-09 神华集团有限责任公司 监测方法及装置
JP6495509B2 (ja) * 2017-05-16 2019-04-03 Ｃａｃｈ株式会社 遠隔状態監視システム及び監視方法
CN107588789B (zh) * 2017-09-05 2019-11-01 华北电力大学（保定） 一种分布式光纤在变压器内部的防护方法
US10712258B2 (en) * 2017-12-06 2020-07-14 California Institute Of Technology Cuvette having high optical transmissivity and method of forming
CN108120656B (zh) * 2017-12-20 2020-05-12 辽宁工程技术大学 测量煤体吸附和解吸煤层气过程温度及热量的系统及方法
JP6951995B2 (ja) * 2018-03-22 2021-10-20 鹿島建設株式会社 グラウト材の注入装置及び注入方法
CN108442725A (zh) * 2018-03-29 2018-08-24 武汉地震工程研究院有限公司 一种多芯光纤复合编制碳纤维布式的结构加固与状态监测装置及方法
GB201808366D0 (en) * 2018-05-22 2018-07-11 Fotech Solutions Ltd Distributed optical fibre vibration sensor
CN109443425B (zh) * 2018-10-23 2024-01-30 南阳理工学院 一种长距离保温输送管道变形及泄漏监测系统
CN109405880A (zh) * 2018-10-24 2019-03-01 东北大学 一种浅部采空区地表与地下相结合监测预警系统及方法
CN110187216B (zh) * 2019-06-28 2024-09-20 中国南方电网有限责任公司超高压输电公司广州局海口分局 一种海底电缆载流量岸上模拟实验装置及方法
JP7376052B2 (ja) * 2020-06-05 2023-11-08 日本電信電話株式会社 形状測定システム及び形状測定方法
EP4266002A4 (en) * 2020-12-21 2024-08-14 Neubrex Co., Ltd. REINFORCED DSS CABLE
KR102648709B1 (ko) * 2021-02-25 2024-03-15 김상욱 가스 배관 모니터링 방법 및 시스템
CN113932944A (zh) * 2021-10-12 2022-01-14 深圳大学 一种用于监测土地内部位移、应变和温度的系统及方法
CN114812667A (zh) * 2022-03-18 2022-07-29 张鹏 一种基于多光纤的输电导线检测方法、结构及装置
CN115655909A (zh) * 2022-09-05 2023-01-31 海南浙江大学研究院 一种co2注入过程中相变地层响应的模拟装置及方法
CN116109033B (zh) * 2022-12-13 2023-08-01 青岛海洋地质研究所 一种海域沉积盆地二氧化碳地质封存适宜性评价方法
AU2024235238A1 (en) * 2023-03-10 2025-10-16 Fiber Sense Limited A fibre optic cable monitoring method and system
WO2024228359A1 (ja) * 2023-05-01 2024-11-07 日本電信電話株式会社 検出システム、検出方法、及び、解析装置
CN117109465B (zh) * 2023-08-31 2024-04-12 交通运输部天津水运工程科学研究所 一种多物理场应变传感信号解耦校准方法
CN117073767B (zh) * 2023-10-13 2024-02-09 新疆智能港环保科技有限公司 一种用于油气水识别的长效示踪监测系统及其方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006001071A1 (ja) * 2004-06-25 2006-01-05 Neubrex Co., Ltd. 分布型光ファイバセンサ
JP2006145465A (ja) * 2004-11-24 2006-06-08 Japan Atomic Energy Agency マルチ情報計測用光ファイバ
JP2010216877A (ja) * 2009-03-13 2010-09-30 Neubrex Co Ltd 分布型光ファイバ圧力センサ
JP2011053146A (ja) * 2009-09-03 2011-03-17 Neubrex Co Ltd 検知用ケーブル及びこれを備えた監視システム

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4427263A (en) * 1981-04-23 1984-01-24 The United States Of America As Represented By The Secretary Of The Navy Pressure insensitive optical fiber
GB9626099D0 (en) * 1996-12-16 1997-02-05 King S College London Distributed strain and temperature measuring system
US6788417B1 (en) * 1999-04-30 2004-09-07 The Regents Of The University Of California Optical fiber infrasound sensor
US6496265B1 (en) * 2000-02-16 2002-12-17 Airak, Inc. Fiber optic sensors and methods therefor
GB2442486B (en) * 2006-10-06 2009-01-07 Schlumberger Holdings Measuring brillouin backscatter from an optical fibre using a tracking signal
JP4441624B2 (ja) 2007-08-07 2010-03-31 日本電信電話株式会社 光ファイバを用いた歪・温度の分布測定方法及び測定装置
US7859654B2 (en) * 2008-07-17 2010-12-28 Schlumberger Technology Corporation Frequency-scanned optical time domain reflectometry
US8699009B2 (en) * 2008-11-27 2014-04-15 Neubrex Co., Ltd. Distributed optical fiber sensor
CN102589714A (zh) * 2012-02-23 2012-07-18 南昌航空大学 一种基于高压气体瑞利-布里渊散射频谱测量温度的装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006001071A1 (ja) * 2004-06-25 2006-01-05 Neubrex Co., Ltd. 分布型光ファイバセンサ
JP2006145465A (ja) * 2004-11-24 2006-06-08 Japan Atomic Energy Agency マルチ情報計測用光ファイバ
JP2010216877A (ja) * 2009-03-13 2010-09-30 Neubrex Co Ltd 分布型光ファイバ圧力センサ
JP2011053146A (ja) * 2009-09-03 2011-03-17 Neubrex Co Ltd 検知用ケーブル及びこれを備えた監視システム

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014181617A1 (ja) * 2013-05-10 2014-11-13 公益財団法人地球環境産業技術研究機構 光ファイバケーブル、光ファイバケーブルの製造方法、および分布型測定システム
JP5980419B2 (ja) * 2013-05-10 2016-08-31 公益財団法人地球環境産業技術研究機構 光ファイバケーブル、光ファイバケーブルの製造方法、および分布型測定システム
US9557196B2 (en) 2013-05-10 2017-01-31 Research Institute Of Innovative Technology For The Earth Optical fiber cable, manufacturing method for the same and distributed measurement system
CN103983309B (zh) * 2014-06-06 2016-05-25 中国工程物理研究院流体物理研究所 一种用于冲击和爆轰实验的三合一光电探头
CN103983309A (zh) * 2014-06-06 2014-08-13 中国工程物理研究院流体物理研究所 一种用于冲击和爆轰实验的三合一光电探头
WO2017033494A1 (ja) * 2015-08-27 2017-03-02 新日鉄住金エンジニアリング株式会社 ケーブル及び光ファイバへの初期引張歪の導入方法
US9958299B2 (en) 2015-08-27 2018-05-01 Nippon Steel & Sumikin Engineering Co., Ltd. Cable and method for introducing initial tensile strain to optical fiber
US10612947B2 (en) 2015-10-06 2020-04-07 Neubrex Co., Ltd. Distributed pressure, temperature, strain sensing cable using metal wires with slot grooves and optical fibers in the slot grooves
WO2017060971A1 (ja) * 2015-10-06 2017-04-13 ニューブレクス株式会社 Ｄｐｔｓｓケーブル
JPWO2017060971A1 (ja) * 2015-10-06 2018-04-05 ニューブレクス株式会社 Ｄｐｔｓｓケーブル
RU2700043C1 (ru) * 2016-06-08 2019-09-12 Ньюбрекс Ко., Лтд. Волоконно-оптический кабель для измерения распределений давления, температуры и деформации
JP2019178976A (ja) * 2018-03-30 2019-10-17 沖電気工業株式会社 光ファイバ歪み及び温度測定装置並びに光ファイバ歪み及び温度測定方法
JP7003807B2 (ja) 2018-03-30 2022-02-10 沖電気工業株式会社 光ファイバ歪み及び温度測定装置並びに光ファイバ歪み及び温度測定方法
WO2020027223A1 (ja) * 2018-07-31 2020-02-06 古河電気工業株式会社 ケーブル、ケーブルの形状センシングシステム、センシングシステム、ケーブル形状のセンシング方法
JPWO2020027223A1 (ja) * 2018-07-31 2021-09-24 古河電気工業株式会社 ケーブル、ケーブルの形状センシングシステム、センシングシステム、ケーブル形状のセンシング方法
US11519759B2 (en) 2018-07-31 2022-12-06 Furukawa Electric Co., Ltd. Cable, cable shape sensing system, sensing system, and cable shape sensing method
JP2024045538A (ja) * 2018-07-31 2024-04-02 古河電気工業株式会社 ケーブル、ケーブルの形状センシングシステム、センシングシステム、ケーブル形状のセンシング方法
JP2022545211A (ja) * 2019-08-16 2022-10-26 ペトロリアム・ナショナル・ブルバド・（ペトロナス） 光ケーブル、光ケーブル監視システム、および坑井の監視方法
JP7673051B2 (ja) 2019-08-16 2025-05-08 ペトロリアム・ナショナル・ブルバド・（ペトロナス） 光ケーブル、光ケーブル監視システム、および坑井の監視方法
US11286773B2 (en) 2020-03-11 2022-03-29 Neubrex Co., Ltd. Using fiber-optic distributed sensing to optimize well spacing and completion designs for unconventional reservoirs
JP2022129025A (ja) * 2021-02-24 2022-09-05 マルイチ エアリアル エンジニア株式会社 評価方法、評価用プログラム、評価システム

Also Published As

Publication number Publication date
CN104583730A (zh) 2015-04-29
CN104583730B (zh) 2016-12-07
JP2014038039A (ja) 2014-02-27
US9829352B2 (en) 2017-11-28
JP5769676B2 (ja) 2015-08-26
US20150211900A1 (en) 2015-07-30

Similar Documents

Publication Publication Date Title
JP5769676B2 (ja) 物質の圧力、温度、ひずみ分布測定システム、これを用いた二酸化炭素地中貯留の監視方法、二酸化炭素注入による地層安定性への影響評価方法、および結氷監視方法
CN112268642B (zh) 基于分布式光纤传感的地下应力测量装置及测量方法
US8800653B2 (en) Systems and methods for monitoring a well
US8776609B2 (en) Use of fiber optics to monitor cement quality
JP4757092B2 (ja) 地下水流動評価方法
Zhang et al. Shear deformation calculation of landslide using distributed strain sensing technology considering the coupling effect
JP5544443B2 (ja) 圧力パルス崩壊試験における不確実性減少技法
JP6351279B2 (ja) 地盤状態監視システムおよび地盤状態監視方法
Zhang et al. Microanchored borehole fiber optics allows strain profiling of the shallow subsurface
Sun et al. Fiber optic distributed sensing technology for real-time monitoring water jet tests: Implications for wellbore integrity diagnostics
CN116124025B (zh) 基于分布式光纤或阵列光栅光纤的形变监测系统及方法
Sasaki et al. Distributed fiber optic strain sensing of bending deformation of a well mockup in the laboratory
Wu et al. Experimental investigation of interfacial behavior of fiber optic cables embedded in frozen soil for in-situ deformation monitoring
Flynn et al. Instrumented concrete pile tests–part 1: a review of instrumentation and procedures
WO2016091972A1 (en) Method for ascertaining characteristics of an underground formation
US11286773B2 (en) Using fiber-optic distributed sensing to optimize well spacing and completion designs for unconventional reservoirs
Leng et al. An integrated investigation of the failure mechanism of loess landslide induced by raining: from field to laboratory
Anjana et al. Development of an optical fibre sensor system for ground displacement and pore water pressure monitoring
CN213543861U (zh) 基于分布式光纤传感的地下应力测量装置
Liu et al. Subsurface multi-physical characterization of mountain excavation and city construction in loess plateau with a fiber-optic sensing system
Viana da Fonseca et al. Improved laboratory techniques for advanced geotechnical characterization towards matching in situ properties
Minardo et al. Soil slope monitoring by use of a Brillouin distributed sensor
Shi et al. DFOS technology-based landslide monitoring: The Majiagou landslide case study (China)
Birid Comparative study of rock mass deformation modulus using different approaches
Ikeda et al. In-situ formation compaction monitoring in deep reservoirs by use of fiber optics

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13879406

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 14417962

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13879406

Country of ref document: EP

Kind code of ref document: A1