Water quality chromaticity rapid detection method with turbidity compensation function

Abstract

The invention relates to the technical field of water quality detection, in particular to a water quality chromaticity rapid detection method with a turbidity compensation function, and the method comprises the following steps: acquiring a transmission signal in a to-be-detected water body of a detection cavity along a physically set optical path direction; determining a transmission light intensity center drift dynamic tendency index of water quality chromaticity detection; determining absorption and scattering coupling indexes of water quality chromaticity detection according to time points where two adjacent wave troughs of the peak value region are located; based on the transmission light intensity center drift dynamic trend index and the absorption and scattering coupling index, determining an equivalent optical path disturbance index of water quality chromaticity detection; and correcting the physically set optical path by using the equivalent optical path disturbance index to obtain an effective optical path, and determining the chromaticity of the water body to be measured according to the effective optical path. Therefore, the stability, the accuracy and the reliability of the chromaticity detection result are improved.

Description

Technical Field

[0001] This invention relates to the field of water quality testing technology, specifically to a rapid method for detecting water color with turbidity compensation function. Background Technology

[0002] Water color is an important sensory and physicochemical indicator reflecting the levels of dissolved organic matter, humus, metal ions, and industrial pollutants in water bodies. It is also a key parameter for drinking water safety assessment, municipal water supply treatment, industrial circulating water management, and surface water environmental monitoring. Actual water bodies often contain suspended particles, colloids, and microorganisms, which can cause turbidity interference, particularly noticeable in the detection of surface water, river water, landscape water, and insufficiently settled raw water.

[0003] In some scenarios, during rapid optical detection of water colorimetry, the dynamic changes in the particle size of micron-sized suspended particles within the water sample can induce multiple scattering and enhanced forward scattering of light, altering the propagation path of light within the detection optical path and causing the effective optical path to deviate from the preset physical optical path. Current colorimetry detection methods largely rely on the Lambert-Beer law, assuming a constant physical optical path and failing to consider optical path drift caused by scattering and the absorption-scattering coupling effect. Therefore, nonlinear systematic errors are easily generated in rapid flow detection, resulting in low stability and accuracy of water colorimetry detection results. Summary of the Invention

[0004] To address the technical problem of low stability and accuracy in water colorimetric detection results, the present invention aims to provide a rapid water colorimetric detection method with turbidity compensation function.

[0005] To solve the above technical problems, the specific technical solution adopted is as follows: In a first aspect, embodiments of the present invention provide a rapid water quality colorimetric detection method with turbidity compensation function, comprising: acquiring a transmission signal along a physically set optical path direction in the water body to be tested within a detection chamber; determining a dynamic trend index of the transmission light intensity center drift of the water quality colorimetric detection based on the peak amplitude of the peak region, the half-width at half-maximum (WHM) of the peak region, the light intensity distribution center of the peak region, the reference WHM of the reference water body, and the reference light intensity distribution center of the reference water body; determining an absorption and scattering coupling index of the water quality colorimetric detection based on the time points of the two adjacent troughs of the peak region; determining an equivalent optical path perturbation index of the water quality colorimetric detection based on the dynamic trend index of the transmission light intensity center drift and the absorption and scattering coupling index; correcting the physically set optical path using the equivalent optical path perturbation index to obtain an effective optical path, and determining the colorimetric properties of the water body to be tested based on the effective optical path.

[0006] Optionally, the dynamic trend index of the transmission intensity center drift of water quality colorimetric detection is determined based on the peak amplitude, half-width at half-maximum (WHM) of the peak region, the light intensity distribution center of the peak region, the reference WHM of the reference water body, and the reference light intensity distribution center of the reference water body. This includes: determining the transient dispersion of the transmission signal based on the peak amplitude and time interval of adjacent peak regions; determining the transmission response expansion of water quality colorimetric detection based on the transient dispersion, the WHM of the peak region, and the reference WHM of the reference water body; determining the light intensity distribution center of the peak region based on each time point and the corresponding light intensity value in the peak region; determining the transmission response expansion, the light intensity distribution center of the peak region, and the reference light intensity distribution center of the reference water body; and determining the dynamic trend index of the transmission intensity center drift of water quality colorimetric detection based on the transmission intensity distribution center shift amplitude and the light intensity distribution centers of adjacent peak regions.

[0007] Optionally, determining the transient dispersion of the transmitted signal based on the peak amplitude and time interval of adjacent peak regions of the transmitted signal includes: calculating a first difference between the peak amplitudes of adjacent peak regions of the transmitted signal in the current observation window, and a second difference between the time intervals between the time points corresponding to the peak amplitudes of adjacent peak regions; calculating a third difference between the second differences of each adjacent time interval; and determining the transient dispersion of the transmitted signal based on the first difference and the third difference.

[0008] Optionally, determining the transmission response extension of water quality colorimetry detection based on transient dispersion, the half-peak width of peak regions, and the reference half-peak width of the reference water body includes: calculating a first ratio between the half-peak width of each peak region during the detection of the water body to be tested and the reference half-peak width of the reference water body; and determining the transmission response extension of water quality colorimetry detection based on each first ratio, the number of all peak regions in the transmission signal during the detection of the water body to be tested, and the transient dispersion.

