Groundwater recharge inversion system of irrigation and drainage dual-purpose lining channel based on water level monitoring

Abstract

This invention discloses a groundwater recharge inversion system for irrigation and drainage lined canals based on water level monitoring, belonging to the field of water conservancy facilities technology. The invention includes a lining state detection unit for acquiring the lining state of the canal lining surface; a groundwater monitoring unit for setting monitoring wells along the canal according to its location distribution; acquiring the water level of the monitoring wells; obtaining the groundwater level distribution in a local area along the canal based on the water level of the monitoring wells; and an inversion calculation unit for obtaining the relationship between the groundwater level distribution along the canal and the lining state of the canal lining surface, and inverting the groundwater level distribution along the canal based on the relationship between the lining state and the groundwater level distribution in the local area along the canal lining surface. This invention provides a comprehensive and accurate understanding of the groundwater state along the canal.

Description

Technical Field

[0001] This invention belongs to the field of water conservancy facilities technology, and in particular relates to a groundwater recharge inversion system for irrigation and drainage lining channels based on water level monitoring. Background Technology

[0002] In the fields of water conservancy engineering and groundwater management, liningd canals for both irrigation and drainage are commonly used facilities. However, the interaction between this type of canal and groundwater is complex, involving both the canal's recharge of groundwater and the groundwater's discharge of water into the canal. The strength of these two effects varies with time, groundwater level, and canal water level.

[0003] However, traditional groundwater recharge inversion methods are mostly based on hydraulic models for mathematical simulation, requiring a large number of geological and channel parameters. Obtaining these parameters usually requires extensive field investigations and experiments, which is not only labor-intensive but also costly. In addition, the accuracy of these methods is often limited due to the complexity of groundwater flow processes.

[0004] Patent CN114935385A discloses a groundwater level surveying device for hydrogeological exploration, including a base frame, a support frame at the top of the base frame, a limiting sleeve at the middle of the top of the base frame, a servo motor at the top of the base frame, the output end of the servo motor being connected to a rotating shaft via a reducer, one end of the rotating shaft extending to the inside of the limiting sleeve and having a transmission gear, a transmission screw movably mounted inside the limiting sleeve, a drill rod measuring device at the bottom of the transmission screw, support legs on both the left and right sides of the bottom of the base frame, a limiting fixing device at the bottom of the base frame, sliding rods movably mounted on the outer walls of the left and right sides of the limiting fixing device, and pressure plates at the bottom of the sliding rods. This solution relies on manual inspection of groundwater levels at each location to detect groundwater, but it depends on the experience of the surveyors for sampling and is difficult to accurately and comprehensively reflect the groundwater condition. Summary of the Invention

[0005] The purpose of this invention is to provide a groundwater recharge inversion system for irrigation and drainage lining channels based on water level monitoring. By sampling and inverting the channel lining status and the groundwater along the channel, the groundwater status along the channel can be obtained comprehensively and accurately.

[0006] To solve the above-mentioned technical problems, the present invention is achieved through the following technical solution:

[0007] This invention provides a groundwater recharge inversion system for irrigation and drainage lined canals based on water level monitoring, comprising:

[0008] The channel information acquisition unit is used to acquire the location distribution of channels;

[0009] A lining status detection unit is used to obtain the lining status of the lining surface of the channel;

[0010] A groundwater monitoring unit is used to install monitoring wells along the channel according to the channel's location distribution;

[0011] Obtain the water level of the monitoring well;

[0012] The distribution of groundwater level in a local area along the canal is obtained based on the water level of the monitoring well;

[0013] The inversion calculation unit is used to obtain the relationship between the groundwater level distribution along the channel and the lining state of the channel based on the lining state of the channel lining surface and the groundwater level distribution in a local area along the channel.

[0014] The groundwater level distribution along the channel is obtained by inverting the relationship between the lining state of the channel lining surface and the groundwater level distribution along the channel with respect to the lining state of the channel lining surface.

[0015] This invention comprehensively and accurately obtains the groundwater status along a canal by sampling and inverting the lining condition and groundwater levels in the surrounding area. In implementation, the location and distribution information of the canal are first acquired. Then, the condition of the canal lining surface is obtained through a lining condition detection unit. Next, monitoring wells are set up along the canal according to its location, and the water levels in these wells are obtained, thus revealing the local groundwater level distribution along the canal. Finally, based on the relationship between the lining surface condition and the groundwater level distribution, an inversion calculation is performed to obtain the groundwater level distribution along the entire canal. Through this process, comprehensive and accurate monitoring of the canal lining condition and the groundwater level in the surrounding area can be achieved.

