CN106646661A

CN106646661A - Comprehensive survey system for hydrogeology of mineral deposit

Info

Publication number: CN106646661A
Application number: CN201610843064.7A
Authority: CN
Inventors: 张春志; 王光栋; 索立涛; 江永
Original assignee: Individual
Current assignee: Individual
Priority date: 2016-09-22
Filing date: 2016-09-22
Publication date: 2017-05-10

Abstract

The present invention discloses a comprehensive survey system for the hydrogeology of a mineral deposit. The system comprises a computer, a detector system, a probe system and a sleeve. The detector system is connected with the computer through a standard video cable. The probe system is connected with the detector system. The probe system is mounted on the sleeve. According to the technical scheme of the invention, by means of the comprehensive survey system for the hydrogeology of the mineral deposit, hydrological data collected in a hydrological monitoring station can be effectively integrated and fused with the data of water information monitoring videos. The automatically collected hydrological data and the water information monitoring videos are displayed in real time on a dedicated display for monitoring videos. Therefore, the intuitiveness of the hydrology and water monitoring process is improved. Meanwhile, data acquired through the measurement of geological loose circles and geophysical mechanical properties are analyzed, so that novel mining parameters in the current mining geological condition are conveniently formulated. According to the technical scheme of the invention, the content neutrality and the real-time performance of image blur degree evaluation are realized. Meanwhile, the ambiguity between any images can be quickly and accurately evaluated and compared.

Description

A kind of mineral deposit hydrogeology Comprehensive Exploration system

Technical field

The invention belongs to geological syntheses reconnoitre apparatus field, more particularly to a kind of mineral deposit hydrogeology Comprehensive Exploration system.

Background technology

Hydrogeology, geology subdiscipline refers to the various change of underground water and the phenomenon of motion in nature.Hydrology ground Matter is the science of Study of The Underground water.It is mainly the distribution of Study of The Underground water and formed rule, the physical property of underground water and Chemical composition, rationally groundwater resources and its utilization, adverse effect and its preventing and treating of the underground water to engineering construction and mining Deng.With the reach of science and the needs of production and construction, hydrogeology is divided into regional hydrogeology, groundwater dynamics again The subdisciplines such as, hydrogeochemistry, water supply hydrogeology, mineral deposit hydrogeology, soil improvement hydrogeology.Closely The research of the aspects such as Nian Lai, hydrogeology and underground heat, earthquake, environmental geology is interpenetrated, and some frontiers are defined again.

Underground water (groundwater)：Preservation and the water migrated in the ground space of underground.Aqueous ground is divided into two Band, top is aeration zone, i.e. unsaturation band, here, than water, also gas；Bottom is zone of saturation, i.e. saturated zone, is satisfied Water is full of water with the space in ground.The underground water of narrow sense refers to the water in zone of saturation.

Underground water：1. it is widely distributed, it is easy in site recovery to use；2. it is clean, be difficult to be contaminated, water quality is generally more excellent；③ It is not take up spatial surface；4. Dynamic comparison is stablized；5. output be affected by climate change it is less, with larger to Regulation capacity Deng.：1. irrational irrigation can cause Secondary Saline；2. excessive exploitation, can cause：In coastal area, seawater invasion, water quality Deteriorate；Surface subsidence, loses building in area stable；Hydraulic connection is induced between different water cut layer, the mixing for producing water is made With making water quality deterioration；Karst region surface collapse；3. other, such as mine pit water discharge, basis and side slope stable problem.

China develops the present situation of groundwater resources：1. the important water source of northern many city domestic waters；2. it is northern Arid, semiarid zone (17 provinces and cities) industrial and agricultural production, unique water source of life；3. southern some areas also begin to utilize underground Water and demand is increasing；4. the construction in big industrial base first has to solve water resources problem.

Develop the future of groundwater resources：1. the Sustainable Exploitation of groundwater resources is realized；2. underground water money is strengthened The scientific management in source；3. the environmental protection relevant with groundwater resources development is strengthened.(the three big problems that the world today faces：People Mouth, resource, environment)

Since the 1950's, social production scale unprecedentedly expands, and science and technology enters new developing period, and just There is new technological revolution, the ability of mankind's nature remodeling strengthens rapidly, and people keeps away via going after profit or gain for ancient times with the relation of water Evil, and modern age reduced levels bring good to and remove all evil, and have developed into the new stage brought good to and remove all evil of modern higher level.This new stage Hydrological science is given with new power and new characteristic.

Firstly, since prominent demand of the mankind to water resource, the research field of hydrological science towards for water resource most The direction of excellent exploitation is developed, to for objective evaluation, reasonable development, make full use of with water conservation provide science according to According to.

Secondly, large-scale mankind's activity is producing many impacts to natural water on natural environment.Grind Study carefully and evaluate the hydrologic effect of mankind's activity and the Significance for Environment of this effect, disclose the rule of hydrology phenomenon under the effect of human activity Rule, and then the new method and new way of hydrological analysis are inquired into, impact of the mankind's activity to hydrologic cycle is prevented towards being unfavorable for people The direction of class living environment is developed, and everything is becoming the new problem that hydrological science faces.

In addition, modern science and technology make the means of acquisition hydrographic information and the method for analysis hydrographic information have considerable to enter Step.For example, the application of remote sensing technology, is possibly realized observation macroscopical hydrology phenomenon interior on a large scale simultaneously；The application of nuclear technology Enable people to obtain microcosmic hydrographic information；Hydrological simulation method, hydrology Stochastic Analysis Method, Hydrologic System, make The ability development of people's research hydrology phenomenon is to new level；The especially application of electronic computer, makes hydrological science from the hydrology The research of basic law is observed, by manpower and mechanically actuated, the automation with electronic computer as core is developed into.

Frontier science between hydrological science and other science constantly rise, and interdisciplinary space is gradually filled out Mend.Meanwhile, people start to see, water becomes the key factor for affecting social development.Water is showing the same of its natural quality When, its social property is also increasingly showed, and is gradually recognized by people.Therefore, hydrological science would be possible to develop into It is a Comprehensive Science with natural science and social science double properties.

Generally, hydrology is from terms of the ambit that it is subordinate to, as a branch of geophysics science, mainly The presence of water, distribution in research earth system, move and circulation change rule, the physics of water, chemical property, and hydrosphere with it is big The correlation in air ring, lithosphere and biosphere；As the important component part of water conservancy subject, the shape of main research water resource Into, spatial and temporal distributions, development and protection, the formation of bloods and droughts, prediction and preventing and treating, and hydraulic engineering and other works Hydrology and water conservancy computing technique in planning, design, construction, the management of Cheng Jianshe.

Although having started in 19th century use the word of hydrogeology one, to 20 beginning of the century scientist Mead this art is just provided The extensive implication of language one：Hydrogeology is the generation and motion for studying the following water of earth's surface.Phase late 1950s is to 80 In this nearly 30 years time of early stage in age, hydrogeology is quickly ripe, becomes the plentiful a member of geoscience assistant. Before nineteen sixty, hydrogeology is mainly the field of geologist, as a natural scientist, for control subsurface flow Dynamic factor and rule, it is with no interest or know little about it, no matter difierence equation goes to be been described by.On the other hand, engineer When the specific capacity and total water yield of well is estimated, merely must calculate, the ash between rock stratum " permeable " and " waterproof " It is at a loss as to what to do among domain.