[0009] Optionally, the offset of the transmitted light intensity distribution center for water quality colorimetric detection is determined based on the degree of transmission response expansion, the light intensity distribution center of the peak region, and the reference light intensity distribution center of the reference water body. This includes: calculating the fourth difference between the light intensity distribution center of the peak region and the reference light intensity distribution center of the reference water body; and determining the offset of the transmitted light intensity distribution center for water quality colorimetric detection based on the number of all peak regions in the transmitted signal during the detection process of the water body to be tested, the fourth difference, and the degree of transmission response expansion.

[0010] Optionally, the dynamic trend index of the transmission light intensity center drift in water quality colorimetry detection is determined based on the offset amplitude of the transmission light intensity distribution center and the light intensity distribution centers of adjacent peak regions. This includes: calculating the fifth difference between the light intensity distribution centers of adjacent peak regions; and determining the dynamic trend index of the transmission light intensity center drift in water quality colorimetry detection based on the number of all peak regions in the transmission signal during the detection process of the water body to be tested, the fifth difference, and the offset amplitude of the transmission light intensity distribution center.

[0011] Optionally, the absorption and scattering coupling index of water quality colorimetry detection is determined based on the time points of the two adjacent troughs of the peak region, including: calculating the sixth difference between the time point of the peak amplitude of the peak region and the time points of the two adjacent troughs of the peak amplitude of the peak region; and determining the absorption and scattering coupling index of water quality colorimetry detection based on the number of all peak regions in the transmitted signal during the detection process of the water body to be tested and the absolute value of the sixth difference.

[0012] Optionally, the equivalent optical path perturbation index for water quality colorimetry detection is determined based on the dynamic trend index of the center drift of the transmitted light intensity and the absorption and scattering coupling index, including: calculating the first product between the dynamic trend index of the center drift of the transmitted light intensity and the absorption and scattering coupling index; and normalizing the first product to obtain the equivalent optical path perturbation index for water quality colorimetry detection.

[0013] Optionally, the effective optical path is obtained by correcting the physical optical path using the equivalent optical path perturbation index, which includes: calculating the sum between the predetermined value and the equivalent optical path perturbation index; and determining the second product between the sum and the physical optical path as the effective optical path.

[0014] Optionally, determining the chromaticity of the water body to be tested based on the effective optical path length includes: calculating a second ratio between the incident light intensity of the incident signal and the transmitted light intensity of the transmitted signal of the water body to be tested, and a third product between the molar absorptivity and the effective optical path length; and determining the chromaticity of the water body to be tested based on the second ratio and the third product.

[0015] This invention offers the following advantages: By extracting key features such as peak amplitude, full width at half maximum (FWHM), light intensity distribution center, and trough time from the transmission signal of the water sample, this invention quantifies the influence of multiple scattering, enhanced forward scattering, and dynamic particle size changes caused by micron-sized suspended particles, reducing the interference of suspended impurities in high-turbidity, high-colloidal water on colorimetric detection. Furthermore, by constructing a dynamic trend index for the drift of the transmitted light intensity center, it tracks the light propagation path shift and effective optical path deviation caused by scattering in real time, eliminating reliance on a fixed physical optical path for calculation and removing systematic errors caused by optical path drift, thus maintaining a stable optical response relationship even in rapid flow-through detection. Secondly, by constructing an absorption and scattering coupling index based on the adjacent trough time points of the peak region, it characterizes the coupling effect of light absorption and scattering in the water sample. Combined with an equivalent optical path perturbation index, it completes optical path correction, solving the nonlinear error caused by neglecting coupling effects in traditional methods. Finally, based on the corrected effective optical path, it calculates the colorimetric value of the water sample, achieving dynamic turbidity compensation and significantly improving the stability, accuracy, and reliability of colorimetric detection results. Attached Figure Description

[0016] To more clearly illustrate the technical solutions and advantages in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0017] Figure 1 A flowchart illustrating a rapid water color detection method with turbidity compensation function provided in an embodiment of the present invention; Figure 2 This is a schematic diagram of a rapid water color detection system with turbidity compensation function provided in an embodiment of the present invention. Detailed Implementation

[0018] To further illustrate the technical means and effects adopted by the present invention to achieve its intended purpose, the following, in conjunction with the accompanying drawings and preferred embodiments, details the specific implementation, structure, features, and effects of a rapid water color detection method with turbidity compensation function proposed according to the present invention. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. Furthermore, specific features, structures, or characteristics in one or more embodiments can be combined in any suitable form.

[0019] The terms "first," "second," etc., used in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such terms are interchangeable where appropriate; this is merely a way of distinguishing objects with the same attributes in the embodiments of this application. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion, so that a process, method, system, product, or apparatus that comprises a series of elements is not necessarily limited to those elements, but may include other elements not explicitly listed or inherent to those processes, methods, products, or apparatuses.

[0020] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains.

[0021] The following describes in detail, with reference to the accompanying drawings, a specific scheme for a rapid water color detection method with turbidity compensation function provided by the present invention.

[0022] Example 1: Please see Figure 1 The diagram illustrates a flowchart of a rapid water color detection method with turbidity compensation function according to an embodiment of the present invention, including: Step S101: Obtain the transmission signal along the physically set optical path direction in the water body to be tested within the detection chamber.