[0016] Of course, any product implementing this invention does not necessarily need to achieve all of the advantages described above at the same time. Attached Figure Description

[0017] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the description of the embodiments 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.

[0018] Figure 1 This is a schematic diagram of the functional units and information interaction of the groundwater recharge inversion system for irrigation and drainage lining channels based on water level monitoring in one embodiment of the present invention.

[0019] Figure 2 This is a schematic diagram of the steps of an embodiment of the groundwater recharge inversion system for irrigation and drainage lining channels based on water level monitoring described in this invention.

[0020] Figure 3 This is a flowchart illustrating step S2 of the present invention in one embodiment;

[0021] Figure 4 This is a flowchart illustrating step S21 of the present invention in one embodiment;

[0022] Figure 5 This is a flowchart illustrating step S212 of the present invention in one embodiment;

[0023] Figure 6 This is a flowchart illustrating step S213 of the present invention in one embodiment;

[0024] Figure 7 This is a flowchart illustrating step S3 of the present invention in one embodiment;

[0025] Figure 8 This is a flowchart illustrating step S32 of the present invention in one embodiment;

[0026] Figure 9 This is a flowchart illustrating step S33 of the present invention in one embodiment;

[0027] Figure 10 This is a flowchart illustrating step S34 of the present invention in one embodiment. Figure 1 ;

[0028] Figure 11 This is a flowchart illustrating step S34 of the present invention in one embodiment. Figure 2 .

[0029] The attached diagram lists the components represented by each number as follows:

[0030] 1-Channel information acquisition unit, 2-Liner status detection unit, 3-Groundwater monitoring unit, 4-Inversion calculation unit. Detailed Implementation

[0031] 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 embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0032] In order to comprehensively and accurately quantify the groundwater replenishment effect of lined channels, the present invention provides the following solution.

[0033] Please see Figures 1 to 2 As shown, this invention provides a groundwater recharge inversion system for irrigation and drainage lined channels based on water level monitoring. Functionally, it can be divided into a channel information acquisition unit 1, a lining status detection unit 2, a groundwater monitoring unit 3, and an inversion calculation unit 4. In practical applications, the channel information acquisition unit 1 can be a data interface for reading input data, the lining status detection unit 2 can be a handheld mobile device with a camera, the groundwater monitoring unit 3 can be a drilling exploration device, and the inversion calculation unit 4 can be a computing server, or of course, a portable computing device.

[0034] In the specific implementation process, the channel information acquisition unit 1 first executes step S1 to obtain the location distribution of the channel. Next, the lining status detection unit 2 executes step S2 to obtain the lining status of the channel's lining surface. Then, the groundwater monitoring unit 3 executes step S3 to set up monitoring wells along the channel according to the channel's location distribution, and then executes step S4 to obtain the water level in the monitoring wells. Next, step S4 is executed to obtain the groundwater level distribution in a local area along the channel based on the water level in the monitoring wells. In practical applications, the obtained groundwater level distribution can also be supplemented and corrected using technologies such as infrared remote sensing to obtain a more accurate groundwater level distribution in a local area along the channel.

[0035] Then, the inversion calculation unit 4 can execute step S5 to obtain the relationship between the groundwater level distribution along the channel and the lining state of the channel lining surface, based on the lining state of the channel lining surface and the groundwater level distribution in local areas along the channel. Finally, step S6 can be executed to invert the groundwater level distribution along the channel lining surface based on the relationship between the lining state of the channel lining surface and the groundwater level distribution along the channel lining surface.

[0036] The above steps, in practice, first require obtaining the geographical distribution of the channel. Then, a lining condition monitoring unit is used to assess the condition of the channel lining surface. Next, monitoring wells are established along the channel's location to obtain water level data, thereby understanding the local groundwater level distribution along the channel. Finally, inversion calculations are used to clarify the correlation between the lining surface condition and the groundwater level distribution, inferring the groundwater level distribution along the channel, thus achieving comprehensive and accurate monitoring of the channel lining condition and the groundwater situation in the surrounding area.

[0037] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0038] # Import the required libraries

[0039] import numpy as np

[0040] #Function Section

[0041] defget_canal_position():

[0042] "″"

[0043] Location distribution of acquisition channels

[0044] "″"

[0045] #Implementation code

[0046] pass defget_canal_status():

[0047] "″"

[0048] Obtain the lining condition of the channel lining surface

[0049] "″"

[0050] #Implementation code

[0051] pass defset_monitoring_wells(canal_position):

[0052] "″"