The essential characteristic of Modern hydrogeology mainly has：1. combine closely with the new theory new disciplines of modern science, than Such as systematology, information theory, control is touched upon corresponding systematic science, environmental science, the information science for producing, to hydrogeology Development generates significant impact；2. the combination of contemporary applied mathematics and hydrogeology, particularly method for numerical simulation obtain general All over application, scale-model investigation becomes the main contents of water resource research, makes hydrogeology develop into quantitative study from qualitative research New stage；3. from the research of ground water regime and natural environment system correlation, expand to and social economic system relation Research.Research to groundwater resources, also develops into the research of administrative model and economic model from Mathematical Modeling；4. it is many new Subdiscipline generation with development, such as regional hydrogeology, karst hydrogeology, Remote Sensing To Hydrology geology, environment Hydrogeology, medical environment geochemistry, pollution hydrogeology and mathematics hydrogeology, water resource hydrogeology Learn；5. new technology, the application of new method, in addition to computer technology, remote sensing technology, isotope technology, automatic monitoring technical are indoor Analogue technique, and high accuracy water quality analytical technology etc., all it is widely used, promote the development of hydrogeology.

As hydrogeological the reach of science, its research contents are more and more extensive, main research can be summarized as six Individual aspect：

(1) the formation of underground water and conversion：Illustrate groundwater origin with the ABC for being formed (including the preservation of underground water Condition), and inquire into atmospheric water, surface water, the soil water and underground water mutually convert, alternate basic law.

(2) the type and feature of underground water：Illustrate the condition of storage and its fundamental type of underground water, including the master of underground water Want physicochemical property.

In zone of saturation and aeration zone moisture and solute motion：The fundamental differential of main Study of The Underground current, bag The flowing of underground water Xiang Jing, canal is included, to disclose the spatial and temporal variation of level of ground water and the water yield.Vadose water is inquired into simultaneously with ground The fundamental equation of lower water solute transfer.

(4) groundwater dynamic and hydrologic(al) budget：Groundwater dynamic under the influence of different natural factors and human factor is discussed Groundwater equilibrium equation under Changing Pattern, and different condition.

(5) Underground water resources calculation and evaluation：Local mining area and regional large area mining area underground water money are discussed respectively The main method that source is evaluated, and specifically introduce the calculating side of relevant aquifer parameter measure and the increment of groundwater and excretion Method.Meanwhile, illustrate the relevant knowledge of Groundwater Quality Evaluation.

(6) water resources system management：ABC in terms of elaboration Groundwater Resource Management and protection, focuses on Water resources system administrative model and its application.

Using technological means：(1) investigation, probing, geophysical exploration and remote sensing technology；(2) it is various to observe and test skill Art (the observation of water level, flow etc.；Bailing test, tracer test and Dispersion Test etc.)；(3) various groundwater simulation technology (numerical value Simulation it is more)；(4) isotope technology etc..

With the continuous improvement of scientific and technological level, hydrogeological computational methods also constantly develop.Hydrogeological calculating side Method substantially has：Analytic solution, physical Modeling Method, numerical solution, systematic analytic method, probabilistic method etc..

Analytic solution

Before the sixties, the head and problems of liquid flow of water-bearing layer underground water is solved, bias toward analytic solution more.Such as " underground water Described in dynamics " course, the either Qiu cotton clothes' formula based on stationary flow, or the Tai Si based on unsteady flow Formula, their derivation has many it is assumed that when hydrogeologic condition meets these hypothesis, having no problem certainly.But solve When certainly large-scale ground water regime is calculated, due to the complexity of hydrogeologic condition, analytic solution is just helpless.

Physical Modeling Method

Physical analogy has electrical analogue, hydraulic analogy, viscous flow field simulation, membrane analogy etc., more with electrical analogue application.It is early In the twenties in this century, Pavlov's Si base of the Soviet Union proposes electrolyte simulation (arn A), and it becomes study at that time water conservancy project Ooze the important means of the problem of catching in building area.The simulation of resistance net is developed into back vent, in the fifties and the sixties, R-C networks Also it is developed with R-R Ah networks simulation.The mid-1960s fork occurs in that the mixer being combined with computer structure.

Numerical solution

The later stage sixties is applied to hydrogeology calculating with the development of electronic computer, people numerical simulation. Because electrical analogue making and parameter testing are all bothered than numerical method, so using being more numerical solution.

5 classes can be substantially summarized as in the numerical method of hydrogeological application in calculation.Finite difference calculus (abbreviation finite difference Method)；Finite Element (abbreviation FInite Element)；Boundary element method (abbreviation boundary element method)；The method of characteristic curves }；Finite Analytic Method.

Finite difference calculus begins to be applied to hydrogeological calculating from the beginning of the sixties.Initially multiplex regular net is solved with lax Method, nineteen sixty-eight introduces alternating direction Henan formula calculus of finite differences, introduces SIP technique again later, is generalized to declension away from situation within 1973, Lan Mate (Lemard) proposes that upstream weighted is limited to block point-score in 1D79.

Finite Element starts to be applied to water history geological calculation from nineteen sixty-eight, bends the parameter Finite Elements such as 1 eight within 1972, Xiu Yankang (Huyakorn) in 1977 and Neil health card (lxlilkuka) etc. propose windward Finite Element.

Finite difference calculus and Finite Element are numerical computation methods the most frequently used in hydrogeological juice grate.

Boundary element method is a kind of new numerical method that the mid-1970s grow up.

Finite Analytic Method is a kind of new numerical computation method for growing up the eighties.It is also a kind of discrete region side Method, it is to carry out discretization by certain parsing approach, obtains prescription journey, then tries to achieve the head approximation of each node and enters One step calculates flow.

Other methods

Systematic analytic method, is a kind of method that combined mathematical module and computer technology are analyzed, in underground water money Developed rapidly in source control.Many countries, call, to reach ground extensive and large-scale river is carried out with the method Lower water is mutually adjusted with river resource watt, unified operation.Systems approach is with natural according to the meteorology of location, geology, landforms etc. The social environment conditions such as the relation and economy, politics of geographical conditions and system, as needed with may, be the system determine- Individual optimal solution.

Stochastic model is also extensively applied in Groundwater Resource Management.Such as time series analysis, also begin to be applied to underground During water is calculated.With the development of computer science, the method for making more renewals is applied in actual production.

Limited by rock mechanics development level, current geologic prospect design there is still a need for by " engineering analogy " with " field monitoring " combines, and could obtain the effect of satisfaction.

As exploitation deep continues to increase, current mining depth has reached nearly 600m, ore bed geological conditions there occurs compared with Whether big change, original parameter does not also know suitable for existing ore deposit geological conditions, is unfavorable for realizing that mine safety is efficient Production.The measure of geology relaxation zone has very important significance to geology physico-mechanical properties.

At present, in the hydrology water regime monitoring scene of China, the technology such as automatic data collection and regimen video monitoring of hydrographic data Means application is wide, but, there is a basic problem therebetween：To the video monitoring data of regimen often through special Method for processing video frequency realize, then play out on special monitor video display；The hydrographic data of automatic data collection is then In depositing in the database in computer server, need to be inquired about and processed by special application software.Two kinds of data Inquiry be separated from each other with display, be unfavorable for carrying out regimen video monitoring and hydrographic data in real time, intuitively monitor.It is existing Multi-antenna, measurement parameter is relatively simple, and hydrologic parameter tends not to have concurrently with water quality parameter, between parameter directly Correlation is difficult to comprehensive analysis.

The content of the invention

It is an object of the invention to provide a kind of mineral deposit hydrogeology Comprehensive Exploration system, it is intended to solve it is existing it is long-range from Dynamic monitoring system, measurement parameter is relatively simple, and hydrologic parameter tends not to have concurrently with water quality parameter, and whether original parameter is suitable for Also do not know in current ore deposit geological conditions, be unfavorable for realizing that mine safety is efficiently produced, the measure of current geology relaxation zone Inaccurate problem is measured to geology physico-mechanical properties.

