CN1174225C - supersonic method of detecting the fixed-point liquid level in container - Google Patents

supersonic method of detecting the fixed-point liquid level in container Download PDF

Info

Publication number
CN1174225C
CN1174225C CNB011336838A CN01133683A CN1174225C CN 1174225 C CN1174225 C CN 1174225C CN B011336838 A CNB011336838 A CN B011336838A CN 01133683 A CN01133683 A CN 01133683A CN 1174225 C CN1174225 C CN 1174225C
Authority
CN
China
Prior art keywords
remained shock
liquid level
eigenwert
container
decay
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CNB011336838A
Other languages
Chinese (zh)
Other versions
CN1348091A (en
Inventor
马志敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CNB011336838A priority Critical patent/CN1174225C/en
Publication of CN1348091A publication Critical patent/CN1348091A/en
Application granted granted Critical
Publication of CN1174225C publication Critical patent/CN1174225C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of Levels Of Liquids Or Fluent Solid Materials (AREA)

Abstract

The present invention relates to a supersonic method of detecting the fixed-point liquid level in a container, which is characterized in that 1.1, high frequency pulse supersonic waves are transmitted into a container through a container wall by an ultrasonic sensor which is coupled on the outer wall of the container; 1.2, the after-vibration energy of the supersonic waves in the container wall is detected by the ultrasonic sensor, and thus, the change of the after-vibration energy can be measured; 1.3, a signal of the after-vibration energy which is measured actually and attenuation characteristics are compared with a preset characteristic value of the after-vibration energy signal to determine whether a liquid level in the container reaches a control point or not; 1.3.1, if the duration of the after-vibration energy signal is short and the attenuation of the after-vibration energy signal is quick, the liquid level reaches or exceeds the measurement control point; 1.3.2, if the duration of the after-vibration energy signal is long and the attenuation of the after-vibration energy signal is slow, the liquid level is lower than the measurement control point; 1.4, according to an actual measurement result, a command or an alarm is sent to relevant equipment or operation persons.

Description

A kind of method of utilizing the ultrasound examination fixed-point liquid level in container
Technical field
The present invention relates to a kind of method of utilizing the ultrasound examination fixed-point liquid level in container.
Background technology
In many industry spot, the high low level indication of the inner liquid level of closed container or pressure vessel has crucial meaning with warning for the control of production run and the safety of production.Because the singularity of container work condition, the non-intervention type detection method, easy for installation because of it, not pressure-bearing, pollution-free, do not have corrosion, easily make blast resistance construction and be subjected to the user and welcome.Especially for the container of splendid attire high pressure, inflammable and explosive and strongly corrosion liquid, often non-intervention type is first-selected method, even in some high pressure vessels, consider and the regulation of relevant safety standard for the security standpoint of structure, do not allow increase equipment mounting hole in these containers.Owing to the pick-up unit that existing similar detection method adopted is installed trouble, poor reliability has influenced applying of this detection method.
Summary of the invention
Technical matters to be solved by this invention provides a kind of method of utilizing the ultrasound examination fixed-point liquid level in container, the capacity volume variance that this method utilizes gas, liquid medium that the difference of acoustic wave transmission rate is caused remaining in the remained shock signal in the chamber wall is carried out the liquid level indication or is reported to the police, thereby can obtain good security and reliability, its device installs also very convenient.The technical scheme that solves the present technique problem is: a kind of method of utilizing the ultrasound examination fixed-point liquid level in container is characterized in that: 1.1 utilize the ultrasonic sensor that is coupling in container outer wall to launch the high-frequency impulse ultrasound wave through chamber wall in container; 1.2 with the ultrasound wave remained shock energy signal in the ultrasonic sensor detection receptacle wall, the variation of promptly measuring the remained shock energy signal; 1.3 the signal of the actual remained shock energy that records and decay characteristics and the remained shock energy signal standard feature value of adjusting are in advance compared, to determine whether liquid level reaches the reference mark in the container; 1.3.1 when recording remained shock energy duration weak point, the decay of remained shock energy signal is fast, promptly liquid level meets or exceeds Surveying Control Point; 1.3.2 when recording remained shock energy longer duration, the decay of remained shock energy signal is slow, promptly liquid level is lower than Surveying Control Point; 1.4 send instruction or warning to the relevant device or the personnel that control according to measured result.And, the length of the ultrasound wave remained shock signal duration in the 2.1 usefulness ultrasonic sensor detection receptacle walls; 2.2 the actual remained shock signal duration that records and the rate of decay and the standard feature value of adjusting are in advance compared, to determine whether liquid level reaches the reference mark in the container; 2.2.1 when recording remained shock duration weak point, remained shock signal attenuation is fast, T≤T F, K≤K F, promptly liquid level meets or exceeds Surveying Control Point; 2.2.2 when recording the remained shock longer duration, remained shock signal attenuation is slow, T 〉=T E, K 〉=K E, promptly liquid level is lower than Surveying Control Point; 2.4 send instruction or warning to the relevant device or the personnel that control according to measured result.In the technical program, when the actuating medium of the internal tank at corresponding sound wave incidence point place not simultaneously, the size of the transmission effect of sound wave also is different, especially transmissivity big or small far different under gas, two kinds of situations of liquid medium differs greatly so remain in the energy of the remained shock signal in the chamber wall under two kinds of operating modes.With two kinds of media of water and air is example, is 9.7 * 10 by steel plate (when chamber wall is steel plate) to the sound transmission rate of water -2, then only be 1.9 * 10 to the sound transmission rate of air by steel plate -5Obviously, for the former, the sound wave in the steel plate is after the reflection of experience several times, and most of energy transmission is in liquid medium, and for the latter, then most of energy is trapped in and comes back reflective in the steel plate, till decay exhausts gradually.Therefore, concerning the former, the remained shock signal attenuation in its steel plate is fast, and the duration is short; Concerning the latter, the remained shock signal attenuation in the steel plate is slow, longer duration.Therefore, the technical program is by the variation of the remained shock energy in the detection receptacle wall (steel plate) or the speed of its signal attenuation, and promptly whether the length of remained shock energy duration and then judgement internal tank liquid level reach Surveying Control Point.Measured test is the result show, the ratio of two kinds of operating mode (liquid level reaches the reference mark and liquid level is lower than the reference mark) following remained shock signal attenuation times can reach 5: 1~10: 1, is enough to differentiate reliably therefore whether liquid level reaches the control limit in the container.
Description of drawings
Fig. 1 is the application synoptic diagram of the method for the invention.
Fig. 2 is remained shock and the remained shock signal waveforms that liquid level reaches or be higher than the reference mark.
Fig. 3 is remained shock and the remained shock signal waveforms that liquid level is lower than the reference mark.
Fig. 4 is the schematic block diagram of a kind of embodiment of equipment therefor of the present invention.
Fig. 5 is the schematic block diagram of second kind of embodiment of equipment therefor of the present invention.
Fig. 6 is the oscillogram that liquid level reaches or the remained shock signal envelope involves comparer output square wave when being higher than the reference mark.
Fig. 7 be liquid level when being lower than the reference mark remained shock signal envelope involve the oscillogram of comparer output square wave.
Fig. 8 be liquid level reach or when being higher than the reference mark integration output (last figure) of remained shock signal envelope ripple export the synoptic diagram that (figure below) concerns with comparer.
Fig. 9 is that the integration output (last figure) of liquid level remained shock signal envelope ripple when being lower than the reference mark is exported the synoptic diagram that (figure below) concerns with comparer.
Embodiment
With reference to Fig. 1~Fig. 9, method of the present invention is carried out according to the following steps: 1.1 utilize the ultrasonic sensor 1 that is coupling in container outer wall to launch the high-frequency impulse ultrasound wave through chamber wall 10 in container; 1.2 with the ultrasound wave remained shock energy signal in the ultrasonic sensor 1 detection receptacle wall 10, the variation of promptly measuring the remained shock energy signal; 1.3 the signal of the actual remained shock energy that records and decay characteristics and the remained shock energy signal standard feature value of adjusting are in advance compared, to determine whether liquid level 11 reaches the reference mark in the container; 1.3.1 when recording remained shock energy duration weak point, the decay of remained shock energy signal is fast, promptly liquid level 11 meets or exceeds Surveying Control Point; 1.3.2 when recording remained shock energy longer duration, the decay of remained shock energy signal is slow, promptly liquid level 11 is lower than Surveying Control Point; 1.4 send instruction or warning to the relevant device or the personnel that control according to measured result.Except that detecting the remained shock energy signal, can also detect the duration of remained shock signal, promptly further technical scheme can be following second kind of embodiment: the length of the ultrasound wave remained shock signal duration in 2.1 usefulness ultrasonic sensors, the 1 detection receptacle wall 10; 2.2 the actual remained shock signal duration that records and the rate of decay and the standard feature value of adjusting are in advance compared, to determine whether liquid level 11 reaches the reference mark in the container; 2.2.1 when recording remained shock duration weak point, remained shock signal attenuation is fast, T≤T f, K≤K f, promptly liquid level 11 meets or exceeds Surveying Control Point; 2.2.2 when recording the remained shock longer duration, remained shock signal attenuation is slow, T 〉=T e, K 〉=K e, promptly liquid level 11 is lower than Surveying Control Point; Then, send instruction or warning according to measured result to the relevant device or the personnel that control.Described chamber wall 10 can be steel plate or other nonmetallic materials; The frequency range of ultrasound wave high-frequency impulse is 1~5MHz.Wherein: T e, K eThe eigenwert of die-away time when being lower than the reference mark for liquid level 11, T f, K fThen for reaching or the eigenwert of die-away time during the high level cadre reference mark.Because different container wall thickness, the die-away time of remained shock can be different, thereby had better not be with the standard of an identical standard feature value as all tank fill level 11 situation judgings, so set up the standard feature value of each container self, specifically can carry out according to the following steps: the eigenwert T during the measuring vessel canful repeatedly FiAnd K FiAnd the eigenwert T during slack tank EiAnd K Ei, wherein, T Fi, K FiFor liquid level 11 reaches or eigenwert die-away time when being higher than the reference mark; T Ei, K EiEigenwert die-away time when being lower than the reference mark for liquid level 11, i=1,2 ..., n;
Average and as the reference standard eigenwert by following formula:
T f=K 0(∑T fi)/n
K f=K 0(∑K fi)/n
T e=K i(∑T ei)/n
K e=K 1(∑K ei)/n
n≤100,K 0=1.10~1.30,K 1=0.75~0.95。
Owing to aging, the volatilization of couplant of sensor performance, the factors such as corrosion of container outer wall steel plate all can change along with the time, and then influence remained shock signal attenuation characteristics value and the device long-term work reliability, for fear of issuable erroneous judgement, suggestion: the eigenwert and the historical standard eigenwert that measure in a period of time are weighted on average, to obtain revised new standard eigenwert, that is:
K eNEW=K 2·K e+K 3·K eOLD
T eNEW=K 2·T e+K 3·T eOLD
K fNEW=K 2·K f+K 3·K fOLD
T fNEW=K 2·T f+K 3·T fOLD
Wherein: T e, K eIndividual features value during for the slack tank that records in the recent period, K f, T fIndividual features value during for the canful that records in the recent period, K 2, K 3Be weighting coefficient, K 2=0.05~0.10, K 3=0.90~0.95, K EOLD, T EOLD, K FOLD, T FOLDBe corresponding historical standard eigenwert.
The actual remained shock signal duration that records and the rate of decay and the new feature value of adjusting are in advance compared, promptly by to parameter from the motion tracking correction, get rid of because of wear out, factors such as volatilization and corrosion are to the influence of eigenwert.The applied device of the present invention (seeing Fig. 1, Fig. 4 and Fig. 5) can comprise: parts such as ultrasonic sensor 1, radiating circuit 2, reception amplifying circuit 3, detecting circuit 4, integrator 4 ' (only embodiment one has), level comparison circuit 5, single chip circuit 6, storer 7, warning indicating circuit 8 and telecommunication circuit 9.Because each several part itself belongs to typical circuit of the prior art, so each circuit is not further given unnecessary details itself.The function of each several part is as follows:
The effect of radiating circuit 2 is under the control of single chip circuit 6, produces high-frequency pulse signal, excitation ultrasonic sensor 1.The effect of ultrasonic sensor 1 is to produce the high-frequency impulse ultrasound wave and be coupled on the outer wall of container under the excitation of emission voltage, meanwhile, sensor 1 also becomes electric signal and feed-in reception amplifier circuit 3 with the ultrasound wave remained shock conversion of signals in the receiving vessel outer wall steel plate.
The effect of amplifying circuit 3 is that faint remained shock signal is amplified, scope when surveying to improve, thus improve reliability when surveying.
The effect of detecting circuit 4 is that the remained shock signal after amplifying is adjusted into unipolar signal, and filtering high fdrequency component, thereby obtain the positive half of envelope ripple (see Fig. 6, Fig. 7, two figure represent detection output, comparer 1 output and comparer 2 outputs from top to bottom respectively) of remained shock signal.
The effect of comparator circuit 5 is according to the comparative level that sets in advance envelope signal to be made into two square-wave signals (seeing Fig. 6, Fig. 7), square wave V 1Width be T 1, be time of 90% of saturated amplitude corresponding to remained shock signal attenuation.Square wave V 2Width be T 2, be time of 10% of saturated amplitude corresponding to remained shock signal attenuation, and go out T by single chip computer measurement 1And T 2Size.
Single-chip microcomputer 6 is as intelligent parts, and its effect is: 1. hyperacoustic emission is regularly carried out in control; 2. measure T die-away time of remained shock signal 1And T 2And calculate T=(T 1+ T 2)/2, K=(T 2-T 1)/T judges the liquid level of container thus.As T 〉=T e, and K 〉=K eThe time, this liquid level of decidable is lower than the reference mark; Otherwise, if T≤T fAnd K≤K fThe time, then decidable liquid level 11 meets or exceeds the reference mark, wherein, and T e, K eThe eigenwert of die-away time when being lower than the reference mark for liquid level, T f, K fThen reach or the eigenwert of die-away time when being higher than the reference mark for liquid level 11.Under other situation, then keep current indication and alarm condition constant.3. control warning indicating circuit 8 and relay output signal; 4. telecommunication administration is so that in time pass to System Control Center with the work information of equipment; 5. survey the automatic adjusting and the adaptive tracing of object remained shock signal characteristic value, and in time the eigenwert that system is new is saved in the storer 7, detects the foundation of judging as the back.
The effect of storer 7 is standard feature values of depositing the remained shock signal attenuation time under the different operating modes (liquid level and no liquid level are arranged).
Integrator 4 (is seen Fig. 4, only embodiment one has) effect be that envelope ripple to the remained shock signal carries out integration, the output of various operating mode lower integral devices 4 and the relation of envelope ripple can be referring to Fig. 6~Fig. 9, integrator 4 ' output give comparer, and compare with the eigenwert of adjusting in advance, to determine whether liquid level reaches the reference mark.The effect of all the other each several parts shown in Figure 4 is identical with Fig. 5 (embodiment two).

Claims (1)

1. method of utilizing the ultrasound examination fixed-point liquid level in container is characterized in that:
In container, launch the high-frequency impulse ultrasound wave 1.1 utilize the ultrasonic sensor that is coupling in container outer wall through chamber wall;
1.2 with the ultrasound wave remained shock energy in the ultrasonic sensor detection receptacle wall, the variation of promptly measuring the remained shock energy, and the respectively standard feature value of this container of adjusting according to the following steps:
1.2.1 the remained shock eigenwert signal duration T during many measuring vessel canfuls FiWith rate of decay eigenwert K FiAnd the remained shock eigenwert signal duration T during slack tank EiWith rate of decay eigenwert K Ei, wherein, T Fi, K FiThe value that reaches or measure when being higher than the reference mark for liquid level; T Ei, K EiThe value of measuring when being lower than the reference mark for liquid level, i=1,2 ..., n; With
1.2.2 average and as the standard feature value by following formula:
T f=K 0(∑T fi)/n
K f=K 0(∑K fi)/n
T e=K 1(∑T ei)/n
K e=K 1(∑K ei)/n
N≤100, K 0=1.10~1.30, K 1=0.75~0.95, T e, K eRemained shock standard feature signal duration value when being respectively the slack tank that records in the recent period and rate of decay standard feature value, T f, K fRemained shock standard feature signal duration value when being respectively the canful that records in the recent period and rate of decay standard feature value;
1.3 the standard feature value and the historical standard eigenwert that record in the recent period in a period of time are weighted on average, to obtain revised new standard eigenwert, that is:
K eNEW=K 2·K e+K 3·K eOLD
T eNEW=K 2·T e+K 3·T eOLD
K fNEW=K 2·K f+K 3·K fOLD
T fNEW=K 2·T f+K 3·T fOLD
Wherein:, K 2, K 3Be weighting coefficient, K 2=0.05-0.10, K 3=0.90-0.95, K EOLD, T EOLD, K FOLD, T FOLDBe corresponding historical standard eigenwert, K ENEW, T ENEW, K FNEW, T FNEWBe corresponding new standard eigenwert;
1.4 replace the historical standard eigenwert with revised new standard eigenwert, with get rid of because of wear out, volatilization and corrosion factor be to the influence of eigenwert;
1.5 actual remained shock energy signal duration T that records and rate of decay K are compared with the remained shock energy signal eigenwert of adjusting in advance respectively, to determine whether liquid level reaches the reference mark in the container;
1.5.1 when recording that remained shock energy signal duration T is shorter than or remained shock standard feature signal duration value T when equaling new canful f, and the rate of decay standard feature value K of remained shock energy signal rate of decay K during faster than new canful f, T≤T f, K≤K f, promptly liquid level meets or exceeds Surveying Control Point; With
1.5.2 when recording that remained shock energy signal duration T is longer than or remained shock standard feature signal duration value T when equaling new slack tank e, and the rate of decay standard feature value K of remained shock energy signal rate of decay K when being slower than new slack tank e, T 〉=T e, K 〉=K e, promptly liquid level is lower than Surveying Control Point;
1.6 send instruction or warning to the relevant device or the personnel that control according to measured result.
CNB011336838A 2001-11-19 2001-11-19 supersonic method of detecting the fixed-point liquid level in container Expired - Fee Related CN1174225C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB011336838A CN1174225C (en) 2001-11-19 2001-11-19 supersonic method of detecting the fixed-point liquid level in container

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB011336838A CN1174225C (en) 2001-11-19 2001-11-19 supersonic method of detecting the fixed-point liquid level in container

Publications (2)

Publication Number Publication Date
CN1348091A CN1348091A (en) 2002-05-08
CN1174225C true CN1174225C (en) 2004-11-03

Family

ID=4672017

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB011336838A Expired - Fee Related CN1174225C (en) 2001-11-19 2001-11-19 supersonic method of detecting the fixed-point liquid level in container

Country Status (1)

Country Link
CN (1) CN1174225C (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101657702B (en) * 2007-03-17 2012-08-22 莫伯蕾有限公司 Ultrasonic level measurement apparatus and method having variable transmit power

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102322920A (en) * 2011-08-24 2012-01-18 四川大学 Method for measuring existence of powdery materials in specific material position of container and implementation material level switch thereof
CN103604483B (en) * 2013-11-29 2017-01-11 国家电网公司 Method for measuring oil level of oil charging insulator of electric power system through ultrasonic waves
CN104501909B (en) * 2014-12-04 2018-05-18 哈尔滨工程大学 A kind of small-range liquid level emasuring device and measuring method based on ultrasonic wave
CN105181069B (en) * 2015-09-01 2017-02-01 深圳麦开网络技术有限公司 Inside-container liquid volume measuring method and device based on acoustic detection
CN106017612A (en) * 2016-03-03 2016-10-12 江苏迅创科技有限公司 Liquid level alarm measuring method not influenced by environment factor
CN105606181B (en) * 2016-03-14 2017-03-29 陕西师范大学 The method of oil level in a kind of oil-filled insulator of ultrasound detection
CN108562349A (en) * 2018-06-01 2018-09-21 中国科学院声学研究所 A kind of ultrasonic wave detecting system and its method
DE102018222339A1 (en) 2018-12-19 2020-06-25 Robert Bosch Gmbh Device for determining a fill level of a container for holding a liquid and a container with such a device
CN110736522A (en) * 2019-08-26 2020-01-31 广西电网有限责任公司电力科学研究院 oil level detection method for sealed oil-filled equipment
CN110487357A (en) * 2019-09-11 2019-11-22 广州秀明环保科技有限公司 A kind of non-contact type ultrasonic liquid level monitor and liquid level inducing method
EP3822660A1 (en) * 2019-11-13 2021-05-19 ABB Schweiz AG Integrity detection system for an ultrasound transducer
CN111974600A (en) * 2020-07-17 2020-11-24 重庆海浦洛自动化科技有限公司 Pipeline cleanable paint conveying and mixing system and pipeline cleaning method thereof
CN112378996A (en) * 2020-11-16 2021-02-19 甘肃黑驴王子生物科技有限公司 Method for rapidly detecting donkey milk adulteration
CN114942042A (en) * 2022-04-15 2022-08-26 深圳市帝拓电子有限公司 Method and circuit for detecting material parameters

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101657702B (en) * 2007-03-17 2012-08-22 莫伯蕾有限公司 Ultrasonic level measurement apparatus and method having variable transmit power

Also Published As

Publication number Publication date
CN1348091A (en) 2002-05-08

Similar Documents

Publication Publication Date Title
CN1174225C (en) supersonic method of detecting the fixed-point liquid level in container
US20070169549A1 (en) Method and apparatus for sensing fuel levels in tanks
US7114373B2 (en) Fluid monitoring
CN101571407B (en) Excitation method of vibrating wire sensor
CA2573518A1 (en) Sonar sand detection
CN102628437B (en) System for indirectly measuring flow and pressure of constant displacement pump driven by permanent magnet servo motor
US7039530B2 (en) Fluid measurement
CN101246213A (en) Double-comparer type ultrasonic distance measurer
US4280205A (en) Ultrasonic fork height measuring apparatus
CN1975348A (en) Ultrasonic wave flow and liquid level measuring method based on multi-pulsion automatic gaining control
CN111323100A (en) Ultrasonic gas meter fault diagnosis system and method
CN1547672A (en) Ultrasonic doppler effect speed measurement
EP3971537A1 (en) Ultrasonic transducer health status monitor
CN104913873B (en) A kind of Ultrasonic-wave Hydraulic measuring method and its system for improving multifactor impact
CN113551741B (en) Ultrasonic liquid level interface quantitative calculation method for closed container containing mixed medium
CN214503464U (en) Wood structure water content measuring device based on wave velocity method
Grabec et al. A comparison of high-performance acoustic emission transducers
CN113808092A (en) Method, system, device and medium for detecting debonding defect of steel pipe concrete interface
CN2508235Y (en) Fixen point liquid level ultrasonic measurer for container
Mai et al. Applicability of radar level gauges in wave monitoring
CN2288425Y (en) Ultrasonic liquid level measurer for oil tank
CN213632325U (en) Ultrasonic liquid level measuring device for underground main drainage pump room
CN1250154A (en) Automatic supersonic oil tank level measurement technology based on half-wave transmission principle
CN113030275B (en) Wood structure water content measuring device and method based on wave velocity method
EP3971536A1 (en) Ultrasonic flowmeter with vibration-resistant operating mode

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee