JPS59108540A - 超音波診断装置 - Google Patents

超音波診断装置

Info

Publication number
JPS59108540A
JPS59108540A JP21693882A JP21693882A JPS59108540A JP S59108540 A JPS59108540 A JP S59108540A JP 21693882 A JP21693882 A JP 21693882A JP 21693882 A JP21693882 A JP 21693882A JP S59108540 A JPS59108540 A JP S59108540A
Authority
JP
Japan
Prior art keywords
echo
frequency spectrum
reflected
subject
echoes
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.)
Pending
Application number
JP21693882A
Other languages
English (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.)
Hitachi Ltd
Original Assignee
Aloka Co Ltd
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 Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP21693882A priority Critical patent/JPS59108540A/ja
Publication of JPS59108540A publication Critical patent/JPS59108540A/ja
Pending legal-status Critical Current

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 本発明は超音波診断装置、特に被検体へ向けて超音波ビ
ームを送波し得られる反射エコーから被検体の診断情報
を得る超音波診断装置に関する。
近年、超音波を利用して診断を行う超音波診断装置が広
く用いられている。この超音波診断装置は、体内の組織
や臓器が異なる音響的特性を有1゛ることを利用してい
る。すなわち、超音波ビームを一定周期で繰り返して体
内に放射すると、超音波は体内を伝搬する途中異なった
組織の境界でその一部が反射されるだめ、この反射エコ
ーを検出し、ブラウン管等に表示すれば、体内組織の音
響的特性の分布を表示させることができる。その際、健
全な組織と腫瘍組織とは音響的特性が異なるので、ブラ
ウン管上に表示されたパターンから、組織内に存在する
異常の有無及びその位置を知ることができる。
このように、超音波診断装置は、被検体組織の診断を非
観血的に、しかも被曝させることなく行うことができ、
更に所望被検部位の診断画像をリアルタイムで表示する
ことができる等の優れた特徴を有している。
第1図には、従来用いられていた超音波診断装置が示さ
れており、この超音波診断装置は、探触子10を被検体
の所定被検部位に向は接触させ超音波ビームの送受波を
行い、得られる受信エコーに基づき被検体所望被検部位
の断層面の診断画像を表示している。ここにおいて、探
触子10を用いた超音波ビームの送波は送信器12から
トリガパルスを探触子10に入力することにより行われ
、この超音波ビームの送波により得られる反射エコーの
受波は探触子10を介して受1言器14にて行われる。
このような超音波ビームの送受波は制御回路16により
一定の同14JJをとりつつ送信器12、受信器I4の
動作を制御することにより繰シ返して行われる。
ところで、反射エコーは探触子1oにおいて電気11τ
号に変換された後、受1言器14に入力されるが、この
ように、受信器14に入力される反射エコーはそのダイ
ナミックレンジが広いだめ、受信器14のダイナミック
レンジも反射エコーのダイナミックレンジに合わすた広
いものとする必要がある。このため、従来の装置におい
ては、受信器14を利得制御器18で制御し、反射エコ
ーの対数増幅を行っている。
そして1、受信器14で受信され対数増幅して出力され
た反射エコーは検波器20にて検波され、更に映像増幅
器22で画像処理され、ブラウン管24上に診断画像と
して画像表示される。
ところで、このような超音波ビームの送受波により得ら
れる受信エコー中には、被検体組織に関する各種の診断
情報、例えば、周波数情報、位相情報、振幅情報、伝搬
時間情報等が含まれている。
しかし、従来の超音波診断装置においては、反射エコー
中に含まれる各種診断情報の有効な活用が図られておら
ず、例えばAモード、Bモード、Mモードの超音波診断
装置では、振幅情報の一部及び伝搬時間情報の一部を利
用しているにすきず、またドプラ型の超音波診断装置で
は、これに加え、周波数及び位相情報の一部を利用して
いるにずき゛なかつた。このため、従来の超音波診断装
置では、得られる反射エコーから必ずしも十分な診断情
報を得ることができず、この面からの改良が必要とされ
ていた。
すなわち、従来の超音波診断装置では、第2図に示すよ
うに、反射エコー100a 、 100b、 100c
中に異なる周波数情報が含まれている場合でも、その波
連長LT、  および振幅が等しい場合には、検波器2
0を介して出力されるこれら反射エコー100a、 1
00b、 100cの電気信号200a、 200b、
 200cd、1rj1図に示すように、波連長LT 