[0053] Monitoring wells will be installed along the canal based on its location.

[0054] "″"

[0055] #Implementation code

[0056] pass defget_well_water_level():

[0057] "″"

[0058] Obtain the water level of the monitoring well

[0059] "″"

[0060] #Implementation code

[0061] pass defget_local_groundwater_distribution(well_water_level):

[0062] "″"

[0063] The distribution of groundwater levels in local areas along the canal was obtained based on the water levels in the monitoring wells.

[0064] "″"

[0065] #Implementation code

[0066] pass

[0067] defget_relationship(canal_status,local_groundwater_distribution):

[0068] "″"

[0069] The relationship between the groundwater level distribution along the canal and the lining condition of the canal is obtained based on the lining condition of the canal lining surface and the groundwater level distribution in local areas along the canal.

[0070] "″"

[0071] #Implementation code

[0072] pass

[0073] definverse_calculation(relationship):

[0074] "″"

[0075] Based on the relationship between the lining condition of the canal and the groundwater level distribution along the canal and the lining condition, the groundwater level distribution along the canal can be derived.

[0076] "″"

[0077] #Implementation code

[0078] pass

[0079] #Main Program Section

[0080] canal_position=get_canal_position()

[0081] canal_status=get_canal_status()

[0082] set_monitoring_wells(canal_position)

[0083] well_water_level=get_well_water_level()

[0084] local_groundwater_distribution=get_local_groundwater_distribution(well_water_level)

[0085] relationship=get_relationship(canal_status,local_groundwater_distribution)

[0086] inverse_calculation(relationship)

[0087] Please see Figure 3 As shown, since the channel lining is constructed using bricks, stones, or other materials, the limitations of the masonry process inevitably result in mortar joints of varying widths, which significantly impacts channel leakage. Therefore, in the implementation of step S2, step S21 first selects the lining inspection area on the channel lining surface. Next, step S22 obtains the length and width of the mortar joints within the inspection area. Then, step S23 calculates the total area of ​​the mortar joints within the inspection area based on their length and width. Next, step S24 obtains the area of ​​the mortar joints per unit area within the inspection area as the lining state of the channel lining surface. Finally, step S25 calculates the lining state at any point along the channel based on the location distribution and lining state of the inspection area. Quantifying the lining state by statistically calculating the total mortar joint area facilitates subsequent calculations.

[0088] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0089] # Import the required libraries

[0090] import numpy as np

[0091] #Function Section

[0092] def select_inspection_area():

[0093] "″"

[0094] Select the lining inspection area on the lining surface of the channel.

[0095] "″"

[0096] #Implementation code

[0097] pass

[0098] defget_crack_length_and_width():

[0099] "″"

[0100] Obtain the length and width of the lining joints within the lining inspection area.

[0101] "″"

[0102] #Implementation code

[0103] pas

[0104] defget_crack_area(length,width):

[0105] "″"

[0106] The total area of ​​the lining joints within the lining inspection area is obtained based on the length and width of the joints.

[0107] "″"

[0108] return length*width

[0109] defget_canal_status(crack_area,unit_area):

[0110] "″"

[0111] The area of ​​the lining joints per unit area within the lining inspection zone is used as the lining condition of the channel lining surface.

[0112] "″"

[0113] return crack_area / unit_area

[0114] defget_status_distribution(status,inspection_area_distribution):

[0115] "″"

[0116] The lining condition of the lining inspection area is obtained at any point along the channel based on the location distribution of the lining inspection area and the lining condition of the lining inspection area.

[0117] "″"

[0118] #Implementation code

[0119] pass

[0120] #Main Program Section

[0121] select_inspection_area()

[0122] length,width=get_crack_length_and_width()

[0123] crack_area=get_crack_area(length,width)

[0124] unit_area = 1.0 # The unit area is 1 hectare

[0125] canal_status=get_canal_status(crack_area,unit_area)

[0126] inspection_area_distribution = [] # The location distribution of the lining inspection area needs to be provided based on the actual situation.

[0127] get_status_distribution(canal_status,inspection_area_distribution)

[0128] Please see Figure 4 As shown, due to the long length of the channel, it is difficult to inspect each lining surface individually. To obtain a more comprehensive lining condition with limited inspection capabilities, it is necessary to select representative lining inspection areas. Therefore, in the implementation of step S21 above, step S211 can be executed first to obtain the channel's transport direction based on its location distribution. Next, step S212 can be executed to select several sampling points along the channel's transport direction. Next, step S213 can be executed to delineate the lining inspection area within the channel's lining surface based on the sampling points and the channel's transport direction. Finally, step S214 can be executed to obtain several lining inspection areas.

