골전도(Bone Conduction, BC)를 야기하는 여러 요인 중 가장 주요한 요소로 중이의 관성, 달팽이관 내부의 유체 관성과 달팽이관 주변 뼈의 압축 효과가 주로 언급된다[1-3]. 이러한 주요 요소가 골전도에 미치는 영향을 정량화하고 골전도 메커니즘을 밝히고자 하는 실험에서는 Temporal Bone과 Cadaver Head가 주로 사용되어 왔다.

예를 들어 Stieger, et al. (2018)은 골전도 메커니즘을 설명하기 위해 Temporal Bone을 사용하여 Promontory 속도와 달팽이관 내압을 측정했다[4]. 그들은 Cochleostomy를 수행하여 Laser Vibrometer로 Promontory 속도를 측정하였다. Temporal Bone을 사용한 연구는 Stenfelt, et al. (2002)에 의해서도 수행되었다. 이들은 Temporal Bone을 Shaker에 고정하여 상하 진동을 이용한 골전도 자극을 가하고 이에 상응하는 Umbo와 Stapes Footplate의 움직임을 측정했다[5]. 그들은 이러한 실험 세팅을 이용하여 중이의 질량 증가에 따른 골전도 청력의 변화, 즉 골전도에 있어 중이의 관성의 영향을 정량화하여 보여주었다. 골전도 메커니즘의 이해를 돕기 위한 실험에는 Temporal Bone 외에도 Cadaver Head가 많이 사용된다.

Cadaver Head를 이용한 실험은 Eeg-Oloffson, et al. (2018)의 연구에서 볼 수 있다. 이들은 가장 큰 달팽이관 반응을 야기하는 골전도 자극 위치를 찾기 위해 Cadaver Head로부터 Mechanical Point Impedance와 Promontory 속도를 측정하였다. Promontory 속도를 측정하기 위해, 고막과 함께 Malleus와 Incus를 제거한 후 Promontory에 작은 유리 반사판을 부착하였다. 일반적인 골전도 보청기(Bone Conduction Devices, BCDs) 부착점과 Ear Canal과 가까운 위치 사이의 여러 위치에서 진동 자극을 가했으며 Laser Doppler Vibrometer로 Promontory 속도를 측정하였다. 이를 통해, 진동 자극을 표준 BCDs 위치보다 달팽이관에 더 가깝게 가하면 Promontory 속도가 증가한다는 것을 밝혔다[6].

이러한 실험 결과들이 골전도 메커니즘을 이해하는데 크게 기여했으나, 골전도의 또 다른 주요 요소인 달팽이관 주변 뼈의 압축 효과에 대한 설명 및 연구는 이에 비해 상대적으로 미흡했다. 뼈의 압축 효과는 뼈의 기계적 물성치에 매우 의존하므로, 달팽이관을 포함하는 두개골의 영률 측정에 관한 연구가 중요하다. 예를 들어, Auperrin, et al. (2014)은 Density Measurements, Geometrical Measurements, Three-Point Bending Tests와 Statistical Analysis를 통해 두개골의 평균적인 영률을 측정하였다[7]. 그 결과 Frontal Bone은 3.81±1.55, Parietal Bones는 5.00±3.12, Temporal Bones는 9.70±5.75 GPa이었다. 또, Rho, et al. (1993)이 측정한 Cortical Bone의 Tensile 영률은 18.6 GPa이었다[8]. 이처럼 두개골의 영률은 여러 값으로 측정되고, 이러한 영률의 차이가 골전도 청력에 미치는 영향, 특히 영률의 차이에 의해 Temporal Bone과 Cadaver의 뼈 압축 효과에 따른 골전도 변화에 대한 연구는 많이 이루어지지 않은 상황이다.

본 연구에서는 청각 기관을 포함하는 전체 Head의 3D 유한 요소(Finite Element, FE) 모델을 사용하여 영률 차이에 따른 Temporal Bone과 Cadaver Head 사이의 골전도 반응에 대한 영향을 조사하였다. 이 모델은 이경화증(Otosclerosis) 환자로부터 얻은 임상 데이터와 청력 상실을 비교하여 검증되었다. 또한, 이경화증 상태를 통해, 골전도에 대한 뼈의 영률에 따른 뼈의 압축 효과는 Temporal Bone과 Cadaver Head 모델 모두에서 조사되었다.

Fig. 3은 정상 상태와 이경화증 상태의 기저막 속도 차이를 통해 계산한 주파수별 청력 손실을 나타낸다. 이 청력 손실은 Figs. 3(a)부터 3(c)에서 흑색 점선으로 표시된 임상 데이터[12]와 비교되었다. 임상 데이터는 이경화증 환자의 골전도 청력을 측정한 Carhart의 데이터를 사용하였다. Fig. 3(a)는 Isolated 모델에 관성을 이용한 시뮬레이션 결과를 나타내며, 흑색과 적색 실선은 각각 뼈의 압축 효과가 없을 때와 있을 때를 나타낸다. Isolated 모델의 시뮬레이션에서 관성뿐만 아니라 뼈의 압축 효과가 있는 경우(적색 실선), 계산된 청력 손실은 약 10 dB의 차이 내로 임상 데이터와 유사함을 나타내었다. 하지만, 골전도 자극에 관성만을 사용한 경우, 600 Hz 미만의 저주파에서 임상 데이터와 비교하여 20-30 dB의 커다란 청력 손실 차이가 관찰되었다. 따라서, Isolated 달팽이관에 Shaker를 사용하여 실험을 할 경우, 추가적인 뼈의 압축을 고려하여 적절한 보상을 감안하지 않는다면 1 kHz 미만의 저주파수에서 얻어질 실험값이 다소 왜곡될 가능성이 있음을 추론할 수 있다.