に相当する時間りの波長をもつ等しい信号波形となって
L斗い、反射エコー100中に含捷れる各種の情報を有
効に利用することかできないという欠点があった。
本発明はこのような従来の課題に鑑みなされたもので、
その目的は反射エコーに含捷れる情報を(i幼に利用1
−1被検体の所定エコー分向における超)l波ビームの
M#R減衰周波数スペクトルを求めることかでき、この
音響減衰周波数スペクトルに基づいて被検体の診断を正
確に行うことを可能とする超音波診断装置を提供するこ
とにある。
この目的達成のだめ、本発明の装置は被検体へ向は超音
波ビームを送波し得られる反射エコーから被検体の診断
情報を得る超音波診断装置において、反射エコーをその
まま記憶する波形メモリと、波形メモリに記憶された反
射エコーから被検体内に存在する所定エコー源の前方及
び後方境界面における反射エコーを読み出しその周波数
スペクトルを求める周波数分析器と、周波数分析器によ
り求められたエコー源の前方及び後方境界面における周
波数スペクトルを差演算しエコー分肉における超音波ビ
ームの音響減衰周波数スペクトルを求める演算器と、を
備え、この演算器により求められる音響減衰周波数スペ
クトルに基づき被検体の診断を行うことを特徴とする。
次に本発明の好適な実施例を図面に基づき説明する。
まず最初に本発明の原理を第3図ないし第5図に基づき
説明する。
第3図には、被検体組織30−a、 30−b  内に
厚さLのエコー源32が混入している音響モデルが示さ
れている。
このような音響モデルに向け、探触子1oがら超音波ビ
ーム300を送波すると、この超音波ビーム300は組
織30−a内を伝搬していき、エコー源32と組織30
−a、 30−b  との境界面A、Bで順次反射され
、これら境界面A、Bで反射きれなかった一部の超音波
ビームは組織30−b内を透過波400として伝搬して
いく。
ここにおいて、各境界面A、Bにて順次反射された反射
工j −100a、 100bは、再度組織30−a内
を伝搬していき、受信用探触子1oにて受波される。
第4図には、このようにしで受波された反射エコー10
0a、 100bの信号波形が示きれており、第5図に
は、これら反射エコー100a、 100bの周M a
 スペクトルが示されている。そして、周波数スペクト
