Method and system for determining flow regime of heterogeneous seepage channel based on solute tracer method

Abstract

The application discloses a kind of non-homogeneous seepage passage flow state determination method and system based on solute tracer method, comprising: under different flow conditions, tracer test is carried out to non-homogeneous seepage passage, obtains solute transport breakthrough curve and extracts characteristic time;Based on the double logarithm relationship between characteristic time and flow, the relative deceleration index is calculated, and through the phenomenon that the relative deceleration index systematically decreases from the theoretical value 1, the critical flow interval of the transition from Darcy flow to non-Darcy flow is identified;Two segment power law fitting is carried out to normalized characteristic time and normalized flow, and the actual critical flow is inversed by solving continuity equation;The inversed actual critical flow is compared with the current actual flow, and the current flow state is judged.The application can realize quantitative and reliable prediction of the flow state transition of non-homogeneous passage without complex non-homogeneous passage geometric parameters or multi-stage differential pressure test, and is especially suitable for seepage safety evaluation and risk early warning of non-homogeneous passage system which is difficult to test directly.

Description

Technical Field

[0001] This invention relates to the field of flow regime determination technology, and specifically to a method and system for determining the flow regime of heterogeneous seepage channels based on solute tracer method. Background Technology

[0002] Heterogeneous seepage channels are critical pathways for fluid transport in projects such as reservoir leakage, geological disposal of nuclear waste, and deep geothermal development. Accurately identifying the seepage state within these channels is crucial for assessing the production efficiency, long-term safety, and environmental risks of engineering systems. When flow velocity increases (e.g., in high-head geothermal extraction) or channel opening becomes large, the inertial effect becomes prominent, and seepage transitions from Darcy to non-Darcy flow. This transition is characterized by a nonlinear relationship between velocity and hydraulic gradient, leading to a sharp increase in flow resistance and changes in solute transport mechanisms. Therefore, accurately predicting the critical point of the Darcy-non-Darcy flow transition is one of the core issues in engineering design and safety assessment. Traditional methods for determining non-Darcy flow rely on empirical criteria such as critical Reynolds number, critical flow rate, or critical hydraulic gradient. However, the critical point values ​​differ significantly for different types of heterogeneous seepage channels (fractures), and these specific critical values ​​are difficult to apply to different media. Recent developments in flow state determination methods are mostly based on nonlinear seepage experiments, obtained through experimental methods. P - Q The flow regime is determined by a pressure drop deviation of 10%. While this method is relatively intuitive, it faces many limitations in practical applications. For well-hidden seepage channels, often only the seepage overflow point can be observed, while the seepage initiation point is difficult to determine, thus making it impossible to obtain the complete upstream and downstream pressure difference of the seepage path. P Meanwhile, the parameterized critical Reynolds number model based on a 10% pressure drop involves a complex parameter calibration process, requiring known geometric parameters of the heterogeneous channel (such as aperture and roughness) to calculate the Reynolds number. In contrast, obtaining solute transport breakthrough curves through single or a small number of tracer tests is a more economical and convenient method for on-site investigation. The characteristic times in the curve (such as initial arrival time, half-concentration time, and mean residence time) comprehensively reflect the connectivity structure, effective flow channels, and velocity distribution information of the seepage channel. Currently, characteristic times are mainly used to invert hydraulic parameters such as dispersion, but there is no technology that can directly establish a reliable predictive relationship between them and the critical conditions of non-Darcy flow where Darcy's law fails. Therefore, there is an urgent need for a new method that does not rely on cumbersome multi-stage differential pressure tests and precise geometric parameters, but only utilizes readily available solute transport data to predict the transition from Darcy to non-Darcy flow in heterogeneous seepage channels. Summary of the Invention

[0003] To overcome the above-mentioned defects of the prior art, the present invention provides a method for determining the flow regime of heterogeneous seepage channels based on solute tracer method. The core of this method is as follows: First, by analyzing the relative rate of change of characteristic time with flow rate, the critical flow range of flow regime transition is identified; then, based on the self-normalized characteristic time data, the specific critical flow value is accurately located by fitting and inverting a two-segment power law model.

