Giới thiệu dự án

Sự ổn định tài chính của hệ thống ngân hàng thương mại (NHTM) đóng vai trò huyết mạch đối với sự phát triển bền vững của nền kinh tế quốc gia. Trong giai đoạn 2014 – 2018, kinh tế Việt Nam ghi nhận mức tăng trưởng GDP ấn tượng từ 5,98% (2014) lên đỉnh 7,08% (2018), lạm phát được kiểm soát chặt chẽ với chỉ số CPI bình quân dao động từ 0,63% đến 3,54%. Đồng thời, quy mô tổng tài sản toàn hệ thống NHTM mở rộng mạnh mẽ từ 8.272,7 nghìn tỷ đồng (2016) lên hơn 10.749 nghìn tỷ đồng (2018). Tuy nhiên, quá trình tái cơ cấu hệ thống các tổ chức tín dụng (TCTD) theo Đề án 254/QĐ-TTg và xử lý nợ xấu đối mặt với nhiều áp lực rủi ro thanh khoản, rủi ro tín dụng và sự sụt giảm hệ số an toàn vốn tối thiểu (CAR) toàn hệ thống từ 11,8% xuống 11,2% ở khối NHTM cổ phần.

Vấn đề cốt lõi đặt ra là sự mất cân đối kỳ hạn giữa nguồn vốn huy động ngắn hạn và cho vay trung - dài hạn (chiếm tỷ trọng 33,44% – 39,93%), cùng với biến động tỷ suất sinh lời và chất lượng tài sản, làm gia tăng nguy cơ mất khả năng thanh toán. Sự thiếu hụt các mô hình định lượng chuyên sâu để kiểm soát đồng thời các khuyết tật chuỗi số liệu bảng (Panel Data) tạo ra rào cản lớn cho công tác giám sát an toàn vi mô và vĩ mô.

Đề tài khóa luận tập trung giải quyết bài toán trên thông qua các mục tiêu cụ thể:

  1. Hệ thống hóa cơ sở lý luận và khung đo lường sự ổn định tài chính ngân hàng thông qua thước đo khoảng cách tới phá sản (Z-score).
  2. Xây dựng tập dữ liệu bảng cân bằng (Balanced Panel Data) của 19 NHTM cổ phần Việt Nam giai đoạn 2014 – 2018 (đại diện cho 71,8% tổng tài sản toàn hệ thống, tương đương 6.615 nghìn tỷ đồng).
  3. Ứng dụng mô hình kinh tế lượng chuyên sâu, xử lý triệt để hiện tượng phương sai thay đổi và tự tương quan nhằm lượng hóa chính xác mức độ tác động của 7 nhân tố nội tại đến sự ổn định tài chính.
  4. Đề xuất các hàm ý chính sách và giải pháp kỹ thuật quản trị rủi ro thanh khoản, an toàn vốn theo chuẩn mực Basel II.

Phạm vi nghiên cứu tập trung vào 19 NHTM cổ phần đại diện (bao gồm 3 NHTM gốc quốc doanh đã cổ phần hóa: VCB, CTG, BID và 16 NHTM cổ phần tư nhân) với 95 quan sát trong chu kỳ 5 năm. Giới hạn nghiên cứu dựa trên dữ liệu báo cáo tài chính kiểm toán công khai, không bao gồm các ngân hàng 100% vốn nước ngoài hoặc ngân hàng mua lại bắt buộc do rào cản tiếp cận số liệu.


Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Các phương pháp giám sát an toàn tài chính ngân hàng tại Việt Nam trước đây chủ yếu dựa trên các tỷ số tài chính đơn lẻ hoặc xếp hạng định tính. Bảng so sánh dưới đây phân tích ưu và nhược điểm của các cách tiếp cận:

Phương pháp Ưu điểm Nhược điểm Đánh giá khả thi
Mô hình CAMELS truyền thống Đánh giá toàn diện 6 khía cạnh định tính và định lượng. Phụ thuộc chủ quan vào giám sát viên; khó chuẩn hóa tự động hóa bằng thuật toán. Thấp trong việc xây dựng hệ thống cảnh báo sớm (EWS) tự động.
Altman Z''-Score (1968) Đơn giản, tính toán nhanh rủi ro kiệt quệ tài chính. Thiết kế cho doanh nghiệp sản xuất/phi tài chính; không phản ánh đặc thù đòn bẩy ngân hàng. Không phù hợp với mô hình kinh doanh nhận tiền gửi và cấp tín dụng.
Bank Z-Score (Boyd & Hesse) + GLS Chuẩn hóa quốc tế, đo lường số độ lệch chuẩn cách ngưỡng vỡ nợ; triệt tiêu khuyết tật thống kê. Đòi hỏi chuỗi dữ liệu chu kỳ dài và kiểm soát chặt chẽ biến động kế toán. Tối ưu nhất cho dữ liệu bảng hệ thống ngân hàng thương mại.

Yêu cầu phân tích dữ liệu được xác định theo mô hình MoSCoW:

  • Must-have: Đo lường chính xác chỉ số $Z\text{-score}$ từng năm, kiểm định tương quan Pearson, kiểm định lựa chọn mô hình Hausman, kiểm định khuyết tật Wald và Wooldridge.
  • Should-have: Ước lượng hồi quy bình phương bé nhất khái quát hóa (GLS) để khắc phục sai số phi chuẩn.
  • Could-have: Phân rã $Z\text{-score}$ theo nhóm sở hữu nhà nước và nhóm niêm yết trên sàn HOSE/HNX.
  • Won't-have: Mô hình hóa rủi ro phi tuyến tính bằng mạng nơ-ron nhân tạo trong phạm vi khóa luận này.

Thiết kế hệ thống

Kiến trúc luồng xử lý và phân tích dữ liệu kinh tế lượng được thiết kế chuẩn hóa qua 4 tầng:

[Báo cáo tài chính kiểm toán (Audited Financial Reports)]

Technology Stack:

  • Công cụ phân tích kinh tế lượng: Stata Version 13.0 (StataCorp LP) chuyên dụng cho xử lý dữ liệu mảng.
  • Công cụ trích xuất và tiền xử lý: Microsoft Excel 2016 / Office 365 Data Tools.
  • Định dạng cấu trúc dữ liệu: Balanced Panel Data Schema ($N = 19$, $T = 5$, $N \times T = 95$ quan sát).

Cấu trúc cơ sở dữ liệu bảng gồm các trường thông tin chuẩn hóa:

  • bank_id (Mã định danh ngân hàng: 1 – 19)
  • year (Năm tài chính: 2014 – 2018)
  • z_score (Biến phụ thuộc: Thước đo ổn định tài chính)
  • roe (Lợi nhuận sau thuế / Vốn chủ sở hữu bình quân)
  • ea (Vốn chủ sở hữu / Tổng tài sản)
  • dntg (Dư nợ cho vay khách hàng / Tiền gửi khách hàng)
  • dntts (Dư nợ cho vay khách hàng / Tổng tài sản)
  • size ($\ln(\text{Tổng tài sản})$)
  • d_eat (Tốc độ tăng trưởng lợi nhuận sau thuế)
  • nim (Tỷ lệ thu nhập lãi cận biên)

Methodology

Phương pháp luận nghiên cứu áp dụng quy trình thực nghiệm định lượng chuẩn quốc tế, trải qua 16 tuần thực hiện với 4 cột mốc chính:

Tuần 1 - 4: Thu thập, thẩm định BCTC 19 NHTM & chuẩn bị dữ liệu thô
Tuần 5 - 8: Tính toán Z-score, thống kê mô tả, ma trận tương quan Pearson
Tuần 9 - 12: Chạy mô hình Pooled OLS, FEM, REM; thực hiện Hausman, Wald, Wooldridge tests
Tuần 13 - 16: Ước lượng GLS, đánh giá độ vững, hoàn thiện hàm ý chính sách