[0023] Specifically, in this embodiment of the invention, a main transmission light channel is set within the detection cavity, wherein the main transmission channel is arranged along a physically predetermined optical path direction. The main transmission light detector collects the transmitted light intensity value in real time, forming a transmission intensity sequence data, i.e., a transmission signal, and records timestamp information. Furthermore, this embodiment of the invention synchronizes and aligns the transmission signal and other monitoring data according to a unified time axis, forming a multi-dimensional time-series monitoring data set.

[0024] Furthermore, the water body to be tested refers to the water body that needs to undergo water color detection.

[0025] Step S102: Based on the peak amplitude of the peak region of the transmitted signal, the half-width of the peak region, the light intensity distribution center of the peak region, the reference half-width of the reference water body, and the reference light intensity distribution center of the reference water body, determine the dynamic trend index of the drift of the transmitted light intensity center for water quality colorimetry detection.

[0026] Specifically, traditional turbidity compensation methods only focus on energy attenuation caused by scattering, while the embodiments of the present invention map the multipath propagation problem caused by scattering into a time structure imbalance problem by observing the stable state of the time structure of the transmitted signal. This transformation enables the separation of absorption changes and path changes in the analytical dimension, providing a physical basis for the subsequent decoupling of absorption-scattering coupling components. Therefore, as an optional embodiment of the present invention, determining the dynamic trend index of the drift of the transmitted light intensity center in water quality colorimetry detection based on the peak amplitude of the peak region of the transmitted signal, the half-width at half-maximum (WHM) of the peak region, the light intensity distribution center of the peak region, the reference WHM of the reference water body, and the reference light intensity distribution center of the reference water body includes: determining the transient dispersion of the transmitted signal based on the peak amplitude and time interval of adjacent peak regions of the transmitted signal; determining the transmission response expansion of water quality colorimetry detection based on the transient dispersion, the WHM of the peak region, and the reference WHM of the reference water body; determining the light intensity distribution center of the peak region based on each time point and the corresponding light intensity value in the peak region; determining the offset amplitude of the transmitted light intensity distribution center in water quality colorimetry detection based on the transmission response expansion amplitude, the light intensity distribution center of the peak region, and the reference light intensity distribution center of the reference water body; and determining the dynamic trend index of the drift of the transmitted light intensity center in water quality colorimetry detection based on the offset amplitude of the transmitted light intensity distribution center and the light intensity distribution centers of adjacent peak regions.

[0027] Specifically, in this embodiment of the invention, a short observation window is first set during the rapid detection of water color, with the window size potentially set to 30 seconds. Then, transmission signal data within the current observation window is obtained using existing detection equipment. Furthermore, the transmission signal data is not simply an absorbance value, but rather a high-temporal-resolution sequence of light intensity values ​​formed by the continuous electrical signal output by the photodetector in the transmission channel after the detection beam emitted by the light source passes through the water sample, and then converted from analog to digital. Therefore, in this embodiment of the invention, the normalized result of the peak amplitude in the peak region is denoted as h in the transmission signal data obtained in the current observation window; then, all peak points in the time series curve are statistically analyzed (all peak points in the time series curve can be directly obtained using a peak point detection algorithm); and the normalized result of the time interval between two adjacent peak points is denoted as h. It is worth noting that the peak amplitudes of the peak regions and the time intervals between peak points used in subsequent calculations are all normalized values, and the transmitted signal data used in subsequent calculations are all transmitted signal data within the observation window. The normalization method described above can be achieved using the min-max normalization method.

[0028] Furthermore, as an optional embodiment of the present invention, determining the transient dispersion of the transmitted signal based on the peak amplitude and time interval of adjacent peak regions of the transmitted signal includes: calculating a first difference between the peak amplitudes of adjacent peak regions of the transmitted signal in the current observation window, and a second difference between the time intervals between the time points corresponding to the peak amplitudes of adjacent peak regions; calculating a third difference between the second differences of each adjacent time interval; and determining the transient dispersion of the transmitted signal based on the first difference and the third difference.

[0029] Specifically, the embodiments of the present invention use the following formula to calculate the transient dispersion of the transmitted signal: In the above formula, It indicates the degree of transient dispersion of the transmitted signal. This indicates the number of peak regions in the current observation window. This represents the peak amplitude of the i-th peak region within the current observation window, and it is a normalized value. Similarly, This represents the peak amplitude of the (i+1)th peak region adjacent to the i-th peak region. This indicates the difference in peak amplitude between adjacent peak regions. This represents the average difference between all time intervals. If the propagation path is stable, the time intervals between peaks are approximately uniform; if there is photon dwell time extension, the difference between time intervals increases.

[0030] The average difference between all time intervals The method for obtaining the values ​​is as follows: Calculate the time interval between the peak amplitudes of each pair of peak regions, calculate the absolute value of the difference between each pair of time intervals, and use this absolute value to represent the difference between each pair of time intervals. The average of the absolute values ​​of the differences between all pair of time intervals is then used as the average of the differences between all time intervals. .