The present invention is achieved in that a kind of mineral deposit hydrogeology Comprehensive Exploration system, and the mineral deposit hydrogeology is comprehensively surveyed System is looked into including computer, detector system, probe system, sleeve pipe；The detector system is by standard video cable connection electricity Brain；The probe system connects detector system；The probe system is arranged on sleeve pipe；

The probe system includes ultrasonic detection device, hydrology sensor；Ultrasonic detection device is sent out including ultrasonic wave Penetrate transducer, ultrasonic wave receive transducer, fluting plastic packaging pipe；The ultrasound transmitting transducer, ultrasonic wave receive transducer are equal In fluting plastic packaging pipe；The hydrology sensor is arranged on sleeve pipe；The fluting plastic packaging pipe is arranged on sleeve pipe side；

Radiating circuit is installed on the ultrasound transmitting transducer；Reception electricity is installed in ultrasonic wave receive transducer Road, control circuit, timing circuit, lithium battery；The receiving circuit, control circuit, timing circuit, lithium battery by wire according to Secondary connection；The control circuit is connected with radiating circuit signal, and radiating circuit is electrically connected with lithium battery by power line；

Described sleeve pipe is provided with more piece, and often section is threaded connection；

The detector system includes that waterlevel data monitoring node, video monitoring node, monitor video display, signal connect Receive module, central processing unit；The monitor video display, signal receiving module connect central processing unit by printed line；Institute State waterlevel data monitoring node, video monitoring node to be connected with central processing unit by signal receiving module；Hydrology sensor It is connected with waterlevel data monitoring node；Waterlevel data detection node and hydrology sensor are connected with video monitoring node.

Further, the waterlevel data detection node constitutes detection node network using distributed water level detection gauge.

Further, the Data Detection node constitutes detection node network using distributed rainfall data detector.

Further, the waterlevel data detection node and hydrology sensor adopt wireless transmission method with video monitoring node Connection；The hydrology sensor is provided with multiple.

Further, it is using the method for ultrasonic detection device test geology wall rock loosening ring：

Wall rock loosening ring method of testing adopts sonic method, sonic method test wall rock loosening ring to propagate in country rock based on sound wave The change of speed, will occur geonetrical attenuation and Physical Attenuation when elastic wave is propagated in rock mass, the different mechanical property in rock mass Elastic wave will be reflected, scattered and thermal losses on the structural plane of matter so that constantly decay causes velocity of wave to reduce to elastic wave energy, By the wave theory of elastic wave, the wave equation in unlimited isotropy medium is：

V in formula_p--- compressional wave；

V_s--- shear wave

Boundary condition and primary condition based on isotropy elastic half-space, obtains corresponding velocity of wave, the elasticity with rock Relational expression between modulus E, Poisson's Ratio σ, density p is：

It is distributed in country rock according to velocity of wave, draws relaxation zone scope.

Further, the acoustic transit time meter of ultrasonic detection device length country rock during velocity of longitudinal wave is by determining drilling Calculate, calculated using single hole method of testing；

Single hole method of testing is：Transmitting transducer F launches in the borehole ultrasonic wave, and along the wall of a borehole propagation is slided, and transmitting is changed Triggering timing circuit timing while energy device F transmitting ultrasonic waves, after receive transducer J receives ultrasound information timing is stopped, Measure propagation time t of the sound wave in F-J length rock mass：

T=Δ t+t₀+Δt

The transmitting transducer that t-instrument shows in formula is to the propagation time between receive transducer；

Δ t-sound wave is in the wall of a borehole and transducer interstitial propagation time；

t₀- special ripple is in transmitting-receive transducer length range along the hole wall propagation time；

φ_D- drilling internal diameter；

φ_d- transducer diameter；

The SVEL of water is coupled in v-drilling.

Spread speed of the sound wave in rock-boring be：

V in formula --- SVEL in drilling；

L --- distance between transducer F and J；

During ultrasonic test, water coincidence Acoustic Wave Propagation is full of in drilling, and drilling is pricked into downwards 3 °～5 °；

When wall-rock crack is more, velocity of wave of the velocity of wave relative to depth completely without crack rock is low, is measured by rock-boring Changes in distribution curve of the sound wave velocity of longitudinal wave in wall rock drill-hole or using time-hole depth curve, judges wall-rock crack scope, Transmitting transducer F launches in the borehole ultrasonic wave, and along the wall of a borehole propagation is slided, and touches while transmitting transducer F transmitting ultrasonic waves Send out timing road and start timing, timing s is stopped after receive transducer I receives ultrasound information, measure sound wave in the propagation of F-I Between, calculate velocity of wave.

Further, the hydrology sensor includes level switch, current amount detector, geological image collector；It is described Level switch, current amount detector are connected with waterlevel data detection node；Geological image collector and video monitoring node It is connected；The monitor video display includes display screen, ambiguity evaluation module, fuzziness adjusting module；The fuzziness is adjusted Mould preparation block is connected with display screen by printed line；The geological image collector is used to obtain the image of geological stratification；

The ambiguity evaluation module is used to obtain the geology tomographic image of geological image collector transmission, and calculates before filtering Image statistics ratio afterwards；

The fuzziness adjusting module is connected with ambiguity evaluation module, and for adjusting original image fuzziness final figure is drawn Picture and image blur evaluation index；

It is to image blur evaluation method using ambiguity evaluation module, fuzziness adjusting module：

Step one, image is obtained, and by geological image collector geology tomographic image to be evaluated is obtained；

Step 2, image gray processing, for convenience of the edge extracting of image, using the R of RGB image in Digital Image Processing, Coloured image is converted into gray level image by the pixel value of G, B each passage and the transformational relation of gray level image pixel value, and formula is such as Under：

Gray=R*0.3+G*0.59+B*0.11；

Step 3, Edge extraction is made using the Roberts operator edge detections technology in digital image processing method For the edge that gray level image obtains image, different detective operators have different edge detection templates, according to concrete template The difference for intersecting pixel is calculated as current pixel value, it is as follows using template：

E (i, j)=| F (i, j)-F (i+1, j+1) |+| F (i+1, j)-F (i, j+1) |；

Step 4, image procossing is filtered process to gray level image to be evaluated to construct using high pass/low pass filter The reference picture of image, using 3*3 mean filters, using Filtering Template traversing graph as each pixel, every time by template center Current pixel is placed in, the mean value of all pixels is newly worth as current pixel using in template, and template is as follows：

Step 5, image border statistical information is calculated, and respective edge half-tone information before and after image filtering, filtering are calculated respectively The image F statistical informations to be evaluated of before processing are sum_orig, and the reference picture F2 statistical informations after filtering process are sum_ Filter, specific formula for calculation is as follows：

Sum_orig=w1 × (| F (i, j)-F (i-1, j) |+| F (i, j)-F (i, j-1) |+| F (i, j)-F (i, j+1) |+ |F(i,j)-F(i+1,j)|)+w2×(|F(i,j)-F(i-1,j-1)|+|F(i,j)-F(i-1,j+1)|+|F(i,j)-F(i+ 1, j-1) |+| F (i, j)-F (i+1, j+1) |),

Sum_filter=w1 × (| F2 (i, j)-F2 (i-1, j) |+| F2 (i, j)-F2 (i, j-1) |+| F2 (i, j)-F2 (i, j+1) |+| F2 (i, j)-F2 (i+1, j) |)+w2 × (| F2 (i, j)-F2 (i-1, j-1) |+| F2 (i, j)-F2 (i-1, j+1) |+| F2 (i, j)-F2 (i+1, j-1) |+| F2 (i, j)-F2 (i+1, j+1) |),

Wherein, w1 and w2 is according to the weights set with a distance from center pixel, w1=1, w2=1/3；

