Method for measuring flow velocity distribution of groundwater using distributed optical fiber with active heating

A method for measuring flow velocity distribution of groundwater using distributed optical fiber with active heating, comprising:

• Step S1:
• Setting and optimizing field test parameters comprising:
• Step S2:
• Conducting a background temperature monitoring test, comprising performing multiple rounds of temperature measurement tests;
• monitoring background temperature of a borehole along a depth profile in each round of temperature measurement tests, with a monitoring time of more than 24 hours;
• and performing a quality check on the temperature measurement test data after each round of temperature measurement tests;
• S2.1 Optimizing distributed temperature sensing (DTS) parameters comprising:
• S2.1.1 Testing circuit design and on-site equipment testing comprising:
• S2.1.2 Laying out DTS optical fiber on-site and testing;
• comprising checking and debugging optical fiber temperature measurement accuracy and parameters, and calibrating a DTS temperature measurement system, including calibrating the length of the optical fiber to be monitored underground, surface temperature calibration duration, and position of the underground temperature probe;
• S2.2 Conducting a formal background temperature measurement test comprising:
• S2.2.1 Conducting a first round of temperature measurement tests comprising:
• Starting DTS to monitor the borehole background temperature-depth profile;
• wherein the monitoring time is more than 24 hours;
• S2.2.2 Performing a Data quality check of the first round of temperature measurement tests to ensure data integrity and quality;
• S2.2.3 Repeating steps S2.2.1 and S2.2.2 to perform background temperature measurement 2-3 times to grasp the overall temperature change characteristics;
• S2.3 Conducting post-experiment processing and data integration analysis comprising:
• S2.3.1 of Calibrating temperature data between a downhole temperature measuring unit and a surface temperature measuring unit comprising:
• calibrating all collected temperature data to ensure that the temperature measurement results at different depths and locations are coordinated with each other;
• S2.3.2 Conducting a wellbore temperature-depth monitoring profile analysis comprising:
• integrating and analyzing all monitoring data, to produce a temperature-depth profile diagram, and to provide basic data for subsequent active heating tests comprising:
• Step S3:
• Performing a line-source active heating distributed temperature measurement test, comprising performing multiple rounds of heating tests;
• comprising:
• in each round of heating tests, controlling a heating device controller to turn on composite heating optical cable for a predetermined time, so that the temperature of groundwater monitoring line rises to a predetermined temperature value, then stopping heating;
• and stopping monitoring when the temperature of the groundwater monitoring line returns to the background temperature;
• and performing a quality check on the temperature measurement test data after each round of heating tests;
• S3.1 Conducting a preliminary test comprising:
• S3.1.1 Repeating the steps in S2.1.2 in addition to circuit connection debugging, and power testing of power supply equipment
• S3.1.2 Conducting a downhole line heating maximum heating power test, a maximum heating temperature test, and optimization and adjustment of heating parameter settings comprising:
• S3.2 Conducting a formal temperature measurement test comprising:
• S3.2.1 Conducting background temperature monitoring, comprising turning on DTS to start monitoring, and selecting the background monitoring duration according to the background temperature measurement results in S2.3;
• S3.2.2 Conducting a first round of heating tests comprising:
• Turning on the intelligent controller of the heating device, and according to the results of S3.1.2, turning on the heating cable for a suitable time to make the temperature of the groundwater monitoring line rise by a certain degree, then turning off the heating cable and stopping DTS monitoring when the temperature of the monitoring line returns to the background temperature;
• S3.2.3 Conducting a quality check of DTS temperature measurement data in the first round of a temperature measurement test;
• S3.2.4 Repeating S3.2.1, S3.2.2 and S3.2.3 and completing 2-3 heating tests;
• S3.3 Conducting post-experiment processing and data integration analysis comprising:
• S3.3.1 Calibrating temperature data of the downhole temperature measuring unit and the surface temperature measuring unit comprising:
• comparing the data collected by the downhole temperature measuring unit and surface temperature measuring unit with standard temperature data, adjusting the deviation, and ensuring that all temperature measurement data meet high accuracy and high reliability standards;
• S3.3.2 Conducting analysis of downhole temperature-depth profiles during the heating and cooling phases comprising:
• analyzing temperature data during the heating phase, focusing on changes in temperature with depth, and the speed and intensity of groundwater response to heating in boreholes of different depths;
• Step S4:
• Denoising thermal plume attenuation signal data of the groundwater obtained in Step 3;
• S4.1 Conducting original signal data trimming and conversion comprising:
• S4.1.1 For the obtained line-source thermal plume attenuation signal data of the groundwater, determine a depth range and depth interval a temperature data recording frequency and a time interval for data recording;
• S4.1.2 Use the determined depth range to segment the original data, and select and obtain the required temperature data according to the determined time interval;
• such that X is a given two-dimensional line source thermal plume attenuation signal data of the groundwater, with a depth range L of [0, M], a time range T of [0, N], then the temperature data X is a two-dimensional array of shape N×M, X={xi,j}, i∈{1, 2, . . . , N};
• j∈{1, 2, . . . , M};
• S4.1.3 Perform a transposition operation on the trimmed temperature data to obtain a two-dimensional array XT of shape M×N, XT={xi,j}, i∈{1, 2, . . . , M};
• j∈{1, 2, . . . , N};
• S4.2 Conducting temperature data signal denoising comprising:
• S4.2.1 Determining a cutoff frequency;
• S4.2.2 Using a Butterworth low-pass filter to perform low-pass filtering on the signal;
• wherein for the Butterworth filter, a relationship between amplitude and frequency of an n Butterworth low-pass filter is expressed by the following formula:
• ❘
• "\[LeftBracketingBar]"
• H
• ⁡
• (
• j
• ⁢
• ω
• )
• ❘
• "\[RightBracketingBar]"
• =
• 1
• 1
• +
• (
• ω
• ω
• C
• )
• 2
• ⁢
• n
• ,
• where |H(jω)| is amplitude of frequency response;
• ω is angular frequency;
• ωc is cutoff frequency, which is a frequency point at which the filter begins to significantly attenuate the output signal;
• and n is the order of the filter;
• Step S5:
• Fitting a groundwater velocity-depth curve in the borehole to obtain a distributed groundwater velocity data set by monitoring and analyzing the temperature difference signal between the internal temperature and the external temperature of the composite heating optical cable in a single depth section of the borehole.