[0004] To address the aforementioned technical problems, the present invention adopts the following technical solution: A method for determining the flow regime of heterogeneous seepage channels based on solute tracer method includes the following steps: Obtain solute transport breakthrough curves of the target channel under different flow conditions, extract feature time from each solute transport breakthrough curve, and form multiple flow-feature time data pairs sorted by flow rate with the corresponding flow rate; Based on the flow-characteristic time data pair, the relative deceleration index, which reflects the degree of solute transport deceleration, is calculated, and the curve of the relative deceleration index versus flow rate is plotted. Based on the starting position of the relative deceleration index continuously decreasing from 1 and falling below a certain set limit in the curve, the initial critical flow rate range for flow regime transition is determined. The actual critical flow rate is inverted within the initial critical flow rate range of the flow regime transition using a continuous condition equation model; The inverted actual critical flow rate is compared with the current actual flow rate to determine the current flow state.

[0005] Furthermore, the characteristic time is the initial arrival time of the tracer, the half-concentration time of the tracer, or the average residence time of the tracer calculated based on the method of moments.

[0006] Furthermore, the relative deceleration index is:

[0007] in, For the first The relative deceleration index under various flow conditions Total number of flow-feature time data pairs; The first species and first Characteristic time of solute transport breakthrough curve under certain flow rate conditions; The first species and first Flow under certain flow conditions.

[0008] Furthermore, the specific criterion for determining whether the relative deceleration index in the change curve continuously decreases from 1 and falls below a certain set limit is: when two consecutive relative deceleration index values ​​satisfy <1- deltaFurthermore, if these relative deceleration index values ​​exhibit a monotonically decreasing or overall decreasing trend, then the first one to satisfy the condition of relative deceleration index < 1- delta The flow rate corresponding to the point is taken as the first non-Darcy flow point, and the point preceding the first non-Darcy flow point is taken as the last Darcy flow point. delta This is a preset threshold.

[0009] Furthermore, the preset threshold delta The value ranges from 0.05 to 0.1.

[0010] Furthermore, the initial critical flow range for the flow regime transition is: Setting the first This point is the last Darcy flow point. If the first non-Darcy flow point is [the initial critical flow range is], then [the range is] ... ,in For the first Traffic per point For the first Traffic flow at each point.

[0011] Furthermore, the inverted actual critical flow rate includes: The flow rate of the last Darcy flow point All previous data points were categorized into Darcy flow segments for fitting, and power-law functions for non-Darcy flow segments were obtained. The flow rate of the first non-Darcy flow point All subsequent data points were classified as non-Darcy flow segments for fitting, and the power-law function of the Darcy flow segment was obtained. The actual critical flow rate is determined by solving the condition equations for the continuity of the power-law function at specific critical points in both the non-Darcy and Darcy sections.

[0012] Furthermore, the power-law function of the non-Darcy flow segment is: ; The Darcy flow segment power-law function is: ; in, 、 These are the power-law function coefficients for the non-Darcy flow section and the Darcy flow section, respectively; For self-normalized flow; For self-normalized characteristic time; 、 These are the power exponents of the power-law functions for the non-Darcy flow segment and the Darcy flow segment, respectively.

[0013] Furthermore, the current actual flow rate through the heterogeneous seepage channel With the inversion of the actual critical flow rate Comparison is used to determine the current flow state of the water in the channel: like If so, the water flow in this channel is currently in Darcy flow mode; like If so, the water flow in the channel is currently in a non-Darcy flow state.

[0014] On the other hand, the present invention provides a system for determining the flow regime of heterogeneous seepage channels based on solute tracer method, comprising: The flow-characteristic time data pair acquisition module is used to acquire solute transport breakthrough curves of the target channel under different flow conditions, extract characteristic time from each solute transport breakthrough curve, and form multiple flow-characteristic time data pairs sorted by flow size with the corresponding flow. The initial critical flow range determination module is used to calculate the relative deceleration index reflecting the degree of solute transport deceleration based on the flow-characteristic time data pair, and to plot the curve of the relative deceleration index changing with flow rate. Based on the starting position of the curve where the relative deceleration index continuously decreases from 1 and falls below a certain set limit, the initial critical flow range of the flow regime change is determined. The actual critical flow inversion module is used to invert the actual critical flow within the initial critical flow range of the flow regime transition using a continuous condition equation model. The flow regime determination module is used to compare the inverted actual critical flow with the current actual flow to determine the current flow regime.