[0031] When the photon propagation path is stable, the transmission signal fluctuation mainly comes from light source noise and fluid smooth fluctuation; when multiple scattering is enhanced, some photon paths are lengthened or folded back, resulting in a wider arrival time distribution, which manifests as an increase in local fluctuation amplitude. The larger the value, the more it indicates that the transient dispersion of the transmitted signal can only be significantly increased when both amplitude and time anomalies exist simultaneously.

[0032] Furthermore, although obtaining the transient dispersion of the transmitted signal can identify whether the photon propagation path distribution has been disturbed, this transient dispersion only characterizes the degree of disruption to the temporal structure stability and does not reflect the overall shift in the photon propagation time distribution. In the above embodiment, all peak regions in the transmitted signal data have been obtained, and then the half-width at half-maximum (WHM) of each peak region is acquired; then, the reference WHM of the reference water body is obtained and denoted as... The reference water body can be clean water. Therefore, as an optional embodiment of the present invention, determining the transmission response extension of water quality colorimetry detection based on transient dispersion, the half-peak width of the peak region, and the reference half-peak width of the reference water body includes: calculating a first ratio between the half-peak width of each peak region during the detection of the water body to be tested and the reference half-peak width of the reference water body; and determining the transmission response extension of water quality colorimetry detection based on each first ratio, the number of all peak regions in the transmission signal during the detection of the water body to be tested, and the transient dispersion.

[0033] Specifically, the embodiments of the present invention use the following formula to calculate the transmission response extension of water colorimetry detection: In the above formula, This indicates the extent of transmission response expansion in water colorimetry detection. It indicates the degree of transient dispersion of the transmitted signal. This represents the half-peak width of the c-th peak region. This indicates the reference half-peak width of the water body. This means comparing the half-peak width of the current water sample's peak region with the baseline half-peak width under clear water conditions. Essentially, it represents how much the time response has been broadened relative to the ideal linear propagation state. c represents one of the peak values. The purpose of calculating the average is that the scattering state inside the water body is not instantaneous, and peaks at different time periods may be affected by different particle distributions.

[0034] in, In physics, this typically represents the coupling strength of two effects, characterizing the degree of light propagation time diffusion modulated by the activity level of particle disturbance. A larger product indicates a stronger coupling between the multiple scattering path extension effect and the dynamic particle disturbance effect in the water body, resulting in a significant expansion and instability of the transmitted light time structure, thus increasing its interference with rapid colorimetric detection results. By calculating the transmission response expansion, the extent to which the light propagation path in the water body is broadened under the background of dynamic suspended particle scattering is quantified, reflecting the intensity of the transient peak shape's temporal structure change. However, analyzing only the peak shape expansion is insufficient to fully describe the interference of light propagation on colorimetric measurements, because changes in particle concentration, particle size distribution, and flow velocity in the water sample can also cause shifts in the overall light attenuation curve.

[0035] Furthermore, after obtaining the extent of transmission response expansion, in order to comprehensively characterize the transient dynamic effects of light propagation in water, it is necessary to further analyze the transient shift characteristics of the transmission attenuation curve. The previous indicator reflects the elongation of the peak width, while the next indicator reflects the drift of the overall light intensity distribution center. Combining the two can comprehensively quantify the influence of multiple scattering and particle disturbance on optical measurements, thereby providing an accurate basis for dynamic turbidity compensation. For the transmission signal time-series data obtained in the above embodiment, all peak regions in the curve are obtained. For each peak region, the light intensity distribution center of a single peak region is obtained. The specific calculation method for the light intensity distribution center is as follows: multiply each time point in the peak region by the corresponding light intensity, then sum and average the products of all time points multiplied by the corresponding light intensity, which gives the light intensity distribution center of the current peak region, denoted as . The reference light intensity distribution center of the clear water was obtained in the same way and denoted as... In the rapid water quality detection process, multiple scattering can cause changes in the length of the light propagation path. An increase in the width of each peak only indicates that the peak shape is elongated, but it does not indicate whether the overall arrival time of the light has shifted. The center of the light intensity distribution directly reflects the average arrival time of the peak.

[0036] Furthermore, as an optional embodiment of the present invention, determining the offset of the transmitted light intensity distribution center of water quality colorimetric detection based on the transmission response expansion degree, the light intensity distribution center of the peak region, and the reference light intensity distribution center of the reference water body includes: calculating the fourth difference between the light intensity distribution center of the peak region and the reference light intensity distribution center of the reference water body; and determining the offset of the transmitted light intensity distribution center of water quality colorimetric detection based on the number of all peak regions in the transmitted signal during the detection process of the water body to be tested, the fourth difference, and the transmission response expansion degree.

[0037] Specifically, the embodiments of the present invention use the following formula to calculate the offset of the center of transmitted light intensity distribution in water quality colorimetry detection: In the above formula, This indicates the magnitude of the center shift in the transmitted light intensity distribution during water colorimetry detection. This indicates the number of all peak regions in the transmitted signal during the detection process of the water body under test. This indicates the extent of transmission response expansion in water colorimetry detection. This represents the time corresponding to the center of the light intensity distribution in the c-th peak region. This indicates the time corresponding to the center of the reference light intensity distribution in the reference water body. This indicates the offset of the light intensity distribution center in the current peak region relative to the ideal clear water. A large offset indicates that the average light path is lengthened or delayed. This represents a sufficiently small non-zero constant, such as 0.0001, to ensure that the calculation is meaningful.

