Abstract: 
We develop a generic numerical scheme to simulate the paired adapted strong solution (i.e., a pair of vector random fields) to a generalized Cauchy (terminalvalue) problem, i.e., a unified system of backward stochastic partial differential equations (BSPDEs including infinitedimensional BSDEs) driven by Brownian motions. This unified system can be used to model many realworld system dynamics such as optimal control and differential game problems. To better understand our numerical scheme and mathematical regime in the wellknown artificial intelligence (AI) convolutional neural network (CNN), we use a CNN to model our numerical process by conditional expectation projection. It consists two convolution parts: W layers of backward networks and L layers of reinforcement iterations. Furthermore, it is a completely discrete and iterative algorithm in terms of both time and space with almost sure (a.s.) convergence supported by both theoretical proof and numerical examples. In doing so, the system is assumed to be highdimensional and vectorvalued and whose drift and diffusion coefficients may involve nonlinear and highorder partial differential operators. Under general local Lipschitz and linear growth conditions, the unique existence of the paired adapted strong solution to the system is proved by constructing a suitable Banach space to handle the difficulty that the partial differential orders on both sides of these equations are different. During the proof, we also develop new techniques of random field Malliavin calculus to show the unique existence of paired adapted strong solutions to two embedded systems of the first and second Malliavin derivative based BSPDEs under random environments. 