Step 6, image blur index is calculated, the image filtering front and rear edges grey-level statistics that step 5 is drawn Ratio as fuzziness index, for convenience of evaluating, take it is larger for denominator, it is less for molecule, keep the value between (0,1) Between；

Step 7, according to the DMOS scopes of the best visual effect draw a corresponding fuzziness indication range [min, Max], specially：

Fuzziness adjusting range is drawn, using the ambiguity evaluation method in above-mentioned steps 174 panel heights in LIVE2 are evaluated This blurred picture, calculates the ambiguity evaluation value of each of which, is then set up using fitting tool plot (value, DMOS) Mapping relations between evaluation of estimate value and DMOS, according to the corresponding DMOS scopes of the best visual effect corresponding one is drawn Fuzzy evaluation value scope [min, max]；

Step 8, image blur adjustment, if image blur index is less than min, according to step 6, judges image filtering Change very big in front and back, original image is excessively sharpened, then be filtered adjustment using low pass filter；If being more than max, the filter of process decision chart picture Vary less after wavefront, original image is excessively obscured, be then filtered adjustment using high-pass filter, to reach more preferably vision effect Really；

Step 9, draws final image and the image blur evaluation index, and shows on a display screen.

Further, filter process image is not that single mode processes view picture evaluation image, but in view of the border of image With center pixel because the difference of position causes difference of the wave filter to its processing mode, according to filter template size correspondingly Ignore recycling filter process image after boundary pixel, then for boundary pixel takes the method that original pixels are filled to carry out Process.

Further, the central processing unit is provided with synchronized orthogonal Frequency Hopping Signal blind source separating module, the orthogonal jump of the step The signal processing method of frequency signal blind source separating includes：

Step one, using the array antenna received containing M array element from multiple synchronized orthogonal frequency hopping sensors Frequency Hopping Signal, to per all the way receive signal sample, the M roads discrete time-domain mixed signal after being sampledThe interaction times of different time piece between collection array antenna node, according to obtaining Data setup time sequence, the interaction times of next timeslice are predicted between node by third index flatness, will hand over Direct trust value of the relative error of mutual number of times predicted value and actual value as node；The concrete calculation procedure of direct trust value For：The interaction times of n timeslice between collection network observations node i and node j：Intervals t are chosen as one Individual observation time piece, it is true to hand over using the interaction times of observer nodes i and tested node j in 1 timeslice as observation index Mutually number of times, is denoted as y_t, the y of n timeslice is recorded successively_n, and save it in the communications records table of node i；Prediction (n+1)th The interaction times of individual timeslice：According to the interaction times setup time sequence of the n timeslice for collecting, put down using three indexes Interaction times between sliding method prediction next one timeslice n+1 interior nodes i and j, predict interaction times, are denoted asCalculate public Formula is as follows：

Predictive coefficient a_n、b_n、c_nValue can be calculated by equation below：

Wherein：Be respectively once, secondary, Three-exponential Smoothing number, calculated by equation below Arrive：

It is the initial value of third index flatness, its value is

α is smoothing factor (0 ＜ α ＜ 1), embodies the time attenuation characteristic trusted, i.e., from predicted value more close to timeslice y_tWeight is bigger, from predicted value more away from timeslice y_tWeight is less；If data fluctuations are larger, and long-term trend change width Degree is larger, and α when substantially rapidly rising or falling trend is presented should take higher value

(0.6～0.8), increases impact of the Recent data to predicting the outcome；When data have fluctuation, but long-term trend change not When big, α can between 0.1～0.4 value；If data fluctuations are steady, α should take smaller value (0.05～0.20)；

Calculate direct trust value：

Direct trust value TD of node j_ijTo predict interaction timesWith true interaction times y_n+ ¹Relative error,

Indirect trust values are calculated using calculating formula obtained from multipath trust recommendation mode；Trusted node is collected to section The direct trust value of point j：Node i to all credible associated nodes for meeting TDik≤φ inquire its direct letter to node j Appoint value, wherein φ for recommended node believability threshold, according to the precision prescribed of confidence level, the span of φ is 0～ 0.4；Calculate indirect trust values：Trust value collected by COMPREHENSIVE CALCULATING, obtains indirect trust values TR of node j_ij,Wherein, Set (i) is interacted and its is direct to have with j nodes in the associated nodes of observer nodes i Trust value meets TD_ikThe node set of≤φ；

Comprehensive trust value is drawn by direct trust value and indirect trust values conformity calculation；Comprehensive trust value (T_ij) calculating it is public Formula is as follows：T_ij=β TD_ij+(1-β)TR_ij, wherein β (0≤β≤1) represent direct trust value weight, when β=0, node i and The calculating that node j does not have direct interaction relation, comprehensive trust value arises directly from indirect trust values, and it is more objective to judge；When β=1 When, node i, all from direct trust value, in this case, judges more subjective to the comprehensive trust value of node j, real Border calculates the value for determining β as needed；

Step 2, carries out overlapping adding window Short Time Fourier Transform to M roads discrete time-domain mixed signal, obtains M mixing letter Number time-frequency domain matrix

P=0,1 ..., P-1, q=0,1 ..., N_fft- 1, wherein P represent total window number, N_fftRepresent FFT length；P, Q) time-frequency index is represented, specific time-frequency value isHere N_fftThe length of FFT is represented, p is represented and added Window number of times, T_sRepresent sampling interval, f_sSample frequency is represented, C is integer, represent the sampling at Short Time Fourier Transform adding window interval Points, C ＜ N_fft, and K_c=N_fft/ C is integer, that is to say, that use the Short Time Fourier Transform for overlapping adding window；

Step 3, to the frequency-hopping mixing signal time-frequency domain square obtained in step 2Enter Row pretreatment；

Step 4, estimates the jumping moment of each jump and respectively jumps corresponding normalized hybrid matrix using clustering algorithm Column vector, Hopping frequencies；It is right at p (p=0,1,2 ... the P-1) momentThe frequency values of expression are clustered, in the cluster for obtaining Heart numberThe carrier frequency number that the expression p moment is present,Individual cluster centre then represents the size of carrier frequency, uses respectivelyRepresent；To each sampling instant p (p=0,1,2 ... P-1), using clustering algorithm pairEnter Row cluster, it is same availableIndividual cluster centre, usesRepresent；To allAverage and round, obtain To the estimation of source signal numberI.e.：

Find outMoment, use p_hRepresent, the p to each section of continuous value_hIntermediate value is sought, is usedTable Show that l sections are connected p_hIntermediate value, thenRepresent the estimation at l-th frequency hopping moment；Obtained according to estimationAnd the 4th estimate that the frequency hopping moment for obtaining estimates in step and each jump correspondingIt is individual Hybrid matrix column vectorSpecifically formula is：

HereRepresent that l is jumped correspondingIndividual mixing Matrix column vector estimate；Estimate the corresponding carrier frequency of each jump, useRepresent that l jumps correspondence 'sIndividual frequency estimation, computing formula is as follows：

Step 5, estimates that the normalization hybrid matrix column vector for obtaining estimates time-frequency domain frequency hopping source signal according to step 4；

Step 6, splices to the time-frequency domain frequency hopping source signal between different frequency hopping points；

Step 7, according to source signal time-frequency domain estimate, recovers time domain frequency hopping source signal；To each sampling instant p (p= 0,1,2 ...) frequency domain data Y_n(p, q), q=0,1,2 ..., N_fft- 1 is N_fftThe IFFT conversion of point, obtains p sampling instants pair The time domain frequency hopping source signal answered, uses y_n(p,q_t)(q_t=0,1,2 ..., N_fft- 1) represent；The time domain that above-mentioned all moment are obtained Frequency hopping source signal y_n(p,q_t) process is merged, final time domain frequency hopping synthesizer Signal estimation is obtained, concrete formula is as follows：