[0015] Compared with the prior art, the present invention has the following beneficial effects: 1. The data requirements of this invention are simple: the flow rate and characteristic time series can be obtained by simply conducting conventional tracer tests, which completely eliminates the dependence on the geometric parameters of the seepage channel and the fluid physical property parameters, and also eliminates the need for complex multi-stage differential pressure tests, which greatly reduces the application threshold and cost.

[0016] 2. The present invention has a clear physical concept and a solid theoretical foundation: it creatively defines the "relative deceleration index" and rigorously demonstrates that its theoretical value is 1 in the Darcy flow region, providing clear mathematical and physical criteria for deviations in non-Darcy flow regions.

[0017] 3. This invention is highly operable and yields intuitive results: the entire process is clearly defined and easy to automate through programming. The final output is a critical flow value that can be directly used for engineering judgment, providing a clear quantitative basis for optimizing engineering parameters and assessing safety risks. Attached Figure Description

[0018] Figure 1 This is a flowchart of a method according to an embodiment of the present invention; Figure 2This is a schematic diagram of the seepage channel tracing test of the dam foundation rock and soil in an embodiment of the present invention; Figure 3 This is a schematic diagram of feature time extraction in an embodiment of the present invention; Figure 4 This is a schematic diagram illustrating the determination of the critical flow range in an embodiment of the present invention; Figure 5 This is the lnq~lnt relationship diagram in an embodiment of the present invention. Detailed Implementation

[0019] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0020] Example 1 like Figure 1 As shown, this embodiment provides a method for determining the flow regime of heterogeneous seepage channels based on solute tracer method, including the following steps: Step S1: Obtain solute transport breakthrough curves of the target channel under different flow conditions, extract the characteristic time from each solute transport breakthrough curve, and form multiple flow-characteristic time data pairs sorted by flow rate with the corresponding flow rate; Step S2: Based on the flow rate-characteristic time data pair, calculate the relative deceleration index, which reflects the degree of solute transport deceleration, and plot the curve of the relative deceleration index as a function of flow rate. Based on the starting position of the relative deceleration index continuously decreasing from 1 and falling below a certain set limit in the curve, determine the initial critical flow rate range for the flow regime transition. Step S3: Invert the actual critical flow rate within the initial critical flow rate range of the flow regime transition using a continuity condition equation model; Step S4: Compare the inverted actual critical flow rate with the current actual flow rate to determine the current flow state.

[0021] In step S1 of this embodiment, a tracer injection experiment is conducted for the target heterogeneous channel based on natural hydrological events. A tracer is released upstream of the recharge point, and data on solute concentration changes over time are collected at downstream monitoring points. Continuous long-term monitoring is performed, covering at least one complete hydrological year or multiple typical rainfall events, to capture the natural flow variation process caused by changes in the natural hydraulic gradient and obtain solute transport breakthrough curves corresponding to each flow stage. .

[0022] In step S2 of this embodiment, feature time extraction is performed from each breakthrough curve. Extracting characteristic time that can represent solute transport velocity This characteristic time can be the initial arrival time (the time when the tracer is first detected), the half-concentration time (the time when the concentration reaches half at the inlet), or the average residence time calculated based on the method of moments, etc. The obtained... n Group data Sort by traffic volume from smallest to largest, denoted as ,in , They are number 1, 2, ... n Group data flow They are number 1, 2, ... n The characteristic time of the group data.

[0023] Based on relative deceleration index Critical range identification: defining the relative deceleration index ,in for and The slope of the numerical difference of the relationship curve characterizes the first... Feature time of group data Traffic Relative sensitivity to change. When the number of data points... n When the value is ≥3, the one-sided difference method can be used at the endpoints, while the three-point central difference method is preferred for calculating the intermediate points.

[0024] In the Darcy flow region, fracture seepage satisfies the following basic physical relationship: Flow rate-velocity relationship: ,in For the first Group data flow The cross-sectional area of ​​the fracture flow is [area]. The average flow rate; the relationship between solute transport characteristic time and flow rate: In tracer experiments, the solute transport characteristic time... (This embodiment uses half-concentration time) and average flow rate Inversely proportional, that is ,in This represents the effective path length for solute transport. Under constant flow conditions... It can be considered a constant. For a specific crack, It is a constant. From the above relationships, we can derive:

[0025] in It is a constant. Therefore, the Darcy flow time can be calculated. for:

[0026] In non-Darcy flow regions, solute transport paths become more complex due to effects such as inertial vortices. Follow The rate of decrease slows down, making Value less than 1. Draw. Follow The changing curve. By observing the curve shape, a reasonable threshold δ (0 < δ < 1, usually taken as 0.05~0.15, the specific value can be determined according to the Darcy flow stage) is set. The fluctuation range of the value is determined, and this stage is taken. (2 to 3 times the standard deviation of the value), when m consecutive (m≥2) values ​​appear. Value less than 1- delta When these values ​​show a decreasing trend, it can be determined that a change in flow regime has occurred. The first condition that satisfies this... <1- delta The flow rate corresponding to the point is denoted as the first non-Darcy flow point. The preceding flow point is the last Darcy flow point. Therefore, the critical flow range is determined as follows: The actual critical flow rate It then lies within this interval.

[0027] In step S3 of this embodiment, a two-segment power-law inversion based on self-normalized characteristic time is performed: for further precise positioning. Perform self-normalization analysis. Using the minimum flow rate... and its characteristic time Calculate the self-normalized flow rate as a reference baseline. and self-normalized characteristic time In a double logarithmic coordinate system ( ln t i ~ ln q i Plot data points in the ) .

[0028] Assuming a critical dimensionless flow rate Divide the data into two segments, when At time (Darcy section), the data should meet the following requirements. t i ∝ q i ^( a 1), i.e., power exponent a 1 = -1. When At time (non-Darcy flow segment), the data satisfy t i ∝ q i ^(a 2), where 0 < a 2<-1; where, 、 These are the power-law function coefficients for the non-Darcy flow section and the Darcy flow section, respectively; For self-normalized flow; For self-normalized characteristic time; 、 These are the power exponents of the power-law functions for the non-Darcy flow segment and the Darcy flow segment, respectively.

[0029] Based on the critical interval identified in step S3 This allows us to determine the corresponding data point partitioning. The last Darcy flow point... All previous data points were grouped into the Darcy flow segment for fitting, and the first non-Darcy flow point was excluded. All subsequent data points were categorized as non-Darcy flow segments for fitting. The least squares method was used to fit a power function to both data segments to obtain the fitting parameters for the Darcy flow segments. , Fitting parameters for non-Darcy flow sections , .

[0030] Based on the physical requirement of continuity, at the critical point At this point, the normalized characteristic times calculated from the two fitting functions should be equal, that is, satisfy the equation:

[0031] Taking the natural logarithm of both sides of the equation and simplifying, we can obtain the precise critical dimensionless flow rate:

[0032] Ultimately, the actual critical flow rate for the inversion is:

[0033] In step S4 of this embodiment, the comprehensive flow regime judgment and prediction involves calculating the precise critical flow rate. This serves as the final prediction. For the target channel that needs to be evaluated, its current workload is obtained. By comparison and : like If so, the current seepage is determined to be in Darcy flow state; like If so, it is determined that it has entered a non-Darcy flow state; If the difference between the two is within the preset engineering tolerance range, it can be determined as a critical transition state.

[0034] like Figure 2 As shown, this embodiment selects a major seepage channel within the dam foundation soil that has been monitored for a long time as the target object. Utilizing the seasonal natural fluctuations in reservoir water level, a tracer experiment is conducted upstream of the seepage channel. The tracer concentration change curve over time is recorded at the stable outlet point of this channel within the dam gallery, thereby simultaneously obtaining different flow rates. Solute transport breakthrough curve in the lower channel .

[0035] like Figure 3 As shown, the breakthrough curves obtained from each experiment In the middle, the extraction half-concentration time T As the characteristic time of this embodiment, all data are sorted in ascending order of flow rate to construct an ordered dataset for subsequent analysis: (2 m³ / h, 780 min), (4 m³ / h, 395 min), (6 m³ / h, 268 min), (8 m³ / h, 202 min), (10 m³ / h, 162 min), (15 m³ / h, 118 min), (20 m³ / h, 95 min).