ルとして示されるこれら各反射エコー100a。
100b  中には以下に述べる各種の情報が含捷れて
いる。
すなわち、組織30−a、 30−b  の音響インピ
ーダンスをそれぞれzl、 Zaとし、エコー源32の
音響インピーダンスをZ2  とすれば、エコー源32
に向は送波された超音波ビーム300の各境界面A。
Bにおける反射率R1+ R2はそれぞれ次式で表わさ
れる。
このため、これら各境界面A、Bから反射される反射エ
コー1.00a、 100b中には、それぞれ各組織3
0−a、 30−b  及びエコー源32における音響
インピーダンスの比が情報として含まれている。
またこのように境界面A、Bにおいて反射される受信エ
コー100a、 100bには、探触子10から境界面
Atでの減衰情報が両者に同様に含まれている。寸だ組
織30−a内に薄膜組織が他のエコー源として存在して
いるような場合でも、これにより受1言エコー100a
、 100bの周波数スペクトル分布に変化が生ずるこ
とはない。
ま/ここのようにして得られる反射エコー100a。
100b  中には、超音波ビームを送受波する超音波
診断装置、特に探触子10の特性に起因する誤差成分が
均等に含まれている。
以上の各種情報及び誤差成分のほかに、境界面J3から
反射される受1gエコー100b  に1は、注目する
エコー源32内における超音波の減衰情報が含まれてい
る。従って、反射エコー100b  からこのエコー源
32に関する減衰情報を取り出すことができれば、この
減衰情報を用いてエコー源32に関する診断を行うこと
が可能となる。
しかし1.受信エコー100a、 100b Kけ、各
境界面A、Hにおける異なる音響インピーダンスの比が
含−まれているため、その絶対値を単に比較してもエコ
ーdぐ32に関する減衰情報を得ることができない。と
ころが、この絶対値の周波数分布に着目してみると、第
5図から明らかなように、各反射エコー100a、 1
00bに関する絶対値は受信エコー100b、すなわち
、遠い位置からの受[言エコーの方が高周波領域で大き
な減衰を受けているのが理解される。そこで、周波数ス
ペクトルによって示されるこのような減衰情報を利用し
て、エコー源32に関する超音波ビーム音響減衰周波数
スペクトルを求めることができる。すなわち、第5図に
示ス受信エコー100a  の周波数スペクトルから受
信エコー100b  の周波数スペクトルの差を求める
ことにより、第5図中に示される音響減衰周波数スペク
トル500を得ることができる。
なお反射エコー100a、 100b中に含まれる探触
子10から境界面Aまでの減衰情報及び探触子10の特
性に基づく誤差情報は各受信エコー1.OOa、。
100b  に均等に加わっているため、前述したよう
に、受信エコー100a  から受信エコー100b 
 を差演算することにより求めた超音波ビームの音響減
衰周波数スペクトル500からは、この組織30−a内
における減衰情報及び探触子10の特性に基づく誤差情
報が除去され、エコー源32内における減衰情報のみが
含まれていることとなる。
ここにおいて、得られる音響減衰周波数スペクトル50
0はエコー源32内における超音波ビームaV減衰が大
きいほどこの傾きが大きくなるため、受信エコー100
a、 100bの受信時間差しからエコー源32の厚さ
Lを目測すれば、この厚さLと音響減衰周波数スペクト
ル500の傾きとからエコー源32内における単位長当
たりの減哀緻周波数分布を得ることができる。
従って、このような原理を用いて、被検体動に生体の診
断を行う場合に、例えば癌細胞、腫瘍組織の超音波ビー
ムに対する音響減衰周波数スペクトルを求めておけは、
得られるエコー源32からの音響減衰周波数スペクトル
と癌細胞あるいは腫瘍組織の音響減衰周波数スペクトル
とを比較することにより、注目するエコー源32の診断
を正確に行うことが5丁能となる。
第6図にIt:l:、以北詳述した原理に基づき形成さ
れた本発明の装置の好適な実施例が示されている。
実施例において、本発明の超音波診断装置は従来の装置
を用いてブラウン管24上眞被検体の超音波断層画像を
描かせ、この断層画像から所望の被検部位を任意に選択
し、当該被検部位の診断情報を得るように構成されてい
る。なお第1図に示す従来装置と対応する部材には、同
一符号を付しその説明を省略する。
まず探触子10から所望の被検部位に向は超音波ビーム
を送波すると、当該被検部位からの反射エコーは探触子
10を介して受[8器14で受信され、更に検波器20
で検波され映像増幅器22で画像処理され、ブラウン管
24上に超音波断層画像として表示される。
またこれと同時に探触子10を介して得られる反射エコ
ーは増幅器40を介1−て波形メモリ42に入力され、
そのま捷の形で記憶される。このように、波形メモリ4
2には、検波器20を介さずに探触子10で得られた反
射エコーがその寸ま高周波信号波形を損われることなく
入力されるので、反射エコーは波形メモリ42にMHz
帯域の原波J杉の捷ま記憶されることになる。
なお、定形メモリ42に反射エコーを入)Jするに当た
って使用する増幅器40には、リニア増!隅2咎を用い
る。
そして、ザンプル個所指示器44を月4いて、フ゛ラウ
ン管24上に表示された超音波診断面1象力・ら所望の
被検部位Qをエコー源32として(旨示するど、その指
示信号は周波数分析器46に入ブ3される。ここにおい
て、被検部位Qの指定は、第3図に示ずように、エコー
源32の前方及び後ブ5境界面A、B17J2点を指示
する必要がある。J月波数分析器46はザンプル個所指
示器44力・ら入ブコされる指示信号に基つき、波形メ
モ1J42力1ら被4黄音昏位Qの境界+fiiA、B
における反射x −r−IQQB。
100b  φ〕読出を行う○ 周波数分析器46は、アナログ捷たはデ・ジタルに関す
る周知の手法を用い、波形メモリ42力・ら1メしみ出
じl’c受信エコー100a、 1()Obをノζワー
に関する周波数スペクトルに変換し、メモリ48に書べ
込み記憶する。
そして、演算器50は、メモリ48内に書き込み記憶さ
れた受信エコーLOOa、 100bに関する周波数ス
ペクトルを読み出し、これを差演算することにより、エ
コー源32内における超音波ビームの音響減衰周波数ス
ペクトル500を求める。ここにおいて、前述したよう
に、エコー源32内における単位長当たりの音響減衰周
波数を求めることができれば、エコー源32の診断の正
確を期する上で好ましい。このため、演′S−器50は
エコー源32の厚さLに関する情報を、距離計側器52
から入力し、エコー源32における単位長当たりの超音
波ビーム音響減衰周波数スペクトルを求め、表示器54
に表示している。ここにおいて、距離計測器52は゛、
エコー源32の厚さLを、サンプル個所指示器44の指
示位置から自動的に演算し、演算器50に入力している
なおこのような距離計測器52を用いずに、エコー源3
2の厚さLを計測する場合には、ブラウン管24−1=
に表示された診断画像から、捷だは、第4図に示すよう
に、得られる受信エコー100a。
100b  の時間差しから求めることもできる。この
場合、エコー源32の厚さが小さくなると、ブラウン管
上から言1測する場合には誤差−が大きくなるので、受
[言エコー100a  と受信エコー100bトの受I
S遅れ時間りを用いて計測するのが好ましい。
本発明に、以上の構成から成り、次にその作用を説明す
る。
捷ず探触子10を用いて被検体に超音波ビームを送受波
すると、ブラウン管24上には被検体断層面の超音波診
断画像が表示されるとともに、波形メモリ42内には探
触子10を介1−て検出された受信エコーがMHz帯域
の原波形のまま記憶される。
ここにおいて、ザンプル個所指示器44を用いてブラウ
ン管24に表示された超音波診断画像から、i9r望の
被検部位Qをエコー源32として指示!I−ると、周波
数分析器46は波形メモリ42に記憶されたMHz帯域
の受1言エコーから、エコー源:32のAil 、)J
及び後方境界面A、Bからの受信エコー100a、 1
00bを読み出し、これら受信エコー100a。
100b  から周波数スペクトルを取り出す。
演算器50はこのように周波数分析器46にて周波数ス
ペクトルに変換された受信エコー100a。
100b  を差演算することにより、エコー源32内
における超音波ビームの音響減衰周波数スペクトル50
0を求め、表示器54上にこれを表示する。
ここにおいて、周波数分析器46にて周波数スペクトル
に変換された受(言エコー1−’OOa、 100bに
は、探触子10とエコー源32の間において超音波ビー
ムが伝搬する際に減衰する減衰誤差情報及び装置特に探
触子10の特性に起因する誤差情報が含まれているが、
これら誤差情報は両スペクトルLOOa、 100bに
均等に含まれているだめ、この両スペクトル100a、
 100bを差演算することによって求められたエコー
源32内における超音波ビームの音響減衰周波数スペク
トル500からは、これら誤差情報が除去され、エコー
源32内における減衰情報のみが含まれることとなる。
従って、表示器54上に表示されるエコー源32内にお
ける超音波ビームの音響減衰周波数スペクトル500は
、エコー源32に関する診断情報を正確に示すこととな
り、例えば、生体各部の診断を行うに当たり、極めて有
効な診断情報を得ることが6丁能となる。
次に本実施例の装置を用いて行った実験結果を示す。
第7図ないし第10図には、探触子10として公称周波
a 2.25 NI−Izの探触子を用い、ゴム板をエ
コー源として行った実験データが示されている0第7図
には、厚さ2.35mmのゴム板1枚をエコー源とし、
このゴム板の表面からの受信エコー100a  と裏面
からの受信エコー100b  の周波数スペクトルがそ
れぞれ示されており、第8図には、この第7図に示す各
受信エコー100a、 100bの周波数スペクトルを
差演算することにより得られた音響減音周波数スペクト
ル500が示されている。
また第9図には、厚さ2.35mm のゴム板を2枚中
ねたものをエコー源と12、このゴム板表面からq〕受
イ1.エコー100a  と、このゴム板裏面からの受
信エコー100b  の周波数スペクトルがそれぞれ示
されてお9、第10図には第9図に示す各受信エコー 
100a、 100bを差演算することにより求めた音
響減衰周波数スペクトル500が示されている。
このように、厚さ2.35mm のゴム板1枚をエコー
源と口て得られた音響減衰周波数スペクトル500の傾
きに比し、その2倍の厚さを有するコ゛ム板をエコー源
として得られた音響減衰周波数スペクトル500の傾き
は約2倍となる。これはエコー源の厚さが約2倍となる
ことにより、そのエコー漂白における超音波ビームの減
衰もデンベル値で約2倍となることに起因するO 以上の実験データからも明らかなように、本発明の超音
波診断装置を用いて被検体へ向は超音波ビームを送受波
することにより、エコー源32内における減衰情報を正
確に表わす音響減衰周波数スペクトル500を求めるこ
とができる。
またと6ようにして求めたエコー源32に関する音響減
衰周波数スペクトル500を、エコー源32の厚さ情報
に基づき単位長当たりの汗響減哀周波数スペクトルとし
て表わすことも可能である。
従って、例えば、本発明の超音波診断装置を用いて生体
内の所望の組織の診断を行う場合には、当該組織が正常
な場合、及びその組織が所定の病気に冒されている場合
等の各種音響減衰周波数スペクトル500をデータ入力
しておけば、本発明の装置によって得られる所望エコー
源からの音響減衰周波数スペクトル500の傾きや曲率
を予めデータ入力されたスペクトルの傾きや曲率と比較
することにより、当該被検部位の診断を極めて正確にし
かも非観血的に行うことが可能となる。
なお前記実施例においては、探触子10を介して得られ
る受信エコーを波形メモリ42内に順次記憶していくも
のを示しだが、本発明はこれに限らず、例えば、サンプ
ル個所指示器44にて指示された所定被検部位Qの各境
界面A、Bにおける反射エコー100a、 100bの
みを記憶する構造とすることも可能であり、またこの2
点A、Bを結ぶ走査線上の受信エコーを記憶する構造と
することも可能である。
また本実施例においては、本発明の1診断対象として生
木を例にとり説明しだが、本発明の診断装置はこれに限
らず、超音波が伝搬する媒体であれば、媒体内に存在す
るエコー源の各種診断を行うことが可能である。
以」二説明したように、本発明によれば、反射エコーに
含まれる情報を有効に利用して、エコー漂白における超
音波ビームの減衰特性を表わす音響減衰周波数スペクト
ルを求め、これに基づき被検体の診断を極めて正確に行
うことが可能となる。
従って、本゛発明の超音波診断装置を生体の診断用に用
いれば、病変の早期発見等の各種診断を確実にかつ正確
に行うことが可能となる。
【図面の簡単な説明】
第1図は従来の超音波診断装置のブロック図、第2図は
第1図の装置において得られる反射エコーの検波の°前
後における波形説明図、第3図は本発明の詳細説明に用
いる音響モデルの説明図、第4図はこの音響モデルから
得られる受信エコーの波形図、第5・図は第4図に示す
受信エコーの周波数スペクトル特性図、第6図は本発明
の超音波診断装置の好適な実施例を示すブロック図、第
7図ないし第10図は第6図に示す装置を用いて行った
実験のデータ図である。 32・・・エコー源    42・・・波形メモリ46
・・・周波数分析器  50・・・演算器100a、 
100b−反射エコー 500・・音響減衰周波数スペクトル A、B・・・境界面。 出願人  アロカ株式会社 第7図 第8図 第9図 第10図

Claims (2)

    【特許請求の範囲】
  1. (1)被検体へ向は超音波ビームを送波し得られる反射
    エコーから被検体の診断情報を得る超音波診断装置にお
    いて、反射エコーをそのまま記憶する波形メモリと、こ
    の波形メモリに記憶された反射エコーから被検体内に存
    在する所定エコー源の前方及び後方境界面における反射
    エコーを読み出しその周波数スペクトルを求める周波数
    分析器と、周波数分析器により求められたエコー源の前
    方及び後方境w面における周波数スペクトルを差演算し
    エコー源内における超音波ビームの音響減衰周波数スペ
    クトルを求める演算器と、を備え、この演”;I器によ
    り求められる音響減衰周波数スペクトルにj、5つき被
    検体の診断を行うことを特徴とする超H)<波計断装置
  2. (2)特i+’l請求の範囲(1)記載の装置において
    、演算器(rl、J−コー源の厚さ情報に基づきエコー
    源内におけ、、:、 llj f、7長当たりの超g波
    ビーム音響減衰周波数スペクトルを求めることを特徴と
    する超音波診断装置。
JP21693882A 1982-12-13 1982-12-13 超音波診断装置 Pending JPS59108540A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP21693882A JPS59108540A (ja) 1982-12-13 1982-12-13 超音波診断装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP21693882A JPS59108540A (ja) 1982-12-13 1982-12-13 超音波診断装置

Publications (1)

Publication Number Publication Date
JPS59108540A true JPS59108540A (ja) 1984-06-23

Family

ID=16696276

Family Applications (1)

Application Number Title Priority Date Filing Date
JP21693882A Pending JPS59108540A (ja) 1982-12-13 1982-12-13 超音波診断装置

Country Status (1)

Country Link
JP (1) JPS59108540A (ja)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS56147082A (en) * 1980-04-17 1981-11-14 Yoshinori Hayakawa Method and device for quantitative ultrasonic inspection using plural frequencies
JPS57209040A (en) * 1981-05-29 1982-12-22 Tokyo Shibaura Electric Co Ultrasonic diagnostic apparatus

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS56147082A (en) * 1980-04-17 1981-11-14 Yoshinori Hayakawa Method and device for quantitative ultrasonic inspection using plural frequencies
JPS57209040A (en) * 1981-05-29 1982-12-22 Tokyo Shibaura Electric Co Ultrasonic diagnostic apparatus

Similar Documents

Publication Publication Date Title
US5361767A (en) Tissue characterization method and apparatus
US4043181A (en) Ultrasonic pulse-echo apparatus
EP0064399B1 (en) Ultrasonic measuring method
US4414850A (en) Measurement method and system utilizing ultrasonic wave
EP0154869B1 (en) Ultrasonic measurement apparatus
EP0146707B1 (en) Ultrasonic measurement method, and apparatus therefor
EP0139235A2 (en) Ultrasonic measurement method and apparatus therefor
US20050080329A1 (en) Ultrasound doppler diagnostic apparatus and image date generating method
JPH0649037B2 (ja) 超音波パルスエコーの皮相周波数決定方法及び皮相周波数計算装置
JPH0713631B2 (ja) 媒体の超音波エコーグラフィック検査方法及び装置
JPH021273B2 (ja)
JPH0715457B2 (ja) 超音波エコーグラフィ検査方法及び装置
JPS59108540A (ja) 超音波診断装置
US4546772A (en) Method and means for determining ultrasonic wave attenuation in tissue using phase locked loop
JPS6083645A (ja) 超音波断層測定方法およびその装置
JPS6053133A (ja) 超音波診断装置
JPS61135640A (ja) 超音波測定方法およびその装置
JPH0478298B2 (ja)
JPS5982836A (ja) 超音波診断装置
JPH0249103B2 (ja)
JPS6272336A (ja) 超音波組織診断装置
JPS61149131A (ja) 超音波診断装置
JPS60173463A (ja) 超音波測定装置
JPS61154546A (ja) 超音波診断装置
JPS60188842A (ja) 超音波探触子