[0038] Furthermore, The product of these two values ​​represents the coupling result of local spread and global drift. That is, when peak spread is strong, if the light intensity center drift is also large, the impact of light propagation perturbation on colorimetry is more significant, meaning the shift in the center of the transmitted light intensity distribution is greater. The shift in the center of the transmitted light intensity distribution mainly reflects the shift of the overall light intensity distribution along the time axis, i.e., the drift of the average arrival time of light relative to clear water. It embodies the static or average drift effect of light propagation, but does not describe the change or fluctuation characteristics of the shift over time.

[0039] Based on this, as an optional embodiment of the present invention, determining the dynamic trend index of the drift of the center of transmitted light intensity in water quality colorimetry detection based on the offset amplitude of the center of transmitted light intensity distribution and the center of light intensity distribution of adjacent peak regions includes: calculating the fifth difference between the centers of light intensity distribution of adjacent peak regions; and determining the dynamic trend index of the drift of the center of transmitted light intensity in water quality colorimetry detection based on the number of all peak regions in the transmitted signal during the detection process of the water body to be tested, the fifth difference, and the offset amplitude of the center of transmitted light intensity distribution.

[0040] Specifically, in this embodiment of the invention, the following formula is used to calculate the dynamic trend index of the center drift of transmitted light intensity in water quality colorimetry detection: In the above formula, This is an index representing the dynamic trend of the center drift of transmitted light intensity in water colorimetry detection. This indicates the magnitude of the center shift in the transmitted light intensity distribution during water colorimetry detection. This represents the time corresponding to the center of the light intensity distribution in the c-th peak region. This represents the time corresponding to the light intensity distribution center of the (c+1)th peak region. This represents a sufficiently small non-zero constant, such as 0.0001, to ensure that the calculation is meaningful.

[0041] Under normal circumstances, the center of light intensity distribution is like a "fixed point", the light propagation path is stable, and the peak shape is clear; This indicates the time interval between the light intensity distribution centers of two adjacent peak regions during the current observation period. This indicates the magnitude of the change in the center of light intensity distribution over time; the larger the value, the more significant the change in the light propagation path over time. It also reflects the overall drift of the light propagation path and its active fluctuations over time. That is, the larger the product, the higher the dynamic trend index of the drift of the center of transmitted light intensity distribution, indicating that the light propagation path is unstable and the colorimetric measurement is easily interfered with.

[0042] Furthermore, the dynamic trend index of the center drift of transmitted light intensity can illustrate the intensity of dynamic disturbances in the light propagation path in the water sample over time. However, changes in transmitted light intensity may originate from either the chromatic absorption of the water sample itself or from multiple scattering disturbances along the light propagation path.

[0043] Step S103: Determine the absorption and scattering coupling index of water colorimetry detection based on the time points of the two adjacent troughs in the peak region.

[0044] Specifically, to achieve accurate analysis of the transmitted signal, it is necessary to further decompose the transmitted signal into an absorption-dominant component and a residence time perturbation component. The absorption-dominant component reflects the true color information of the water sample, while the residence time perturbation component reflects the perturbation effect caused by the dynamic changes in the light propagation path due to particles. For all peak regions in the transmitted signal data obtained within the current observation window as described in the preceding embodiments, the time points where the peak amplitude of the peak region is located are marked and denoted as... Then, the time points where the two troughs on both sides of the peak amplitude in the peak region are located are respectively recorded as follows: Then, calculate the interval between the time point of the peak amplitude of a single peak region and the time points of the troughs located to the left and right of the peak amplitude of the single peak region. That is, subtract the time point of the peak amplitude of the single peak region from the time points of the troughs on the left and right sides, and then take the absolute value, which is recorded as follows: .

[0045] Furthermore, as an optional embodiment of the present invention, determining the absorption and scattering coupling index of water quality colorimetry detection based on the time points of the two adjacent troughs of the peak region includes: calculating the sixth difference between the time point of the peak amplitude of the peak region and the time points of the two adjacent troughs of the peak amplitude of the peak region; and determining the absorption and scattering coupling index of water quality colorimetry detection based on the number of all peak regions in the transmitted signal during the detection process of the water body to be tested and the absolute value of the sixth difference.

[0046] Specifically, the embodiments of the present invention use the following formula to calculate the absorption-scattering coupling index for water colorimetry detection: In the above formula, This represents the absorption-scattering coupling index for water colorimetry detection. It represents the absolute value of the difference between the time point where the peak amplitude of a single peak region is located and the time point where the trough is located to its left. It represents the absolute value of the difference between the time point where the peak amplitude of a single peak region is located and the time point where the trough is located to its right. This represents a sufficiently small non-zero constant, such as 0.0001, to ensure that the calculation is meaningful.

[0047] in, It is used to indicate whether the peak center is biased to the left or right in a single peak region. When there is particle scattering or dynamic disturbance in the water sample, the light propagation path is lengthened and the peak shape is shifted to one side, that is, the larger the difference is. This indicates one of the peak regions; It is used to calculate the average, taking the average of n peak values ​​to obtain the overall coupling strength within the observation window during the detection process. The larger the value, the more likely the change in light intensity is caused by the combined effects of absorption and scattering. In other words, the larger the absorption-scattering coupling index, the stronger the coupling effect.