[0129] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0130] # Import the required libraries

[0131] import numpy as np

[0132] #Function Section

[0133] defget_canal_direction(canal_position):

[0134] "″"

[0135] Based on the location distribution of the channels, obtain the delivery direction of the channels.

[0136] "″"

[0137] #Implementation code

[0138] pass

[0139] defget_sample_points(canal_direction):

[0140] "″"

[0141] Select several sampling points along the transport direction of the channel.

[0142] "″"

[0143] #Implementation code

[0144] pass

[0145] defdefine_inspection_area(sample_points,canal_direction):

[0146] "″"

[0147] Based on the sampling points and the transport direction of the channel, the lining inspection area is delineated within the lining surface of the channel.

[0148] "″"

[0149] #Implementation code

[0150] pass

[0151] defget_inspection_areas(inspection_area):

[0152] "″"

[0153] Several lining inspection areas were obtained.

[0154] "″"

[0155] #Implementation code

[0156] pass

[0157] #Main Program Section

[0158] canal_position = [] # The location distribution of the channel; data needs to be provided based on the actual situation.

[0159] canal_direction=get_canal_direction(canal_position)

[0160] sample_points=get_sample_points(canal_direction)

[0161] inspection_area=define_inspection_area(sample_points,canal_direction)

[0162] inspection_areas=get_inspection_areas(inspection_area)

[0163] Please refer to Figure #. To obtain a representative lining inspection area, representative sampling points are needed. Therefore, in implementing step S212 above, step S2121 can be executed first, continuously counting the number of lining blocks used for lining the channel within a unit area along the channel's transport direction. Next, step S2122 can be executed to determine if the rate of change of the number of lining blocks within a unit area exceeds a set quantity. If so, step S2123 can be executed to set sampling points at locations where the rate of change of the number of lining blocks within a unit area exceeds the set quantity. Otherwise, the process can return to step S2121.

[0164] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0165] # Import the required libraries

[0166] import numpy as np

[0167] #Function Section

[0168] defcount_blocks(canal_direction):

[0169] "″"

[0170] Along the transport direction of the channel, continuously count the number of blocks per unit area used for lining the channel.

[0171] "″"

[0172] #Implementation code

[0173] pass

[0174] defcheck_rate_change(block_count,threshold):

[0175] "″"

[0176] Determine whether the rate of change of the number of lining blocks per unit area exceeds the set quantity.

[0177] "″"

[0178] #Calculate the rate of change

[0179] rate_change=np.diff(block_count) / block_count[:-1]

[0180] return np.abs(rate_change)>threshold

[0181] defset_sample_points(change_rate):

[0182] "″"

[0183] Sampling points were set at locations where the rate of change of the number of lining blocks per unit area exceeded a predetermined number.

[0184] "″"

[0185] #Implementation code

[0186] pass

[0187] #Main Program Section

[0188] canal_direction = [] # The direction of channel delivery; data needs to be provided based on the actual situation.

[0189] block_count=count_blocks(canal_direction)

[0190] threshold=0.1 #Sets the threshold for the rate of change of quantity

[0191] change_rate=check_rate_change(block_count,threshold)

[0192] ifnp.any(change_rate):

[0193] sample_points=set_sample_points(change_rate)

[0194] else:

[0195] block_count=count_blocks(canal_direction)

[0196] Please see Figure 6 As shown, in order to obtain the corresponding lining detection area based on the sampling point, the above step S213 can first be implemented by step S2131 to obtain the channel cross-section that passes through the sampling point and is perpendicular to the channel's conveying direction. Finally, step S2132 can be implemented to translate the channel cross-section along the channel's conveying direction by a set length, obtaining the range of the channel cross-section sweeping across the channel's lining surface as the lining detection area corresponding to the sampling point.

[0197] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0198] # Import the required libraries

[0199] import numpy as np

[0200] #Function Section

[0201] defget_canal_cross_section(sample_point,canal_direction):

[0202] "″"

[0203] Obtain the channel cross-section with oversampling points perpendicular to the channel's transport direction.

[0204] "″"

[0205] #Implementation code

[0206] pass

[0207] deftranslate_cross_section(cross_section,canal_direction,length):

[0208] "″"

[0209] The channel cross-section is translated by a set length along the conveying direction of the channel, and the range of the channel lining surface swept by the channel cross-section is taken as the lining detection area corresponding to the sampling point.

[0210] "″"

[0211] #Implementation code

[0212] pass

