SCHIESSL SABRINA (DE)
| US20050072226A1 | 2005-04-07 | |||
| US20110288795A1 | 2011-11-24 | |||
| US6587796B2 | 2003-07-01 |
| Verfahren zum Betreiben einer akustischen Messvorrichtung (10), bei dem mindestens ein Schallmesswandler (12) der akustischen Messvorrichtung (10) bereitgestellt wird, wobei der mindestens eine Schallmesswandler (12) ausgebildet ist zum Aussenden und/ oder zum Empfangen von Schallsignalen, ein Ausgangssignal (40) von dem mindestens einen Schall¬ messwandler (12) abgegeben und ein Echosignal (42) des Ausgangssignals (40) von dem mindestens einen Schall¬ messwandler (12) empfangen wird, ein zeitlicher Verlauf des Echosignals (42) erfasst wird, aus dem zeitlichen Verlauf des Echosignals (42) ein erster Schätzwert (T_sig_l) einer Laufzeit des Schallsignals vom Beginn des Ausgangssignals (40) bis zum Beginn des Echosignals (42) bestimmt wird, mittels einer Kreuzkorrelation zwischen dem Echosignal (42) und einem Referenzsinussignal (44) ein Verlauf einer Kreuzkorrelationsfunktion (R) abhängig von einer Frequenz und einer Phase des Referenzsinussignals (44) ermittelt wird, ein Wert (f_0) der Frequenz des Referenzsinussignals (44) und ein Wert (phi_0) der Phase des Referenzsi¬ nussignals (44) ermittelt werden, für die ein Maxi¬ malwert der Kreuzkorrelationsfunktion (R) erreicht wird und, abhängig von dem Wert (f_0) der Frequenz des Referenzsinussignals (44) und dem Wert (phi_0) der Phase des Referenzsinussignals (44), für die der Maximalwert der Kreuzkorrelationsfunktion (R) erreicht wird, und dem ersten Schätzwert (T_sig_l) der Laufzeit des Schall¬ signals vom Beginn des Ausgangssignals (40) bis zum Beginn des Echosignals (42) ein aktualisierter Schätzwert (T_sig_2) der Laufzeit des Schallsignals vom Beginn des Ausgangssignals (40) bis zum Beginn des Echosignals (42) bestimmt wird. 2. Verfahren nach Anspruch 1, wobei die Kreuzkorrelations¬ funktion (R) der Kreuzkorrelation zwischen dem Echosignal (42) und dem Referenzsinussignal (44) abhängig von der Frequenz und für Phasenwerte 0 und Π/2 des Referenzsi¬ nussignals (44) ermittelt wird. 3. Verfahren nach Anspruch 1 oder 2, wobei das Bestimmen des aktualisierten Schätzwerts (T_sig_2) der Laufzeit des Schallsignals vom Beginn des Ausgangssignals (40) bis zum Beginn des Echosignals (42) mittels eines Iterations¬ verfahrens erfolgt. 4. Verfahren nach einem der vorhergehenden Ansprüche, wobei das Echosignal (42) auf einen Wert eines maximalen Betrags einer Amplitude des Echosignals (42) normiert wird. 5. Verfahren nach Anspruch 4, wobei der erste Schätzwert (T_sig_l) der Laufzeit des Schallsignals vom Beginn des Ausgangssignals (40) bis zum Beginn des Echosignals (42) bestimmt wird abhängig von einem Zeitpunkt (tl) eines ersten Nulldurchgangs des normierten Echosignals (42). 6. Verfahren nach Anspruch 4 oder 5, wobei der aktualisierte Schätzwert (T_sig_2) der Laufzeit des Schallsignals vom Beginn des Ausgangssignals (40) bis zum Beginn des Echo¬ signals (42) bestimmt wird abhängig von einem Zeitpunkt (t2) eines Nulldurchgangs des Referenzsinussignals (44). 7. Verfahren nach einem der vorhergehenden Ansprüche, wobei der mindestens eine Schallmesswandler (12) ein Ultraschallmesswandler ist. 8. Vorrichtung zum Betreiben einer akustischen Messvorrichtung, die dazu ausgebildet ist mindestens ein Schallmesswandlers (12) der akustischen Messvorrichtung (10) bereitzustellen, wobei der mindestens eine Schallmesswandler (12) ausgebildet ist zum Aussenden und/oder zum Empfangen von Schallsignalen, - ein Ausgangssignal (40) von dem mindestens einen Schallmesswandler (12) abzugeben und ein Echosignal (42) des Ausgangssignals (40) von dem mindestens einen Schallmesswandler (12) zu empfangen, einen zeitlichen Verlaufs des Echosignals (42) zu erfassen, aus dem zeitlichen Verlauf des Echosignals (42) einen ersten Schätzwert (T_sig) einer Laufzeit des Schall¬ signals vom Beginn des Ausgangssignals (40) bis zum Beginn des Echosignals (42) zu bestimmen, - mittels einer Kreuzkorrelation zwischen dem Echosignal (42) und einem Referenzsinussignal (44) einen Verlauf einer Kreuzkorrelationsfunktion (R) abhängig von einer Frequenz und einer Phase des Referenzsinussignals (44) zu ermitteln, - einen Wert (f_0) der Frequenz des Referenzsinussignals (44) und einen Wert (phi_0) der Phase des Referenz¬ sinussignals (44) zu ermitteln, für die ein Maximalwert der Kreuzkorrelationsfunktion (R) erreicht wird, abhängig von dem Wert (f_0) der Frequenz des Refe- renzsinussignals (44) und dem Wert (phi_0) der Phase des Referenzsinussignals (44), für die der Maximalwert der Kreuzkorrelationsfunktion (R) erreicht wird, und dem ersten Schätzwert (T_sig_l) der Laufzeit des Schall¬ signals vom Beginn des Ausgangssignals (40) bis zum Beginn des Echosignals (42) einen aktualisierten Schätzwert (T_sig_2) der Laufzeit des Schallsignals vom Beginn des Ausgangssignals (40) bis zum Beginn des Echosignals (42) zu bestimmen. |
Verfahren und Vorrichtung zum Betreiben einer akustischen Messvorrichtung
Die Erfindung betrifft ein Verfahren und eine entsprechende Vorrichtung zum Betreiben einer akustischen Messvorrichtung.
Eine akustische Messvorrichtung kann insbesondere zur Durch- flussmessung von Fluiden in Fluidleitungen sowie zur Bestimmung der Zusammensetzung des Fluids verwendet werden. Ein Schallmesswandler der akustischen Messvorrichtung kann sowohl als Schallerzeuger als auch als Schallempfänger arbeiten. Aus der Laufzeit der Schallimpulse können Rückschlüsse auf die Ge- schwindigkeit des Fluids in der Fluidleitung gezogen werden. Des Weiteren kann die akustische Messvorrichtung insbesondere zum Bestimmen einer Höhe eines Fluidniveaus in einem Fluidbehälter verwendet werden. Die Aufgabe der Erfindung ist, ein Verfahren und eine korrespondierende Vorrichtung zum Betreiben einer akustischen Messvorrichtung zu schaffen, das beziehungsweise die eine hohe Messgenauigkeit der akustischen Messvorrichtung ermöglicht. Die Aufgabe wird gelöst durch die Merkmale der unabhängigen Patentansprüche. Vorteilhafte Weiterbildungen der Erfindung sind in den Unteransprüchen gekennzeichnet.
Die Erfindung zeichnet sich durch ein Verfahren und eine korrespondierende Vorrichtung zum Betreiben einer akustischen Messvorrichtung aus, bei dem mindestens ein Schallmesswandler der akustischen Messvorrichtung bereitgestellt wird, wobei der mindestens eine Schallmesswandler ausgebildet ist zum Aussenden und/oder zum Empfangen von Schallsignalen. Ein Ausgangssignal wird von dem mindestens einen Schallmesswandler abgegeben, und ein Echosignal des Ausgangssignals wird von dem mindestens einen Schallmesswandler empfangen. Ein zeitlicher Verlauf des Echo- Signals wird erfasst. Aus dem zeitlichen Verlauf des Echosignals wird ein erster Schätzwert einer Laufzeit des Schallsignals vom Beginn des Ausgangssignals bis zum Beginn des Echosignals bestimmt . Mittels einer Kreuzkorrelation zwischen dem Echosignal und einem Referenzsinussignal wird ein Verlauf einer Kreuz ¬ korrelationsfunktion abhängig von einer Frequenz und einer Phase des Referenzsinussignals ermittelt. Ein Wert der Frequenz des Referenzsinussignals und ein Wert der Phase des Referenzsi ¬ nussignals werden ermittelt, für die ein Maximalwert der Kreuz- korrelationsfunktion erreicht wird. Abhängig von dem Wert der Frequenz des Referenzsinussignals und dem Wert der Phase des Referenzsinussignals, für die der Maximalwert der Kreuzkor ¬ relationsfunktion erreicht wird, und dem ersten Schätzwert der Laufzeit des Schallsignals vom Beginn des Ausgangssignals bis zum Beginn des Echosignals wird ein aktualisierter Schätzwert der Laufzeit des Schallsignals vom Beginn des Ausgangssignals bis zum Beginn des Echosignals bestimmt.
Der Wert der Frequenz des Referenzsinussignals, für den der Maximalwert der Kreuzkorrelationsfunktion erreicht wird, ist repräsentativ für eine maximale Übereinstimmung der Frequenz des Referenzsinussignals mit der Frequenz des Echosignals. Der Wert der Phase des Referenzsinussignals, für den der Maximalwert der Kreuzkorrelationsfunktion erreicht wird, ist repräsentativ für eine maximale Übereinstimmung der Phase des Referenzsinussignals mit der Phase des Echosignals.
Dies hat den Vorteil, dass der aktualisierte Schätzwert der Laufzeit des Schallsignals von einem Beginn des Ausgangssignals bis zu einem Beginn des Echosignals mittels der Kreuzkorrelation sehr genau bestimmt werden kann. Damit ist es des Weiteren möglich, das Echosignal mit einer niedrigen Abtastrate abzu ¬ tasten, und dabei dennoch eine genaue Bestimmung der Laufzeit des Schallsignals vom Beginn des Ausgangssignals bis zum Beginn des Echosignals zu erreichen. Aufgrund der Möglichkeit des Abtastens des Echosignals mit einer niedrigen Abtastrate kann insbesondere auch die Notwendigkeit des Einsatzes eines schnellen Analog- Digital-Wandlers zur Bestimmung der Laufzeit des Schallsignals von einem Beginn des Ausgangssignals bis zu einem Beginn des Echosignals vermieden werden. In einer vorteilhaften Ausgestaltung wird die Kreuzkorrelationsfunktion der Kreuzkorrelation zwischen dem Echosignal und dem Referenzsinussignal abhängig von der Frequenz und für Phasenwerte 0 und Π/2 des Referenzsinussignals ermittelt. Dies hat den Vorteil, dass aufgrund der Periodizität der Kreuz- korrelationsfunktion bei der Berechnung der Laufzeit des Schallsignals auf eine umfassende Phasenvariation verzichtet werden kann, und die Berechnung auf die Benutzung von zwei Phasenwerten beschränkt werden kann. Damit kann die Berechnung der Laufzeit des Schallsignals mit sehr geringem Rechenaufwand erfolgen.
In einer weiteren vorteilhaften Ausgestaltung erfolgt das Bestimmen des aktualisierten Schätzwerts der Laufzeit des Schallsignals vom Beginn des Ausgangssignals bis zum Beginn des Echosignals mittels eines Iterationsverfahrens. Dies hat den Vorteil, dass die Berechnung der Laufzeit des Schallsignals vom Beginn des Ausgangssignals bis zum Beginn des Echosignals sehr effektiv erfolgen kann.
In einer weiteren vorteilhaften Ausgestaltung wird das Echo- signal auf einen Wert eines maximalen Betrags einer Amplitude des Echosignals normiert. Dies hat den Vorteil, dass eine Bestimmung des ersten Schätzwerts der Laufzeit des Schallsignals vom Beginn des Ausgangssignals bis zum Beginn des Echosignals in einfacher Weise erfolgen kann.
In einer weiteren vorteilhaften Ausgestaltung wird der erste Schätzwert der Laufzeit des Schallsignals vom Beginn des Ausgangssignals bis zum Beginn des Echosignals bestimmt abhängig von einem Zeitpunkt eines ersten Nulldurchgangs des normierten Echosignals. Dies hat den Vorteil, dass der erste Nulldurchgang des normierten Echosignals gut zur Bestimmung des Beginns des Echosignals herangezogen werden kann. In einer weiteren vorteilhaften Ausgestaltung wird der aktualisierte Schätzwert der Laufzeit des Schallsignals vom Beginn des Ausgangssignals bis zum Beginn des Echosignals bestimmt, abhängig von einem Zeitpunkt eines Nulldurchgangs des Refe ¬ renzsinussignals. Dies hat den Vorteil, dass Nulldurchgänge des Referenzsinussignals gut zur Bestimmung des Beginns des Echo ¬ signals herangezogen werden können.
In einer weiteren vorteilhaften Ausgestaltung ist der mindestens eine Schallmesswandler ein Ultraschallmesswandler. Dies hat den Vorteil, dass Ultraschall für viele Anwendungen mit Fluiden, beispielsweise das Ermitteln einer Strömungsgeschwindigkeit des Fluids in einer Fluidleitung, sehr gut und vielfältig einsetzbar ist .
Ausführungsbeispiele der Erfindung sind im Folgenden anhand der schematischen Zeichnungen erläutert.
Es zeigen:
Figur 1 eine schematische Ansicht einer akustischen Mess ¬ vorrichtung,
Figur 2 ein zeitlicher Verlauf eines Ausgangssignals und von
Echosignalen,
Figur 3 ein weiterer zeitlicher Verlauf eines normierten
Echosignals ,
Figur 4 ein weiterer zeitlicher Verlauf eines normierten
Echosignals und eines Referenzsinussignals und
Figur 5 ein Ablaufdiagramm eines Programms zum Betreiben einer akustischen Messvorrichtung. Elemente gleicher Konstruktion oder Funktion sind figurenübergreifend mit den gleichen Bezugszeichen versehen.
Figur 1 zeigt eine akustische Messvorrichtung 10, die mindestens einen Schallmesswandler 12 hat. Der mindestens eine Schallmesswandler 12 ist vorzugsweise als Ultraschallmesswandler ausgebildet .
In einem Messrohr 16 mit einer Längsachse L strömt ein Fluid 20. Das Fluid 20 hat eine Strömungsgeschwindigkeit v. Das Fluid 20 kann eine Flüssigkeit oder ein gasförmiges Medium sein. Ins ¬ besondere kann das Fluid 20 die Frischluft oder das Abgas einer Brennkraftmaschine sein. Der mindestens eine Schallmesswandler 12 kann entweder in Richtung oder in jedem beliebigen Winkel zu der Längsachse L Schallsignale aussenden und empfangen. Aus einer Laufzeit der Schallsignale kann die Strömungsgeschwindigkeit v des Fluids 20 berechnet werden. Ein möglicher Weg der Schall ¬ signale von dem mindestens einen Schallmesswandler 12 zu dem Messrohr 16 und wieder zurück oder zu einem zweiten Schallwandler erfolgt beispielsweise auf einem Pfad 24 durch das Fluid 20. In einer alternativen Ausführungsform weist die akustische Messvorrichtung 10 zwei Schallmesswandler 12 auf, von denen der eine dazu ausgebildet ist, Schallsignale auszusenden und der andere dazu ausgebildet ist, Schallsignale zu empfangen.
In Figur 1 ist weiter eine Vorrichtung 30 zum Betreiben der akustischen Messvorrichtung 10 dargestellt. Der Vorrichtung 30 sind Sensoren zugeordnet, die verschiedene Messgrößen erfassen und jeweils den Wert der Messgrößen ermitteln können. Die Vorrichtung 30 ist dazu ausgebildet, in Abhängigkeit von min ¬ destens einer der Messgrößen Werte zu ermitteln, mittels der ein Betrieb der akustischen Messvorrichtung 10 ermöglicht wird.
Im Folgenden soll die Funktionsweise der akustischen Mess- Vorrichtung 10 kurz dargestellt werden: Schallsignale werden von dem mindestens einen Schallmesswandler 12 ausgesendet und gelangen in das Fluid 20. Je nach der Strömungsgeschwindigkeit v des Fluids 20 werden unterschiedliche Laufzeiten der Schallsignale beginnend vom Aussenden des Schallsignals durch das Sende- und Empfangselement des min ¬ destens einen Schallmesswandlers 12 bis zum Empfang des Schall ¬ signals an dem Sende- und Empfangselement des mindestens einen Schallmesswandlers 12 erreicht. Dadurch ist eine Berechnung der Strömungsgeschwindigkeit v des Fluids möglich. Des Weiteren kann daraus auch ein Massenstrom des Fluids berechnet werden. Zu ¬ sätzlich kann durch die Auswertung der Schallgeschwindigkeit auf die Zusammensetzung des Fluids geschlossen werden.
Figur 2 zeigt einen zeitlichen Spannungsverlauf eines von dem mindestens einen Schallmesswandler 12 abgegebenen Ausgangssignals 40. Des Weiteren zeigt Figur 2 zeitliche Verläufe mehrerer Echosignale 42, die von dem mindestens einen Schall ¬ messwandler 12 empfangen werden können. Für eines der Echosignale 42 ist des Weiteren eine Laufzeit T_sig des Schallsignals vom Beginn des Ausgangssignals 40 bis zum Beginn des Echosignals 42 dargestellt .
In Figur 3 ist ein zeitlicher Verlauf des Echosignals 42 mit einer hohen zeitlichen Auflösung gezeigt. Das Echosignal 42 ist hier dargestellt normiert auf einen Betrag eines Werts einer maximalen Amplitude des Echosignals 42.
Figur 4 zeigt das normierte Echosignal 42 entsprechend der Figur 3 und den Verlauf eines Referenzsinussignals 44. Die Bedeutung des Referenzsinussignals 44 wird im Rahmen der Beschreibung zu der Figur 5 näher erläutert.
Ein in Figur 5 schematisch dargestelltes Programm, mittels dem die Vorrichtung 30 zum Betreiben der akustischen Messvorrichtung 10 betrieben wird, ist vorzugsweise auf einem Speichermedium der Vorrichtung 30 gespeichert. Das Programm wird vorzugsweise in einem Schritt S10 gestartet, in dem gegebenenfalls Variablen initialisiert werden. Dies geschieht vorzugsweise bei Inbetriebnahme der akustischen Mess ¬ vorrichtung 10.
In einem Schritt S12 wird das Ausgangssignal 40 von dem Schall ¬ messwandler 12 abgegeben und das Echosignal 42 des Ausgangs ¬ signals 40 von dem Schallmesswandler 12 empfangen. Der zeitliche Verlauf des Echosignals 42 wird erfasst. Die Bearbeitung des Programms wird anschließend in einem Schritt S14 fortgesetzt.
In dem Schritt S14 wird eine Normierung N des zeitlichen Verlaufs des Echosignals 42 auf einen Wert eines Betrags einer maximalen Amplitude des Echosignals 42 durchgeführt. Der derart normierte zeitliche Verlauf des Echosignals 42 ist in den Figuren 3 und 4 gezeigt . Die Bearbeitung des Programms wird anschließend in einem Schritt S16 fortgesetzt.
In dem Schritt S16 wird ein Wert tl eines Zeitpunkts eines ersten Nulldurchgangs des zeitlichen Verlaufs des normierten Echo ¬ signals 42 ermittelt (siehe auch Figur 3) . Ein erster Schätzwert T_sig_l einer Laufzeit des Schallsignals vom Beginn des Aus ¬ gangssignals 40 bis zum Beginn des Echosignals 42 wird bestimmt abhängig von dem Wert tl des Zeitpunkts des ersten Nulldurchgangs des zeitlichen Verlaufs des normierten Echosignals 42. Der erste Schätzwert T_sig_l repräsentiert einen ersten Näherungswert für die Laufzeit des Schallsignals vom Beginn des Ausgangssignals 40 bis zum Beginn des Echosignals 42. Die Bearbeitung des Programms wird anschließend in einem Schritt S18 fortgesetzt.
In dem Schritt S18 wird mittels einer Kreuzkorrelation zwischen dem Echosignal 42 und dem Referenzsinussignal 44 ein Verlauf einer Kreuzkorrelationsfunktion R, abhängig von einer Frequenz und einer Phase des Referenzsinussignals 44, ermittelt. Ins- besondere wird die Frequenz des Referenzsinussignals 44 derart variiert, dass sie in einem vorgegebenen Bereich um eine Frequenz des Ausgangssignals 40 liegt. Beispielsweise wird die Frequenz des Referenzsinussignals 44 derart variiert, dass sie in einem Bereich von +/- 5% um die Frequenz des Ausgangssignals 40 liegt. Vorzugsweise wird die Kreuzkorrelationsfunktion R für Phasenwerte 0 und Π/2 des Referenzsinussignals 44 ermittelt. Die Frequenz und die Phase des Referenzsinussignals 44 werden vorzugsweise gemeinsam variiert. Alternativ können die Frequenz und die Phase des Referenzsinussignals 44 in getrennten Schritten variiert werden. Beispielsweise kann zunächst die Frequenz des Referenzsinussignals 44 variiert werden, bis eine maximale Übereinstimmung der Frequenzen des Echosignals 42 und des
Referenzsinussignals 44 erreicht ist. In einem weiteren Schritt wird dann die Phase des Referenzsinussignals 44 variiert, bis eine möglichst gute Übereinstimmung zwischen dem Echosignal 42 und dem Referenzsinussignal 44 erhalten wird. Die Bearbeitung des Programms wird anschließend in einem Schritt S20 fortgesetzt.
In dem Schritt S20 wird ein Wert f_0 der Frequenz des Refe ¬ renzsinussignals 44, für den die Kreuzkorrelationsfunktion R einen Maximalwert erreicht, ermittelt. Der Wert f_0 der Frequenz des Referenzsinussignals 44, für den die Kreuzkorrelations ¬ funktion R einen Maximalwert erreicht, ist repräsentativ für eine maximale Übereinstimmung der Frequenz des Referenzsinussignals 44 mit der Frequenz des Echosignals 42. Des Weiteren wird ein Wert phi_0 der Phase des Referenzsinussignals 44, für den die Kreuz- korrelationsfunktion R den Maximalwert erreicht, ermittelt. Der Wert phi_0 der Phase des Referenzsinussignals 44, für den die Kreuzkorrelationsfunktion R den Maximalwert erreicht, ist repräsentativ für eine maximale Übereinstimmung der Phase des Referenzsinussignals 44 mit der Phase des Echosignals 42 (siehe auch Figur 4) . Die Bearbeitung des Programms wird anschließend in dem Schritt S22 fortgesetzt.
In dem Schritt S22 wird abhängig von dem Wert f_0 der Frequenz und dem Wert phi_0 der Phase des Referenzsinussignals 44, für den die Kreuzkorrelationsfunktion R einen Maximalwert erreicht, und dem ersten Schätzwert T_sig_l der Laufzeit des Schallsignals vom Beginn des Ausgangssignals 40 bis zum Beginn des Echosignals 42 ein Wert t2 des Zeitpunkts des Nulldurchgangs des Referenz ¬ sinussignals 44 bestimmt. Abhängig von dem Wert t2 des Zeitpunkts des Nulldurchgangs des Referenzsinussignals 44 wird ein ak ¬ tualisierter Schätzwert T_sig_2 der Laufzeit des Schallsignals vom Beginn des Ausgangssignals 40 bis zum Beginn des Echosignals 42 bestimmt. Vorzugsweise wird der aktualisierte Schätzwert T_sig_2 der Laufzeit des Schallsignals vom Beginn des Ausgangs ¬ signals 40 bis zum Beginn des Echosignals 42 mittels eines Iterationsverfahrens bestimmt. Die Bearbeitung des Programms wird anschließend in dem Schritt S24 fortgesetzt. In dem Schritt S24 wird das Programm beendet. Vorzugsweise wird das Programm jedoch während des Betriebs des Schallmesswandlers 20 regelmäßig abgearbeitet . Ein Vorteil besteht darin, dass der aktualisierte Schätzwert T_sig_2 der Laufzeit des Schallsignals vom Beginn des Ausgangs ¬ signals 40 bis zum Beginn des Echosignals 42 mittels der Kreuz ¬ korrelation des Echosignals 42 und des Referenzsinussignals 44 sehr genau bestimmt werden kann. Damit kann insbesondere darauf verzichtet werden, digitale Signalprozessoren mit hoch auflösender Abtastung und sehr schnelle Analog-Digital-Wandler für eine Bestimmung der Laufzeit des Schallsignals vom Beginn des AusgangsSignals 40 bis zum Beginn des Echosignals 42 einzusetzen. Vielmehr können preiswerte digitale Signalprozessoren und Analog-Digital-Wandler zum Einsatz kommen.
Next Patent: DEVICE FOR PROCESSING FISH