[0048] Step S104: Based on the dynamic trend index of the center drift of transmitted light intensity and the absorption and scattering coupling index, determine the equivalent optical path disturbance index for water quality colorimetry detection.

[0049] Specifically, embodiments of the present invention reduce nonlinear errors caused by scattering from suspended particles and microparticles by obtaining equivalent optical path parameters, effectively reducing the impact of short-term fluctuations on colorimetric measurements and improving detection accuracy. Therefore, as an optional embodiment of the present invention, determining the equivalent optical path perturbation index for water quality colorimetric detection based on the transmitted light intensity center drift dynamic trend index and the absorption and scattering coupling index includes: calculating the first product between the transmitted light intensity center drift dynamic trend index and the absorption and scattering coupling index; and normalizing the first product to obtain the equivalent optical path perturbation index for water quality colorimetric detection.

[0050] Specifically, the embodiments of the present invention use the following formula to calculate the equivalent optical path disturbance index for water colorimetry detection: In the above formula, This represents the equivalent optical path perturbation index for water colorimetry testing. This represents the absorption-scattering coupling index for water colorimetry detection. The index represents the dynamic trend of the center drift of transmitted light intensity in water colorimetry detection, reflecting the instability of the light propagation path. This demonstrates the enhancing effect of two-factor coupling on optical path perturbation. This represents the normalization function, used to quantize and constrain the calculation results to the range of 0-1, ensuring that the exponent can be directly used for equivalent optical path correction or turbidity compensation; that is, the larger the product, the larger the equivalent optical path perturbation exponent. Among these, The maximum-minimum normalization method can be used.

[0051] Step S105: Correct the physical optical path using the equivalent optical path perturbation index to obtain the effective optical path, and determine the color of the water body to be tested based on the effective optical path.

[0052] Specifically, embodiments of the present invention utilize the equivalent optical path perturbation index to physically set the optical path. The correction is performed. In one optional embodiment of the invention, the effective optical path is obtained by correcting the physical optical path using an equivalent optical path perturbation index, which includes: calculating the sum of a predetermined value and the equivalent optical path perturbation index; and determining the second product between the sum and the physical optical path as the effective optical path.

[0053] Specifically, the predetermined value in the embodiments of the present invention can be 1. The embodiments of the present invention specifically use the following formula to calculate the effective optical path: In the above formula, Indicates the effective optical path. This indicates the physical setting of the optical path. This represents the equivalent optical path disturbance index for water colorimetry detection. When The smaller the value, the closer the equivalent optical path perturbation index is to the physical setting, and the smaller the impact of turbidity; conversely, when... The larger the value, the more the equivalent optical path perturbation index deviates from the physical setting, and the greater the turbidity effect, meaning a larger optical path correction range.

[0054] Furthermore, after obtaining the effective optical path, the chromaticity is solved using existing calculation methods (such as the Lambert-Beer law). After optical path correction, the chromaticity calculation takes into account the optical path drift caused by turbidity, making the measurement more accurate. Finally, the calculated chromaticity is output as a rapid detection result. At the same time, the equivalent optical path perturbation index can be combined with other transient indicators as a quality control or data reliability assessment indicator.

[0055] Furthermore, as an optional embodiment of the present invention, determining the chromaticity of the water body to be tested based on the effective optical path length includes: calculating a second ratio between the incident light intensity of the incident signal and the transmitted light intensity of the transmitted signal of the water body to be tested, and a third product between the molar absorptivity and the effective optical path length; and determining the chromaticity of the water body to be tested based on the second ratio and the third product.

[0056] Specifically, the embodiments of the present invention use the following formula to calculate the color of the water body to be tested: In the above formula, This indicates the color of the water body being tested. Indicates the effective optical path. The molar absorptivity is an inherent property of the water body being tested, representing its ability to absorb light at a specific wavelength of the incident signal. Its unit is molar absorptivity. It is related to the water body being tested, the wavelength of the incident signal, and the temperature. This represents the intensity of the incident light that caused the incident signal. This indicates the intensity of the transmitted light in the transmitted signal.

[0057] This invention extracts key features from the transmission signal of the water sample, such as peak amplitude, full width at half maximum (FWHM), light intensity distribution center, and trough time, to quantify the effects of multiple scattering, enhanced forward scattering, and dynamic particle size changes caused by micron-sized suspended particles. This reduces the interference of suspended impurities in high-turbidity, high-colloidal water on colorimetric detection. Furthermore, by constructing a dynamic trend index for the drift of the transmitted light intensity center, it tracks the light propagation path shift and effective optical path deviation caused by scattering in real time, eliminating the reliance on a fixed physical optical path for calculation and removing systematic errors caused by optical path drift. This ensures that rapid flow-through detection maintains a stable optical response. Secondly, based on the time points of adjacent troughs in the peak region, an absorption and scattering coupling index is constructed to characterize the coupling effect of light absorption and scattering in the water sample. Combined with an equivalent optical path perturbation index, optical path correction is completed, solving the nonlinear errors caused by neglecting coupling effects in traditional methods. Finally, the colorimetric properties of the water sample are calculated based on the corrected effective optical path, achieving dynamic turbidity compensation and significantly improving the stability, accuracy, and reliability of colorimetric detection results.