[0213] #Main Program Section

[0214] sample_point = [] # Sampling point, data needs to be provided according to the actual situation.

[0215] canal_direction = [] # The direction of channel delivery; data needs to be provided based on the actual situation.

[0216] length = 0.1 # Translation length, which needs to be set according to the actual situation.

[0217] cross_section=get_canal_cross_section(sample_point,canal_direction)

[0218] inspection_area=translate_cross_section(cross_section,canal_direction,length)

[0219] Please see Figure 7 As shown, the monitoring wells can include exploration wells and supplementary wells. First, exploration wells can be used to conduct preliminary detection of the water level along the canal, and then supplementary wells can be used for further detection. Therefore, in the implementation of step S3 above, step S31 can be executed first to obtain the set number of exploration wells. Next, step S32 can be executed to allocate the set number of exploration wells according to the lining state of the canal lining surface, obtaining the location distribution of the exploration wells along the canal. Next, step S33 can be executed to obtain the groundwater level change rate along the canal based on the water level of the exploration wells. Finally, step S34 can be executed to set up supplementary wells according to the groundwater level change rate along the canal.

[0220] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0221] # Import the required libraries

[0222] import numpy as np

[0223] #Function Section

[0224] defallocate_survey_wells(well_count,lining_state):

[0225] "″"

[0226] The number of logging wells is allocated based on the lining condition of the channel, resulting in the distribution of logging well locations along the channel.

[0227] "″"

[0228] #Implementation code

[0229] pass

[0230] defget_groundwater_rate_change(survey_well_levels):

[0231] "″"

[0232] The rate of change of groundwater level along the canal was obtained based on the water level of the survey wells.

[0233] "″"

[0234] #Calculate the rate of change

[0235] rate_change=np.diff(survey_well_levels) / survey_well_levels[:-1]

[0236] return rate_change

[0237] defset_supplementary_wells(rate_change):

[0238] "″"

[0239] Replenishment wells are set up according to the groundwater level change rate along the canal.

[0240] "″"

[0241] #Implementation code

[0242] pass

[0243] #Main Program Section

[0244] well_count = 10 # The set number of exploration wells

[0245] lining_state = [] # The lining state of the channel lining surface, which needs to be provided according to the actual situation. survey_well_distribution = allocate_survey_wells(well_count,lining_state)

[0246] survey_well_levels = [] # Water level of the survey well. Data needs to be provided based on the actual situation.

[0247] rate_change=get_groundwater_rate_change(survey_well_levels)

[0248] supplementary_wells=set_supplementary_wells(rate_change)

[0249] Please see Figure 8 As shown, to improve the comprehensiveness of well detection with a limited number of wells, step S32 can be implemented by first executing step S321 to evenly divide the channel into multiple segments. Next, step S322 can be executed to obtain the ratio of the total area of ​​the joints between the multiple channel segments based on the lining condition of the channel lining. Next, step S323 can be executed to distribute the set number of wells along the channel segments based on the ratio of the total area of ​​the joints, obtaining the number of wells allocated to each channel segment's area. Next, step S324 can be executed to obtain the set survey distance. Next, step S325 can be executed to define the area within the set distance of the channel as the survey area, obtaining the survey sub-area corresponding to each channel segment. Next, step S326 can be executed to evenly distribute the wells allocated to each channel segment's area within the corresponding survey sub-area, obtaining the location of the wells within each survey sub-area. Finally, step S327 can be executed to summarize the locations of all wells in the survey sub-regions, thus obtaining the distribution of well locations along the channel.

[0250] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0251] defdivide_canal_into_segments(canal):

[0252] "″"

[0253] The channel is evenly divided into multiple sections.

[0254] "″"

[0255] #Implementation code

[0256] pass

[0257] defcalculate_area_ratio(segments):

[0258] "″"

[0259] The ratio of the total area of ​​the lining joints between multiple channel sections is obtained based on the lining condition of the channel lining surface.

[0260] "″"

[0261] #Implementation code

[0262] pass

[0263] defdistribute_wells(well_count,area_ratio):

[0264] "″"

[0265] The number of exploration wells is allocated among the areas along the canal sections based on the ratio of the total area of ​​the joints between multiple canal sections.

[0266] The number of survey wells allocated along the route of each canal section was obtained.

[0267] "″"

[0268] #Implementation code

[0269] pass

[0270] defdefine_survey_subareas(segments,survey_distance):

[0271] "″"

[0272] The area within a set distance of the canal is designated as the survey area, resulting in a survey sub-area for each canal segment.

[0273] "″"

[0274] #Implementation code