Here K_c=N_fftThe sampling number that/C, C are spaced for Short Time Fourier Transform adding window, N_fftFor the length of FFT.

The hydrology that the mineral deposit hydrogeology Comprehensive Exploration system that the present invention is provided can will be gathered effectively in hydrologic monitoring website Data carry out data integration fusion with regimen monitor video, and the hydrographic data of automatic data collection and regimen monitor video can be shown in real time Show and be exclusively used on the display of monitor video, improve the sex chromosome mosaicism directly perceived of hydrology regimen monitoring；

The present invention to geology relaxation zone according to determining and geology physico-mechanical properties determines obtained data and carried out point Analysis, is conducive to and has formulated operational parameter new under current coal mining geological condition；

The picture appraisal of the present invention is different from traditional evaluation method, and the present invention sets up special in image self structure to be evaluated On the basis of point, from the angle of relative evaluation, the reference picture of image to be evaluated is constructed using wave filter, before and after calculating change The ratio of image border statistical information is used as evaluation index；The principle of the present invention is simple, realizes the interior of image blur evaluation Hold independence and real-time, fuzziness that can quick and precisely between any image of evaluation comparison.

The present invention under conditions of any channel information is not known, believe by the mixing according only to the multiple Frequency Hopping Signals for receiving Number, frequency hopping source signal is estimated, multiple Frequency Hopping Signals can be carried out under conditions of reception antenna number is less than source signal number Blind estimate, with only Short Time Fourier Transform, and amount of calculation is little, easily realize, the method is carrying out blind point to Frequency Hopping Signal From while, moreover it is possible to partial parameters are estimated, it is practical, with stronger popularization and using value.

Description of the drawings

Fig. 1 is mineral deposit hydrogeology Comprehensive Exploration system schematic provided in an embodiment of the present invention；

In figure：1st, computer；2nd, detector system；3rd, probe system；4th, sleeve pipe.

Fig. 2 be it is provided in an embodiment of the present invention be single hole test job principle schematic.

Fig. 3 is ultrasonic detection device schematic diagram provided in an embodiment of the present invention.

Fig. 4 is image blur evaluation method general flow schematic diagram provided in an embodiment of the present invention.

Specific embodiment

In order that the objects, technical solutions and advantages of the present invention become more apparent, with reference to embodiments, to the present invention It is further elaborated.It should be appreciated that specific embodiment described herein is not used to only to explain the present invention Limit the present invention.

Below in conjunction with the accompanying drawings and specific embodiment to the present invention application principle be further described.

As shown in figure 1, mineral deposit hydrogeology Comprehensive Exploration system provided in an embodiment of the present invention, the mineral deposit hydrogeology is comprehensive Closing investigation system includes computer 1, detector system 2, probe system 3, sleeve pipe 4；The detector system 2 is by normal video electricity Cable connects computer 1；The probe system 3 connects detector system 2；The probe system 3 is arranged on sleeve pipe 4；

V in formula_p--- compressional wave；

V_s--- shear wave

T=Δ t+t₀+Δt

φ_D- drilling internal diameter；

φ_d- transducer diameter；

The SVEL of water is coupled in v-drilling.

Spread speed of the sound wave in rock-boring be：

V in formula --- SVEL in drilling；

L --- distance between transducer F and J；

Gray=R*0.3+G*0.59+B*0.11；

E (i, j)=| F (i, j)-F (i+1, j+1) |+| F (i+1, j)-F (i, j+1) |；

Sum_filter=w1 × (| F2 (i, j)-F2 (i-1, j) |+| F2 (i, j)-F2 (i, j-1) |+| F2 (i, j)-F2 (i,j+1)|+|F2(i,j)-F2(i+1,j)|)+w2×(|F2(i,j)-F2(i-1,j-1)|+|F2(i,j)-F2(i-1,j+1) |+| F2 (i, j)-F2 (i+1, j-1) |+| F2 (i, j)-F2 (i+1, j+1) |),

Step one, using the array antenna received containing M array element from multiple synchronized orthogonal frequency hopping sensors Frequency Hopping Signal, to per all the way receive signal sample, the M roads discrete time-domain mixed signal after being sampledThe interaction times of different time piece between collection array antenna node, according to obtaining Data setup time sequence, the interaction times of next timeslice are predicted between node by third index flatness, will hand over Direct trust value of the relative error of mutual number of times predicted value and actual value as node；The concrete calculation procedure of direct trust value For：The interaction times of n timeslice between collection network observations node i and node j：Intervals t are chosen as one Individual observation time piece, it is true to hand over using the interaction times of observer nodes i and tested node j in 1 timeslice as observation index Mutually number of times, is denoted as y_t, the y of n timeslice is recorded successively_n, and save it in the communications records table of node i；Prediction (n+1)th The interaction times of individual timeslice：According to the interaction times setup time sequence of the n timeslice for collecting, put down using three indexes Interaction times between sliding method prediction next one timeslice n+1 interior nodes i and j, predict interaction times, are denoted asComputing formula It is as follows：

It is the initial value of third index flatness, its value is

α is smoothing factor (0 ＜ α ＜ 1), embodies the time attenuation characteristic trusted, i.e., from predicted value more close to timeslice y_tWeight is bigger, from predicted value more away from timeslice y_tWeight is less；If data fluctuations are larger, and long-term trend change width Degree is larger, and α when substantially rapidly rising or falling trend is presented should take higher value (0.6～0.8), increase Recent data to prediction As a result impact；When data have a fluctuation, but long-term trend change it is little when, α can between 0.1～0.4 value；If data wave Dynamic steady, α should take smaller value (0.05～0.20)；

Calculate direct trust value：

Indirect trust values are calculated using calculating formula obtained from multipath trust recommendation mode；Collect trusted node pair The direct trust value of node j：Node i meets TD to all_ikThe credible associated nodes of≤φ inquire its direct letter to node j Appoint value, wherein φ for recommended node believability threshold, according to the precision prescribed of confidence level, the span of φ is 0～ 0.4；Calculate indirect trust values：Trust value collected by COMPREHENSIVE CALCULATING, obtains indirect trust values TR of node j_ij,Wherein, Set (i) is interacted and its is direct to have with j nodes in the associated nodes of observer nodes i Trust value meets TD_ikThe node set of≤φ；

The detection of loosening area of surrounding rocks, is application technology that elastic wave testing technology is surveyed as boundary condition.Due to vertical Ripple has the characteristics such as spread speed is fast, energy is high, live easy realization, and merely with compressional wave the detection of loosening area of surrounding rocks is engaged in. In general, acoustic wave propagation velocity is relevant with following factors：

1) impact of the crack to elastic wave in country rock.When the parallel crack in elastic wave propagation direction, velocity of wave is unchanged.And work as When elastic wave propagation direction is vertical with crack, velocity of wave then will be reduced, and the number of reduction and the width in crack, the size in crack, be filled Fill out physical property matter and lithology is relevant.

2) impact of the mechanical property suffered by rock mass to elastic wave.Different elastic fluids, acoustic wave propagation velocity is also different.

3) impact of the stress suffered by rock mass to SVEL.With stress increase, increase in all directions velocity of wave, parallel to plus Carrying direction velocity of wave increases maximum, and it is advanced the speed gradually less, and perpendicular to loading direction, velocity of wave increases little.

So, it is distributed in country rock according to velocity of wave, so that it may draw relaxation zone scope.