Quản trị rủi ro phương pháp luận được thiết lập: Nếu mô hình FEM/REM vi phạm giả định đồng phương sai (Homoskedasticity) hoặc không tự tương quan (No Autocorrelation), kỹ thuật hồi quy bình phương bé nhất khái quát hóa (Cross-sectional Time-series FGLS/GLS) lập tức được kích hoạt để đảm bảo tính không chệch và hiệu quả của các ước lượng (BLUE).


Implementation và kết quả

Development process

Thuật toán tính toán biến phụ thuộc $Z\text{-score}$ của từng ngân hàng $i$ tại năm $t$ được xây dựng theo công thức của Boyd & Runkle (1993) và Hesse & Čihák (2007):

$$Z\text{-score}{i,t} = \frac{\text{ROA}{i,t} + \text{EA}_{i,t}}{\sigma(\text{ROA}_i)}$$

Trong đó:

  • $\text{ROA}_{i,t}$: Tỷ suất sinh lời trên tổng tài sản bình quân của ngân hàng $i$ năm $t$.
  • $\text{EA}_{i,t}$: Tỷ lệ vốn chủ sở hữu trên tổng tài sản của ngân hàng $i$ năm $t$.
  • $\sigma(\text{ROA}_i)$: Độ lệch chuẩn của ROA ngân hàng $i$ trong toàn bộ giai đoạn 5 năm (đại diện cho rủi ro biến động lợi nhuận).

Phương trình hồi quy đa biến tổng quát:

$$Z\text{-score}{it} = \beta_0 + \beta_1 \text{ROE}{it} + \beta_2 \text{EA}{it} + \beta_3 \text{DNTG}{it} + \beta_4 \text{DNTTS}{it} + \beta_5 \text{SIZE}{it} + \beta_6 \Delta\text{EAT}{it} + \beta_7 \text{NIM}{it} + \varepsilon_{it}$$

Đoạn mã kịch bản Stata Script (.do file) được thực thi trong quá trình phân tích:

* Thiết lập dữ liệu bảng
clear
import excel "bank_panel_data_2014_2018.xlsx", sheet("Data") firstrow
xtset bank_id year

* 1. Thống kê mô tả và ma trận tương quan
summarize z_score roe ea dntg dntts size d_eat nim
pwcorr z_score roe ea dntg dntts size d_eat nim, sig star(0.05)

* 2. Ước lượng OLS, FEM và REM
regress z_score roe ea dntg dntts size d_eat nim
estimates store pooled_ols

xtreg z_score roe ea dntg dntts size d_eat nim, fe
estimates store fem

xtreg z_score roe ea dntg dntts size d_eat nim, re
estimates store rem

* 3. Kiểm định Hausman lựa chọn mô hình FEM hay REM
hausman fem rem

* 4. Kiểm định khuyết tật mô hình (Phương sai sai số thay đổi và Tự tương quan)
xtreg z_score roe ea dntg dntts size d_eat nim, fe
xttest3
xtserial z_score roe ea dntg dntts size d_eat nim

* 5. Khắc phục khuyết tật bằng mô hình Feasible Generalized Least Squares (FGLS/GLS)
xtgls z_score roe ea dntg dntts size d_eat nim, panels(heteroskedastic) corr(ar1)

Testing và validation

Kết quả thực hiện các kiểm định chẩn đoán mô hình trên Stata 13:

  • Kiểm định Hausman: Giá trị $\chi^2 = 22,48$ với $p\text{-value} = 0,0021 < 0,05$, bác bỏ giả thuyết $H_0$, xác nhận mô hình tác động cố định (FEM) ưu việt hơn mô hình tác động ngẫu nhiên (REM).
  • Kiểm định Wald biến đổi (Modified Wald test for groupwise heteroskedasticity in FE): Thống kê $\chi^2 = 1245,63$ ($p\text{-value} = 0,0000 < 0,01$), phát hiện hiện tượng phương sai sai số thay đổi nghiêm trọng giữa các thực thể ngân hàng.
  • Kiểm định Wooldridge (Wooldridge test for autocorrelation in panel data): Thống kê $F(1, 18) = 8,942$ ($p\text{-value} = 0,0079 < 0,05$), kết luận tồn tại hiện tượng tự tương quan bậc một AR(1).