[0275] pass

[0276] defdistribute_wells_in_subareas(wells_per_segment,survey_subareas):

[0277] "″"

[0278] The survey wells allocated along the canal section are evenly distributed within their corresponding survey sub-regions to obtain the location of the survey wells within each sub-region.

[0279] "″"

[0280] #Implementation code

[0281] pass

[0282] defaggregate_well_positions(well_positions_subareas):

[0283] "″"

[0284] The locations of all wells within the survey sub-regions are summarized to obtain the distribution of well locations along the canal.

[0285] "″"

[0286] #Implementation code

[0287] pass

[0288] Please see Figure 9 As shown, in order to quantify the rate of change of groundwater level along the canal, step S33 above can first be implemented by dividing the area along the canal into multiple measurement zones, where each measurement zone includes multiple survey wells. Next, step S332 can be implemented to obtain the water levels of the multiple survey wells within each measurement zone. Finally, step S333 can be implemented to calculate the standard deviation of the water levels of the multiple survey wells within the measurement zone as the rate of change of groundwater level in the measurement zone.

[0289] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0290] import numpy as np

[0291] defdivide_into_measurement_areas(canal_line):

[0292] "″"

[0293] The area along the canal was divided into multiple measurement zones, each of which included multiple logging wells.

[0294] "″"

[0295] #Implementation code

[0296] pass

[0297] defget_well_water_levels(measurement_areas):

[0298] "″"

[0299] Obtain the water level of multiple survey wells within each measurement area.

[0300] "″"

[0301] #Implementation code

[0302] pass

[0303] defcalculate_water_level_variance(water_levels):

[0304] "″"

[0305] The standard deviation of water levels in multiple wells within the measurement area is calculated as the rate of change of water level in the measurement area.

[0306] "″"

[0307] return{area:np.std(levels)for area,levels in water_levels.items()}

[0308] Please see Figure 10 As shown, to improve the detection effect of supplementary wells, the above-mentioned step S34 can first be implemented by executing step S341 to obtain the set number of supplementary wells. Next, step S342 can be executed to distribute the supplementary wells among multiple measurement areas according to the ratio of water level change rates between measurement areas, obtaining the number of supplementary wells in each measurement area. Next, step S343 can be executed to distribute the supplementary wells evenly or randomly within the measurement area according to the number of supplementary wells in each measurement area, obtaining the location of the supplementary wells in each measurement area. Finally, step S344 can be executed to summarize the locations of the supplementary wells in each measurement area to obtain the location distribution of the supplementary wells along the channel.

[0309] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0310] defget_supplement_well_number():

[0311] "″"

[0312] Get the set number of replenishment wells

[0313] "″"

[0314] #Implementation code

[0315] pass

[0316] def distribute_supplement_wells(measurement_areas,water_level_variances,supplement_well_number):

[0317] "″"

[0318] Based on the ratio of water level change rates between measurement areas, supplementary wells are allocated among multiple measurement areas to obtain the number of supplementary wells in each measurement area.

[0319] "″"

[0320] #Implementation code

[0321] pass

[0322] defdistribute_wells_in_area(supplement_well_distribution):

[0323] "″"

[0324] Based on the number of supplementary wells in each measurement area, the supplementary wells are distributed evenly or randomly within the measurement area to obtain the location of the supplementary wells in each measurement area.

[0325] "″"

[0326] #Implementation code

[0327] pass

[0328] defaggregate_well_positions(supplement_well_positions):

[0329] "″"

[0330] The locations of supplementary wells within each measurement area are summarized to obtain the distribution of supplementary well locations along the canal.

[0331] "″"

[0332] #Implementation code

[0333] pass

[0334] Please see Figure 11As shown, water seeping from the channel may also replenish groundwater in areas beyond the survey distance. To obtain the groundwater status beyond the survey distance, the monitoring well also includes a secondary replenishment well. In the implementation of step S34, step S345 can first be executed to obtain the measurement area whose straight-line distance from the channel is equal to the survey distance and mark it as the edge measurement area. Next, step S346 can be executed to determine whether the water level change rate of the edge measurement area is less than a set value. If so, step S347 can be executed without taking any action; otherwise, step S348 can be executed to set a new edge measurement area outside the existing edge measurement area and install a replenishment well in the newly set edge measurement area.

[0335] To provide supplementary explanations of the above steps, due to space limitations, only the source code of some functional modules is provided, with explanations in the comments.

[0336] defidentify_edge_measurement_area(distance_to_canal,measurement_area):

[0337] "″"