Velocity of longitudinal wave is calculated by the acoustic transit time of certain distance (probe length) country rock in measure drilling, There are " diplopore is to surveying " and " single hole test " two methods, diplopore needs a pair of parallel drilling to surveying, wherein transmitting is installed in a hole passing Sensor, another hole is installed in respective depth and receives sensor, and what it was reflected is Radial crack feature, and diplopore is flat to drilling to surveying Row degree requires higher, operation inconvenience, and at present application is gradually few.

Single hole test reflection is annular crack feature, and Fig. 2 is single hole test job principle schematic.Transmitting transducer F Launch ultrasonic wave in the borehole, along the wall of a borehole propagation is slided.

Triggering timing circuit timing while transmitting transducer F transmitting ultrasonic waves, when receive transducer J receives ultrasonic wave letter Stop timing after breath, measure propagation time t of the sound wave in F-J length rock mass：

T=Δ t+t₀+Δt

φ_D- drilling internal diameter；

φ_d- transducer diameter；

The SVEL of water is coupled in v-drilling.

Spread speed of the sound wave in rock-boring be：

V in formula --- SVEL in drilling；

L --- distance between transducer F and J.

As shown in Figure 3：Ultrasonic detection device is by transmitting transducer and receive transducer via fluting plastic packaging pipe connection group Into two transducers of transmitting-receiving are used interchangeably.

Sleeve pipe is made up of copper pipe, has a yardstick groove per 10cm, totally 20 section, thread connection.Transducer couples plastic tube, is facing upward Tiltedly adhesive tape is needed to swathe during test, with preventing water leakage.

During ultrasonic test, need, full of water coincidence Acoustic Wave Propagation, generally testing bore holes should be laid in drilling Roadway's sides, and drilling is slightly pricked downwards 3 °～5 ° and tests more convenient so as to water, so, when drilling is faced upward tiltedly or upwards When, it is to ensure to fill water in drilling, hole packer need to be used.

Spread speed of the ultrasonic wave in rock mass is relevant with rock mass stress and crack degree, when wall-rock crack (rupture Meet) it is many when, velocity of wave of the velocity of wave relative to depth completely without crack (non-loose damage) rock mass is low.By rock-boring (diameter 40 ～45mm) measure changes in distribution curve or " time --- hole depth " curve of the sound wave velocity of longitudinal wave in wall rock drill-hole, you can sentence Determine wall-rock crack (loosening) scope.Velocity of longitudinal wave is by determining the Acoustic Wave Propagation for boring aerial certain distance (probe length) country rock Time Calculation is out.Transmitting transducer F launches in the borehole ultrasonic wave, and along the wall of a borehole propagation is slided.Transmitting transducer F launches Triggering timing road starts timing while ultrasonic wave, and timing s is stopped after receive transducer I receives ultrasound information, measures sound Thus ripple calculates velocity of wave in the propagation time of F-I.

Continuous moving ultrasonic detection device in the borehole, you can measure on whole drillable length " velocity of wave-hole depth " or " time-hole depth " curve.The maximum hole depth of velocity of wave or time change, the as size of wall rock loosening ring in curve.

Sonic method test relaxation zone test basic principle

Affect relaxation zone factor it is more, by non-stressed concentrate or non-stressed overlap-add region in, the loosening number of turns for stabilizing Value is defined as a reference value of relaxation zone, and what it reflected is the essential characteristic of rock strength and the stress of primary rock.Undergo to adopt when geology, When the factors such as tomography affect, relaxation zone numerical value is higher than a reference value, should supplement this measure.Relaxation zone numerical value to measuring is carried out Analysis, after excluding abnormal data, using larger numerical value as the ultimate depth of relaxation zone.

Presently preferred embodiments of the present invention is the foregoing is only, not to limit the present invention, all essences in the present invention Any modification, equivalent and improvement made within god and principle etc., should be included within the scope of the present invention.

Claims

1. a kind of mineral deposit hydrogeology Comprehensive Exploration system, it is characterised in that the mineral deposit hydrogeology Comprehensive Exploration system includes Computer, detector system, probe system, sleeve pipe；The detector system connects computer by standard video cable；The probe System connects detector system；The probe system is arranged on sleeve pipe；

The probe system includes ultrasonic detection device, hydrology sensor；Ultrasonic detection device includes that ultrasonic wave transmitting is changed Can device, ultrasonic wave receive transducer, fluting plastic packaging pipe；The ultrasound transmitting transducer, ultrasonic wave receive transducer are respectively mounted In fluting plastic packaging pipe；The hydrology sensor is arranged on sleeve pipe；The fluting plastic packaging pipe is arranged on sleeve pipe side；

Radiating circuit is installed on the ultrasound transmitting transducer；Receiving circuit, control are installed in ultrasonic wave receive transducer Circuit processed, timing circuit, lithium battery；The receiving circuit, control circuit, timing circuit, lithium battery are connected successively by wire Connect；The control circuit is connected with radiating circuit signal, and radiating circuit is electrically connected with lithium battery by power line；

The detector system includes that waterlevel data monitoring node, video monitoring node, monitor video display, signal receive mould Block, central processing unit；The monitor video display, signal receiving module connect central processing unit by printed line；The water Position data monitoring node, video monitoring node are connected by signal receiving module with central processing unit；Hydrology sensor and water Position data monitoring node is connected；Waterlevel data detection node and hydrology sensor are connected with video monitoring node.

2. mineral deposit hydrogeology Comprehensive Exploration system as claimed in claim 1, it is characterised in that the waterlevel data detection section Point constitutes detection node network using distributed water level detection gauge.

3. mineral deposit hydrogeology Comprehensive Exploration system as claimed in claim 1, it is characterised in that the Data Detection node is adopted Detection node network is constituted with distributed rainfall data detector.

4. mineral deposit hydrogeology Comprehensive Exploration system as claimed in claim 1, it is characterised in that the waterlevel data detection section Point and hydrology sensor are connected with video monitoring node using wireless transmission method；The hydrology sensor is provided with multiple.

5. mineral deposit hydrogeology Comprehensive Exploration system as claimed in claim 1, it is characterised in that utilize ultrasonic detection device Test geology wall rock loosening ring method be：

Wall rock loosening ring method of testing adopts sonic method, sonic method test wall rock loosening ring to be based on sound wave spread speed in country rock Change, there is geonetrical attenuation and Physical Attenuation when elastic wave is propagated in rock mass, the different mechanical properties in rock mass Elastic wave will be reflected, scattered and thermal losses on structural plane so that constantly decay causes velocity of wave to reduce to elastic wave energy, by bullet The wave theory of property ripple, the wave equation in unlimited isotropy medium is：

V in formula_p--- compressional wave；

V_s--- shear wave

Boundary condition and primary condition based on isotropy elastic half-space, obtains corresponding velocity of wave, the elastic modelling quantity with rock Relational expression between E, Poisson's Ratio σ, density p is：

V_{p} = \sqrt{\frac{E}{ρ} \cdot \frac{1 - σ}{(1 - 2 σ) (1 + σ)}}

V_{s} = \sqrt{\frac{E}{ρ} \cdot \frac{1}{2 (1 + σ)}}

6. mineral deposit hydrogeology Comprehensive Exploration system as claimed in claim 5, it is characterised in that velocity of longitudinal wave is bored by determining The acoustic transit time of ultrasonic detection device length country rock is calculated in hole, is calculated using single hole method of testing；

Single hole method of testing is：Transmitting transducer F launches in the borehole ultrasonic wave, and along the wall of a borehole propagation, transmitting transducer F are slided Triggering timing circuit timing while transmitting ultrasonic wave, after receive transducer J receives ultrasound information timing is stopped, and measures sound Propagation time t of the ripple in F-J length rock mass：

T=Δ t+t₀+Δt

Δ t = \frac{φ_{D} - φ_{d}}{v}

φ_D- drilling internal diameter；