Fig. 3 
Comparison of the hearing loss level calculated in the simulation with clinical data. The used FE models are (a) The isolated cochlear model, and (b)-(c) The full head model. In (a), the black-solid line represents the only inertia input whereas the red-solid line indicates the inertia and bone compression inputs. In (b) and (c), the blue-solid line means the case when the Young’s modulus for cortical bone and cochlear bone was 10 GPa. On the other hand, the green-solid line describes the case when the Young’s modulus was 20 GPa. In all the figures, the black-dotted line represents the clinical data obtained by Carhart. The BC input is the dynamic force in (b) whereas the input is the sinusoidal displacement in (c). The displacement in (c) was applied on a few nodes on the surface of skull

Figs. 3(b)와 3(c)의 청색과 녹색 실선은 각각 Cortical Bone과 Cochlear Bone의 영률이 10과 20 GPa인 경우를 나타낸다. Fig. 3(b)는 골전도 입력값으로 Sinusoidal Force를 가했을 경우 이경화증에 의한 청력 손실을 나타낸다. 뼈의 영률이 10 GPa일 경우, 1 kHz 이상에서는 이경화증 환자의 골전도 청력 손실이 크지 않아 임상 데이터와 비교했을 때 상대적으로 이경화증 환자로부터 청력 증가가 나타났다. 반면, 영률이 20 GPa의 경우에는, 1 kHz 이상에서 임상데이터와 10 dB 이내의 차이를 보이며 약 1.3 kHz 부근에서 Notch (Carhart’s Notch)가 나타났다.

Fig. 3(c)는 Head 모델을 Shaker에 고정하여 Sinusoidal Displacement 조건을 할당한 실험[14]과 같이 골전도 입력값으로 머리 유한요소 모델에 Sinusoidal Displacement를 적용하여 시뮬레이션을 수행했을 때 청력 손실을 나타낸다. 영률이 10 GPa(청색)일 때보다 20 GPa(녹색)일 때, 1 kHz 이하의 저주파에서 10-15 dB 정도의 보다 큰 청력 손실이 발생했다. 이는 주변 뼈의 영률이 클 때 뼈의 압축 효과가 더 적게 발생함을 보여준다. 다시 말해, 골전도 입력값이 BCDs와 같은 Sinusoidal Force일 경우 두개골의 영률이 고주파에서 골전도 청력에 큰 영향을 미치는 것을 알 수 있다.

이 결과를 실험 조건에 적용해 보았을 때, Cadaver Head를 사용한 골전도 실험은 두개골의 특성상 시간이 지날수록 굳어지기 때문에 영률 차이에 따른 골전도 청력에 대한 뼈의 압축 효과를 보다 신중하게 고려할 필요가 있다. 또한 이전 연구[7,8]에 따르면, 두개골의 영률은 약 2-19 GPa 범위로 보고되고 있다. 두개골 영률에 있어 존재하는 개인차는 골전도 청력에 영향을 끼칠 것으로 예측된다. 따라서, 정확한 실험값 도출이나 환자에 따른 올바른 BCDs 사용을 위해서는 각 개인의 두개골 영률을 고려하는 것이 필요할 것으로 사료된다.

본 연구는 머리뼈의 영률에 따른 골전도 반응 차이 분석을 위해 두개골에 직접 자극을 주는 Implantable Device를 가정하여 시뮬레이션을 진행하였다. 하지만, 최근 부착형 골전도 보청기도 많이 보급되는 것을 상기할 필요가 있다. 기본적인 에너지 전달 과정 및 머리뼈가 골전도에 미치는 영향은 삽입형과 부착형 서로 간에 큰 차이가 없을지라도 부착형 골전도 보청기에 대한 연구가 추후 진행될 필요가 있다.

본 연구에서는 두개골의 영률이 골전도 청력에 미치는 영향을 두 종류의 3차원 청각기관 유한요소 모델(Isolated Cochlear Model과 Head Model)을 이용하여 연구하였다. 이경화증으로 인한 청력 손실을 임상 데이터와 비교하여 모델을 검증하였고, 이를 통해 두개골의 영률이 골전도 청력에 미치는 영향, 특히 뼈 압축력이 골전도 청력에 미치는 영향을 살펴보았다. Temporal Bone에 Shaker가 사용된 골전도 실험의 경우, 저주파에서 뼈의 압축력이 더 고려되어야 하는 반면, Cadaver Head를 이용한 실험의 경우 고주파에서 뼈의 압축력이 지나치게 강조되지 않도록 주의를 기울여야 할 것이다.