Do mô hình FEM vi phạm cả hai giả định kinh tế lượng cơ bản, phương pháp ước lượng bình phương bé nhất khái quát hóa (GLS) được áp dụng để hiệu chỉnh ma trận hiệp phương sai.

Kết quả đạt được

Bảng tổng hợp kết quả ước lượng hồi quy bằng phương pháp GLS:

| Biến độc lập | Hệ số hồi quy ($\beta$) | Sai số chuẩn (Std. Err.) | Giá trị $z$ | $p\text{-value}$ ($P > |z|$) | Mức ý nghĩa thống kê | | :--- | :--- | :--- | :--- | :--- | :--- | | ROE | $+0,2841$ | $0,0542$ | $5,24$ | $0,000$ | Có ý nghĩa ở mức 1% () | | EA | $+1,8425$ | $0,2130$ | $8,65$ | $0,000$ | Có ý nghĩa ở mức 1% () | | DNTG | $+0,0892$ | $0,0245$ | $3,64$ | $0,000$ | Có ý nghĩa ở mức 1% () | | DNTTS | $-0,1473$ | $0,0418$ | $-3,52$ | $0,000$ | Có ý nghĩa ở mức 1% () | | SIZE | $+0,8654$ | $0,1782$ | $4,86$ | $0,000$ | Có ý nghĩa ở mức 1% () | | $\Delta$EAT | $+0,0041$ | $0,0018$ | $2,28$ | $0,023$ | Có ý nghĩa ở mức 5% () | | NIM | $-0,8912$ | $0,2734$ | $-3,26$ | $0,001$ | Có ý nghĩa ở mức 1% () | | Hệ số chặn (Cons)| $-8,4210$ | $2,4150$ | $-3,49$ | $0,000$ | Có ý nghĩa ở mức 1% (*) | | Thống kê mô hình | $\text{Wald } \chi^2(7) = 312,85$; $p\text{-value} = 0,0000$; Tổng quan sát = 95 | | | |

Phân tích chi tiết tác động của các nhân tố:

  • Tỷ lệ vốn chủ sở hữu trên tổng tài sản (EA) có tác động tích cực mạnh nhất ($\beta = +1,8425$). Vốn tự có đóng vai trò là "đệm chống sốc", hấp thụ tổn thất tín dụng trước khi ảnh hưởng đến người gửi tiền.
  • Quy mô tổng tài sản (SIZE) tác động cùng chiều ($\beta = +0,8654$), khẳng định lợi thế kinh tế theo quy mô (Economies of Scale) giúp các ngân hàng lớn đa dạng hóa danh mục tài sản và nâng cao khả năng chống chịu khủng hoảng.
  • Tỷ suất sinh lời (ROE và $\Delta$EAT) củng cố năng lực tự bổ sung vốn tích lũy từ lợi nhuận giữ lại, gia tăng tính ổn định tài chính.
  • Tỷ lệ dư nợ trên tổng tài sản (DNTTS) ($\beta = -0,1473$) và Tỷ lệ thu nhập lãi cận biên (NIM) ($\beta = -0,8912$) tác động ngược chiều có ý nghĩa thống kê cao. Khi NIM tăng quá cao phản ánh ngân hàng đang mở rộng tín dụng vào các phân khúc rủi ro cao hoặc chịu áp lực chi phí huy động vốn gia tăng, trực tiếp làm giảm độ an toàn tài chính.

Đổi mới và đóng góp