[0338] The area where the straight-line distance from the channel equals the survey distance is marked as the edge measurement area.

[0339] "″"

[0340] #Implementation code

[0341] pass

[0342] defcheck_water_level_variance(edge_measurement_area,threshold):

[0343] "″"

[0344] Determine whether the rate of change of water level in the edge measurement area is less than the set value.

[0345] "″"

[0346] #Implementation code

[0347] pass

[0348] defadd_new_edge_measurement_area(edge_measurement_area):

[0349] "″"

[0350] Set up a new edge measurement region outside the existing edge measurement region.

[0351] "″"

[0352] #Implementation code

[0353] pass

[0354] defset_supplement_wells_in_new_area(new_edge_measurement_area):

[0355] "″"

[0356] Set up supplementary wells in the newly established edge measurement area

[0357] "″"

[0358] #Implementation code

[0359] pass

[0360] #Setting settings

[0361] threshold = 0.05

[0362] #The area where the straight-line distance from the channel equals the survey distance is marked as the edge of the measurement area.

[0363] edge_measurement_area=identify_edge_measurement_area(distance_to_canal,measurement_area)

[0364] #Determine if the rate of change of water level in the edge measurement area is less than the set value.

[0365] ifcheck_water_level_variance(edge_measurement_area,threshold):

[0366] #If so, no action will be taken.

[0367] pass

[0368] else:

[0369] #If not, then set a new edge measurement region outside the existing edge measurement region.

[0370] new_edge_measurement_area=add_new_edge_measurement_area(edge_measurement_area)

[0371] #And set up supplementary wells in the newly established edge measurement area.

[0372] set_supplement_wells_in_new_area(new_edge_measurement_area)

[0373] In summary, this scheme can comprehensively and accurately grasp the groundwater conditions along the canal. During implementation, the first step is to collect information on the canal's location and distribution. Next, a lining condition monitoring unit is used to investigate the condition of the canal lining surface. Then, based on the canal's location and distribution, monitoring wells are deployed along the canal to obtain water level information, thereby understanding the local groundwater level distribution along the canal. Finally, inversion calculations are used to extrapolate the groundwater level distribution along the canal based on the relationship between the lining surface condition and the groundwater level distribution. Through these methods, the canal lining condition and the groundwater in the surrounding area can be monitored comprehensively and accurately.

[0374] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of apparatus, systems, methods, and computer program products according to various embodiments of this application. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of an instruction containing one or more executable instructions for implementing a specified logical function. In some alternative implementations, the functions marked in the blocks may occur in a different order than those shown in the drawings. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved.

[0375] It should also be noted that each block in the block diagram and / or flowchart, as well as combinations of blocks in the block diagram and / or flowchart, can be implemented using hardware that performs the corresponding function or action, such as circuits or ASICs (Application Specific Integrated Circuits), or using a combination of hardware and software, such as firmware.

[0376] Although the invention has been described herein in conjunction with various embodiments, those skilled in the art will understand and implement other variations of the disclosed embodiments by reviewing the accompanying drawings, disclosure, and appended claims in carrying out the claimed invention. In the claims, the word "comprising" does not exclude other components or steps, and "a" or "an" does not exclude a plurality. A single processor or other unit can implement several functions listed in the claims. While different dependent claims may recite certain measures, this does not mean that these measures cannot be combined to produce good results.

[0377] The various embodiments of this application have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope of the described embodiments. The terminology used herein is chosen to best explain the principles, practical application, or improvement of the technology in the market, or to enable others skilled in the art to understand the embodiments disclosed herein.

Claims

1. A groundwater recharge inversion system for irrigation and drainage lined canals based on water level monitoring, characterized in that, include, The channel information acquisition unit is used to acquire the location distribution of channels; The lining status detection unit is used to obtain the lining status of the lining surface of the channel and to quantitatively express the lining status by statistically calculating the total area of ​​the lining joints. A groundwater monitoring unit is used to install monitoring wells along the channel according to the channel's location distribution; Obtain the water level of the monitoring well; The distribution of groundwater level in a local area along the canal is obtained based on the water level of the monitoring well; The inversion calculation unit is used to obtain the relationship between the groundwater level distribution along the channel and the lining state of the channel based on the lining state of the channel lining surface and the groundwater level distribution in a local area along the channel. The groundwater level distribution along the channel is obtained by inverting the relationship between the lining state of the channel lining surface and the groundwater level distribution along the channel with respect to the lining state of the channel lining surface.