φ_d- transducer diameter；

The SVEL of water is coupled in v-drilling；

Spread speed of the sound wave in rock-boring be：

v = \frac{L}{t - 2 Δ t}

V in formula --- SVEL in drilling；

L --- distance between transducer F and J；

When wall-rock crack is more, velocity of wave of the velocity of wave relative to depth completely without crack rock is low, and by rock-boring sound wave is measured Changes in distribution curve of the velocity of longitudinal wave in wall rock drill-hole or using time-hole depth curve, judges wall-rock crack scope, transmitting Transducer F launches in the borehole ultrasonic wave, and along the wall of a borehole propagation, triggering meter while transmitting transducer F transmitting ultrasonic waves are slided When road start timing, after receive transducer I receives ultrasound information stop timing s, measure propagation time of the sound wave in F-I, Calculate velocity of wave.

7. mineral deposit hydrogeology Comprehensive Exploration system as claimed in claim 1, it is characterised in that the hydrology sensor includes Level switch, current amount detector, geological image collector；The level switch, current amount detector with water level number It is connected according to detection node；Geological image collector is connected with video monitoring node；The monitor video display include display screen, Ambiguity evaluation module, fuzziness adjusting module；The fuzziness adjusting module is connected with display screen by printed line；The geology Image collection device is used to obtain the image of geological stratification；

The ambiguity evaluation module is used to obtain the geology tomographic image of geological image collector transmission, and calculates figure before and after filtering As statistical information ratio；

The fuzziness adjusting module is connected with ambiguity evaluation module, for adjust original image fuzziness draw final image and Image blur evaluation index；

Step 2, image gray processing, for convenience of the edge extracting of image, R, G, the B using RGB image in Digital Image Processing is each Coloured image is converted into gray level image by the pixel value of individual passage with the transformational relation of gray level image pixel value, and formula is as follows：

Gray=R*0.3+G*0.59+B*0.11；

Step 3, Edge extraction, using the Roberts operator edge detections technical role in digital image processing method in Gray level image obtains the edge of image, and different detective operators have different edge detection templates, according to concrete formwork calculation Intersect the difference of pixel as current pixel value, it is as follows using template：

E (i, j)=| F (i, j)-F (i+1, j+1) |+| F (i+1, j)-F (i, j+1) |；

Step 4, image procossing is filtered process to gray level image to construct image to be evaluated using high pass/low pass filter Reference picture, using 3*3 mean filters, using Filtering Template traversing graph as each pixel, template center is placed in every time Current pixel, the mean value of all pixels is newly worth as current pixel using in template, and template is as follows：

\frac{1}{9} \times [\begin{matrix} 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 1 \end{matrix}];

Step 5, image border statistical information is calculated, and respective edge half-tone information, filtering process before and after image filtering are calculated respectively Front image F statistical informations to be evaluated are sum_orig, and the reference picture F2 statistical informations after filtering process are sum_filter, Specific formula for calculation is as follows：

\begin{matrix} s u m_o r i g = w 1 \times (| F (i, j) - F (i - 1, j) | + | F (i, j) - F (i, j - 1) | + | F (i, j) - F (i, j + 1) | + | F (i, j) - F (i + 1, j) |) \\ + w 2 \times (| F (i, j) - F (i - 1, j - 1) | + | F (i, j) - F (i - 1, j + 1) | + | F (i, j) - F (i + 1, j - 1) | + | F (i, j) - F (i + 1, j + 1) |) \end{matrix},

\begin{matrix} s u m_f i l t e r = w 1 \times (| F 2 (i, j) - F 2 (i - 1, j) | + | F 2 (i, j) - F 2 (i, j - 1) | + | F 2 (i, j) - F 2 (i, j + 1) | + | F 2 (i, j) - F 2 (i + 1, j) |) \\ + w 1 \times (| F 2 (i, j) - F 2 (i - 1, j - 1) | + | F 2 (i, j) - F 2 (i - 1, j + 1) | + | F 2 (i, j) - F 2 (i + 1, j - 1) | + | F 2 (i, j) - F 2 (i + 1, j + 1) |) \end{matrix},

Step 6, image blur index is calculated, the ratio of the image filtering front and rear edges grey-level statistics that step 5 is drawn Value as fuzziness index, for convenience of evaluating, take it is larger for denominator, it is less for molecule, keep the value between (0,1) it Between；

Step 7, according to the DMOS scopes of the best visual effect a corresponding fuzziness indication range [min, max] is drawn, tool Body is：

Fuzziness adjusting range is drawn, using the ambiguity evaluation method in above-mentioned steps 174 width Gaussian modes in LIVE2 are evaluated Paste image, calculates the ambiguity evaluation value of each of which, is then set up using fitting tool plot (value, DMOS) and is evaluated Mapping relations between value value and DMOS, show that corresponding one obscures according to the corresponding DMOS scopes of the best visual effect Evaluation of estimate scope [min, max]；

Step 8, image blur adjustment, if image blur index is less than min, according to step 6, before and after judging image filtering Change is very big, and original image is excessively sharpened, then be filtered adjustment using low pass filter；If being more than max, before judging image filtering After vary less, original image is excessively obscured, then be filtered adjustment using high-pass filter, to reach more preferably visual effect；

8. mineral deposit hydrogeology Comprehensive Exploration system as claimed in claim 7, it is characterised in that filter process image is not Single mode processes view picture evaluation image, but in view of the border of image and center pixel cause wave filter due to the difference of position Difference to its processing mode, according to filter template size recycling filter process figure after boundary pixel is correspondingly ignored Picture, then for boundary pixel takes the method that original pixels are filled to process.

9. mineral deposit hydrogeology Comprehensive Exploration system as claimed in claim 1, it is characterised in that the central processing unit is arranged The signal processing method for having synchronized orthogonal Frequency Hopping Signal blind source separating module, the step quadrature frequency hopping signal blind source separating includes：

Step one, utilizes the array antenna received containing M array element from the Frequency Hopping Signal of multiple synchronized orthogonal frequency hopping sensors, To sampling per reception signal all the way, the M roads discrete time-domain mixed signal after being sampled(k=1,2 ... .) m= 1,2,…,M；The interaction times of different time piece between collection array antenna node, according to the data setup time sequence for obtaining, lead to Cross third index flatness to predict the interaction times of next timeslice between node, by interaction times predicted value and actual value Direct trust value of the relative error as node；The concrete calculation procedure of direct trust value is：Collection network observations node i and section The interaction times of n timeslice between point j：Intervals t are chosen as an observation time piece, with observer nodes i With interaction times of the tested node j in 1 timeslice as observation index, true interaction times are denoted as y_t, n is recorded successively The y of individual timeslice_n, and save it in the communications records table of node i；The interaction times of (n+1)th timeslice of prediction：According to The interaction times setup time sequence of the n timeslice for collecting, using third index flatness next timeslice n+1 is predicted Interaction times between interior nodes i and j, predict interaction times, are denoted asComputing formula is as follows：

{\hat{y}}_{n + 1} = a_{n} + b_{n} + c_{n}

a_{n} = 3 {\hat{y}}_{n + 1}^{(1)} - 3 {\hat{y}}_{n + 1}^{(2)} + {\hat{y}}_{n + 1}^{(3)}

b_{n} = \frac{α}{2 {(1 - α)}^{2}} [(6 - 5 α) {\hat{y}}_{n + 1}^{(1)} - 2 (5 - 4 α) {\hat{y}}_{n + 1}^{(2)} + (4 - 3 α) {\hat{y}}_{n + 1}^{(3)}]

c_{n} = \frac{α^{2}}{2 {(1 - α)}^{2}} [{\hat{y}}_{n + 1}^{(1)} - 2 {\hat{y}}_{n + 1}^{(2)} + {\hat{y}}_{n + 1}^{(3)}]