[0058] Furthermore, this invention acquires transmission signal time-series data through the current observation window, then calculates the transient dispersion of the transmission signal based on the amplitude differences and peak distribution state in the transmission signal time-series data. It then assesses the transmission response expansion during rapid water quality detection by combining the difference between the width of the peak region and the standard control group. Next, it analyzes the difference between the light intensity distribution center and the control group to obtain the shift amplitude of the transmission light intensity distribution center. Then, it combines the change in the light intensity distribution center of each peak region over time during the observation period to obtain the dynamic trend of the transmission light intensity center drift. Next, it analyzes the symmetry of the peak region to resolve the dominant absorption component and the residence time perturbation component, thereby obtaining the absorption and scattering coupling index during rapid water quality color detection. Finally, it combines the dynamic trend of the transmission light intensity center drift and the absorption and scattering coupling index to obtain the equivalent optical path perturbation index during the detection stage. Using the equivalent optical path perturbation index, the physical optical path is corrected to the effective optical path, and color is calculated based on the corrected optical path to achieve dynamic turbidity compensation, thereby improving the stability and accuracy of rapid detection.

[0059] Example 2: Corresponding to the rapid water color detection method with turbidity compensation function provided in the above embodiments, based on the same technical concept, this embodiment of the invention also provides a rapid water color detection system with turbidity compensation function. This rapid water color detection system with turbidity compensation function is used to execute the above-described rapid water color detection method with turbidity compensation function. Figure 2 This is a schematic diagram of a rapid water color detection system with turbidity compensation function provided in one embodiment of the present invention, as shown below. Figure 2 As shown. Water quality color rapid detection systems with turbidity compensation function can vary considerably depending on configuration or performance. They may include one or more processors 201 and memory 202. The memory 202 stores computer programs that can run on the processor 201. The processor 201 executes the programs stored in the memory 202 to achieve the above... Figure 1 The various steps in the Chinese method embodiment. The memory 202 can be temporary or persistent storage. The application stored in the memory 202 may include one or more modules (not shown in the figures), each module may include a series of computer-executable instructions for the rapid water color detection system with turbidity compensation function.

[0060] Furthermore, the processor 201 can be configured to communicate with the memory 202 and execute a series of computer-executable instructions stored in the memory 202 on the water quality color rapid detection system with turbidity compensation function. The water quality color rapid detection system with turbidity compensation function may also include one or more power supplies 203, one or more wired or wireless network interfaces 204, one or more input / output interfaces 205, and one or more keyboards 206.

[0061] Specifically, in this embodiment, the rapid water color detection system with turbidity compensation function includes a processor, a communication interface, a memory, and a communication bus; wherein, the processor, communication interface, and memory communicate with each other via the bus; the memory is used to store computer programs; the processor is used to execute the programs stored in the memory to achieve the above... Figure 1 The various steps in the method embodiments are the same as those in the above method embodiments, and have the same beneficial effects. To avoid repetition, the embodiments of the present invention will not be described again here.

[0062] It should be noted that the water quality color rapid detection system with turbidity compensation function provided in this embodiment of the invention and the water quality color rapid detection method with turbidity compensation function provided in this embodiment of the invention are based on the same application concept. Therefore, the specific implementation of this embodiment can refer to the implementation of the aforementioned water quality color rapid detection method with turbidity compensation function, and has the same or similar beneficial effects. Repeated parts will not be described again.

[0063] It should be noted that the order of the above embodiments of the present invention is merely for descriptive purposes and does not represent the superiority or inferiority of the embodiments. The processes depicted in the accompanying drawings do not necessarily require a specific or sequential order to achieve the desired result. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0064] The various embodiments in this specification are described in a progressive manner. The same or similar parts between the various embodiments can be referred to each other. Each embodiment focuses on describing the differences from other embodiments.

[0065] This invention also proposes a computer-readable storage medium storing one or more programs. When executed by a rapid water quality color detection system with turbidity compensation function, which includes multiple application programs, the one or more programs cause the rapid water quality color detection system with turbidity compensation function to perform... Figure 1 The methods disclosed in the embodiments shown achieve the functions and beneficial effects of the methods in the preceding method embodiments, and will not be repeated here.

[0066] Computer-readable storage media include read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

Claims

1. A rapid water color detection method with turbidity compensation function, characterized in that, The rapid water color detection method with turbidity compensation function includes: Acquire the transmission signal along the physically set optical path direction in the water body to be tested within the detection chamber; Based on the peak amplitude of the peak region of the transmitted signal, the half-width at half-peak of the peak region, the light intensity distribution center of the peak region, the reference half-width at half-peak of the reference water body, and the reference light intensity distribution center of the reference water body, the dynamic trend index of the drift of the transmitted light intensity center in water quality colorimetry detection is determined. The absorption-scattering coupling index for water colorimetry detection is determined based on the time points of the two adjacent troughs in the peak region. Based on the transmitted light intensity center drift dynamic trend index and the absorption and scattering coupling index, the equivalent optical path perturbation index for water quality colorimetry detection is determined. The physical optical path is corrected using the equivalent optical path perturbation index to obtain the effective optical path, and the color of the water body to be tested is determined based on the effective optical path.