[0036] Once the critical flow range is determined, the sum corresponding to each flow point is calculated. The relative deceleration index for each flow point is calculated using the three-point central difference method (with one-sided difference at the endpoints):

[0037] draw Follow Changing curves, such as Figure 4 As shown, when 3 consecutive The value satisfies: <1- delta ( delta = delta 0.15), and these Value varies with flow If it exhibits a monotonically decreasing trend, then the first one to be satisfied will be... <1- ln t The flow corresponding to the point is used as the first non-Darcy point. The point preceding it is the last Darcy point. Thus, it was initially determined that the critical flow rate was between [8 m³ / h, 10 m³ / h].

[0038] The critical flow rate is inverted using a continuity condition equation model, with the minimum flow rate as the minimum. and the corresponding characteristic time For reference, calculate the normalized flow for all data points. qi = Q i / Q 1 and normalized time t i = T i / T 1. Plot in a double logarithmic coordinate system ( ln q ~ Figure 5 Relationship curve, such as ​ As shown, the data points are divided into two groups based on the normalization critical interval [4,5] initially determined by S3: Darcy section: contains data points corresponding to (1,1)~(4,0.26).

[0039] Non-Darcy flow segment: contains data points corresponding to (5,0.21)~(10,0.12).

[0040] Apply the power function to each of the two sets of data. Perform fitting: To fit the Darcy section, the fitting function can be obtained by analyzing the first four data points. t =1.008× q ^(-0.975). Analysis of the non-Darcy flow section yields the fitting function for the last three data points. t =0.524× q ^(-0.528). At the critical flow rate At this point, the two piecewise function values ​​should be equal, that is, satisfy the equation:

[0041] Solving this equation yields the accurate critical normalized flow rate. :

[0042] Then, the predicted critical seepage flow rate for the flow regime transition in this seepage channel was calculated:

[0043] Based on the calculation results and the currently measured actual flow rate Compare the states: like If the flow rate is ≤9.07 m³ / h, then the current flow in the seepage channel is in Darcy flow mode.

[0044] like > If the flow rate is 9.07 m³ / h, then the current seepage flow in this seepage channel is in a non-Darcy flow state.

[0045] Example 2 This embodiment provides a system for determining the flow regime of heterogeneous seepage channels based on solute tracer method, including: The flow-characteristic time data pair acquisition module is used to acquire solute transport breakthrough curves of the target channel under different flow conditions, extract characteristic time from each solute transport breakthrough curve, and form multiple flow-characteristic time data pairs sorted by flow size with the corresponding flow. The initial critical flow range determination module is used to calculate the relative deceleration index reflecting the degree of solute transport deceleration based on the flow-characteristic time data pair, and to plot the curve of the relative deceleration index changing with flow rate. Based on the starting position of the curve where the relative deceleration index continuously decreases from 1 and falls below a certain set limit, the initial critical flow range of the flow regime change is determined. The actual critical flow inversion module is used to invert the actual critical flow within the initial critical flow range of the flow regime transition using a continuous condition equation model. The flow regime determination module is used to compare the inverted actual critical flow with the current actual flow to determine the current flow regime.

[0046] Although preferred embodiments of the invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including both the preferred embodiments and all changes and modifications falling within the scope of the invention.

[0047] Obviously, those skilled in the art can make various modifications and variations to the embodiments of the present invention without departing from the spirit and scope of the embodiments of the present invention. Thus, if these modifications and variations to the embodiments of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention also intends to include these modifications and variations.

[0048] All other parts not described in detail are existing technologies.

Claims

1. A method for determining the flow regime of heterogeneous seepage channels based on solute tracer method, characterized in that, Includes the following steps: Obtain solute transport breakthrough curves of the target channel under different flow conditions, extract feature time from each solute transport breakthrough curve, and form multiple flow-feature time data pairs sorted by flow rate with the corresponding flow rate; Based on the flow-characteristic time data pairs, a relative deceleration index reflecting the degree of solute transport deceleration is calculated, and a curve showing the relative deceleration index changing with flow rate is plotted. Based on the starting position of the relative deceleration index continuously decreasing from 1 and falling below a certain set limit in the curve, the initial critical flow rate range for flow regime transition is determined; the relative deceleration index is: in, For the first The relative deceleration index under various flow conditions Total number of flow-feature time data pairs; The first species and first Characteristic time of solute transport breakthrough curve under certain flow rate conditions; The first species and first Flow under certain flow conditions; The initial critical flow range for the flow regime transition is: Setting the first This point is the last Darcy flow point. If the first non-Darcy flow point is [the initial critical flow range is], then [the range is] ... ,in For the first Traffic per point For the first Traffic flow at each point; The actual critical flow rate is inverted within the initial critical flow rate range of the flow regime transition using a continuous condition equation model; The inverted actual critical flow rate is compared with the current actual flow rate to determine the current flow state.

2. The method for determining the flow regime of heterogeneous seepage channels based on solute tracer method according to claim 1, characterized in that, The characteristic time is the initial arrival time of the tracer, the half-concentration time of the tracer, or the average residence time of the tracer calculated based on the method of moments.

3. The method for determining the flow regime of heterogeneous seepage channels based on solute tracer method according to claim 1, characterized in that, The specific criterion for determining whether the relative deceleration index in the change curve continuously decreases from 1 and falls below a certain set limit is: when two consecutive relative deceleration index values ​​satisfy <1- δ Furthermore, if these relative deceleration index values ​​exhibit a monotonically decreasing or overall decreasing trend, then the first one to satisfy the condition of relative deceleration index < 1- δ The flow rate corresponding to the point is taken as the first non-Darcy flow point, and the point preceding the first non-Darcy flow point is taken as the last Darcy flow point. δ This is a preset threshold.

4. The method for determining the flow regime of heterogeneous seepage channels based on solute tracer method according to claim 3, characterized in that, The preset threshold δ The value ranges from 0.05 to 0.

1.

5. The method for determining the flow regime of heterogeneous seepage channels based on solute tracer method according to claim 1, characterized in that, The inverted actual critical flow rate includes: The flow rate of the last Darcy flow point All previous data points were grouped into the Darcy flow segment for fitting, and the power-law function of the Darcy flow segment was obtained. The flow rate of the first non-Darcy flow point All subsequent data points were categorized into non-Darcy flow segments for fitting, and the power-law function of the non-Darcy flow segments was obtained. The actual critical flow rate is determined by solving the condition equations for the continuity of the power-law function at specific critical points in both the non-Darcy and Darcy sections.

6. The method for determining the flow regime of heterogeneous seepage channels based on solute tracer method according to claim 5, characterized in that, The power-law function of the non-Darcy flow segment is: ; The Darcy flow segment power-law function is: ; in, 、 These are the power-law function coefficients for the non-Darcy flow section and the Darcy flow section, respectively; For self-normalized flow; For self-normalized characteristic time; 、 These are the power exponents of the power-law functions for the non-Darcy flow segment and the Darcy flow segment, respectively.

7. The method for determining the flow regime of heterogeneous seepage channels based on solute tracer method according to claim 1, characterized in that, Current actual flow rate through the heterogeneous seepage channel With the inversion of the actual critical flow rate Comparison is used to determine the current flow state of the water in the channel: like If so, the water flow in this channel is currently in Darcy flow mode; like If so, the water flow in the channel is currently in a non-Darcy flow state.

8. A system for determining the flow regime of a heterogeneous seepage channel based on solute tracer method, characterized in that, include: The flow-characteristic time data pair acquisition module is used to acquire solute transport breakthrough curves of the target channel under different flow conditions, extract characteristic time from each solute transport breakthrough curve, and form multiple flow-characteristic time data pairs sorted by flow size with the corresponding flow. The initial critical flow range determination module is used to calculate the relative deceleration index reflecting the degree of solute transport deceleration based on the flow-characteristic time data pair, and to plot the curve of the relative deceleration index changing with flow rate. Based on the starting position of the curve where the relative deceleration index continuously decreases from 1 and falls below a certain set limit, the initial critical flow range of the flow regime change is determined. The actual critical flow inversion module is used to invert the actual critical flow within the initial critical flow range of the flow regime transition using a continuous condition equation model. The flow regime determination module is used to compare the inverted actual critical flow with the current actual flow to determine the current flow regime. The heterogeneous seepage channel flow regime determination system based on solute tracer method is used to perform the steps in the heterogeneous seepage channel flow regime determination method based on solute tracer method as described in any one of claims 1-7.

Images