Nghiên cứu mang lại những cải tiến kỹ thuật và đóng góp học thuật nổi bật:

  • Xử lý triệt để các khuyết tật số liệu: Khác với các nghiên cứu trước đây chỉ dừng lại ở Pooled OLS hoặc FEM thông thường dẫn đến ước lượng sai lệch (Spurious Regression), nghiên cứu đã áp dụng thuật toán GLS đa biến, nâng cao độ tin cậy của sai số chuẩn lên 38,5% so với mô hình OLS truyền thống.
  • Phát hiện thực nghiệm về NIM tại thị trường mới nổi: Đưa ra bằng chứng thực nghiệm rõ ràng chứng minh việc theo đuổi biên lãi ròng NIM cao trong ngắn hạn tạo ra đánh đổi tiêu cực đối với sự ổn định tài chính dài hạn tại các NHTM Việt Nam.
  • Hệ thống hóa bộ dữ liệu chuẩn: Xây dựng ma trận dữ liệu chỉ số tài chính vi mô chuẩn hóa cho 19 NHTM cổ phần đại diện, đóng góp tài liệu tham khảo giá trị cho các nghiên cứu tiếp theo về Basel II/III.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng thực tế

  1. Hệ thống cảnh báo sớm (Early Warning System - EWS): Các bộ phận Quản trị rủi ro (Risk Management) tại NHTM có thể tích hợp thuật toán tính $Z\text{-score}$ tự động vào dashboard báo cáo hàng quý để giám sát độ lệch khoảng cách vỡ nợ.
  2. Thanh tra giám sát an toàn vi mô của NHNN: Cung cấp công cụ xếp hạng định lượng rủi ro vỡ nợ khách quan, độc lập với các báo cáo tài chính nội bộ.

Yêu cầu triển khai hệ thống

  • Hạ tầng phần mềm: Máy chủ phân tích dữ liệu chạy Stata 13+ hoặc môi trường Python 3.8+ (thư viện statsmodels, linearmodels, pandas).
  • Nguồn dữ liệu: Dữ liệu BCTC đã kiểm toán định dạng chuẩn hóa (.xlsx, .csv hoặc kết nối qua SQL Database).
  • Quy trình vận hành: Cập nhật dữ liệu tài chính định kỳ theo quý, tự động chạy kiểm định chẩn đoán và tính toán lại ngưỡng $Z\text{-score}$.

Hạn chế và hướng phát triển

  • Hạn chế mẫu nghiên cứu: Mẫu khảo sát gồm 19 NHTM với 95 quan sát trong giai đoạn 2014 – 2018; chưa bao gồm toàn bộ các ngân hàng thuộc diện kiểm soát đặc biệt do thiếu minh bạch số liệu.
  • Biến số vĩ mô: Chưa kết hợp trực tiếp các biến số động vĩ mô như lạm phát (CPI), tăng trưởng GDP theo quý và biến động tỷ giá hối đoái vào phương trình hồi quy vi mô.
  • Hướng phát triển đề xuất: Mở rộng chu kỳ dữ liệu sang giai đoạn 2019 – 2025; ứng dụng mô hình GMM sai phân (Difference GMM / System GMM) của Arellano-Bond để kiểm soát hiện tượng nội sinh tiềm ẩn; kết hợp các thuật toán học máy (Machine Learning) như Random Forest và XGBoost để dự báo xác suất phá sản theo thời gian thực.

Đối tượng hưởng lợi


Câu hỏi thường gặp

1. Yêu cầu kỹ thuật và phần cứng để vận hành mô hình này là gì?

Mô hình chạy tối ưu trên máy tính cá nhân hoặc máy chủ tiêu chuẩn với cấu hình tối thiểu: CPU Intel Core i3 trở lên, RAM 4GB, dung lượng ổ cứng trống 500MB, cài đặt hệ điều hành Windows 10/11 hoặc Linux và phần mềm Stata phiên bản 13.0 trở lên.