Wherein：Be respectively once, secondary, Three-exponential Smoothing number, be calculated by equation below：

{\hat{y}}_{n + 1}^{(1)} = α \times y_{n} + (1 - α) \times {\hat{y}}_{n}^{(1)}

{\hat{y}}_{n + 1}^{(2)} = α \times {\hat{y}}_{n + 1}^{(1)} + (1 - α) \times {\hat{y}}_{n}^{(2)}

{\hat{y}}_{n + 1}^{(3)} = α \times {\hat{y}}_{n + 1}^{(2)} + (1 - α) \times {\hat{y}}_{n}^{(3)}

It is the initial value of third index flatness, its value is

{\hat{y}}_{0}^{(1)} = {\hat{y}}_{0}^{(2)} = {\hat{y}}_{0}^{(3)} = \frac{y_{1} + y_{2} + y_{3}}{3}

α is smoothing factor (0 ＜ α ＜ 1), embody trust time attenuation characteristic, i.e., from predicted value more close to timeslice y_tWeight It is bigger, from predicted value more away from timeslice y_tWeight is less；If data fluctuations are larger, and long-term trend amplitude of variation compared with Greatly, α when substantially rapidly rising or falling trend is presented should take higher value (0.6～0.8), increase Recent data to predicting the outcome Impact；When data have a fluctuation, but long-term trend change it is little when, α can between 0.1～0.4 value；If data fluctuations are put down Surely, α should take smaller value (0.05～0.20)；

Calculate direct trust value：

Direct trust value TD of node j_ijTo predict interaction timesWith true interaction times y_n+1Relative error,

Indirect trust values are calculated using calculating formula obtained from multipath trust recommendation mode；Trusted node is collected to node j's Direct trust value：Node i meets TD to all_ikThe credible associated nodes of≤φ inquire its direct trust value to node j, wherein φ is the believability threshold of recommended node, and according to the precision prescribed of confidence level, the span of φ is 0～0.4；Calculate letter indirectly Appoint value：Trust value collected by COMPREHENSIVE CALCULATING, obtains indirect trust values TR of node j_ij,Its In, Set (i) is interacted and its direct trust value meets TD to have with j nodes in the associated nodes of observer nodes i_ikThe section of≤φ Point set；

Comprehensive trust value is drawn by direct trust value and indirect trust values conformity calculation；Comprehensive trust value (T_ij) computing formula such as Under：T_ij=β TD_ij+(1-β)TR_ij, wherein β (0≤β≤1) represents the weight of direct trust value, when β=0, node i and node The calculating that j does not have direct interaction relation, comprehensive trust value arises directly from indirect trust values, and it is more objective to judge；When β=1, section Point i, all from direct trust value, in this case, judges more subjective, Practical Calculation to the comprehensive trust value of node j The value of β is determined as needed；

Step 2, carries out overlapping adding window Short Time Fourier Transform to M roads discrete time-domain mixed signal, obtains M mixed signal Time-frequency domain matrixP=0,1 ..., P-1, q=0,1 ..., N_fft- 1, wherein P represents total Window number, N_fftRepresent FFT length；P, q) time-frequency index is represented, specific time-frequency value isHere N_fftThe length of FFT is represented, p represents adding window number of times, T_sRepresent sampling interval, f_sSample frequency is represented, C is integer, represented The sampling number at Short Time Fourier Transform adding window interval, C ＜ N_fft, and K_c=N_fft/ C is integer, that is to say, that use weight The Short Time Fourier Transform of superposition window；

Step 3, to the frequency-hopping mixing signal time-frequency domain square obtained in step 2Carry out pre- Process；

Step 4, using clustering algorithm estimate each jump jumping moment and respectively jump corresponding normalized mixed moment array to Amount, Hopping frequencies；It is right at p (p=0,1,2 ... the P-1) momentThe frequency values of expression are clustered, a cluster centre for obtaining NumberThe carrier frequency number that the expression p moment is present,Individual cluster centre then represents the size of carrier frequency, uses respectivelyRepresent；To each sampling instant p (p=0,1,2 ... P-1), using clustering algorithm pairEnter Row cluster, it is same availableIndividual cluster centre, usesRepresent；To allAverage and round, obtain To the estimation of source signal numberI.e.：

\hat{N} = r o u n d (\frac{1}{p} Σ_{p = 0}^{P - 1} {\hat{N}}_{p});

Find outMoment, use p_hRepresent, the p to each section of continuous value_hIntermediate value is sought, is usedL=1,2 ... represent L sections are connected p_hIntermediate value, thenRepresent the estimation at l-th frequency hopping moment；Obtained according to estimationp≠p_hAnd the 4th estimate that the frequency hopping moment for obtaining estimates in step and each jump correspondingIt is individual mixed Close matrix column vectorSpecifically formula is：

{\hat{a}}_{n} (l) = \{\begin{matrix} \frac{1}{{\overset{&OverBar;}{p}}_{h} (1)} \cdot Σ_{p = 1, p &NotEqual; p_{h}}^{{\overset{&OverBar;}{p}}_{h} (1)} b_{n, p}^{0} & l = 1, \\ \frac{1}{{\overset{&OverBar;}{p}}_{h} (l) - {\overset{&OverBar;}{p}}_{h} (l - 1)} \cdot Σ_{p = {\overset{&OverBar;}{p}}_{h} (l - 1) + 1, p &NotEqual; p_{h}}^{{\overset{&OverBar;}{p}}_{h} (l)} b_{n, p}^{0} & l > 1, \end{matrix}, n = 1, 2, ..., \hat{N};

HereRepresent that l is jumped correspondingIndividual hybrid matrix Column vector estimate；Estimate the corresponding carrier frequency of each jump, useRepresent that l is jumped corresponding Individual frequency estimation, computing formula is as follows：

{\hat{f}}_{c, n} (l) = \{\begin{matrix} \frac{1}{{\overset{&OverBar;}{p}}_{h} (1)} \cdot Σ_{p = 1, p &NotEqual; p_{h}}^{{\overset{&OverBar;}{p}}_{h} (1)} f_{o}^{n} (p) & l = 1, \\ \frac{1}{{\overset{&OverBar;}{p}}_{h} (l) - {\overset{&OverBar;}{p}}_{h} (l - 1)} \cdot Σ_{p = {\overset{&OverBar;}{p}}_{h} (l - 1) + 1, p &NotEqual; p_{h}}^{{\overset{&OverBar;}{p}}_{h} (l)} f_{o}^{n} (p) & l > 1, \end{matrix}, n = 1, 2, ..., \hat{N};

Step 7, according to source signal time-frequency domain estimate, recovers time domain frequency hopping source signal；To each sampling instant p (p=0,1, 2 ...) frequency domain data Y_n(p, q), q=0,1,2 ..., N_fft- 1 is N_fftThe IFFT conversion of point, obtains p sampling instants corresponding Time domain frequency hopping source signal, uses y_n(p,q_t)(q_t=0,1,2 ..., N_fft- 1) represent；The time domain frequency hopping that above-mentioned all moment are obtained Source signal y_n(p,q_t) process is merged, final time domain frequency hopping synthesizer Signal estimation is obtained, concrete formula is as follows：

s_{n} [k C : (k + 1) C - 1] = \{\begin{matrix} Σ_{m = 0}^{k} y_{n} [m, (k - m) C : (k - m + 1) C - 1] & k < K_{c} \\ Σ_{m = k - K_{c} + 1}^{k} y_{n} [m, (k - m) C : (k - m + 1) C - 1] & k &GreaterEqual; K_{c} \end{matrix}, k = 0, 1, 2, ...