2. The rapid water color detection method with turbidity compensation function according to claim 1, characterized in that, The determination of the dynamic trend index of transmitted light intensity center drift in water quality colorimetry detection based on the peak amplitude of the peak region of the transmitted signal, the half-width at half-maximum of the peak region, the light intensity distribution center of the peak region, the reference half-width at half-maximum of the reference water body, and the reference light intensity distribution center of the reference water body includes: The transient dispersion of the transmitted signal is determined based on the peak amplitude and time interval of adjacent peak regions of the transmitted signal. Based on the transient dispersion, the half-peak width of the peak region, and the reference half-peak width of the reference water body, the transmission response spread of water quality colorimetry detection is determined. Based on each time point and the corresponding light intensity value in the peak region, the light intensity distribution center of the peak region is determined; Based on the degree of transmission response expansion, the light intensity distribution center of the peak region, and the reference light intensity distribution center of the reference water body, the offset amplitude of the transmission light intensity distribution center of water quality colorimetry detection is determined. Based on the shift amplitude of the center of transmitted light intensity distribution and the light intensity distribution centers of adjacent peak regions, the dynamic trend index of the drift of the center of transmitted light intensity in water quality colorimetry detection is determined.

3. The rapid water color detection method with turbidity compensation function according to claim 2, characterized in that, Determining the transient dispersion of the transmitted signal based on the peak amplitude and time interval of adjacent peak regions includes: Calculate the first difference between the peak amplitudes of adjacent peak regions of the transmitted signal in the current observation window, and the second difference between the time intervals between the time points corresponding to the peak amplitudes of adjacent peak regions; Calculate the third difference between the second differences of each adjacent time interval; The transient dispersion of the transmitted signal is determined based on the first difference and the third difference.

4. The rapid water color detection method with turbidity compensation function according to claim 2, characterized in that, The determination of the transmission response spread of water quality colorimetry detection based on the transient dispersion, the half-maximum width of the peak region, and the reference half-maximum width of the reference water body includes: Calculate the first ratio between the half-peak width of each peak region during the detection process of the water body to be tested and the reference half-peak width of the reference water body; Based on each of the first ratios, the number of all peak regions in the transmission signal during the detection of the water body to be tested, and the transient dispersion, the transmission response expansion of water quality colorimetry detection is determined.

5. The rapid water color detection method with turbidity compensation function according to claim 2, characterized in that, The step of determining the shift amplitude of the transmitted light intensity distribution center in water quality colorimetry detection based on the transmission response expansion, the light intensity distribution center of the peak region, and the reference light intensity distribution center of the reference water body includes: Calculate the fourth difference between the light intensity distribution center of the peak region and the reference light intensity distribution center of the reference water body; The offset amplitude of the center of transmission light intensity distribution in water quality colorimetry detection is determined based on the number of all peak regions in the transmission signal during the detection process of the water body to be tested, the fourth difference, and the degree of transmission response expansion.

6. The rapid water color detection method with turbidity compensation function according to claim 2, characterized in that, The method for determining the dynamic trend index of the drift of the transmitted light intensity center in water quality colorimetry detection based on the shift amplitude of the transmitted light intensity distribution center and the light intensity distribution centers of adjacent peak regions includes: Calculate the fifth difference between the light intensity distribution centers of adjacent peak regions; Based on the number of all peak regions in the transmitted signal during the detection process of the water body to be tested, the fifth difference, and the offset amplitude of the center of transmitted light intensity distribution, the dynamic trend index of the center drift of transmitted light intensity in water quality colorimetry detection is determined.

7. The rapid water color detection method with turbidity compensation function according to claim 1, characterized in that, The step of determining the absorption-scattering coupling index of water colorimetry based on the time points of the two adjacent troughs in the peak region includes: Calculate the sixth difference between the time point where the peak amplitude of the peak region is located and the time points where the two troughs adjacent to the time point where the peak amplitude of the peak region is located; The absorption and scattering coupling index for water quality colorimetry detection is determined based on the number of all peak regions in the transmitted signal during the detection process of the water body to be tested and the absolute value of the sixth difference.

8. The rapid water color detection method with turbidity compensation function according to claim 1, characterized in that, The determination of the equivalent optical path disturbance index for water quality colorimetry based on the transmitted light intensity center drift dynamic trend index and the absorption and scattering coupling index includes: Calculate the first product between the transmitted light intensity center drift dynamic trend index and the absorption and scattering coupling index; The first product is normalized to obtain the equivalent optical path perturbation index for the water colorimetric detection.

9. The rapid water color detection method with turbidity compensation function according to claim 1, characterized in that, The step of correcting the physically set optical path using the equivalent optical path perturbation index to obtain the effective optical path includes: Calculate the sum between the predetermined value and the equivalent optical path perturbation index; The second product between the sum and the physical set optical path is determined as the effective optical path.

10. The rapid water color detection method with turbidity compensation function according to claim 1, characterized in that, The step of determining the color of the water body to be tested based on the effective optical path includes: Calculate the second ratio between the incident light intensity of the incident signal and the transmitted light intensity of the transmitted signal of the water body to be tested, and the third product between the molar absorptivity and the effective optical path. The color of the water body to be tested is determined based on the second ratio and the third product.

Images