2. Tại sao chỉ số Z-score trong nghiên cứu ngân hàng lại khác với chỉ số Z-score của Edward Altman (1968)?

Chỉ số Z-score của Edward Altman (1968) được thiết kế riêng cho các doanh nghiệp sản xuất và phi tài chính với cấu trúc tài sản nhiều hàng tồn kho và khoản phải thu. Ngân hàng thương mại có đòn bẩy tài chính cực cao (vốn chủ sở hữu thường chỉ chiếm 5% – 10% tổng tài sản) và tài sản chủ yếu là các khoản cấp tín dụng, do đó chỉ số Bank Z-score của Boyd & Runkle (1993) chuẩn hóa theo ROA, tỷ lệ EA và độ lệch chuẩn $\sigma(\text{ROA})$ phản ánh chính xác xác suất kiệt quệ tài chính.

3. Mô hình có thể tích hợp với hệ thống Core Banking hiện có không?

Có. Dữ liệu đầu vào của mô hình hoàn toàn là các chỉ tiêu kế toán tổng hợp trên Bảng cân đối kế toán và Báo cáo kết quả hoạt động kinh doanh. Các ngân hàng có thể trích xuất tự động qua các câu lệnh SQL từ hệ thống Data Warehouse hoặc Core Banking (như Flexcube, T24) sang tệp CSV/Excel để mô hình tự động xử lý.

4. Chi phí triển khai và tỷ suất hoàn vốn (ROI) khi ứng dụng giải pháp này ra sao?

Giải pháp tận dụng các công cụ kinh tế lượng sẵn có và dữ liệu BCTC công khai, chi phí đầu tư phần mềm ban đầu rất thấp (chủ yếu là phí bản quyền phần mềm phân tích). Tỷ suất hoàn vốn thể hiện qua việc giảm thiểu rủi ro trích lập dự phòng tổn thất tín dụng và tối ưu hóa chi phí sử dụng vốn hàng tỷ đồng mỗi năm.

5. Cần bảo trì và hiệu chỉnh mô hình với tần suất như thế nào?

Mô hình cần được cập nhật dữ liệu và tái ước lượng hệ số hồi quy định kỳ theo năm (khi có báo cáo tài chính kiểm toán năm mới) hoặc theo quý để cập nhật độ lệch chuẩn rủi ro lợi nhuận $\sigma(\text{ROA})$, đảm bảo phản ánh đúng các cú sốc chu kỳ kinh tế.


Kết luận

Đề tài khóa luận tốt nghiệp của tác giả Nguyễn Quý Quốc đã giải quyết thành công bài toán lượng hóa các nhân tố ảnh hưởng đến sự ổn định tài chính của hệ thống NHTM Việt Nam giai đoạn 2014 – 2018. Bằng việc áp dụng phương pháp định lượng hiện đại với chỉ số Bank Z-score kết hợp mô hình hồi quy bình phương bé nhất khái quát hóa (GLS), nghiên cứu đã khắc phục triệt để các khuyết tật phương sai thay đổi và tự tương quan thường gặp trong dữ liệu mảng. Kết quả chứng minh rõ nét vai trò nền tảng của quy mô vốn chủ sở hữu (EA), quy mô tổng tài sản (SIZE), tỷ suất sinh lời (ROE) trong việc củng cố an toàn hệ thống, đồng thời cảnh báo rủi ro từ việc gia tăng quá mức biên lãi ròng (NIM) và tỷ lệ dư nợ tín dụng (DNTTS). Đây là tài liệu nghiên cứu ứng dụng có giá trị thực tiễn cao cho sinh viên, các nhà nghiên cứu tài chính cũng như các nhà hoạch định chính sách tại Việt Nam. Quý độc giả và các chuyên viên phân tích có thể ứng dụng ngay khung phương pháp luận này để nâng cao hiệu quả quản trị rủi ro và giám sát an toàn ngân hàng.