2. The system according to claim 1, characterized in that, The step of obtaining the lining state of the lining surface of the channel. include, Select a lining inspection area on the lining surface of the channel; Obtain the length and width of the lining joints within the lining detection area; The total area of ​​the joints within the lining detection area is obtained based on the length and width of the joints within the lining detection area; The area of ​​the lining joints per unit area in the lining detection area is obtained as the lining state of the channel lining surface. The lining state of the lining surface at any point along the channel is obtained based on the location distribution of the lining detection area and the lining state of the lining detection area.

3. The system according to claim 2, characterized in that, The step of selecting the lining inspection area on the lining surface of the channel includes, Based on the location distribution of the channels, the conveying direction of the channels is obtained; Several sampling points are selected along the transport direction of the channel; Based on the sampling points and the transport direction of the channel, a lining detection area is delineated within the lining surface of the channel; Several lining inspection areas were obtained.

4. The system according to claim 3, characterized in that, The step of selecting several sampling points along the transport direction of the channel, include, Along the transport direction of the channel, the number of blocks used to lay the lining surface of the channel per unit area is continuously counted; Determine whether the rate of change of the number of lining blocks per unit area exceeds the set quantity; If so, sampling points are set at locations where the rate of change of the number of lining blocks per unit area exceeds the set number; If not, the number of blocks used for lining the channel per unit area will be continuously counted.

5. The system according to claim 3, characterized in that, The step of defining a lining detection area within the lining surface of the channel based on the sampling points and the channel's transport direction includes, Obtain the channel cross-section that passes through the sampling point and is perpendicular to the transport direction of the channel; The channel cross section is shifted by a set length along the conveying direction of the channel to obtain the range of the channel cross section sweeping across the lining surface of the channel, which is taken as the lining detection area corresponding to the sampling point.

6. The system according to claim 2, characterized in that, The step of setting monitoring wells along the channel according to the channel's location distribution includes: The monitoring wells include exploration wells and supplementary wells, and the set number of exploration wells is obtained; The number of exploration wells is allocated according to the lining state of the channel lining surface to obtain the location distribution of the exploration wells along the channel. The rate of change of groundwater level along the canal is obtained based on the water level of the survey well. Supplementary wells are set up according to the rate of change of groundwater level along the channel.

7. The system according to claim 6, characterized in that, The step of allocating a predetermined number of logging wells based on the lining condition of the channel lining surface to obtain the location distribution of the logging wells along the channel includes, The channel is evenly divided into multiple channel segments; The ratio of the total area of ​​the lining joints between multiple channel sections is obtained based on the lining state of the channel lining surface; The predetermined number of survey wells is allocated among the areas along the canal sections based on the ratio of the total area of ​​the joints between the multiple canal sections, so as to obtain the number of survey wells allocated to the area along the canal section for each canal section. Obtain the set survey distance; The area within a set length of the channel is designated as the survey area, thus obtaining the survey sub-area corresponding to each channel segment; The survey wells allocated along the route of each canal section are evenly distributed within the corresponding survey sub-regions to obtain the location of the survey wells within each survey sub-region; The locations of the wells in all the survey sub-regions are summarized to obtain the distribution of the wells along the canal.

8. The system according to claim 6, characterized in that, The step of obtaining the rate of change of groundwater level along the canal based on the water level of the survey well includes, The area along the channel is divided into multiple measurement areas, wherein each measurement area includes multiple survey wells; Obtain the water level of multiple survey wells within each of the measurement areas; The standard deviation of the water level of multiple wells within the measurement area is calculated as the water level change rate of the measurement area.

9. The system according to claim 8, characterized in that, The step of setting up replenishment wells based on the groundwater level change rate along the canal includes, Obtain the set number of the replenishment wells; The replenishment wells are distributed among multiple measurement areas according to the ratio of the water level change rates between the measurement areas, so as to obtain the number of replenishment wells in each measurement area. Based on the number of supplementary wells in each of the measurement areas, the supplementary wells are evenly or randomly distributed in the measurement areas to obtain the location of the supplementary wells in each of the measurement areas. The locations of the supplementary wells in each of the measurement areas are summarized to obtain the distribution of the supplementary wells along the channel.

10. The system according to claim 8, characterized in that, The step of setting up replenishment wells based on the groundwater level change rate along the canal also includes, The monitoring well also includes a secondary replenishment well; The measured area whose straight-line distance from the channel is equal to the survey distance is marked as the edge measured area; Determine whether the rate of change of water level in the edge measurement area is less than a set value; If so, no action will be taken; If not, a new edge measurement area is set outside the existing edge measurement area, and the supplementary well is set in the newly set edge measurement area.

Images