Purpose – The purpose of this paper is to develop a new method based on operational matrices of Bernoulli wavelet for solving linear stochastic Itô-Volterra integral equations, numerically. Design/methodology/approach – For this aim, Bernoulli polynomials and Bernoulli wavelet are introduced, and their properties are expressed. Then, the operational matrix and the stochastic operational matrix of integration based on Bernoulli wavelet are calculated for the first time. Findings – By applying these matrices, the main problem would be transformed into a linear system of algebraic equations which can be solved by using a suitable numerical method. Also, a few results related to error estimate and convergence analysis of the proposed scheme are investigated. Originality/value – Two numerical examples are included to demonstrate the accuracy and efficiency of the proposed method. All of the numerical calculation is performed on a personal computer by running some codes written in MATLAB software.